Transcriptomics: Open Access

Transcriptomics: Open Access
Open Access

ISSN: 2329-8936

Research Article - (2014) Volume 2, Issue 1

The Role of TLR2 in Cigarette Smoke-Induced Gene Induction

Wright WR1, Mazi SI1,3, Parzych K1, Mataragka S1, Mein C2, Galloway-Phillipps N1, Mitchell JA1 and Paul-Clark MJ1*
1Department of Cardiothoracic Pharmacology, Pharmacology and Toxicology, National Heart and Lung Institute, Imperial College London, Guy Scadding Building, Dove house Street, London, UK
2Genome Centre, Barts and The London School of Medicine and Dentistry, John Vane Centre, Charterhouse Sq., London, UK
3King Fahad Cardiac Center, King Khaled University Hospital , College of Medicine, King Saud University, Riyadh 11472, Saudi Arabia
*Corresponding Author: Paul-Clark MJ, Department of Cardiothoracic Pharmacology, Pharmacology and Toxicology, National Heart and Lung Institute, Imperial College London, Guy Scadding Building, Dove House Street, UK, Tel: +44 (0)29 2087 5107 Email:

Abstract

Background: Previous work form our group has shown that innate immune cells perceive cigarette smoke as an inflammatory and xenobiotic stimulus, which alters the immune response to invading pathogens. It is unclear as to the specific cellular and molecular mechanisms driving cigarette smoke-induced inflammation, but we have demonstrated that this response is, in part, driven by Toll-Like Receptor (TLR) 2 in innate immune cells in vitro and in mice in vivo . Methodology/Principle findings: To address the impact TLR2/6, has on cigarette smoke-induced gene induction, HEK293 cells were permanently transfected with TLR2/6 or Null vector, and were stimulated with 10% cigarette smoke extract (CSE) for 8 h. Total RNA was extracted and the transcriptome analysed using Illumina Bead Chip arrays. In HEK293 Null cells, CSE induced 33 genes and down-regulated 41. In HEK293 TLR2/6 cells, CSE induced 23 genes and down-regulated 44. Further analysis revealed that 42 genes were regulated in a TLR2/6-dependent manner. Comparison of these genes with those induced by smoke in human primary monocytes revealed that 5 were mutually regulated. The major pathways affected were those associated with anti-oxidant pathways, tumorgenesis and cell survival. Conclusions: Our data suggest that the innate immune receptor TLR2/6 has a critical role for the expression of a particular cassette of genes induced by cigarette smoke. Pathway analysis indicates these are related to functions in both cell survival and tumorgenesis. Future validation of the relative importance of these pathways using more complex models is required a may lead to improved understanding of the pathology of cigarette smoke-induced diseases such as COPD, cardiovascular disease and cancer.

<

Introduction

Cigarette smoke contains in excess of 4,000 compounds [1], many of which are highly toxic to the cell and include many oxidants. The outcomes of cigarette smoke on immunity are broad and complex, with both pro-inflammatory and suppressive effects being reported [2]. Our group and others have shown that exposure of cells to cigarette smoke results in cellular oxidative stress and can result in the damage of proteins, DNA, and lipids [3,4]. It has also been shown that cigarette smoke promotes autoimmunity and the modification of antigen presentation [5], which, in blood from smokers, can manifest itself in the sensitisation of blood to pathogens [6]. The mechanism by which cigarette smoke causes inflammation has been attributed to activation of the transcription factors NF-κβ and AP-1 [7]. Binding of these transcription factors to response elements in DNA leads to de novo synthesis and release of pro-inflammatory mediators such as CXCL8 and TNF-α [8-11].

Toll-like receptors (TLRs) are innate immune receptors that are part of a larger family of pattern recognition receptors (PRRs) that recognise conserved sequences known as pathogen-associated molecular patterns (PAMPs), and the more recently identified dangerassociated molecular patterns (DAMPs) [12]. DAMPs consist of a growing list of molecules that include hyaluronan, fibrinogen and oxidants that contribute to what is now termed “sterile inflammation” [12]. The inflammation induced by oxidants contained in cigarette smoke has been shown both in vivo and in innate immune cells in vitro to be sensed, in part, through TLR2 and TLR4 [13,14]. There is also increasing evidence that cigarette smoke detrimentally affects viral TLR signaling on innate immune and lung epithelial cells [15-17].

There have been a number of studies assessing gene changes in smokers compared with healthy non-smokers, typically identifying changes in inflammatory, oxidative stress and xenobiotic genes [18-20]. Moreover, when evaluating changes in gene expression of macrophages and airway epithelial cells in response to cigarette smoke extract (CSE), the major genes altered were also associated with oxidative stress and xenobiotic pathways [21-23]. Previous work in our group used a transcriptomics approach to assess gene expression changes induced by acute cigarette smoke exposure in the human monocytic cell line, THP-1, and in primary human peripheral blood mononuclear cells (PBMCs). It revealed that THP-1 cells and PBMCs respond to cigarette smoke extract by up-regulating oxidative stress response and inflammatory pathways [24]. We have also shown that oxidants present in cigarette smoke induce an inflammatory response in vivo, and this is, in part, regulated by TLR2 and TLR4 [14]. Furthermore, in this study we showed that smoke-induced CXCL8 release from human monocytes was TLR2/6 dependent. However, the role of TLR2/6 in smoke-induced transcriptomic changes is not known.

