Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genomic and transcriptomic comparison of allergen and silver nanoparticle-induced mast cell degranulation reveals novel non-immunoglobulin E mediated mechanisms

  • Monica Johnson,

    Roles Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Nasser Alsaleh,

    Roles Investigation

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Ryan P. Mendoza,

    Roles Investigation

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Indushekhar Persaud,

    Roles Investigation

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Alison K. Bauer,

    Roles Conceptualization, Funding acquisition, Writing – review & editing

    Affiliation Department of Environmental and Occupational Health, Colorado School of Public Health, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Laura Saba,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯
  • Jared M. Brown

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Writing – original draft, Writing – review & editing

    jared.brown@ucdenver.edu

    Affiliation Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of Colorado, Aurora, CO, United States of America

    ⨯

Abstract

Mast cells represent a crucial cell type in host defense; however, maladaptive responses are contributing factors in the pathogenesis of allergic diseases. Previous work in our laboratory has shown that exposure to silver nanoparticles (AgNPs) results in mast cell degranulation via a non-immunoglobulin E (IgE) mechanism. In this study, we utilized a systems biology approach to identify novel genetic factors playing a role in AgNP-induced mast cell degranulation compared to the classical activation by antigen-mediated FcεRI crosslinking. Mast cell degranulation was assessed in bone marrow-derived mast cells isolated from 23 strains of mice following exposure to AgNPs or FcεRI crosslinking with dinitrophenyl (DNP). Utilizing strain-dependent mast cell degranulation, an association mapping study identified 3 chromosomal regions that were significantly associated with mast cell degranulation by AgNP and one non-overlapping region associated with DNP-mediated degranulation. Two of the AgNP-associated regions correspond to genes previously reported to be associated with allergic disorders (Trac2 on chromosome 1 and Traf6 on chromosome 2) and an uncharacterized gene identified on chromosome 1 (Fam126b). In conjunction, RNA-sequencing performed on mast cells from the high and low responder strains revealed 3754 and 34 differentially expressed genes that were unique to DNP and AgNP exposures, respectively. Select candidate genes include Ptger4, a gene encoding a G-protein coupled receptor in addition to a multifunctional adaptor protein, Txnip, that may be driving mast cell degranulation by AgNP. Taken together, we identified novel genes that have not been previously shown to play a role in nanoparticle-mediated mast cell activation. With further functional evaluation in the future, these genes may be potential therapeutic targets in the treatment of non-IgE mediated mast cell-linked disorders.

Citation: Johnson M, Alsaleh N, Mendoza RP, Persaud I, Bauer AK, Saba L, et al. (2018) Genomic and transcriptomic comparison of allergen and silver nanoparticle-induced mast cell degranulation reveals novel non-immunoglobulin E mediated mechanisms. PLoS ONE 13(3): e0193499. https://doi.org/10.1371/journal.pone.0193499

Editor: Amitava Mukherjee, VIT University, INDIA

Received: October 17, 2017; Accepted: February 12, 2018; Published: March 22, 2018

Copyright: © 2018 Johnson 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.

Data Availability: All data will be available in the manuscript and supporting information.

Funding: This study was supported by National Institute of Environmental Health Science (NIEHS) grant R01 ES019311 (JMB).

Competing interests: The authors have declared that no competing interests exist.

Introduction

In the last decade, development of new nanotechnologies has contributed to the advancement of many fields. The unique properties of engineered nanomaterials (ENMs) gives them novel electrical, magnetic, mechanical, thermal, or imaging properties that are highly desirable for applications in commercial, medical, and environmental sectors [1]. Although ENMs provide many advantages, increased manufacturing and lack of safety testing raises health concerns due to increased human and environmental exposures [2, 3]. Specifically, adverse immune outcomes by ENM exposures have been reported [4].

Mast cells, a crucial cell type involved in allergic inflammation, have been shown to be activated in response to ENM exposures [5]. Specifically, exposure to a variety of physicochemically distinct ENMs including carbon-based and metal-based nanoparticles induce mast cell degranulation in vivo and in vitro [6–9]. Recently, we demonstrated a robust mast cell degranulation response following silver nanoparticle (AgNP) exposure, which was dependent on physicochemical properties such a size, shape and surface chemistry [10, 11]. AgNPs are one of the most abundantly manufactured ENM due to their antimicrobial/fungal properties and are currently utilized in more than 400 consumer products including wound dressings, food storage containers, and athletic apparel [1]. Understanding the role of AgNP exposure in immunomodulation is critical for evaluating ENM safety for consumer health, which is especially concerning for the population susceptible to mast cell mediated-diseases such as asthma, allergic dermatitis and hyperactive mast cell activation disorders [12–14].

To improve our understanding about the role of AgNPs in mast cell activation, it is crucial to determine the underlying mechanism, which is poorly defined. Conventional mast cell activation is a concerted event initiated by immunoglobulin E (IgE)-dependent cross-linking of high affinity IgE-bound receptors (i.e. FcεRI) on the surface of mast cells [15, 16]. FcεRI crosslinking by IgE and an allergen results in downstream signaling events that accumulate in the release of preformed as well as newly synthesized inflammatory mediators (i.e. histamine, proteases, leukotrienes, IL-4, IL-9, TNF-α, etc.) [13, 15]. Importantly, previous studies have shown that AgNP initiates mast cell activation via a non-IgE mediated pathway, independent of IgE sensitization [17]. Mast cell activation by non-IgE mediated mechanisms is not uncommon, as other factors have been identified to activate mast cells via alternate receptors include physical factors (pH, temperature), toxins, and endogenous signaling components [18–27]. However, the degranulation response by AgNPs has yet to be fully elucidated.

Little is known regarding the mechanism resulting in non-IgE mediated mast cell activation, however prior evidence supports a strong genetic component to allergic diseases [28, 29]. This remains particularly true in the helper T cell type 2 (TH2) responses (i.e. asthma, allergic inflammatory diseases) in which a number of studies have used linkage and association analysis to identify genes responsible for disease pathogenesis [30–32]. However, even less is understood about the direct genetic factors playing a role in nanoparticle-induced inflammation. Several studies utilized various strains of mice to evaluate a genetic component to this response. For example, susceptibility of quantum dot-induced lung inflammation, which resulted in neutrophil infiltration and increases in cytokines, was strain-dependent and heritable [33]. Jones et al. observed a strain-dependent effect on nanoparticle clearance, including a slower rate of clearance in strains that are prone to T cell helper type II (TH1; i.e. C57BL/6, B10D2) compared to TH2-prone (i.e. BALB/c, DBA/2) mice [34]. Overall, previous findings provide evidence that a complex set of genes regulates allergic diseases, with the potential for discovering mechanisms regulating nanoparticle-induced mast cell degranulation.

Therefore, to achieve this broader understanding of AgNP-induced mast cell activation, the current study utilized a modified hybrid mouse diversity panel consisting of 23 strains of recombinant and inbred strains of mice to determine strain-dependent susceptibility to mast cell degranulation by both AgNP and antigen-mediated FcεRI crosslinking by IgE. Using the phenotype data across multiple strains, genetic factors were elucidated using a systems biology approach. First, advanced association mapping methods were utilized to identify genetic loci associated with mast cell degranulation following AgNP exposure as well as identify novel genes mediating the classical IgE-mediated response. Lastly, transcriptomic analysis (via RNA-seq) measuring differential gene expression patterns was performed on high and low responder strains to identify novel pathways and gene targets involved in non-IgE mediated mast cell responses following AgNP exposures. To our knowledge, no gene expression data exists for murine mast cells, especially in multiple strains. Using this novel data, the current study uncovers pathways and receptors potentially playing a role in non-IgE mediated mast cell activation.

Materials and methods

Characterization of silver nanoparticles (AgNP)

20 nm spherical AgNPs suspended in citrate buffer were purchased from NanoComposix (San Diego, CA) at a concentration of 1 mg/ml. Primary size determined by transmission electron microscopy (TEM, Hitach H7600) & hydrodynamic diameter and zeta potential in solution were determined using ZetaSizer Nano dynamic light scattering (DLS, Malvern). All measurements were performed with 3 technical replicates at concentrations of 25 μg/ml, a dose chosen based off previously reported data [11]. It is important to note that this dose is probably higher than the average human acute exposure and is not designed to equate to human risk [35]. However, several studies have evaluated AgNP concentrations in consumer products and have demonstrated variability ranging from 1.3 μg/ml to 270,000 μg Ag /g product [35]. Dissolution of Ag+ was determined by ICP-MS to be 0.14 μg/ml in a solution of 25 μg/ml AgNPs in HEPES buffer for 1 hr. Importantly, we previously reported dissolution of Ag+ and from these particles and the ionic Ag did not contribute to mast cell degranulation (10).