In the current study, therefore, we used basic human cell models to investigate the role of TLR2/6 receptors in transcriptomic effects induced by CSE. To elucidate which genes were controlled in a TLR2/6- dependent and -independent manner, we examined the transcriptome of HEK293 cells that were transfected with TLR2/6, and compared these with control vector-transfected null cells, after stimulation with CSE. These cells offer a stable and well-established method for assessing the effects of agents of interest on various responses. Moreover, to contextualize and assess using HEK293 cells as a model for this purpose, we compared our results to previously analyzed data in THP-1 monocytes and PBMCs. This may help us to identify target genes that are TLR2/6-dependent, and can affect inflammation caused by cigarette smoke.

Materials and Methods

Materials

All cell culture plastics and general disposables were obtained from Fisher Scientific (Loughborough, Leicestershire, UK). Invitrogen (Paisley, Renfrewshire, UK) supplied all cell culture reagents, unless stated otherwise. General laboratory reagents were purchase from Sigma-Aldrich (Poole, Dorset, UK).

Cell Culture

Stably transfected HEK293 cells with pUNO-mcs (Null control), or pDUO-hTLR2/6, were purchased from In vivo Gen (CA, USA) and maintained according to the manufacturer’s instructions.

Preparation of Cigarette Smoke Extract (CSE)

To prepare CSE, four full strength Marlboro cigarettes (filters removed) were combusted through a modified 60ml syringe apparatus, and the smoke passed through 100ml of RPMI 1640. Each cigarette yielded 5 draws of the syringe (to 60 ml mark), with each individual draw taking approximately 10 seconds to complete. CSE was then passed through a 0.25 μm filter to sterilise and remove particulate matter and was used immediately at a 10% concentration diluted in control media [25]. As described previously [26] smoke extract “strength” was evaluated by measuring nitrite using the Griess reaction to ensure continuity between batches. In all experiments nitrite levels in 100% cigarette smoke extract was between 12 and 16μM. CSE (1- 100%) made using this method does not contain detectable levels of bacterial antigens for TLR4 (endotoxin; measured by the E-toxate™ kit, Sigma, UK) or TLR2 (LTA; measured by in house ELISA).

Treatment of Cells

TLR2/6 cells and nulls were plated at a density of 1×106 cells/ml onto 6-well plates and allowed to equilibrate for 16h before being treated for 8h with fresh 10% CSE or DMEM [14,41]. After 8 h, RNA was then extracted using RNeasy kit (Qiagen, Sussex, UK) as described in the manufacturers’ handbook. A control and a CSE sample set were prepared on each day for three consecutive days using fresh reagents on each day.

Gene Array

Total RNA was subject to standard microarray procedures. Samples were converted to cDNA, labelled, and hybridized to the Illumina HumanRef-8v3 Bead Chip array (Illumina, UK) at St Bartholomew’s and The London Genome Centre (BLGC), Queen Mary, University of London. Quality control and basic interpretation of data was also performed at BLGC before datasets were received for further analysis.

PCR verification of target genes regulated by CSE in a TLR2/6-dependant manner

THP-1 cells and PBMCs were treated with 10% CSE for 8 and 24h. Total RNA from these cells was extracted using a Qiagen RNeasy™ mini kit with a DNAase treatment step as instructed by the manufacturer. The purity of RNA and concentration of RNA was determined using a Nanodrop ND-1000 spectrophotometer (Thermo Scientific, Basingstoke, UK). RNA was converted to cDNA using a qScript™cDNA synthesis kit (Quanta Biosciences, Gaithersburg, MD, UK) and probed for PPM1F (Protein Phosphatase, Mg2+/Mn2+ Dependent, 1F), SRGN (Serglycin), EIF2C2 (Eukaryotic Translation Initiation Factor 2C, 2) and GAPDH using off the shelf TaqMan® probes under conditions stated by the manufacturer (Invitrogen, CA, USA). All levels of genes assessed were normalised against GAPDH and reported as fold change compared with control.

Data analysis

Analysis of datasets was performed using GeneSpring GX 11.3 (Agilent Technologies). Raw data were pre-processed to remove variability across and within array samples. To minimize non-biological variability across arrays, raw data were first log2 transformed and then quantile normalized [27], which is the recommended normalization algorithm for Illumina Bead Chip array analyses [28]. Normalization at the level of genes was performed on GeneSpring GX 11.3.

Samples were sorted into conditions based on the treatment applied and the cell type: HEK293 null controls (1-3) and HEK293 null 10% CSE (1-3) and HEK293 TLR2/6 controls (1-3) and HEK293 TLR2/6 10% CSE (1-3). A very stringent filtering of the dataset was performed, by selecting only the genes that had detectable signal intensity value in all samples (filter by flags present in all samples). Data from the THP-1 cells and PBMCs were previously published [24], but were reanalysed for this study using the moderated t-test. Fold change differences between control- and CSE-treated samples were calculated using cut-offs of 1.5-fold for statistically significant genes. These genes were identified using moderated unpaired t-test (p<0.05) with Benjamini-Hochberg False Discovery Rate (FDR) correction for single group comparisons [29]. Where more than one group were analysed, a one-way analysis of variance followed by a Tukey’s post-hoc test and Benjamini-Hochberg FDR correction was used. Significantly differentially expressed genes at each fold-change cut-off were used to generate hierarchical clustering plots using Pearson’s centred correlation and Ward’s linkage rule, and were displayed as heat maps.

Data sets consisting of genes significantly altered by ≥1.5-fold (compiled on GeneSpring GX11.3) were uploaded onto Ingenuity Pathway Analysis (IPA; Ingenuity® Systems, www.ingenuity.com) and mapped to Ingenuity’s Knowledge Base. The significance of the association between the dataset and the canonical pathways were measured using a ratio of number of genes from the dataset that mapped to the pathway divided by the total number of genes in that pathway, and p-value was generated using Fisher’s exact test. A threshold of 0.10 was used to indicate canonical pathways that are significantly represented by genes in a gene list.

Database searches using PubMed, GeneCards and other sources were performed to identify differentially expressed genes from the current study that were also found to be differentially expressed or important in the response to cigarette smoke in previous studies.

Results

Pre-Processing and analysis of raw and normalised global gene expression data

A single time point of 8 h was chosen based on previously published data showing this time point as appropriate to demonstrate both transcriptomic changes [24] and expression of inflammatory and antioxidant proteins of interest [7]. HEK293 Null and TLR2/6 positive cells were treated for 8 h with RPMI medium (n=3) ± 10% CSE-conditioned medium. A 10% CSE concentration was chosen because a 3% solution failed to induce CXCL8 protein release from cells, and a 20% solution resulted in >50% cell death. Extracted RNA was hybridised onto Illumina HumanRef-8v3 Bead Chip arrays and analysed for global changes in the expression of genes in the transcriptome.

Box-whisker plots (not shown) of raw signal intensity values for all samples revealed that the overall distribution of intensities was similar in both HEK293 Null and HEK293 TLR2/6 cells. Quantile normalization of the dataset made the distribution profiles of all samples identical (data not shown). Stringent filtering was used to select only probe sets that were measured as “present” or “marginal” in all samples. This resulted in 22837 probe sets for HEK293 Nulls and 21372 probe sets for HEK293 TLR2/6 cells, which were suitable for further analysis. All data are MIAME compliant and all results have been deposited on the Array Express database (http://www.ebi.ac.uk/cgi-bin/microarray/magetab.cgi) under the experimental title “inflammatory transcriptome profiling of human monocytes exposed acutely to cigarette smoke” as an additional data set named “HEK293 data”.

Gene expression patterns in HEK293 Nulls and HEK293 cells treated with 10% CSE or control media

Unsupervised principal component analysis was performed; using 3 principal components selected using the elbow method and is represented as a 3D scatter plot (Figure 1).The samples clearly stratify between cell type and treatment, indicating detectable differences in the total gene transcription patterns of these two groups. The influence of individual components on global transcriptome gene expression is as follows: X-axis, component 1 (cell treatment) is 59.5% Y-axis, component 2 (cell type) is 18.5% Z-axis, component 3 (n number) is 12.1%.

transcriptomics-Distinct-clustering

Figure 1: HEK293 cells treated for 8h with control media or 10% CSE cluster according to treatment conditions. Principal Component Analysis (PCA) plot on HEK293 Null (blue) and TLR2/6 (red) samples treated for 8h with control media (squares) or 10% CSE (triangles). Distinct clustering by cell type/treatment group indicates that global gene expression is altered by transfection with TLR2/6 and by stimulation with cigarette smoke. Data are n=3, and percentage influence of each component was 59.5% (PC1), 18.5% (PC2), 12.1% PC3.

Effect of 10% CSE on the expression of genes in HEK293 Null cells

Fold change analysis followed by a moderated student t-test and Benjamini-Hochberg FDR correction identified 74 genes that were significantly altered by ≥ 1.5-fold in HEK293 Null cells treated with CSE versus control media. 41 of these genes were downregulated and 33 up regulated (Figure 2A). Gene and conditions are displayed in order of how they cluster using Pearson’s centred rank correlation distance metric and Ward’s linkage rule on the entire 22837 probes gene sets. Using this analysis, there was a clear distinction between controls and smoke-treated samples. Gene symbols displayed more than once represent transcript variants of the same gene.

transcriptomics-cigarette-smoke-extract

Figure 2: Expression of genes that were significantly altered in HEK293 Null and TLR2/6 cells after treatment with cigarette smoke extract (CSE) for 8h. (A) Heat map representation of normalized signal intensity values for genes altered by ≥ 1.5-fold by 8h exposure to 10% CSE in HEK293 Null cells. Red denotes high expression and turquoise denotes low expression. Order of samples was dictated by hierarchical clustering and genes represented are listed in the accompanying table in the order that they appear. (B) Heat map representation of normalized signal intensity values for genes altered by ≥ 1.5-fold by 8h exposure to 10% CSE in HEK293 TLR2/6 cells. Red denotes high expression and turquoise denotes low expression. Order of samples was dictated by hierarchical clustering and genes represented are listed in the accompanying table in the order that they appear. Statistical significance (p<0.05) was calculated using moderated student’s t-test followed by Benjamini-Hochberg false discovery rate correction on GeneSpring GX11.0.2 software. Fold change represents a comparison between mean normalised signal intensity between treatment groups.

Effect of 10% CSE on the expression of genes in HEK293 TLR2/6 cells

Analysis of the data identified 67 genes that were significantly altered by ≥1.5-fold in HEK293 TLR2/6 cells treated with CSE versus controls. 43 of these genes were downregulated and 24 up regulated (Figure 2B). As was observed in the HEK293 Null cells, the HEK293 TLR2/6 cells also showed a clear distinction in their hierarchical clustering between controls and smoke-treated samples. Gene symbols displayed more than once represent transcript variants of the gene.

Genes altered by 10% CSE in a TLR2/6-dependent manner in HEK293 cells

A Venn diagram was constructed to assess which genes were differentially altered by 10% CSE in a TLR2/6-dependent manner (Figure 3). This consisted of a comparison of the mutual genes that were: 1) significantly altered by ≥1.5-fold in HEK293 TLR2/6 cells after treatment with 10% CSE (67 entities), 2) not altered in these cells under control conditions (18377 entities), and 3) not altered by 10% CSE in HEK293 Null cells (18556 entities). This resulted in the identification of 42 TLR2/6-dependent genes altered by cigarette smoke. 29 of these genes were downregulated and 13 upregulated (Figure 3). Ingenuity pathway analysis of these 42 TLR2/6-dependent CSEregulated genes highlighted a significant role in cancer, cardiovascular disease, developmental disorders, hereditary disorder and infectious diseases. These genes play a significant role in the following canonical pathways retinoate biosynthesis, aryl hydrocarbon receptor signaling, diacylglycerol biosynthesis, the visual cycle and phosphatidyl glycerol biosynthesis (Figure 4A). When looking at the top network functions, 34 of the 42 genes that were reliant on TLR2/6, NOTCH 1 and the micro RNA activator EIF2C2 and NURP1 seemed to form central hubs (Figure 4B). Network analysis of these genes showed a significant association with connective tissue disorder, inflammatory disease and the inflammatory response.