Mice

Twenty-three strains of mice were utilized in the study. An inbred panel of 18 strains of mice were chosen based on their respective genetic diversity using the ‘Hybrid Mouse Diversity Panel’ as a reference (C57BL/6J, BALB/cJ, FVB/NJ, C3H/HeJ, RIIIS/J, SWR/J, DBA/2J, 129X1/SvJ, NOD/ShiLtJ, NZW/LacJ, PL/J, CBA/J, CAST/EiJ, BALB/cByJ, LP/J, SJL/J, PWD/PhJ, I/LnJ) [36]. In addition, five recombinant inbred (RI) strains derived from parental lines C57BL/6J and DBA/2J (BXD50/RwwJ, BXD61/RwwJ, BXD73/RwwJ, BXD90/RwwJ, BXD100/RwwJ) were chosen based on the significant mast cell degranulation phenotype observed between the parental strains. Femurs from all strains were obtained from Jackson Laboratories (Bar Harbor, ME).

Cell culture

Femoral bones isolated from 6–8 week old male mice (23 different strains as mentioned above) were used for progenitor cell isolation. Bone marrow progenitor cells (pooled from 2–3 mice) were cultured in supplemented media containing 300 ng/ml purified recombinant mouse interleukin-3 (PeproTech, Rocky Hill, NJ) for 4–6 weeks until fully differentiated into bone marrow derived mast cells (BMMCs; 37°C incubator, 5% CO2). The RPMI media also contained 10% FBS, 100 μg/ml streptomycin, 100 U/ml penicillin, 100 μg/ml primocin (Invivogen, San Diego, CA), 25 mM N-2-hydroxyethylpiperazine-N-2-ethane sulfonic acid (HEPES), 1.0 mM sodium pyruvate, nonessential amino acids (BioSource International, Camarillo, CA), and 0.00035% 2-mercapto-ethanol. All animal procedures were approved and conducted in accordance with the National Institutes of Health guidelines and approved by the University of Colorado Denver Institutional Animal Care and Use Committee. Cells were maintained in non-adherent flasks until treatment. To confirm mast cell responses in a human model, we utilized LUVA cells, a human mast cell line, which was obtained from Kerafast (Boston, MA). Cells were cultured at 37°C and 5% CO2 in StemPro®-34 SFM Complete Medium (Thermo Fisher Scientific Waltham, MA) containing 200 mM L-glutamine, 100 U/ml penicillin, 100 μg/ml streptomycin, 100 μg/ ml PrimocinTM. Cells were plated at 75,000 cells per well in HEPES, exposed to AgNPs at 25 or 50 μg/ml for 1 hour, and then cell degranulation was assessed as described below.

Flow cytometric analysis

To assess surface receptor expression, BMMCs were incubated with the following antibodies (Thermo Fisher Scientific) for 30 min at room temperature: PE-conjugated anti-FcεRI (MAR-1) and FITC-conjugated anti-CD117 (c-Kit). Flow cytometric analysis was performed on 1 x 105 cells performed in triplicate from 3 individual batches of BMMCs using an Accuri C6 Flow Cytometer (BD Biosciences, San Jose, CA) with FCS Express 4 software (De Novo software).

Mast cell degranulation assay

BMMCs cultured for 4–6 weeks were plated at confluency at 5 x 104 cells per well in non-adherent plates in biological buffer (HEPES, pH 7.4). β-hexosaminidase (β-hex) enzyme release was analyzed in BMMC and LUVA human mast cells following exposure of dinitrophenyl (DNP, 100 ng/ml, 30 min), AgNP (25 μg/ml, 1 h) with or without pretreatment (30 min) of the selective competitive antagonist of the EP4 receptor (10 μM, GW627368X, Cayman Chemicals), prostaglandin E (PGE2, 10 μM, Sigma) or non-selective inhibitor of cyclooxygenase (COX)- 1 and -2 (indomethacin, 10 μM, Sigma) as previously described [10, 37]. For DNP-mediated degranulation, BMMCs were sensitized overnight with IgE anti-DNP (100 ng/ml, Sigma- Aldrich, St. Louis, MO), followed by treatment with DNP-HSA for 30 min. After 30 min incubation of DNP-HSA or 1h incubation of AgNPs (25 μg/ml), p-nitrophenyl-N-acetyl-ß-D-glucopyranoside (PNAG; Sigma–Aldrich, St. Louis, MO), a chromogenic substrate of ß-hexosaminidase, was added to cell supernatants and lysates and incubated for 90 min at 37°C. The reaction was stopped with glycine (0.4M) and optical density was read at 405 nm using a SynergyTM HT Multi-Mode Microplate Reader (BioTek Instruments Inc, Winooski, VT). ß-hex release was calculated as percent total cell release after subtracting background release from untreated cells and normalized to non-treated control cells. Basal levels of ß-hex were similar across all strains examined. All experiments were performed in triplicate from 3 individual batches of mature mast cells grown in IL-3 supplemented media.

Statistical analyses

Graphs for mast cell maturation and degranulation were designed and analyzed using GraphPad Prism 7 software (San Diego, CA). Mast cell degranulation measurements are represented as mean technical replicates of all biological replicates (n = 2–3) ± standard error of the mean (SEM). Statistical comparisons were performed by Student’s t-test for two groups and analysis of variance (ANOVA) for more than two groups, with differences between groups assessed using a Bonferroni post hoc test.

The heritability of mast cell degranulation by AgNP and DNP was estimated from the variance components of the linear mixed model that accounts for population structure. The pylmmGWAS multiPhHeri.py program distributed with the MultiTrans software was used to estimate these variance components [38].

Inbred and RI genotyping