transcriptomics-Benjamini-Hochberg-false

Figure 3: Expression of genes that were significantly altered in a TLR2/6-dependent manner after treatment with 10% CSE for 8h. (A) Venn diagram constituting mutual genes differentially expressed between; 1) genes not altered by 10% CSE in HEK293 Null cells (green circle), 2) genes not altered by transfection with TLR2/6 (blue circle) and 3) genes significantly upregulated by ≥ 1.5-fold in HEK293 TLR2/6 cells after treatment with 10% CSE for 8 hours (red circle). (B) Heat map representation of normalized signal intensity values for genes altered by 10% CSE in a TLR2/6-dependent manner. Red denotes high expression and turquoise denotes low expression. Order of samples was dictated by hierarchical clustering, and genes represented are listed in the accompanying table in the order that they appear. Data represent n=3 and statistical significance (p<0.05) was calculated using moderated student’s t-test followed by Benjamini-Hochberg false discovery rate correction on GeneSpring GX11.0.2 software.

transcriptomics-Ingenuity-pathway-analysis

Figure 4: Ingenuity pathway analysis of genes that were significantly altered in a TLR2/6-dependent manner after treatment with 10% CSE for 8 h. (A) Shows the top five canonical pathways that the TLR2/6 dependent genes are significantly involved in. (B) Top associated network functions of the TLR2/6-dependent genes and their association with other key genes in these networks. Genes that were upregulated are in green and downregulated are in red. Genes added in due to pathway association are represented in white.

Comparison of genes altered by 10% CSE in HEK293 cells and THP-1 cells

To find TLR2/6 independent genes, a Venn diagram comparison was made between HEK293 Null cells (74 entities) and THP-1 cells (282 entities),after both sets of cells were exposed to 10% CSE for 8 hours, to identify mutual genes that were significantly altered by ≥ 1.5-fold. 26 genes were identified of which 14 were down regulated and 12 up regulated (Figure 5A). When a similar comparison was made between THP-1 cells (282 entities) and HEK293 TLR2/6 cells (67 entities) treated with 10% CSE, to identify TLR2/6 dependent genes, 18 genes were mutually expressed, 7 were downregulated and 11 upregulated (Figure 5B). When genes that were significantly modulated by 10% CSE in THP-1 cells (282 entities) were compared with genes significantly regulated by 10% CSE in a TLR2/6 dependent manner from HEK293 cells (42 entities), only 1 gene was up regulated and this was SRGN (Figure 5C).

transcriptomics-Venn-diagram-gene

Figure 5: Expression of genes that were mutually altered by 10% CSE in HEK293 cells and THP-1 monocytes after 8h. (A) Venn Venn diagram and gene list of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 Null cells (red circle) and THP-1 monocytes (green circle) after treatment with 10% CSE for 8 h. (B) Venn diagram and gene list of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 TLR2/6 cells (red circle) and THP-1 monocytes (green circle) after treatment with 10% CSE for 8 hours. (C) Venn diagram and gene list of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 cells in a TLR2/6 dependent manner (red circle) and THP-1 monocytes (green circle) after treatment with 10% CSE for 8 h. Statistical significance (p<0.05) was calculated using moderated student’s t-test followed by Benjamini-Hochberg false discovery rate correction on GeneSpring GX11.0.2 software.HEK293 Null and TLR2/6 cells have an n=3 and THP-1 cells have an n=6.

Comparison of genes altered by 10% CSE in HEK293 cells and PBMCs

A Venn diagram comparison was made between HEK293 Null cells (74 entities) and PBMCs (683 entities), after cells were exposed to 10% CSE for 8 h, to identify mutual genes that were significantly altered by ≥ 1.5-fold. 13 genes were identified; 1 was downregulated and 12 upregulated (Figure 6A). When a similar comparison was made between PBMCs (683 entities) and HEK293 TLR2/6 cells (67 entities) treated with 10% CSE, 12 genes were mutually expressed. Of these, 2 were downregulated and 10 upregulated (Figure 6B). When genes that were significantly modulated by 10% CSE in PBMCs (683 entities) were compared with genes significantly regulated by 10% CSE in a TLR2/6 dependent manner from HEK293 cells, 5 gene were identified (Figure 6C).

transcriptomics-Benjamini-Hochberg-false

Figure 6: Expression of genes that were mutually altered by CSE in HEK293 cells and PBMCs after 8h. (A) Venn diagram of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 Null cells (red circle) and PBMCs (green circle) after treatment with 10% CSE for 8 hours. (B) Venn diagram of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 TLR2/6 cells (red circle) and PBMCs (green circle) after treatment with 10% CSE for 8 hours. (C) Venn diagram of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 cells in a TLR2/6 dependent manner (red circle) and PBMCs (green circle) after treatment with 10% CSE for 8 h. Statistical significance (p<0.05) was calculated using moderated student’s t-test followed by Benjamini-Hochberg false discovery rate correction on GeneSpring GX11.0.2 software. HEK293 Null, HEK293 TLR2/6 and PBMCs cells all have an n=3.

Verification of gene modified by CSE in THP-1 cells and PBMCs

A number of genes that were modified by CSE were confirmed as being altered using qPCR. The genes chosen were taken from the list generated in Figure 6C, from genes altered by CSE in PBMCs in a TLR2/6 manner. PPM1F, SGRN and EIF2C2 were chosen as PPMIF is a known negative regulator of cellular stress and SGRN and EIF2C2 are regulators of the immune function. In THP-1 cells no increase in these genes was observed at 8h, however by 24h there was an observed fold change of 1.48 ± 0.22 (PPM1F), 1.96 ± 0.08 (SRGN) and 1.41 ± 0.06 (EIF2C2). By contrast these genes were significantly regulated in PBMCs at 8h, with a reduction in mRNA levels for PPM1F (0.62 ± 0.23) and an induction in SRGN (2.50 ± 0.27) and EIF2C2 (2.09 ± 0.69). Similar results were obtained in PBMCs treated with CSE for 24 h.

Comparison of genes altered by 10% CSE in HEK293 cells, THP-1 cells and PBMCs

A Venn diagram comparison was made between HEK293 Null cells (74 entities), THP-1 cells (282 entities) and PBMCs (683 entities), after all cells were treated with 10% CSE for 8 h, to identify mutual genes that were significantly altered by ≥1.5-fold. 10 genes were identified of which 1 was downregulated and 9 upregulated (Figure 7A). When a similar comparison was made between PBMCs (683 entities), THP- 1 cells (282 entities) modulated by 10% CSE and TLR2/6-dependent CSE-regulated genes (42 entities) only one gene was mutually expressed and this once again was an upregulation of SRGN (Figure 7B). Ingenuity pathway analysis of the 10 CSE-regulated genes in all three-cell types once again highlighted a significant role in very similar bio functions; these were cancer, cardiovascular disease, developmental disorder, hematological disease and hereditary disorder. Although the canonical pathways that these genes were involved with were slightly different: Nrf2 signaling, antioxidant actions of vitamin C, glutathione biosynthesis, hemedegradation and thioredoxin pathway (Figure 8A). When looking at the top network functions, all 10 of the genes analyzed, formed central hubs around HMOX1 and TNFα (Figure 8B). Network analysis of these genes showed a significant association with cell death and survival, cellular movement, hematological system development and function, free radical scavenging, small molecule biochemistry (Figure 9).

transcriptomics-separate-experiments

Figure 7: Conformation of the expression of TLR2/6-dependant that were altered by CSE in THP-1 monocyte and PBMCs after 8 and 24h. Cells were treated with cigarette smoke for 8 and 24h, mRNA extracted and the expression of PPMIF, SRGN and EIF2C2 were assessed by qPCR in (A) THP-1 monocytes and (B) PBMCs. THP-1 data represented is from n=3 separate experiments and PBMC data is from n=3 individual donors.

transcriptomics-Expression-genes

Figure 8: Expression of genes that were mutually altered by CSE in HEK293 cells, THP-1 monocytes and PBMCsafter 8 h. (A) Venn diagram of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 Null cells (red circle), THP-1 monocytes (blue circle) and PBMCs (green circle) after treatment with 10% CSE for 8 h. (B) Venn diagram of mutual genes that were significantly modulated by ≥ 1.5-fold in HEK293 cells in a TLR2/6 dependent manner (red circle), THP-1 monocytes (blue circle) and PBMCs (blue circle) after treatment with 10% CSE for 8 h. Statistical significance (p<0.05) was calculated using moderated student’s t-test followed by Benjamini-Hochberg false discovery rate correction on GeneSpring GX11.0.2 software.HEK293 Null, HEK293 TLR2/6 and PBMCs cells all have an n=3 and THP-1 cells have an n=6.

transcriptomics-Ingenuity-pathway-analysis

Figure 9: Ingenuity pathway analysis of common genes that were significantly altered by CSE in HEK293 cells, THP-1 monocytes and PBMCs after 8 h. (A) Shows the top seven canonical pathways that were
altered in all three cell types. (B) Top associated network functions of common genes in all three cell types and their association with other key genes in the networks. Genes that were upregulated are in green and downregulated are in red. Genes added in due to pathway association are represented in white.

Discussion