For the association mapping, single nucleotide polymorphism (SNP) information on the 23 strains of mice was downloaded from the Mouse Phenome Database (Jackson laboratories, http://phenome.jax.org, [39, 40]. Of the 415,431 fully informative SNPs retrieved after removing heterozygous calls or missing genotyping data, 198,345 informative SNPs with a minor allele frequency greater than 10% and positions within the NCBI-build-38 version of the mouse genome were used for association mapping.

Association mapping

Genome-wide association mapping was performed using phenotype data from mast cell degranulation assays following either AgNP or DNP treatment as described above. Association mapping was executed using the linear efficient mixed-model association (emma) package (version 1.1.2) in R statistical software (version 3.3.0), taking into account population structure among strains and variance between biological replicates [41]. One strain, BALB/cByJ was not included in the mapping study due to its high degree of genetic similarity to BALB/cJ. The genome-wide significance threshold was calculated using a Bonferroni adjustment for the number of unique strain distribution patterns within the SNP data. Upper and lower bounds of an associated region were defined as the location of the most proximal and most distal genotyped SNP nominally association with the genotype (unadjusted p-value<0.01) within 10 Mb of a SNP that reached genome-wide statistical significance (Bonferroni adjusted p-value<0.05). P-values from association mapping were–log10 transformed for graphic visualization.

The UCSC Genome Browser (https://genome.ucsc.edu) was used to identify genes residing in the associated regions [42]. The NHGRI-EBI GWAS catalog (https://www.ebi.ac.uk/gwas) was utilized to find any candidate genes that had been previously identified in human GWAS for allergic immune response traits (i.e. asthma, allergies, allergic dermatitis, general immunological responses) [43].

RNA-sequencing

High and low responder strains were selected based on their differences in mast cell degranulation response following AgNP (i.e. C57BL/6, LP/J) or DNP (i.e. RIIIS/J, C57BL/6J) treatment. For each strain, two biological replicates of BMMCs were produced. Within a biological replicate, cells were separated into 2 groups (3 groups for C57BL/6J), one group remained untreated and the other group was treated with either AgNP (25 μg/ml, 1h) or DNP (100 ng/ml, 30 min) as described above. Following treatment, total RNA from BMMCs was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and purification using a Direct-zol RNA miniprep kit (Zymo Research) following the manufacturers protocol. RNA quality and quantity was estimated by QC (Agilent 2100 Bioanalyzer, Santa Clara, CA). The detailed mRNA library preparation was performed following the manufacturers protocol using TruSeq Stranded mRNA Sample Preparation kit (Illumina, San Diego, California). The library products were sequenced on the Illumina HiSeqTM4000 (Illumina; 50 cycles, single-read sequencing, 7 samples per lane). Raw sequencing reads were trimmed for adapters and for quality using the default settings of trim_galore (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). Trimmed reads were aligned to the Ensembl mouse transcriptome (version GRCm38.84) using the RSEM (RNA-Seq by Expectation-Maximization) package (v1.2.31) [44]. Gene-level expected counts were used for differential expression analysis with the DESeq2 R package (version 1.12.2). Prior to differential expression analysis, genes were removed if their average expected read count across the 14 samples was less than 5. Two separate models were used to examine the effects of AgNP and DNP. For each treatment type, we examined the effect of the specific treatment compared to control animals by initially producing a single p-value that simultaneously tested for either common treatment effect across strains or a strain-specific treatment effect within each gene in a model that also included an effect for the animal from which the mast cells were derived. We also specifically tested for difference in treatment effects between inbred mouse strains. The differences in expression levels were estimated from regularized log transformed values of the original count data (rlog function in the DESeq2 package). A false discovery rate (FDR; [45]) was used to adjust for multiple comparisons across genes. Statistical overrepresentation tests of gene ontology (GO) categories enriched in the differentially expressed gene sets were analyzed using Panther Classification System (http://.pantherdb.org/, [46]). The binomial test compares an input (test) gene list to a reference gene list and determines whether there is statistical overrepresentation or underrepresentation of categories in the input gene set. ‘Expected’ values are the number of genes that would be expected in a gene list given the number of genes mapped to a specific GO category and ‘actual’ values are the input genes. Significance was set at p-value cut off of 0.05 with a fold enrichment cutoff of 5.

Quantitative real-time polymerase chain reaction (qPCR)

Following 1 h treatment of AgNP (25 μg/ml), total RNA from BMMCs was collected as described above. RNA was reverse-transcribed using the Quantitect reverse transcription kit (Bio-Rad, Hercules, CA). Quantitative real-time polymerase chain reaction (qPCR) was performed using a Quantitech SYBR Green PCR kit (Bio-Rad) and the StepOne Plus real-time PCR system (Applied Biosciences) to obtain cycle threshold (Ct) values for target and internal reference cDNA levels. Gene specific primers for Ptger4, Txnip, Trib2, Fam43a, and Gadd45g were obtained from Invitrogen. Gene expression levels were calculated relative to Gapdh (glyceraldehyde 2-phosphate dehydrogenease) as an internal endogenous control according to the ΔCT method where ΔCt is defined as Ct target- Ct internal reference. Values are represented as the mean ± SEM from 3 individual batches of mature mast cells grown in IL-3 supplemented media.

Western blotting

Protein expression analysis by western blotting was performed similar to previously described methods [17]. Briefly, BMMCs isolated from high responder (C57BL/6) and low responder (LP/J) strains of mice were exposed to AgNP at 25 μg/ml for 1 hour. Following isolation of cell supernatant by lysis, protein concentration was determined using Pierce BCA protein assay kit (Thermo Scientific #23225). 20 μg of total protein was loaded and separated using a 12% Tris-Glycine gel and transferred onto nitrocellulose membrane. Primary goat anti-mouse antibodies for β-actin and TXNIP (Cell Signaling Technologies), were incubated overnight, followed by incubation of horseradish-peroxidase conjugated secondary antibodies. Membranes were developed using Pierce Chemiluminescent Substrate (Thermo/Fisher Scientific, Waltham, MA). Relative densitometry was evaluated using Image J software (NIH, Bethesda, MD, USA). Immunoblots are representative images from three independent experiments.

Results

Characterization of silver nanoparticles (AgNP)

Primary size of silver nanoparticles (AgNP) at 25 μg/ml as determined by transmission electron microscopy was 20 ± 4.2 nm. In solution, the hydrodynamic diameter measured by dynamic light scattering measured 16 ± 1.1 nm. The zeta potential was -38 mV, indicating high stability in solution.

Genetic effects of mast cell degranulation caused by silver nanoparticles (AgNP) or FcεRI crosslinking by dinitrophenyl (DNP)

To characterize the genetic contribution to mast cell responses, primary bone marrow derived mast cells (BMMCs) were isolated from 23 strains of mice and evaluated for degranulation by measuring the release of ß -hexosaminidase following FcεRI crosslinking by DNP or AgNP (25 μg/ml) treatment (Fig 1). Prior to degranulation experiments, BMMC maturation was confirmed by expression of CD117 (c-Kit) and the high affinity immunoglobulin E (IgE) receptor (FcεRI)). All strains except I/LnJ, PWD/PhJ, DBA/2J, and CAST/EiJ fully expressed both receptors at 5 weeks following culture in IL-3 (S1 Fig). Strain differences in mast cell degranulation following DNP stimulation revealed a strong genetic influence (heritability = 96%; Fig 1B). For example, the highest responder strain, RIIIS/J (74%) had a 3.6-fold increase in mast cell degranulation compared to C57BL/6J (16%), which is a commonly used strain in mast cell studies. Four non-responsive strains: I/LnJ, PWD/PhJ, DBA/2J and CAST/EiJ had nominal responses below 2% degranulation as expected due to low FcεRI and c-Kit receptor expression.

thumbnail
Download:
Fig 1. Mast cell degranulation by AgNP and DNP across 23 strains of mice.

Bone marrow-derived mast cells (BMMCs) were isolated from 23 inbred and recombinant inbred (BxD) strains of mice. Degranulation was assessed by measuring ß-hexosaminidase release following (A) AgNP treatment at 25 μg/ml for 1 h in un-sensitized BMMCs, or (B) DNP treatment at 100 ng/ml for 0.5 h after BMMC overnight sensitization with IgE anti-DNP. Arrows represent low and high responder strains for each treatment used in RNA-seq (AgNP: green-LP/J, red-C57BL/6J; DNP: red-C57BL/6J, blue-RIIIS/J). Values are expressed as mean ± SEM (n = 3/strain) normalized to strain-specific non-treated control values.

https://doi.org/10.1371/journal.pone.0193499.g001

Interestingly, mast cell degranulation by AgNP also exhibited a strong genetic influence (heritability = 83%) but followed a strikingly different strain distribution pattern compared to DNP (Fig 1B). C3H/HeJ (36%) and C57BL/6J (35%) were the highest responder strains with an 8-fold increase in response compared to the lowest responder strain, LP/J (4%). The four strains previously identified with low surface receptor expression and subsequently non-responsive to DNP-induced mast cell degranulation (I/LnJ, PWD/PhJ, DBA/2J, CAST/EiJ), displayed very different responses following AgNP treatment. Specifically, AgNP treatment caused significant degranulation in the two strains PWD/PhJ (22%) and DBA/2J (15%). These results further demonstrate that AgNPs are not activating mast cells via an FcεRI mechanism but likely a novel receptor or cell membrane-driven mechanism. Lastly, we confirmed the mast degranulation response using LUVA human mast cells. As shown in S2 Fig, a dose dependent increase in mast cell degranulation was observed following AgNP treatment, consistent with murine mast cell responses.

Based on the unique strain-dependent differences in mast cell degranulation following FcεRI crosslinking by DNP as compared to non-IgE-mediated AgNP exposure, we utilized a systems biology approach (outlined in Fig 2) to:

  1. Map genetic regions/loci and report “susceptibility” genes that are associated with differential mast cell degranulation via a non-IgE mediated (AgNP) mechanism
  2. Identify “response” genes that are differentially expressed following treatment to better understand the unique transcriptomic effect of AgNP compared to DNP
  3. Elucidate novel mechanisms of mast cell activation, which may be important in identifying potential targets involved in non-IgE mediated mast cell activation commonly observed in idiopathic disorders
thumbnail
Download:
Fig 2. Systems biology workflow.

Association mapping study utilizing 22 strains of mice were performed in conjunction with RNA-sequencing on high and low responder strains. ‘Susceptibility’ and ‘response’ genes were reported.

https://doi.org/10.1371/journal.pone.0193499.g002

Association mapping identifies unique susceptibility genes involved in mast cell degranulation

An association mapping study was performed on 22 strains of mice surveyed for mast cell degranulation. Using an efficient mixed-model approach (emma) on 198,278 informative SNPs, we identified chromosomal loci associated with quantitative differences in mast cell degranulation by AgNP and DNP (Fig 3). The mapping study identified three loci associated with AgNP-induced mast cell degranulation, which were mapped on chromosomes 1, 2, and 18 (Bonferroni corrected p-value: 4.3 × 10−3, 8.5 × 10−3, 4.6 × 10−6, respectively) (Fig 3A). One locus was associated with DNP-mediated mast cell degranulation on chromosome 9 (Bonferroni corrected p-value: 1.3 × 10−2) (Fig 3B). All associated loci exceeded a significance threshold of Bonferroni-corrected p-value <0.05. Genes residing in any of the 4 associated regions between both treatments were narrowed by only including genes that were expressed in the transcriptomic data (Table 1). The ≤20 closest genes to the SNP peak marker for each locus are outlined in Table 1. Utilizing the NHGRI-EBI catalog of previously reported GWAS, candidate genes were further narrowed to include QTLs and genetic variants previously associated with allergic immune response traits. Genes of interest are trafficking protein, kinesin binding 2 (Trak2) on chromosome 1 and TNF receptor–associated factor 6 (Traf6) on chromosome 2 because they were discovered together in several previously reported GWAS studies associated with atopic dermatitis, psoriasis, and rheumatoid arthritis [47–49].

thumbnail
Download:
Fig 3. Association mapping of AgNP or DNP treatment with mast cell degranulation.

Manhattan plot of chromosomal loci associated with (A) AgNP or (B) DNP treatment across 22 strains of mice. The figure inserts include further information on the peak marker, chromosome, p-values, and associated regions. Loci with p-values ≤ 7×10−5 (red line) and adjusted p-values ≤ 0.05 were considered significant.

https://doi.org/10.1371/journal.pone.0193499.g003

thumbnail
Download:
Table 1. Chromosomal loci significantly associated with mast cell degranulation by AgNP or DNP and genes residing in each locus.

https://doi.org/10.1371/journal.pone.0193499.t001

Transcriptomic analysis of high and low responder strains distinguishes novel response genes in AgNP-induced mast cell degranulation

RNA sequencing (RNA-seq) was performed on BMMCs grown from the high and low responder strains following DNP (RIIIS/J, C57BL/6J) or AgNP (C57BL/6J, LP/J) mediated degranulation. The filtered data set comprised approximately 750 million—50 base pair reads, almost evenly distributed between analyzed samples. Of the 16,504 genes from Mus musculus reference transcriptome that were expressed in mast cells, 3805 and 85 protein-coding genes were differentially expressed (DE) following DNP or AgNP treatment, respectively (FDR<0.05). Of the total 3890 DE genes, 2702 (DNP) and 21 (AgNP) genes displayed DE patterns that differed between the high and low responder strains (strain-dependent treatment effect; FDR<0.1). A complete list of all DE genes following AgNP and DNP treatment are outlined in the supplemental material (S1 Table).

Quantitative comparison of the transcriptomic effects due to AgNP and DNP are reported in Fig 4. A large majority of the total DE genes following both treatments were unique to DNP treatment (3754 out of 3890) (Fig 4A). Interestingly, 34 genes were DE following AgNP treatment, which were further evaluated as novel genetic targets involved in non-IgE-mediated mast cell activation. We examined global expression patterns for the DE genes following AgNP treatment using cluster analysis (Fig 4B). Significant variability in expression of individual genes within treatment groups was apparent. The heat map clearly outlines the abundance of DE genes observed in the low responder stain (LP/J) with marginal expression changes observed in the high responder strain (C57BL/6J), potentially indicating that many DE genes were protective in response to AgNP treatment (Fig 4B). Taken together, these initial analyses revealed that AgNP treatment has specific and different transcriptomic effects compared to DNP treatment, and the AgNP effects were further differentiated between the low (LP/J) and high responder (C57BL/6J) strains.

thumbnail
Download:
Fig 4. Quantitative comparison of differentially expressed transcripts and the biological responses that are significantly enriched in mast cells by AgNPs.

(A) Venn diagram showing the number of differentially expressed (DE) genes produced in response to DNP treatment (black) and AgNP treatment (grey) with a combined total of 3890 genes. (B) Heat map of the relative expression of 85 DE genes (FDR≤0.05) following AgNP exposure in the low responder (LP/J) and high responder strain (C57BL/6J). Expression patterns for individual genes are in rows. The color represents the direction of difference in expression with red indicating increased expression and green indicated decreased expression after treatment; The intensity of the color represents the magnitude of difference in expression (log2 scale) between AgNP treated cells and non-treated cells. (C) Gene ontology (GO) categories overrepresented in the 34 AgNP-responsive gene set were assessed using Pantherdb.org. GO categories are ranked by number of genes, listing the number of AgNP-responsive genes in each category (white bar) verses the number of genes expected in each category (black bar). Significance was set at p≤0.05 with a fold enrichment cutoff of ≥5. RTK, receptor tyrosine kinase; MAPK, mitogen-activated protein kinase.

https://doi.org/10.1371/journal.pone.0193499.g004

We focused further studies on the candidate gene list comprised of 34 AgNP-responsive genes (Table 2). Of specific interest are the 11 DE genes in which AgNP-induced contrasting responses in the high versus the low responder strains (Table 2, bold). We tested for gene ontology (GO) categories overrepresented in the gene set to identify unique biological responses of mast cells to AgNP (Fig 4C). Treatment of BMMCs with AgNP had a strong, positive effect on a number of fundamental and interconnected biological processes involved in stress response including histone modifications and signal transduction pathways. Specifically, several histone components (Hist1h1d, H3f3b, Hist3h2a), signal transduction components including the MAPK signaling cascade (Trib2, Ptger4, Fam43a), and stress response genes involved in kinase inhibitor activity and ubiquitin protein ligase activity (Txnip, Arrdc3, Arrdc4, Gadd45g, Neurl3) were significantly enriched in the AgNP-responsive genes. These results largely confirmed that AgNP-induced mast cell degranulation is a complex process that is likely mediated through a novel, receptor-driven mechanism leading to activation of cellular stress response pathways ultimately culminating in degranulation.

thumbnail
Download:
Table 2. Candidate genes identified by RNA-seq that are differentially expressed following AgNP treatment.

https://doi.org/10.1371/journal.pone.0193499.t002

Validation of select candidate genes

Expression of several AgNP-responsive genes were further validated to assess the reliability of our RNA-seq data set and identify potential gene biomarkers. Genes were chosen for a validation study based on their inclusion in the set of susceptible genes identified by association mapping in addition to the 34 AgNP-responsive genes (Fig 5). These genes include: Fam43a, Trib2, Dusp5, Sqstm1, and Traf6. qPCR analysis confirmed a significant difference (p<0.05) in effects of AgNP on expression between low and high responder strains for Txnip (Fig 5A). The remaining 6 genes demonstrated similar trends observed in the RNA-seq analysis, however the differences were not significant (S3 Fig). Consistent with the gene expression analysis, a significant decline in protein expression of TXNIP was observed following AgNP treatment in the high responder strain (C57BL/6J), where no observable effect was demonstrated in the low responder strain (Fig 5B).

thumbnail
Download:
Fig 5. Expression of Txnip following AgNP treatment.

Bone marrow-derived mast cells (BMMCs) were treated with AgNP for 1 h at 25 μg/ml and A) mRNA expression of Txnip was quantified by real-time quantitative polymerase chain reaction (qPCR). Values are expressed as fold change compared to non-treated cells (n = 3/group) normalized to Gapdh. B) Representative immunoblots of TXNIP expression in both high (C57BL/6J) and low (LP/J) responder strains following AgNP treatment (1 h at 25 μg/ml). Values are expressed as mean ± SEM of at least 3 independent experiments. * Indicates significant difference from controlled group (p≤ 0.05). # Indicates significant difference between strains (p≤ 0.05).

https://doi.org/10.1371/journal.pone.0193499.g005

Functional studies were performed on the prostaglandin E receptor (EP4) which was identified in the RNA-seq analysis (Ptger4), and demonstrated a 25% decrease in gene expression following AgNP exposure in the low responder strain with no significant expression effect in the high responder strain (Table 2). To evaluate a plausible role for EP4 receptor function in mast cell degranulation, BMMCs isolated from C57BL/6J were assessed for degranulation following pre-treatment with EP4 antagonist and agonists. BMMCs pre-treated with an EP4 selective antagonist caused a significant reduction in AgNP-induced mast cell degranulation (Fig 6). In addition, no significant effect was observed following pre-treatment with PGE2, an EP4 ligand. Similar results were observed following pre-treatment with indomethacin, a nonselective inhibitor for COX-1 and -2, which prevents PGE2 synthesis (Fig 6). This suggests that AgNP may possibly be activating mast cells through an EP4-mediated mechanism, which is independent from PGE2 stimulation.

thumbnail
Download:
Fig 6. Prostaglandin receptor (EP4) expression and inhibition of AgNP-directed mast cell degranulation.

Bone marrow-derived mast cell (BMMCs) isolated from C57BL/6J mice were evaluated for degranulation by measuring ß-hexosaminidase release following 1 h exposure of AgNP at 25 μg/ml ± pre-treatment with EP4 selective inhibitor (10 μM), prostaglandin (PGE2, 10 μM), or Cox-1 and -2 inhibitors (indomethacin, 5 μM). Values are expressed as mean ± SEM (n = 3). * Indicates significant difference compared to non-treated control and # indicates significant difference compared to AgNP treatment alone (p ≤ 0.05).

https://doi.org/10.1371/journal.pone.0193499.g006

Discussion

Using a systems biology approach, we aimed to identify the genetic factors and potential mechanisms of activation that are contributing to non-IgE-mediated mast cell degranulation by silver nanoparticles (AgNPs). In our study, utilizing mast cells from 23 strains of mice, we evaluated mast cell degranulation induced by AgNP exposure and compared the response to classical activation induced by crosslinking of the high affinity IgE receptor (e.g. FcεRI) by dinitrophenyl (DNP). Interestingly, we observed contrasting treatment-induced responses across strains indicating a distinct genetic component to AgNP directed mast cell degranulation. Following FcεRI stimulation, the highest responding strain was RIIIS/J, while the lower responder strain was C57BL/6J, a well-characterized strain commonly used in biomedical research. On the contrary, strain distribution was quite different following AgNP treatment, in which C57BL/6J was one of the highest responding strains. The dissimilar strain distribution patterns between the two treatments clearly indicates that AgNPs are activating mast cells via a novel mechanism compared to FcεRI crosslinking. Previous electron microscopy studies in our laboratory have demonstrated minimal uptake of AgNPs by mast cells, suggesting a potential mechanism mediated via surface interaction of the particle with a membrane receptor [11, 17]. Following bioinformatics analyses, putative risk loci associated with AgNP and DNP responses and represented susceptibility genes were reported. In addition, RNA sequencing (RNA-seq) of high and low responder strains identified a subset of AgNP-responsive genes that displayed a unique profile compared to DNP stimulation.

No previous genetic studies have identified mechanisms of mast cell degranulation by nanoparticle exposure; therefore, our interests were in the identification of novel genes of biological relevance to non-IgE mediated degranulation. An association mapping study of approximately 200,000 SNPs characterized across 22 strains (one strain was excluded) of mice identified 3 loci (located on chromosome 1, 2, and 18) associated with AgNP-induced mast cell degranulation compared to 1 non-overlapping locus (chromosome 9) associated with DNP stimulation (Bonferroni corrected p-value <0.05). The differences in associated loci between treatments, AgNP and DNP, further suggest that AgNP is likely mediated through a receptor-driven pathway distinct from IgE crosslinking leading to mast cell degranulation.

Several of the identified AgNP-associated loci (Chromosome 2 and 18) in the current study have previously been implicated in other association mapping studies concomitant with circulating toxic metals [50] and immune-mediated diseases [51–54]. Interestingly, the AgNP-associated loci identified Trak2 and Traf6 on chromosomes 1 and 2, previously reported together as genetic variants associated with several common inflammatory diseases including atopic dermatitis, psoriasis, and rheumatoid arthritis [47–49]. Traf6 is an essential signal transducer downstream from tumor necrosis factor (TNF) receptor and toll-like receptor (TLR) superfamily members and has been shown to contribute to mast cell activation [55, 56]. Interestingly, Traf6-deficient mast cells demonstrate normal IgE-mediated degranulation, raising the possibility that this protein may regulate alternative pathways involved in non-IgE-mediated degranulation. For example, Madera-Salcedo et. al., demonstated a role for a Traf6/ß-arrestin complex in suppressing mast cell activation by modulating opioid receptors, a class of G protein-coupled receptors (GPCRs) that recruit ß-arrestins and mediate signaling pathways leading to inhibition of Ca2+ channels [57, 58]. Taken together, we postulate that GPCR activation may be contributing to AgNP-mediated mast cell degranulation.

Additionally, the loci on chromosome 1 contains a cohort of biologically interesting genes (Cflar, Clk1, Ndfufb3, Nif3l1, Bzw1, Ppil3, Casp8, Fam126b) recently reported in another association mapping study to be predictive markers in cellular toxicity to triptolide, an anti-inflammatory and immunosuppressive compound commonly used in traditional Chinese medicine to treat several immune complex diseases such as rheumatoid arthritis and systemic lupus erythematous [59, 60]. Studies have shown that triptolide can alter multiple signaling pathways leading to the down-regulation of c-Kit and inhibition of nuclear factor-kappa B (NF-κB) and phosphoinositide 3-kinase (PI3K) signaling, key cell signaling molecules in mast cell activation [61, 62]. These findings may be important, as we have recently reported that PI3K phosphorylation is required to elicit mast cell degranulation by AgNPs [17], however no functional data exists for many of the identified genes in evaluating their respective roles in mast cell degranulation and therefore warrants further study.

Using RNA-seq for the first time in murine mast cells exposed to DNP or AgNP, we have uncovered genetic targets involved in non-IgE mediated mast cell responses compared to the traditional activation by FcεRI cross-linking. Transcriptomic analysis identified a combined 3890 DE genes following exposure to either DNP or AgNP. Interestingly, many of the genes were DE following DNP treatment (3805), indicating the abundant gene expression changes occurring following stimulation with antigen crosslinking of the IgE receptor. Although interesting, we chose to focus our current study on the DE genes following AgNP treatment in the aim of better understanding mechanisms of activation by non-IgE triggers. In so forth, the 34 DE genes in response to AgNP were the major focus of the discussion, with special attention on the 11 genes displaying DE patterns that differed between the high and low responder strains.

Statistical enrichment analysis of the 34 AgNP-responsive genes demonstrated that many of the genes were involved in stress response and signal transduction (Fig 4C). Central to these responses are Txnip, Arrdc3 and Arrdc4, all three are members of the arrestin superfamily which have been shown to both negatively and positively regulate GPCR signaling [63, 64]. In addition to their previously mentioned interaction with Traf6, they are also involved in complement receptor C3aR desensitization and NF-κB activation in mast cell degranulation [65]. A member of the α-arrestin family, thioredoxin interacting protein (Txnip), (also called vitamin D3 up-regulated protein (VDUP1) and thioredoxin-binding protein-2 (TB-2)) may play a protective role in AgNP-induced mast cell degranulation. We observed a 300% increase in expression in the low responder strain (LP/J) in comparison to minimal expression changes in the high responder strain (C57BL/6J) in response to AgNP. Txnip encodes a multifunctional protein that has emerged as a key component of cellular redox balance by inhibiting the reducing activity of thioredoxin (TRX) [66]. It has been established that oxidative stress can contribute to IgE-mediated mast cell degranulation, as observed by a significant suppression of histamine release by FcεRI crosslinking in a TRX-transgenic mouse model compared to wild-type [67]. In addition to redox-dependent functions, Txnip also functions as an arrestin protein modulating expression and function of key receptor(s) and downstream signal transduction pathways [68].

We have recently reported that AgNPs are not readily taken up in mast cells, rather electron microscopy imaging suggests that initiation of mast cell degranulation appears to be membrane or receptor-mediated [17]. Therefore, we hypothesize that Txnip may be regulating a key receptor in the alternative, non-IgE mediated mast cell degranulation pathway. In the transcriptomic study, we identified the prostaglandin E receptor 4 (EP4) as an AgNP-responsive gene (Ptger4). This gene warrants further analysis through knockout or knockdown studies as it was the only receptor identified in the transcriptomic analysis. EP receptor subtypes 1–4 respond to prostaglandin (PGE2) and activate various signaling pathways [69]. An arachidonic acid metabolite, PGE2 functions as both a pro- and anti-inflammatory lipid mediator that balances TH1/TH2 responses via stimulation of different EP receptors. In relation to mast cell degranulation, different receptor subtypes vary the biological response. PGE2-mediated activation of EP3 and EP4 receptors can induce mast cell degranulation by increasing intracellular Ca2+, while EP2 stimulation has the opposite effect as shown by inhibition of IgE-mediated histamine release in human lung mast cells [70–72]. Additionally, Nguyen et. al. observed a genetic effect in PGE2-induced inflammation by mast cells [73]. Currently, we demonstrated that BMMCs pre-treated with a selective EP4 antagonist partially, but significantly inhibited AgNP-induced mast cell degranulation in the high responder strain (C57BL/6J; Fig 6). Additionally, there was no significant effect following pre-treatment with PGE2 or indomethacin. This suggests that AgNPs may be indirectly mediating a response via EP4 activation however the exact mechanism is not known and should be further investigated with additional functional studies using knockout or knockdown models. We have previously reported a role for scavenger receptor B1 (SR-B1) in mast cell degranulation by AgNPs, however this receptor was not identified in our transcriptomics analysis. Our statistical analyses were utilized to identify a subset of genes with expression patterns that significantly differed between strains, however a shortcoming in the study may have missed genes that were significantly expressed in both strains. Also, the lack of identification could be due to the short time frame of the study (1 hr) [10].

Taken together, we conclude that AgNP-induced mast cell degranulation is a complex process that is likely dependent on many factors, potentially mediated through the regulation of GPCRs upstream of key signaling cascades. Most research efforts are focused on mast cell signaling through the high affinity IgE receptor (FcεRI), although an alternative pathway may contribute in relation to AgNP exposure which may involve GPCR-mediated activation of PI3K and MAPK signaling cascades [62]. Still to date, no studies have linked GPCR activation to nanoparticle-induced allergic inflammation. Therefore, it may be of interest to target GPCR-regulators, such as the α/β-ARR family members in future studies. Txnip may be a plausible therapeutic target because of its dual role in redox-dependent and redox-independent regulation of GPCR-mediated signaling pathways. For example, treatment with an anti-asthmatic drug, eugenol, led to the up-regulation of Txnip and down-regulation of NF-κB, ameliorating the progression of asthma in an OVA-induced mouse model [74]. In addition, we observed a novel connection between EP4 activation and mast cell degranulation by a non-IgE trigger. However, future functional studies need to be performed to validate the genes identified, which is beyond the scope of this study. By our knowledge, this is the first study performed combining bioinformatics analyses to study the association of genetic variants by association mapping and transcriptional changes by RNA-sequencing to determine a mechanism of mast cell activation by a non-IgE trigger. Association mapping has become a powerful tool in reproducibly identifying genetic regions associated with common disease traits due to recent advances in genotyping and publicly available SNP databases [75]. RNA-sequencing is gaining traction as means to characterize genomic expression and has been utilized in numerous studies to acquire knowledge about mechanisms of complex diseases [76]. These complementary analyses offer additional insights into novel disease pathology that one test alone could not infer. This study has provided novel gene targets and proposed mechanisms to better understand non-IgE mediated mast cell degranulation in addition to addressing the emerging concern regarding nanoparticle exposures and human health and safety.

Supporting information

S1 Fig. Surface expression of mature bone marrow derived mast cells (BMMCs).

FcεRI (white bar) and c-Kit (grey bar) were analyzed on the surface of BMMCs at week 5 via flow cytometry. Values are expressed as mean ± SEM (n = 2-3/group).

https://doi.org/10.1371/journal.pone.0193499.s001

(TIF)

S2 Fig. Dose-dependent increase in human mast cell degranulation following AgNP exposure.

Human mast cells (LUVA; Kerafast) were evaluated for degranulation by measuring ß-hexosaminidase release following 1 h exposure of AgNP at 25 and 50 μg/ml. Values are expressed as mean ± SEM of at least 3 independent experiments. * Indicates significant difference from controlled group (p≤ 0.05).

https://doi.org/10.1371/journal.pone.0193499.s002

(TIFF)

S3 Fig. Expression of selected AgNP-responsive genes.

Bone marrow-derived mast cells (BMMCs) were stimulated with AgNP for 1 h at 25 μg/ml and the mRNA expression of select genes was quantified by real-time quantitative polymerase chain reaction (qPCR). Values are expressed as fold change compared to non-treated cells (n = 3/group) normalized to Gapdh. Values are expressed as mean ± SEM of at least 3 independent experiments.

https://doi.org/10.1371/journal.pone.0193499.s003

(TIF)

S1 Table. RNA-seq data for genes differentially expressed following AgNP and DNP treatment.

Ensembl gene IDs and MGI gene symbols are reported. *Indicates p-value of treatment. **Indicates p-value of strain-dependent treatment effect. Chr., chromosome. Start, chromosomal start.

https://doi.org/10.1371/journal.pone.0193499.s004

(XLSX)

References

  1. 1. Vance ME, Kuiken T, Vejerano EP, McGinnis SP, Hochella MF Jr., Rejeski D, et al. Nanotechnology in the real world: Redeveloping the nanomaterial consumer products inventory. Beilstein J Nanotechnol. 2015;6:1769–80. Epub 2015/10/02. pmid:26425429; PubMed Central PMCID: PMCPMC4578396.
  2. 2. Cohen J, Deloid G, Pyrgiotakis G, Demokritou P. Interactions of engineered nanomaterials in physiological media and implications for in vitro dosimetry. Nanotoxicology. 2013;7(4):417–31. pmid:22393878; PubMed Central PMCID: PMCPMC4479135.
  3. 3. Podila R, Brown JM. Toxicity of engineered nanomaterials: a physicochemical perspective. Journal of biochemical and molecular toxicology. 2013;27(1):50–5. pmid:23129019; PubMed Central PMCID: PMC3778677.
  4. 4. Wiesner MR, Lowry GV, Alvarez P, Dionysiou D, Biswas P. Assessing the risks of manufactured nanomaterials. Environ Sci Technol. 2006;40(14):4336–45. pmid:16903268.
  5. 5. Wang X, Reece SP, Brown JM. Immunotoxicological impact of engineered nanomaterial exposure: mechanisms of immune cell modulation. Toxicol Mech Methods. 2013;23(3):168–77. pmid:23256453; PubMed Central PMCID: PMC3773497.
  6. 6. Chen EY, Garnica M, Wang YC, Mintz AJ, Chen CS, Chin WC. A mixture of anatase and rutile TiO(2) nanoparticles induces histamine secretion in mast cells. Part Fibre Toxicol. 2012;9:2. pmid:22260553; PubMed Central PMCID: PMCPMC3275461.
  7. 7. Huang YF, Liu H, Xiong X, Chen Y, Tan W. Nanoparticle-mediated IgE-receptor aggregation and signaling in RBL mast cells. J Am Chem Soc. 2009;131(47):17328–34. pmid:19929020; PubMed Central PMCID: PMCPMC2786779.
  8. 8. Maurer-Jones MA, Lin YS, Haynes CL. Functional assessment of metal oxide nanoparticle toxicity in immune cells. ACS Nano. 2010;4(6):3363–73. pmid:20481555.
  9. 9. Yang W, Lee S, Lee J, Bae Y, Kim D. Silver nanoparticle-induced degranulation observed with quantitative phase microscopy. J Biomed Opt. 2010;15(4):045005. pmid:20799800.
  10. 10. Aldossari AA, Shannahan JH, Podila R, Brown JM. Influence of physicochemical properties of silver nanoparticles on mast cell activation and degranulation. Toxicol In Vitro. 2015;29(1):195–203. pmid:25458489; PubMed Central PMCID: PMCPMC4294974.
  11. 11. Johnson MM, Mendoza R, Raghavendra AJ, Podila R, Brown JM. Contribution of engineered nanomaterials physicochemical properties to mast cell degranulation. Sci Rep. 2017;7:43570. Epub 2017/03/07. pmid:28262689; PubMed Central PMCID: PMCPMC5337938.
  12. 12. Brown JM, Wilson TM, Metcalfe DD. The mast cell and allergic diseases: role in pathogenesis and implications for therapy. Clin Exp Allergy. 2008;38(1):4–18. pmid:18031566.
  13. 13. Galli SJ, Grimbaldeston M, Tsai M. Immunomodulatory mast cells: negative, as well as positive, regulators of immunity. Nat Rev Immunol. 2008;8(6):478–86. pmid:18483499; PubMed Central PMCID: PMCPMC2855166.
  14. 14. Akin C. Mast cell activation disorders. J Allergy Clin Immunol Pract. 2014;2(3):252–7 e1; quiz 8. pmid:24811013.
  15. 15. Metcalfe DD, Baram D, Mekori YA. Mast cells. Physiol Rev. 1997;77(4):1033–79. pmid:9354811.
  16. 16. Tkaczyk C, Okayama Y, Woolhiser MR, Hagaman DD, Gilfillan AM, Metcalfe DD. Activation of human mast cells through the high affinity IgG receptor. Mol Immunol. 2002;38(16–18):1289–93. pmid:12217397.
  17. 17. Alsaleh NB, Persaud I, Brown JM. Silver Nanoparticle-Directed Mast Cell Degranulation Is Mediated through Calcium and PI3K Signaling Independent of the High Affinity IgE Receptor. PLoS One. 2016;11(12):e0167366. pmid:27907088; PubMed Central PMCID: PMCPMC5131952.
  18. 18. Woolhiser MR, Brockow K, Metcalfe DD. Activation of human mast cells by aggregated IgG through FcgammaRI: additive effects of C3a. Clin Immunol. 2004;110(2):172–80. pmid:15003814.
  19. 19. McCurdy JD, Olynych TJ, Maher LH, Marshall JS. Cutting edge: distinct Toll-like receptor 2 activators selectively induce different classes of mediator production from human mast cells. J Immunol. 2003;170(4):1625–9. pmid:12574323.
  20. 20. Medic N, Vita F, Abbate R, Soranzo MR, Pacor S, Fabbretti E, et al. Mast cell activation by myelin through scavenger receptor. J Neuroimmunol. 2008;200(1–2):27–40. pmid:18657868.
  21. 21. Brown JM, Swindle EJ, Kushnir-Sukhov NM, Holian A, Metcalfe DD. Silica-directed mast cell activation is enhanced by scavenger receptors. Am J Respir Cell Mol Biol. 2007;36(1):43–52. pmid:16902192; PubMed Central PMCID: PMCPMC1899302.
  22. 22. Zhang D, Spielmann A, Wang L, Ding G, Huang F, Gu Q, et al. Mast-cell degranulation induced by physical stimuli involves the activation of transient-receptor-potential channel TRPV2. Physiol Res. 2012;61(1):113–24. pmid:21574765.
  23. 23. Nakamura Y, Oscherwitz J, Cease KB, Chan SM, Munoz-Planillo R, Hasegawa M, et al. Staphylococcus delta-toxin induces allergic skin disease by activating mast cells. Nature. 2013;503(7476):397–401. pmid:24172897; PubMed Central PMCID: PMCPMC4090780.
  24. 24. Miyazaki D, Nakamura T, Toda M, Cheung-Chau KW, Richardson RM, Ono SJ. Macrophage inflammatory protein-1alpha as a costimulatory signal for mast cell-mediated immediate hypersensitivity reactions. J Clin Invest. 2005;115(2):434–42. pmid:15650768; PubMed Central PMCID: PMCPMC544033.
  25. 25. Jolly PS, Bektas M, Olivera A, Gonzalez-Espinosa C, Proia RL, Rivera J, et al. Transactivation of sphingosine-1-phosphate receptors by FcepsilonRI triggering is required for normal mast cell degranulation and chemotaxis. J Exp Med. 2004;199(7):959–70. pmid:15067032; PubMed Central PMCID: PMCPMC2211871.
  26. 26. Schafer B, Piliponsky AM, Oka T, Song CH, Gerard NP, Gerard C, et al. Mast cell anaphylatoxin receptor expression can enhance IgE-dependent skin inflammation in mice. J Allergy Clin Immunol. 2013;131(2):541–8 e1-9. pmid:22728083; PubMed Central PMCID: PMCPMC3597773.
  27. 27. Gaudenzio N, Sibilano R, Marichal T, Starkl P, Reber LL, Cenac N, et al. Different activation signals induce distinct mast cell degranulation strategies. J Clin Invest. 2016;126(10):3981–98. Epub 2016/09/20. pmid:27643442; PubMed Central PMCID: PMCPMC5096814.
  28. 28. Cho HY, Jedlicka AE, Reddy SP, Zhang LY, Kensler TW, Kleeberger SR. Linkage analysis of susceptibility to hyperoxia. Nrf2 is a candidate gene. Am J Respir Cell Mol Biol. 2002;26(1):42–51. pmid:11751202.
  29. 29. Nichols JL, Gladwell W, Verhein KC, Cho HY, Wess J, Suzuki O, et al. Genome-wide association mapping of acute lung injury in neonatal inbred mice. FASEB J. 2014;28(6):2538–50. pmid:24571919; PubMed Central PMCID: PMCPMC4021442.
  30. 30. Brewer JP, Kisselgof AB, Martin TR. Genetic variability in pulmonary physiological, cellular, and antibody responses to antigen in mice. Am J Respir Crit Care Med. 1999;160(4):1150–6. pmid:10508801.
  31. 31. Kelada SN, Wilson MS, Tavarez U, Kubalanza K, Borate B, Whitehead GS, et al. Strain-dependent genomic factors affect allergen-induced airway hyperresponsiveness in mice. Am J Respir Cell Mol Biol. 2011;45(4):817–24. pmid:21378263; PubMed Central PMCID: PMCPMC3208613.
  32. 32. Whitehead GS, Walker JK, Berman KG, Foster WM, Schwartz DA. Allergen-induced airway disease is mouse strain dependent. Am J Physiol Lung Cell Mol Physiol. 2003;285(1):L32–42. pmid:12626335.
  33. 33. Scoville DK, White CC, Botta D, McConnachie LA, Zadworny ME, Schmuck SC, et al. Susceptibility to quantum dot induced lung inflammation differs widely among the Collaborative Cross founder mouse strains. Toxicol Appl Pharmacol. 2015;289(2):240–50. pmid:26476918; PubMed Central PMCID: PMCPMC4651801.
  34. 34. Jones SW, Roberts RA, Robbins GR, Perry JL, Kai MP, Chen K, et al. Nanoparticle clearance is governed by Th1/Th2 immunity and strain background. J Clin Invest. 2013;123(7):3061–73. pmid:23778144; PubMed Central PMCID: PMCPMC3696555.
  35. 35. Yang Y, Westerhoff P. Presence in, and release of, nanomaterials from consumer products. Adv Exp Med Biol. 2014;811:1–17. Epub 2014/04/01. pmid:24683024.
  36. 36. Ghazalpour A, Rau CD, Farber CR, Bennett BJ, Orozco LD, van Nas A, et al. Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits. Mamm Genome. 2012;23(9–10):680–92. pmid:22892838; PubMed Central PMCID: PMC3586763.
  37. 37. Iwaki S, Tkaczyk C, Satterthwaite AB, Halcomb K, Beaven MA, Metcalfe DD, et al. Btk plays a crucial role in the amplification of Fc epsilonRI-mediated mast cell activation by kit. J Biol Chem. 2005;280(48):40261–70. pmid:16176929.
  38. 38. Sul JH, Raj T, de Jong S, de Bakker PI, Raychaudhuri S, Ophoff RA, et al. Accurate and fast multiple-testing correction in eQTL studies. Am J Hum Genet. 2015;96(6):857–68. pmid:26027500; PubMed Central PMCID: PMCPMC4457958.
  39. 39. Bogue MA, Peters LL, Paigen B, Korstanje R, Yuan R, Ackert-Bicknell C, et al. Accessing Data Resources in the Mouse Phenome Database for Genetic Analysis of Murine Life Span and Health Span. J Gerontol A Biol Sci Med Sci. 2016;71(2):170–7. pmid:25533306; PubMed Central PMCID: PMCPMC4707687.
  40. 40. Yang H, Wang JR, Didion JP, Buus RJ, Bell TA, Welsh CE, et al. Subspecific origin and haplotype diversity in the laboratory mouse. Nat Genet. 2011;43(7):648–55. pmid:21623374; PubMed Central PMCID: PMCPMC3125408.
  41. 41. Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, et al. Efficient control of population structure in model organism association mapping. Genetics. 2008;178(3):1709–23. pmid:18385116; PubMed Central PMCID: PMCPMC2278096.
  42. 42. Rosenbloom KR, Armstrong J, Barber GP, Casper J, Clawson H, Diekhans M, et al. The UCSC Genome Browser database: 2015 update. Nucleic Acids Res. 2015;43(Database issue):D670–81. pmid:25428374; PubMed Central PMCID: PMCPMC4383971.
  43. 43. Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42(Database issue):D1001–6. pmid:24316577; PubMed Central PMCID: PMCPMC3965119.
  44. 44. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. pmid:21816040; PubMed Central PMCID: PMCPMC3163565.
  45. 45. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;57(1):289–300.
  46. 46. Mi H, Poudel S, Muruganujan A, Casagrande JT, Thomas PD. PANTHER version 10: expanded protein families and functions, and analysis tools. Nucleic Acids Res. 2016;44(D1):D336–42. pmid:26578592; PubMed Central PMCID: PMCPMC4702852.
  47. 47. Paternoster L, Standl M, Waage J, Baurecht H, Hotze M, Strachan DP, et al. Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis. Nat Genet. 2015;47(12):1449–56. pmid:26482879; PubMed Central PMCID: PMCPMC4753676.
  48. 48. Baurecht H, Hotze M, Brand S, Buning C, Cormican P, Corvin A, et al. Genome-wide comparative analysis of atopic dermatitis and psoriasis gives insight into opposing genetic mechanisms. Am J Hum Genet. 2015;96(1):104–20. pmid:25574825; PubMed Central PMCID: PMCPMC4289690.
  49. 49. Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014;506(7488):376–81. pmid:24390342; PubMed Central PMCID: PMCPMC3944098.
  50. 50. Ng E, Lind PM, Lindgren C, Ingelsson E, Mahajan A, Morris A, et al. Genome-wide association study of toxic metals and trace elements reveals novel associations. Hum Mol Genet. 2015;24(16):4739–45. pmid:26025379; PubMed Central PMCID: PMCPMC4512629.
  51. 51. Wan YI, Shrine NR, Soler Artigas M, Wain LV, Blakey JD, Moffatt MF, et al. Genome-wide association study to identify genetic determinants of severe asthma. Thorax. 2012;67(9):762–8. pmid:22561531.
  52. 52. Yucesoy B, Kaufman KM, Lummus ZL, Weirauch MT, Zhang G, Cartier A, et al. Genome-Wide Association Study Identifies Novel Loci Associated With Diisocyanate-Induced Occupational Asthma. Toxicol Sci. 2015;146(1):192–201. pmid:25918132; PubMed Central PMCID: PMCPMC4560052.
  53. 53. Duan QL, Lasky-Su J, Himes BE, Qiu W, Litonjua AA, Damask A, et al. A genome-wide association study of bronchodilator response in asthmatics. Pharmacogenomics J. 2014;14(1):41–7. pmid:23508266; PubMed Central PMCID: PMCPMC3706515.
  54. 54. Vonk JM, Scholtens S, Postma DS, Moffatt MF, Jarvis D, Ramasamy A, et al. Adult onset asthma and interaction between genes and active tobacco smoking: The GABRIEL consortium. PLoS One. 2017;12(3):e0172716. Epub 2017/03/03. pmid:28253294; PubMed Central PMCID: PMCPMC5333809.
  55. 55. Rothe M, Wong SC, Henzel WJ, Goeddel DV. A novel family of putative signal transducers associated with the cytoplasmic domain of the 75 kDa tumor necrosis factor receptor. Cell. 1994;78(4):681–92. pmid:8069916.
  56. 56. Yang YJ, Chen W, Carrigan SO, Chen WM, Roth K, Akiyama T, et al. TRAF6 specifically contributes to FcepsilonRI-mediated cytokine production but not mast cell degranulation. J Biol Chem. 2008;283(46):32110–8. pmid:18772140.
  57. 57. Madera-Salcedo IK, Cruz SL, Gonzalez-Espinosa C. Morphine prevents lipopolysaccharide-induced TNF secretion in mast cells blocking IkappaB kinase activation and SNAP-23 phosphorylation: correlation with the formation of a beta-arrestin/TRAF6 complex. J Immunol. 2013;191(6):3400–9. pmid:23960234.
  58. 58. Rozenfeld R, Devi LA. Receptor heterodimerization leads to a switch in signaling: beta-arrestin2-mediated ERK activation by mu-delta opioid receptor heterodimers. FASEB J. 2007;21(10):2455–65. pmid:17384143; PubMed Central PMCID: PMCPMC3131006.
  59. 59. Chauhan L, Jenkins GD, Bhise N, Feldberg T, Mitra-Ghosh T, Fridley BL, et al. Genome-wide association analysis identified splicing single nucleotide polymorphism in CFLAR predictive of triptolide chemo-sensitivity. BMC Genomics. 2015;16:483. pmid:26121980; PubMed Central PMCID: PMCPMC4485866.
  60. 60. Pan J. RNA polymerase—an important molecular target of triptolide in cancer cells. Cancer Lett. 2010;292(2):149–52. pmid:20045594.
  61. 61. Zhou GS, Hu Z, Fang HT, Zhang FX, Pan XF, Chen XQ, et al. Biologic activity of triptolide in t(8;21) acute myeloid leukemia cells. Leuk Res. 2011;35(2):214–8. pmid:20691473.
  62. 62. Gilfillan AM, Tkaczyk C. Integrated signalling pathways for mast-cell activation. Nat Rev Immunol. 2006;6(3):218–30. pmid:16470226.
  63. 63. Luttrell LM, Lefkowitz RJ. The role of beta-arrestins in the termination and transduction of G-protein-coupled receptor signals. J Cell Sci. 2002;115(Pt 3):455–65. pmid:11861753.
  64. 64. Lefkowitz RJ. G protein-coupled receptors. III. New roles for receptor kinases and beta-arrestins in receptor signaling and desensitization. J Biol Chem. 1998;273(30):18677–80. pmid:9668034.
  65. 65. Vibhuti A, Gupta K, Subramanian H, Guo Q, Ali H. Distinct and shared roles of beta-arrestin-1 and beta-arrestin-2 on the regulation of C3a receptor signaling in human mast cells. PLoS One. 2011;6(5):e19585. pmid:21589858; PubMed Central PMCID: PMCPMC3093384.
  66. 66. Junn E, Han SH, Im JY, Yang Y, Cho EW, Um HD, et al. Vitamin D3 up-regulated protein 1 mediates oxidative stress via suppressing the thioredoxin function. J Immunol. 2000;164(12):6287–95. pmid:10843682.
  67. 67. Son A, Nakamura H, Kondo N, Matsuo Y, Liu W, Oka S, et al. Redox regulation of mast cell histamine release in thioredoxin-1 (TRX) transgenic mice. Cell Res. 2006;16(2):230–9. pmid:16474438.
  68. 68. Ferguson SS, Barak LS, Zhang J, Caron MG. G-protein-coupled receptor regulation: role of G-protein-coupled receptor kinases and arrestins. Can J Physiol Pharmacol. 1996;74(10):1095–110. pmid:9022829.
  69. 69. Feng C, Beller EM, Bagga S, Boyce JA. Human mast cells express multiple EP receptors for prostaglandin E2 that differentially modulate activation responses. Blood. 2006;107(8):3243–50. pmid:16357326; PubMed Central PMCID: PMCPMC1895755.
  70. 70. Peachell PT, MacGlashan DW Jr., Lichtenstein LM, Schleimer RP. Regulation of human basophil and lung mast cell function by cyclic adenosine monophosphate. J Immunol. 1988;140(2):571–9. pmid:2447182.
  71. 71. Kay LJ, Gilbert M, Pullen N, Skerratt S, Farrington J, Seward EP, et al. Characterization of the EP receptor subtype that mediates the inhibitory effects of prostaglandin E2 on IgE-dependent secretion from human lung mast cells. Clin Exp Allergy. 2013;43(7):741–51. pmid:23786281.
  72. 72. Nguyen M, Solle M, Audoly LP, Tilley SL, Stock JL, McNeish JD, et al. Receptors and signaling mechanisms required for prostaglandin E2-mediated regulation of mast cell degranulation and IL-6 production. J Immunol. 2002;169(8):4586–93. pmid:12370397.
  73. 73. Nguyen M, Pace AJ, Koller BH. Age-induced reprogramming of mast cell degranulation. J Immunol. 2005;175(9):5701–7. pmid:16237060.
  74. 74. Pan C, Dong Z. Antiasthmatic Effects of Eugenol in a Mouse Model of Allergic Asthma by Regulation of Vitamin D3 Upregulated Protein 1/NF-kappaB Pathway. Inflammation. 2015;38(4):1385–93. pmid:25588851.
  75. 75. Manolio TA, Collins FS. The HapMap and genome-wide association studies in diagnosis and therapy. Annu Rev Med. 2009;60:443–56. pmid:19630580; PubMed Central PMCID: PMCPMC2717504.
  76. 76. Costa V, Aprile M, Esposito R, Ciccodicola A. RNA-Seq and human complex diseases: recent accomplishments and future perspectives. Eur J Hum Genet. 2013;21(2):134–42. pmid:22739340; PubMed Central PMCID: PMCPMC3548270.