There is a clear and unquestionable link between smoking cigarettes and the pathogenesis of many diseases, in particular COPD and cardiovascular disease. We have previously shown that, in an acute exposure model using monocytes, cigarette smoke activates gene networks involved with oxidative stress, xenobiotic metabolism and inflammation [24]. We have also described the dependency of pro-inflammatory gene and protein production on TLR2 and its heterodimer complexes in monocytes [7]. In the current study, a transcriptomic approach was adopted to give an insight into TLR2/6’s function in cigarette smoke-induced cellular responses at the level of the transcriptome.

CSE significantly altered 42 genes by greater than 1.5-fold in a TLR2/6-dependent manner in HEK293 cells. Further analysis of these 42 genes identified pathways relating to retinoate biosynthesis, aryl hydrocarbon receptor signaling and CDP-diacylglycerol biosynthesis, with genes encoding for the proteins argonaute-2 (EIF2C2) and p53 (TP53) being central to the gene networks involved. Only 1 of these 42 genes was altered by cigarette smoke in THP-1 monocytes. When the 42 TLR2/6-dependent genes were compared with the genes changed by cigarette smoke in PBMCs, 5 were induced: ACN9, PCYOX1, SRGN, EIF2C2 and PPM1F.

The identification of genes altered by cigarette smoke in a TLR2/6- dependent manner is supported by our previous findings, which proposed a role for TLR2/6 in the acute cellular response to cigarette smoke. One of the genes, EIF2C2, is central to the above described gene networks and encodes for the protein argonaute-2. This protein has endonuclease activity required for RNA-mediated gene silencing, and is found in high concentrations in the processing bodies present in the cytoplasm [15,30]. Data directly linking EIF2C2 to cellular processes involved in inflammation caused by cigarette smoke is extremely limited [31]. However, a recent review highlighted the potential importance of smoke-induced alterations in non-coding RNA, for which argonaute activity is important due to its role in micro-RNA-associated gene silencing [32]. Although the direct link of this gene to acute inflammation caused by cigarette smoke is sparse, the over expression of EIF2C2 was found in precursor lesions of lung adenocarcinomas [33]. Our data implicates an increased expression of EIF2C2 in cigarette smoke-induced inflammation, and highlights a possible mechanism for its link to lung adenocarcinoma development.

Aryl hydrocarbon receptor (AhR) signaling was highlighted as an important pathway affected by TLR 2/6-dependent CSE-induced cellular responses and has been implicated in adenocarcinoma formation moted by cigarette smoke [34,35]. An overexpression of AhR is found lung cancer cells and this can promote cancer cell invasion by increasing the release of matrix metalloproteinase from these cells [36]. Elevated levels of AhR signaling also promotes DNA strand breaks [37], and stimulates clonogenicity of cancer cells [38].The oxidative stress sensitive transcription factor Nrf2, is critically involved in protecting the liver and gastrointestinal tract against disease by regulating a multifaceted cellular antioxidant defence. Its major target, NQO1, was shown to be insensitive to AhR expression level. However, the induction of AhR expression caused a concurrent increase in mRNA expression of xenobiotic-metabolizing gene, CYP1B1 [39] a gene that we identified as being TLR2/6-dependent and others have shown to be involved in cytogenetic damage [34].

CYP1B1 itself has been implicated in the carcinogenicity of cigarette smoke. It has been directly linked to NNAL (4-(methylnitrosamino)- 1-(3-pyridyl)-1-butanol and its glucuronides) serum levels, which is a prognostic factor for lung cancer risk. In fact, a single nucleotide polymorphism in CY1PB1 was the only polymorphism out of 11 common polymorphisms on four genes known to be sensitive to tobacco carcinogenesis [40].

Four more TLR2/6-dependent genes were also altered by CSE in PBMCs. As mentioned, EIF2C2 may have a link to smoke-related disease pathogenesis, while PCYOX1 is associated with the release of bound cysteine, which is essential in antioxidant activity. PPM1F is associated with cell recovery, and it is of interest to note that breast cancer cell migration is repressed by microRNA-200c targeted repression of PPM1F [41]. This may be a direct mechanism behind the positive association of cigarette smoke and breast cancer, and may also be a novel mechanism by which smoke-induced PPM1F could lead to cancer progression [42].

When we compared TLR2/6-dependent CSE-regulated genes with CSE-regulated genes in both THP-1 monocytes and PBMCs, there was just one gene found as a match. SRGN, the gene encoding serglycin was upregulated by CSE in all three cell types. Serglycin is essential for the maturation of mast cell secretory granules, and is also the major proteoglycan secreted by peritoneal macrophages [43] and control TNF-α release from peritoneal macrophages [43]. Thus, in the setting of cigarette smoking, up regulated SRGN expression in PBMCs and monocytes may provide an explanation for increased cytokine release, which may, if the observation translates to people, impact on susceptibility to infection in smokers.

In the current, using HEK293 cells, we identify a list of genes that are altered by CSE in a TLR2/6-dependent and –independent manner. We have previously used HEK293 cells to show that oxidants can be perceived through a TLR2-dependent mechanism, in these studies HEK293 cells had a similar response to oxidants as THP-1 cells and monocytes [44]. To compensate for the HEK293 cells being embryonic cell line, we have cross-referenced these genes with CSE-induced genes present in monocytes. The potential importance of our findings have been given further validity, since a recent review has highlighted the importance of the proteoglycan serglycin on immune cell function [45]. The fact that this gene may be controlled through a TLR2-dependent mechanism is potentially important to immunologists.

In this study we have identified both TLR2/6 dependent and TLRindependent genes altered by CSE. This study provides a unifying answer for others in the literature, where either TLR-dependent or -independent smoke induced inflammation has been reported. The TLR2/6-dependent genes identified seem to be involved in cancer cells development and therefore warrant further investigation to facilitate the development of novel therapies for smoking-related diseases.

Acknowledgement

MJPC and JAM are funded by the Wellcome Trust Foundation.

References

  1. Hecht SS (1999) Tobacco smoke carcinogens and lung cancer. J Natl Cancer Inst 91: 1194-1210.
  2. Gonçalves RB, Coletta RD, Silvério KG, Benevides L, Casati MZ, et al. (2011) Impact of smoking on inflammation: overview of molecular mechanisms. Inflamm Res 60: 409-424.
  3. Barreiro E, Peinado VI, Galdiz JB, Ferrer E, Marin-Corral J, et al. (2010) Cigarette smoke-induced oxidative stress: A role in chronic obstructive pulmonary disease skeletal muscle dysfunction. Am J Respir Crit Care Med 182: 477-488.
  4. Walters MJ, Paul-Clark MJ, McMaster SK, Ito K, Adcock IM, et al. (2005) Cigarette smoke activates human monocytes by an oxidant-AP-1 signaling pathway: implications for steroid resistance. Mol Pharmacol 68: 1343-1353.
  5. Sopori M (2002) Effects of cigarette smoke on the immune system. Nat Rev Immunol 2: 372-377.
  6. Paul-Clark MJ, Sorrentino R, Bailey LK, Sriskandan S, Mitchell JA (2008) Gram-positive and Gram-negative bacteria synergize with oxidants to release CXCL8 from innate immune cells. Mol Med 14: 238-246.
  7. Walters MJ, Paul-Clark MJ, McMaster SK, Ito K, Adcock IM, et al. (2005) Cigarette smoke activates human monocytes by an oxidant-AP-1 signaling pathway: implications for steroid resistance. Mol Pharmacol 68: 1343-1353.
  8. Albensi BC, Mattson MP (2000) Evidence for the involvement of TNF and NF-kappaB in hippocampal synaptic plasticity. Synapse 35: 151-159.
  9. Baeuerle PA, Baichwal VR (1997) NF-kappa B as a frequent target for immunosuppressive and anti-inflammatory molecules. Adv Immunol 65: 111-137.
  10. Fitzgerald DC, Meade KG, Mc Evoy AN, Lillis L, Murphy EP et al. (2007) Tumour necrosis factor-alpha (TNF-alpha) increases nuclear factor kappaB (NFkappaB) activity inandinterleukin8 (IL8) release from bovine mammary epithelial cells. Vet Immunol Immunopathol. 116:59-68.
  11. Rahman I (2003) Oxidative stress, chromatin remodeling and gene transcription in inflammation and chronic lung diseases. J Biochem Mol Biol 36: 95-109.
  12. Paul-Clark MJ, George PM, Gatheral T, Parzych K, Wright WR, et al. (2012) Pharmacology and therapeutic potential of pattern recognition receptors. Pharmacol Ther 135: 200-215.
  13. Doz E, Noulin N, Boichot E, Guénon I, Fick L, et al. (2008) Cigarette smoke-induced pulmonary inflammation is TLR4/MyD88 and IL-1R1/MyD88 signaling dependent. J Immunol 180: 1169-1178.
  14. Paul-Clark MJ, McMaster SK, Sorrentino R, Sriskandan S, Bailey LK, et al. (2009) Toll-like receptor 2 is essential for the sensing of oxidants during inflammation. Am J Respir Crit Care Med 179: 299-306.
  15. Castro SM, Chakraborty K, Guerrero-Plata A (2011) Cigarette smoke suppresses TLR-7 stimulation in response to virus infection in plasmacytoid dendritic cells. Toxicol In Vitro 25: 1106-1113.
  16. KoaraiA, YanagisawaS, SugiuraH, IchikawaT, AkamatsuK, et al. (2012) Cigarette smoke augments the expression and responses of toll-likereceptor3inhumanmacrophages. Respirology 17: 1018-1025.
  17. Mortaz E, Adcock IM, ItoK, Kraneveld AD, Nijkamp FP et al. (2010) Cigarette smoke induces CXCL8 production by human neutrophils via activation of TLR9 receptor. Eur Respir J 36:1143-1154.
  18. Halliwell B (1987) Oxidants and human disease: some new concepts. FASEB J 1: 358-364.
  19. HeguyA, O'ConnorTP, LuettichK, Worgall S, Cieciuch A, et al. (2006) Gene expression profiling of human a lveolarmacrophages of phenotypically normal smokers and non-smokers reveals a previously un recognized sub set of genes modulated by cigarette smoking. J Mol Med. 84: 318-328.
  20. Hutcheson R, Rocic P (2012) The metabolic syndrome, oxidative stress, environment, and cardiovascular disease: the great exploration. Exp Diabetes Res 2012: 271028.
  21. Gebel S, Gerstmayer B, Bosio A, Haussmann HJ, Van Miert E, et al. (2004) Gene expression profiling in respiratory tissues from rats exposed to mainstream cigarette smoke. Carcinogenesis 25: 169-178.
  22. Gümüş ZH, Du B, Kacker A, Boyle JO, Bocker JM, et al. (2008) Effects of tobacco smoke on gene expression and cellular pathways in a cellular model of oral leukoplakia. Cancer Prev Res (Phila) 1: 100-111.
  23. Kent L, Smyth L, Clayton C, Scott L, Cook T, et al. (2008) Cigarette smoke extract induced cytokine and chemokine gene expression changes in COPD macrophages. Cytokine 42: 205-216.
  24. Wright WR, Parzych K, Crawford D, Mein C, Mitchell JA, et al. (2012) Inflammatory transcriptome profiling of human monocytes exposed acutely to cigarette smoke. PLoS One 7: e30120.
  25. Walters MJ, Paul-Clark MJ, McMaster SK, Ito K, Adcock IM, et al. (2005) Cigarette smoke activates human monocytes by an oxidant-AP-1 signaling pathway: implications for steroid resistance. Mol Pharmacol 68: 1343-1353.
  26. Rahman I, Bel A, Mulier B, Donaldson K, MacNee W (1998) Differential regulation of glutathione by oxidants and dexamethasone in alveolar epithelial cells. Am J Physiol 275: L80-86.
  27. Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185-193.
  28. Schmid R, Baum P, Ittrich C, Fundel-Clemens K, Huber W, et al. (2010) Comparison of normalization methods for Illumina BeadChip HumanHT-12 v3. BMC Genomics 11: 349.
  29. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B 57: 289-300.
  30. Ghildiyal M, Zamore PD (2009) Small silencing RNAs: an expanding universe. Nat Rev Genet 10: 94-108.
  31. Xi S, Xu H, Shan J, Tao Y, Hong JA, et al. (2013) Cigarette smoke mediates epigenetic repression of miR-487b during pulmonary carcinogenesis. J Clin Invest 123: 1241-1261.
  32. Maccani MA, Knopik VS (2012) Cigarette smoke exposure-associated alterations to non-coding RNA. Front Genet 3: 53.
  33. Chiosea S, Jelezcova E, Chandran U, Luo J, Mantha G, et al. (2007) Overexpression of Dicer in precursor lesions of lung adenocarcinoma. Cancer Res 67: 2345-2350.
  34. Dertinger SD, Nazarenko DA, Silverstone AE, Gasiewicz TA (2001) Aryl hydrocarbon receptor signaling plays a significant role in mediating benzo[a]pyrene- and cigarette smoke condensate-induced cytogenetic damage in vivo. Carcinogenesis 22: 171-177.
  35. Lin P, Chang H, Tsai WT, Wu MH, Liao YS, et al. (2003) Overexpression of aryl hydrocarbon receptor in human lung carcinomas. Toxicol Pathol 31: 22-30.
  36. Ishida M, Mikami S, Kikuchi E, Kosaka T, Miyajima A, et al. (2010) Activation of the aryl hydrocarbon receptor pathway enhances cancer cell invasion by upregulating the MMP expression and is associated with poor prognosis in upper urinary tract urothelial cancer. Carcinogenesis 31: 287-295.
  37. Park JH, Mangal D, Frey AJ, Harvey RG, Blair IA, et al. (2009) Aryl hydrocarbon receptor facilitates DNA strand breaks and 8-oxo-2'-deoxyguanosine formation by the aldo-keto reductase product benzo[a]pyrene-7,8-dione. J Biol Chem 284: 29725-29734.
  38. Opitz CA, Litzenburger UM, Sahm F, Ott M, Tritschler I, et al. (2011) An endogenous tumour-promoting ligand of the human aryl hydrocarbon receptor. Nature 478: 197-203.
  39. Yang X, Solomon S, Fraser LR, Trombino AF, Liu D, et al. (2008) Constitutive regulation of CYP1B1 by the aryl hydrocarbon receptor (AhR) in pre-malignant and malignant mammary tissue. J Cell Biochem 104: 402-417.
  40. Church TR, Haznadar M, Geisser MS, Anderson KE, Caporaso NE, et al. (2010) Interaction of CYP1B1, cigarette-smoke carcinogen metabolism, and lung cancer risk. Int J Mol Epidemiol Genet 1: 295-309.
  41. Jurmeister S, Baumann M, Balwierz A, Keklikoglou I, Ward A, et al. (2012) MicroRNA-200c represses migration and invasion of breast cancer cells by targeting actin-regulatory proteins FHOD1 and PPM1F. Mol Cell Biol 32: 633-651.
  42. Terry PD, Rohan TE (2002) Cigarette smoking and the risk of breast cancer in women: are views of the literature. Cancer Epidemiol Biomarkers Prev. 11:953-971.
  43. Zernichow L, Abrink M, Hallgren J, Grujic M, Pejler G et al. (2006) Serglycin is the major secreted proteoglycan in macrophages and has a role in the regulation of macrophage tumor necrosis factor-alpha secretionin response to lipopolysaccharide. TheJournal of biological chemistry 281: 26792-26801.
  44. Paul-Clark MJ, McMaster SK, Sorrentino R, Sriskandan S, Bailey LK, et al. (2009) Toll-like receptor 2 is essential for the sensing of oxidants during inflammation. Am J Respir Crit Care Med 179: 299-306.
  45. Kolset SO, Pejler G (2011) Serglycin: a structural and functional chameleon with wide impact on immune cells. J Immunol 187: 4927-4933.
Citation: Wright WR, Mazi SI, Parzych K, Mataragka S, Mein C, et al. (2013) The Role of TLR2 in Cigarette Smoke-Induced Gene Induction. Transcriptomics 1:102.

Copyright: © 2013 Wright WR, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Top