- Open Access
Toxicogenomic analysis of susceptibility to inhaled urban particulate matter in mice with chronic lung inflammation
© Thomson et al; licensee BioMed Central Ltd. 2009
- Received: 09 December 2008
- Accepted: 11 March 2009
- Published: 11 March 2009
Individuals with chronic lung disease are at increased risk of adverse health effects from airborne particulate matter. Characterization of underlying pollutant-phenotype interactions may require comprehensive strategies. Here, a toxicogenomic approach was used to investigate how inflammation modifies the pulmonary response to urban particulate matter.
Transgenic mice with constitutive pulmonary overexpression of tumour necrosis factor (TNF)-α under the control of the surfactant protein C promoter and wildtype littermates (C57BL/6 background) were exposed by inhalation for 4 h to particulate matter (0 or 42 mg/m3 EHC-6802) and euthanized 0 or 24 h post-exposure. The low alveolar dose of particles (16 μg) did not provoke an inflammatory response in the lungs of wildtype mice, nor exacerbate the chronic inflammation in TNF animals. Real-time PCR confirmed particle-dependent increases of CYP1A1 (30–100%), endothelin-1 (20–40%), and metallothionein-II (20–40%) mRNA in wildtype and TNF mice (p < 0.05), validating delivery of a biologically-effective dose. Despite detection of striking genotype-related differences, including activation of immune and inflammatory pathways consistent with the TNF-induced pathology, and time-related effects attributable to stress from nose-only exposure, microarray analysis failed to identify effects of the inhaled particles. Remarkably, the presence of chronic inflammation did not measurably amplify the transcriptional response to particulate matter.
Our data support the hypothesis that health effects of acute exposure to urban particles are dominated by activation of specific physiological response cascades rather than widespread changes in gene expression.
- Intratracheal Instillation
- Internal Dose
- Urban Particle
- Alveolar Surface Area
- Genotype Main Effect
Increased levels of particulate matter (PM) are associated with respiratory and cardiovascular morbidity and mortality [1–3]. Acute adverse health effects associated with PM are generally attributed to susceptible subgroups of the population, including individuals with respiratory diseases such as chronic obstructive pulmonary disease (COPD) [4–7]. Exacerbation of existing inflammation and increased generation of oxidative stress provoking an acute phase response is one mechanism through which PM may cause adverse cardiovascular effects . Inhaled particles aggravate existing epithelial lesions in the lungs , and impairment of epithelial barrier function can result in increased translocation of particles to the interstitium . Once in the interstitium, particles are less likely to be cleared via macrophage phagocytosis, and can cause interstitial inflammation, induce direct effects on local cell populations (including macrophages, fibroblasts, endothelial cells, and neutrophils) and drain to local lymph nodes [11, 12]. Chronic inflammation and oxidative stress may also prime the lungs to respond to particles with increased production of reactive oxygen species and inflammatory mediators, exacerbating the existing disease state . In addition, inflammation and tissue damage may facilitate direct interaction of PM factors with the pulmonary endothelium, disturbing vascular function [14, 15] and increasing the risk of acute cardiac events from plaque instability or reduced myocardial perfusion. Investigating the interaction of inhaled particles and host factors governing susceptibility may be critical to understanding acute health effects related to PM exposure.
Urban PM is a heterogeneous mixture of organic and inorganic compounds likely to provoke a number of transcriptional and physiological responses in lung tissue [16, 17]. We have previously established that mRNA and plasma levels of the potent vasoconstrictor endothelin (ET)-1 are increased 25–50% after particle inhalation, a response that does not require acute lung injury [9, 14, 18–20]. Increases of plasma ET-1 in this range have high predictive value for chronic heart failure , indicating that the magnitude of ET-1 change need not be large to have pathophysiological relevance. Substantiating these animal studies, elevated plasma ET-1 and increased mean pulmonary arterial pressure have been measured in children exposed to high pollutant levels in Mexico City . While lung injury may not be a prerequisite for pathway activation, existing injury and inflammation can exacerbate or alter the response to inhaled contaminants [9, 19]. Interaction of the physiological state of the lungs and the physicochemical properties of inhaled particles likely involves multiple pathways, and such complexity is difficult to investigate by conventional means. Genomic approaches can contribute to the elucidation of mechanisms by providing context for defined toxicant-related changes and by identifying novel toxicant-responsive pathways.
Few studies have examined responses of the pulmonary transcriptome to urban PM. In an intratracheal instillation study, lungs of spontaneously hypertensive male rats exposed to the urban PM preparation EHC-93 exhibited differential expression of 132 genes with fold-changes of at least 1.5-fold 2–6 h post-exposure, including genes involved in oxidative stress, inflammatory response, and transcription, with fewer genes affected 15–21 h and 24–40 h post-exposure . As samples were pooled for each time point it was not possible to assess the statistical significance of these effects. A subchronic inhalation exposure study of ApoE/low-density-lipoprotein double knockout mice exposed to concentrated ambient PM 6 h/day, 5 d/wk, for 4 mo did not identify any differentially expressed genes in lung tissue 3–4 d after the final exposure according to statistical analysis of results, although 95 probes were up- or down-regulated by at least 1.5-fold . A macroarray study with 1176 probes found little differential expression (9 genes) 3 d after intratracheal instillation of 10 mg Cardiff PM, and there was poor agreement with PCR data . To date, no microarray study has examined immediate or early effects of inhaled urban PM on pulmonary gene expression. Because of the rapid onset of cardiovascular effects detected in epidemiologic work [26, 27], and the early response of endothelin system genes in previous targeted gene expression studies [19, 28], this is an important window for investigation.
In the present study we used a toxicogenomic approach to investigate how chronic inflammation modifies the transcriptional response of the lungs to inhaled urban PM. We hypothesized that a transcriptome-wide approach would allow detection of effects not investigated by conventional means. Furthermore, we hypothesized that chronic inflammation would amplify the pulmonary transcriptional response to inhaled PM. Transgenic mice with constitutive pulmonary expression of tumour necrosis factor (TNF)-α under the control of the surfactant protein (SP)-C promoter develop lung inflammation, characterized by influx of macrophages, lymphocytes, and neutrophils, alveolar disruption and airspace enlargement [29, 30]. The resulting lung pathology has been well-characterized through imaging, histology, and functional analysis [29–32], and the mice have been used for assessment of effects of repeated exposure to PM and ozone . Using commercial microarrays we assessed expression of roughly 21,000 transcripts in TNF mice and their wildtype (WT) littermates exposed by inhalation to urban PM or air. In this study, the combination of a relatively large sample size (n = 5/group, 40 arrays total), paired with a standard reference design controlling for confounding variables, provided a robust approach to examine transcriptional profiles associated with PM inhalation in normal and inflamed lungs. Our objectives were: 1) to evaluate the utility of microarrays to investigate pollutant-phenotype interactions at a relevant internal dose of PM; 2) to identify genes and pathways that respond to inhaled PM; 3) to investigate how inflammation modifies this response; and 4) to identify biomarkers of exposure that are robust to the physiological status of the lungs.
Lung histology and differential cell counts
Confirmation of PM effects on key biological pathways
Microarray processing and quality control
Number of differentially expressed probes for all pair-wise analyses as determined by MAANOVA or fold-change cut-off.
Significant probes (# above background)a
Probes with > 1.5 FCb
WT,0 h (Air vs. EHC)
TNF, 0 h (Air vs. EHC)
WT, 24 h (Air vs. EHC)
TNF, 24 h (Air vs. EHC)
WT, Air (0 vs. 24 h)
TNF, Air (0 vs. 24 h)
WT, EHC (0 vs. 24 h)
TNF, EHC (0 vs. 24 h)
0 h, Air (WT vs. TNF)
24 h, Air (WT vs. TNF)
0 h, EHC (WT vs. TNF)
24 h, EHC (WT vs. TNF)
Functional analysis of Genotype effects
Functional analysis of probes that differ according to factor Genotype.*
Number of genes
% of total
immune system process
4.19 × 10-27
antigen processing and presentation
1.61 × 10-15
4.48 × 10-9
response to wounding
4.55 × 10-9
6.90 × 10-9
3.19 × 10-9
response to external stimulus
3.48 × 10-8
cytoskeleton organization and biogenesis
6.72 × 10-8
hemopoietic or lymphoid organ development
5.14 × 10-7
1.28 × 10-6
4.19 × 10-6
mast cell activation
1.45 × 10-5
lysosome organization and biogenesis
2.83 × 10-5
8.12 × 10-5
Functional analysis of Time effects
Functional analysis of probes that differ according to factor Time.*
Number of genes
% of total
response to biotic stimulus
3.86 × 10-6
5.78 × 10-6
2.23 × 10-5
5.59 × 10-4
intracellular signalling cascade
5.94 × 10-4
protein amino acid methylation
nitrogen compound catabolism
response to stress
PCR validation of microarray data
Interestingly, although metallothionein-II was selected for PCR validation based on the Time effect, with no significant PM effect identified by the array analysis, there was a statistically significant increase in the levels of this transcript after particle inhalation when measured by PCR (Treatment main effect, p = 0.03; Figure 6a). In contrast, of the 4 probes identified by MAANOVA as exhibiting PM-dependent changes, real-time PCR failed to confirm significant differential expression in the three probes for which primers were available (RIKEN cDNA 4921531P07 gene, zinc finger protein 217, CD9 antigen), suggesting that these results were false-positives (data not shown).
Gene set analysis
Gene set analysis of probes selected with relaxed statistical stringency may detect subtle effects not identified by stringent statistical analysis, and so this approach was used to search for PM effects. Screening for gene ontology term enrichment in gene lists of the top 50 genes according to unadjusted p-value for all pairwise comparisons failed to identify coherent functional grouping of differentially expressed genes for factor Treatment [see Additional file 2]. In contrast, functional analysis of Time effects revealed enrichment of terms such as immune response, apoptosis, signalling pathways, blood vessel development and lipid metabolism [see Additional file 3], consistent with the analysis using all probes with FDR-adjusted P < 0.05 (Table 3), indicating that the lists were sufficiently large and representative of effects to enable identification of biological processes.
Fold change (FC)
A fold-change cut-off approach can be used as a non-statistical method to determine a list of candidate genes, with the most likely candidates being probes appearing in multiple lists across the treatment groups. Probe lists generated by filtering on FC = 1.5 resulted in the selection of 1327–3024 (Genotype), 46–1055 (Time), and 0–572 (Treatment) probes (intersection-union; Table 1). Again, no probes were common to all air vs. PM comparisons. However, two probes were increased by at least 1.5-fold immediately after PM exposure in both WT and TNF mice, and were therefore chosen as top candidates for PCR validation: A_51_P279693 (NM_009992; CYP1A1; FC = 1.5 in WT, 1.9 in TNF; previously validated by real-time PCR) and A_51_P442894 (AK048310; RIKEN cDNA 1700055N04 gene, aldehyde dehydrogenase 3 family, member B2; Aldh3b2; FC = 1.5 in WT, 1.6 in TNF). While altered CYP1A1 expression was confirmed (Figure 3c), real-time PCR analysis indicated that the Treatment effect on Aldh3b2 mRNA was not significant (data not shown).
A reference design was used in the present study for its simplicity, flexibility, and ease of data analysis for multi-factorial experiments . It is, however, well-recognized that a reference design may add variability as a result of dividing by a potentially uncertain measurement, which could reduce the power to detect treatment effects or contribute to the generation of spurious results. Focussing on PM effects, a one-colour analysis was conducted within each genotype to verify that use of a common reference did not impede detection of effects. Only three unidentified proteins in WT mice (NM_025901, AK029835, AK036646), and one probe in TNF mice (D83146, Mus musculus mRNA for Six5_ partial cds; FC = 1.25) were differentially expressed according to Treatment. None of these genes were common to both one and two-colour analyses. The lack of PM effects in the one-colour analysis is in line with the two-colour analysis and indicates that removal of the common reference did not improve detection of statistically significant PM effects.
Despite evidence of transcription factor activation in lung cells after PM exposure , few studies have examined effects of PM on pulmonary gene expression (reviewed in ). We hypothesized that by evaluating differential gene expression at the global transcript level, a microarray approach would enable identification of PM-activated pathways not investigated by conventional approaches. In the present work, the impact of PM inhalation on transcript profiles in healthy mice and in mice with chronic lung inflammation was assessed immediately and 24 h post-exposure to detect effects that could be related to the rapid response indicated by epidemiologic studies [26, 27]. Using real-time PCR, we confirmed positive effects of PM inhalation on three distinct pathways: 1) up-regulation of CYP1A1, relevant to the presence of combustion by-products in the inhaled urban PM and activation of a procarcinogenic pathway; 2) up-regulation of metallothionein-II, relevant to the presence of toxic metals in the inhaled particles and activation of pathways to protect against metal toxicity and oxidative stress; and 3) up-regulation of ET-1, relevant to endothelial dysfunction generated by deposited toxicants and downstream cardiovascular effects. The activation of CYP1A1 and metallothionein-II is in line with our previous report relating chemical composition with the biological potency of particles . These responses establish delivery of a biologically-effective internal dose in our study despite the absence of visible impacts on lung lavage cytology, and are consistent with previous work showing that pronounced inflammation and acute lung injury are not prerequisites for induction of physiologically-relevant pathways by PM [9, 14]. The early time points assessed reduce the likelihood that interpretation of molecular data would be confounded by infiltrating inflammatory cells, which often peak 24–48 h after acute injury. Remarkably, although microarray analysis identified significant differential gene expression attributable to genotype and to duration of recovery post-exposure, no significant PM effects were found. Since array quality was internally validated by detection of positive effects, the data indicate that any effects of inhaled particles were below the threshold of detection.
As the site of gas exchange, the lungs are continuously exposed to airborne contaminants, and possess innate defence mechanisms that facilitate clearance of inhaled particles without provoking excessive inflammation that could lead to tissue damage. Chronic inflammation can impair lung defences by increasing the permeability of the epithelial barrier and slowing mucocilliary clearance, facilitating passage of inhaled particles into the tissue . Depletion of anti-oxidant defences, activation of redox-sensitive transcription factors, and chromatin remodelling in the inflamed lungs could all alter the transcriptional response to PM . Microarray data confirmed that the pulmonary transcriptome of TNF mice was substantially different from that of WT littermates, with activation of inflammatory and immune response pathways consistent with the lung phenotype [29, 30]. Clearly, differential expression detected in this model reflects import of message due to changes in cellularity, such as influx of inflammatory cells, in addition to strict transcriptional activation by TNF-α. Surprisingly, this persistent inflammatory state did not measurably amplify or alter the response to inhaled particles. It is unclear whether the lack of response is due to functional transcription-independent defences in TNF animals, or to insufficient sensitivity to detect subtle PM effects by microarray at the whole lung level. Despite the higher numbers of inflammatory cells in the lungs of TNF mice, levels of ρ-tyrosine, a product of hydroxylation of phenylalanine indicative of oxidative stress, are not significantly higher compared to WT littermates, nor increased in either genotype after repeated exposure to PM and ozone . This suggests that anti-oxidant defences may still be functional in the TNF model. After adjusting for the decreased alveolar surface area in TNF mice , PM effects on CYP1A1 and ET-1 were remarkably similar to those observed in WT animals. These results indicate that the altered expression of specific genes relevant to adverse PM effects was robust to the physiological status of the lungs, despite significant genotype differences in the basal levels of the transcripts. The lower basal metallothionein-II mRNA level and apparent reduced response in TNF mice could reduce scavenging of toxicologically-relevant metals [34, 39].
The data indicate that effects of urban particles on the lung transcriptome following a single acute exposure by inhalation are subtle. The EHC-6802 urban particle preparation used in this study is an environmental material collected in Ottawa, Canada, and is equivalent to the well-characterized EHC-93 particle preparation [9, 18, 34]. The inherent toxicity of EHC-93 particles has been demonstrated in vivo using intratracheal instillation [39, 40], and inhalation studies have revealed endothelin system activation in the absence of pronounced inflammatory changes [14, 18, 19]. Activation of biologically-relevant pathways by EHC-6802 is confirmed here by real-time PCR. Thus, the lack of overt effects detected by microarray analysis is not likely due to an absence of factors associated with health effects of particles in this particle preparation. It remains possible that repeated exposure to elevated PM levels could induce measurable effects that were not apparent after the single exposure used here. In the only other microarray study published in the literature that has, to our knowledge, examined pulmonary effects of urban PM in an inhalation model, subchronic exposure to concentrated ambient PM yielded little evidence of differential expression in the whole lung . Note that in this previous study expression was measured 3–4 days after the last exposure, and may have missed subsiding transcriptional events. In the present study, tissues were recovered immediately and 24 h post-exposure to capture rapid effects of exposure on gene expression. Statistical analysis of array data was conducted at varying levels of stringency such that the absence of detected effects was unlikely to be due to overly stringent filters. Moreover, sample and array quality were confirmed by detection of Time and Genotype effects in line with expected responses. Compared with the coherence and PCR validation of Time and Genotype effects, the lack of consistent Treatment effects suggests that the few probes identified as responsive to PM were false-positives, and that the absence of PM effects was not due to the analytical approach chosen. Rather, the data suggest that overt, coordinated transcriptional responses, such as would be expected in an inflammatory signalling cascade, may not be a major part of the pulmonary response to inhaled urban particles in this model.
Internal dose and method of delivery should be considered when comparing results across studies. Previous micro- and macroarray studies examining PM effects on pulmonary gene expression include exposure to saline suspension of urban PM [23, 25], concentrated ambient particles , diesel exhaust particles , residual oil fly ash , and ultrafine carbon particles . Studies using an intratracheal instillation approach employed particle doses ranging from 1.25–10 mg, delivered as a bolus in 0.5 mL saline [23, 25, 41, 42]. In the present study, the dose of EHC-6802 within the respiratory compartment was estimated to be 16 μg (32 ng/cm2 for WT mice and 64 ng/cm2 for TNF mice). The estimated doses were only 25-fold (WT mice) or 50-fold (TNF mice) higher than the modelled deposition for a plausible 24 h human exposure scenario (1.3 ng/cm2 ), but more importantly, were 100–1000 times lower than the dose of particles delivered in the intratracheal instillation studies cited above. Exposure by inhalation, a lower internal dose of PM, and a slower dose rate are less likely to provoke an inflammatory response that could impact transcript levels, and indeed did not cause a significant influx of inflammatory cells, in agreement with our previous work in rats [9, 14]. Our inhalation approach should therefore be more sensitive to changes in parenchymal gene expression and avoid artifacts of an overly high internal dose of particles. For ethical reasons, nose-only exposures should be kept to a minimum duration, and therefore the dose-rate in our study was considerably higher than for a 24 h environmental exposure. Nevertheless, the pulmonary deposition of PM in our inhalation study is directly relevant to the human experience, once a number of reasonable uncertainty factors are considered. These include possible decreased potency of EHC-6802 by comparison to fresh particles, interspecies differences in sensitivity to air pollutants, and the heightened sensitivity within a subset of the human population, such as the known increased adverse risk of individuals with COPD, congestive heart failure, and atherosclerosis [4, 5].
While training in nose-only exposure tubes in the days preceding exposure is aimed at reducing animal stress, exposure altered the expression of a number of genes in the lungs, even in animals exposed to filtered air alone. These changes show some similarity with effects of stress on other organs, such as the spleen , and likely relate to systemic effects of stress. Since both air and PM-exposed animals were similarly exposed to the inhalation system in the present study, such changes do not confound interpretation of results. Moreover, within each genotype, air and particle exposed animals exhibited similar Time effects. It is, however, not known to what degree these changes may impede detection of PM effects in the subset of genes affected. It is noteworthy that although metallothionein-II was selected for PCR validation based on its response to Time, with no significant PM effect according to array analysis, expression was significantly increased by PM when measured by the more sensitive real-time PCR method. The higher levels of metallothionein-II mRNA in PM-exposed animals indicate that, at least for this gene, effects of exposure to PM and stress were additive, consistent with work showing that the mechanism of metallothionein-II induction by restraint stress is distinct from that by chemical stress .
Toxicogenomic analysis in vivo at doses relevant to human health can be problematic when subtle changes are expected to arise. The present study shows that inhaled particles at an internal dose relevant to the human experience may not cause significant changes in gene expression detectable at the level of the whole lung using a moderate number of microarrays. Moreover, small but physiologically-relevant changes in gene expression, such as activation of endothelin and xenobiotic-metabolism pathways, may not be identified by microarray analysis. To detect statistically significant differential expression of a few probes from among 21,000 requires relatively large effect sizes, or a very high number of arrays. Using our data and taking a conservative approach of 90% confidence, we calculated that 60–70 arrays would be required to detect 1.5-fold PM effects in animals of a given genotype, assuming 100 differentially expressed genes with FDR-adjusted p < 0.05 (Figure 7). Experiments requiring this number of arrays may not be practical for many laboratories. Using fewer arrays and an arbitrary fold cut-off rather than statistical analysis risks inclusion of false-positive results and/or exclusion of small (< 1.5-fold) but physiologically relevant changes. This type of approach would require extensive RT-PCR on multiple biological replicates to validate any findings. Alternatively, taking a more traditional hypothesis-driven approach through selection of probe subsets diminishes the likelihood of discovering unexpected effects. A third possibility, increasing dose to generate larger effect sizes, may create artifacts such as overload of clearance pathways. In light of our observations, there are several approaches that might be considered for the assessment of pulmonary effects of contaminant exposure. Since experiments examining expression in whole lung tissue may be insensitive to focal responses in target cells, an anatomically-biased or cell-type specific approach such as laser-capture microdissection might uncover local impacts of PM inhalation using fewer arrays. This should be fairly straightforward for the study of effects in WT animals. However, complex structural changes and intensive inflammatory infiltration will complicate analysis in models such as the TNF mice. Secondly, appropriate pooling of animals may reduce the impact of biological variability, provided a sufficient number of arrays hybridized with independently-pooled RNA are used to permit statistical analysis [46, 47]. However, it is likely that a large number of arrays would still be required to detect subtle effects. Lastly, a combination of transcriptome screening and conventional toxicologic analyses may be appropriate to examine subtle effects of exposure.
The present study verified increased expression of ET-1, CYP1A1, and metallothionein-II in both wildtype and TNF mice exposed by the nose-only inhalation route, consistent with in vitro  and in vivo [19, 35] analyses, and in line with our understanding of the impact of PM on cardiovascular, toxicant metabolism, and metal-responsive pathways. These effects were robust to the physiological status of the lungs, and remarkably, did not differ greatly between genotypes. Although the particles clearly induce pathways relevant to endothelial dysfunction (ET-1), carcinogenicity (CYP1A1), and metal toxicity (metallothionein-II), the exposures did not measurably initiate inflammation (wildtype) or exacerbate existing inflammation (TNF) in the lungs as measured by lavage cytology and global pulmonary gene expression. Given the increased susceptibility of individuals with respiratory disease [4–7], it may be that susceptibility arises as a result of the interaction of primary effects with host factors, in addition to any enhancement or alteration of these primary effects. For example, individuals with existing endothelial dysfunction and an ineffective compensation for the vasopressor effect of ET-1 may respond adversely to a PM-induced increase of circulating ET-1, while healthy individuals may compensate with release of nitric oxide and prostacyclin and not exhibit hemodynamic changes in response to the same relative increase of circulating ET-1. While microarray analysis is a valuable approach for genome-wide screening of pronounced effects, such as characterization of transcriptional changes resulting from TNF over-expression, it may not be particularly useful for dissecting pathways at the level of the whole lung after delivery of an environmentally-relevant internal dose of contaminants. Increased power, site-specific analysis, or more sensitive animal models may be required to investigate effects of acute particle exposure on pulmonary gene expression. In addition, focussed approaches targeting key endpoints should be considered to complement genome-wide transcript profiling. Our data support the hypothesis that adverse health effects of acute exposure to urban PM may be dominated by physiological response cascades rather than widespread changes in the expression of genes escaping homeostatic controls. Nevertheless, transcriptional activation of certain key pathways, such as endothelin synthesis and xenobiotic transformation, may be relevant to these health effects.
SP-C/TNF-α mice were provided by Dr. R. J. Mason (National Jewish Medical and Research Centre, Denver, CO, USA) and crossed with C57BL/6 mice (Charles River Laboratories, St. Constant, QC, Canada). Note that the SP-C/TNF-α line has been maintained as a heterozygous line by repeated backcrossing to C57BL/6 mice for over a decade [29, 30, 33]. Male transgenic TNF mice and their WT littermates were genotyped by PCR analysis of genomic DNA . Animals were exposed as three cohorts, each approximately 3 months apart, and were 131 ± 5 days old and 28.9 ± 2.8 g (WT: 30.8 ± 2.5 g; TNF: 27.1 ± 1.7) at the time of exposure. Mice were housed in individual plexiglass cages on wood-chip bedding under HEPA-filtered air and held to a 12 h dark/light cycle. Food and water were provided ad libitum. All experimental protocols were reviewed and approved by the Animal Care Committee of Health Canada as set forth in the Guidelines of the Canadian Council on Animal Care.
The urban particles EHC-6802 consist of a blend of total suspended PM recovered from filters of the single-pass air-purification system at the Environmental Health Centre (Tunney's Pasture, Ottawa, ON, Canada) in 1996, 1998, 2000, and 2002. Particles were mechanically sieved using a 36 μm mesh filter, and combined in equal proportions. EHC-6802 was recovered at the same site and in the same manner as the urban particles EHC-93 [9, 18, 34], and chemical characterization of the EHC-6802 material has confirmed that it is equivalent to the EHC-93 material (data not shown).
Activation of the aryl hydrocarbon receptor (AhR) by EHC-6802
The AhR-based luciferase reporter cell line H1L1.1c2 (kindly provided by Dr. M. S. Denison, University of California, Davis, CA, USA) was used to assess AhR-activation by the EHC-6802 particles. H1L1.1c2 cells were exposed in triplicate to EHC-6802 (1, 10, 100 μg/cm2), 10 nM benzo(a)pyrene (positive control), and TiO2 (negative control) in 24-well plates (2 cm2/well) for 4 h, and luciferase activity was measured as previously described .
Mice were trained in nose-only exposure tubes over 4 consecutive days and then exposed for 4 h to clean air or 42 mg/m3 EHC-6802 by nose-only exposure as described previously [9, 19]. Particle concentration was monitored during each exposure at the inhalation ports by isokinetic sampling using 0.2 μm Teflon filters. Filter weight was divided by the sampling volume to provide a direct estimate of the time-weighted average particle concentration. Particle counts and size measurements were performed at the inhalation ports (optical size range of 0.3–10 μm; Lasair Model 301; Particle Measuring Systems, Boulder, CO, USA). Aerodynamic size characteristics were determined by gravimetric cascade impactor analysis of isokinetic samples at the inhalation ports (seven-stage Mercer cascade impactor, 1 L/min, 0.2–5.1 μm effective cut-off diameter; Intox, Albuquerque, NM, USA) or the chamber exhaust (seven-stage Mercer cascade impactor, 10 L/min, 0.2 to 9.8 μm effective cut-off diameter; Intox). Multimodal particle size distribution analysis  and pulmonary deposition modelling were performed as previously described . Deposition model assumptions for mice were strict nasal breathing, a minute ventilation of 34 mL/min, an alveolar surface area of 500 cm2, and a tracheobronchial surface area of 3.5 cm2 (default values for the Regional Deposited Dose Ratio RDDR2 modelling software, US EPA).
Collection of biological samples
Mice were anaesthetized by administration of sodium pentobarbital (60 mg/kg, i.p.) and euthanized by exsanguination immediately or 24 h after termination of exposure (n = 5 animals of each genotype/treatment/time). Lungs were washed by bronchoalveolar lavage with warm saline (37°C) at 30 mL/kg body weight, then flash frozen in liquid nitrogen and stored at -80°C until isolation of RNA. Cells recovered by centrifugation (1500 rpm for 10 min at 4°C) were counted in a Coulter Multisizer II (Coulter Canada, Ville St. Laurent, QC, Canada) and differential cell counts were obtained from cytospin preparations. Lungs were homogenized in TRIzol and RNA was isolated according to manufacturer's instructions (Invitrogen Canada Inc., Burlington, ON, Canada), and further purified using RNeasy Mini Kits (Qiagen Inc., Mississauga, ON, Canada). RNA was quantified using the RiboGreen RNA Quantitation Reagent and Kit (Molecular Probes, Eugene, OR, USA), and quality was confirmed using the Agilent 2100 Bioanalyzer and RNA 6000 NanoLab Chip Kit (Agilent Technologies Canada Inc., Mississauga, ON, Canada). For histology, lungs of unexposed animals (n = 3/genotype) were inflated in chest at 25 cm H20 static pressure by intratracheal instillation of 4% paraformaldehyde in phosphate-buffered saline. The lungs were excised, immersed in fixative, and stored at 4°C. Tissue blocks were dehydrated in ethanol and embedded in paraffin. Sections (0.75 μm) were stained with hematoxylin and eosin.
Real-time reverse transcription-polymerase chain reaction (RT-PCR)
Total RNA was reverse transcribed using MuLV reverse transcriptase and random hexamers according to manufacturer's instructions (Applied Biosystems, Mississauga, ON, Canada). Twenty ng of cDNA was combined with Quantitect primer assays and real-time PCR supermix (Qiagen) and run according to the product protocol on the iCycler thermal cycler (Bio-Rad Laboratories Canada Ltd., Mississauga, ON, Canada) using 50 μL volumes . Post-run melt curves were routinely inspected to verify product purity. Expression was calculated relative to reference genes (β-actin and TATA binding protein) using the comparative threshold cycle method . Statistical significance was assessed by three-way ANOVA with Treatment (Air, EHC), Genotype (WT, TNF), and Time (0, 24 h) as factors, followed by the Holm-Sidak multiple comparison procedure to elucidate the pattern of significant effects (α = 0.05; Sigma Stat 3.0, SPSS Inc., Chicago, IL, USA). PM effects were represented graphically as fold-change relative to the time and genotype-matched air-exposed controls.
An unbalanced block factorial design  was used for the factors Treatment (0, 50 mg/m3), Time (0, 24 h), and Genotype (WT, TNF) and blocked for the nuisance factors Date of exposure and Day of hybridization . Five biological replicates per condition were used for a total of 40 microarrays. Individual 2.5 μg aliquots of RNA from each sample were amplified and labelled using the Low RNA Input Fluorescent Linear Amplification Kit (Agilent). Agilent Mouse G4121A Microarrays (containing approximately 21,000 probes) were hybridized with 5 μg Cy5-labelled lung RNA and 5 μg Cy3-labelled Universal Mouse Reference RNA (Stratagene, CA, USA), used as a common reference on all arrays . Arrays were incubated overnight at 60°C in Agilent hybridization solution and washed according to manufacturer's instructions. Arrays were scanned using a ScanArray Express (Perkin-Elmer Life Sciences, Woodbridge, ON, Canada), and data were acquired with ImaGene 5.5 (BioDiscovery, CA, USA).
Statistical analysis of microarray data
The background signal for each array was determined using the negative control (-)3xSLv1 probe. Spots with median signal intensities within three standard deviations of the (-)3xSLv1 probe mean were flagged as absent. Lowess normalization  was performed using the SAS/STAT software, Version 8.2 of the SAS System for Windows (1999–2001 SAS Institute Inc., Cary, NC, USA). Data displays produced in R, including ratio intensity plots for the raw and normalized data, comparison boxplots, heatmaps, dendrograms, and volcano plots, were used for the assessment of data quality and for visual comparison of the impact of each factor . The logarithm base 2 relative intensities were used for subsequent analyses.
The present experiment was designed to evaluate pollutant-phenotype interactions, and as such required the analysis of multiple factors and interactions among factors. A factorial design is a suitable statistical framework to maximize resources, and permits analysis of multiple factors with the same precision as if the experiment were designed to examine one factor, assuming common variance. The Microarray Analysis of Variance (MAANOVA) library  in R was used to identify differentially expressed genes. The statistical model included the main effects Treatment, Time, and Genotype, the three-way interaction, and all two-way interactions. The Fs statistic , a shrinkage estimator for the gene-specific variance components, was used, and p-values for all statistical tests were estimated using the permutation method (1000 permutations with residual shuffling). P-values were then adjusted for multiple comparisons using the false discovery rate (FDR) approach . Group means for the fold-change calculation were based on the adjusted relative intensity for each gene after subtracting estimated Date of exposure and Day of hybridization effects from the normalized ratio. Genes with FDR-adjusted p < 0.05 and present on at least 4 of 5 arrays in the group showing a significant effect were considered differentially expressed. Gene ontology enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) , available at http://david.abcc.ncifcrf.gov/home.jsp. Agilent GeneSpring GX was used for filtering on fold-change. Microarray data were deposited in the NCBI Gene Expression Omnibus (Accession # GSE11037; http://www.ncbi.nlm.nih.gov/projects/geo).
Power scenarios were investigated using the samr.assess.samplesize function in the SAMR library  to determine the sample size (n animals, 1 array/animal) required to detect 1.25–2.5 fold changes in expression when comparing control and exposed animals for each genotype at each time point. The required sample size was identified through inspection of plots generated from this application, and corresponded to the value at which the 90th percentile for FDR was less than 0.05. Results for the four analyses were averaged for each fold change, and the number of arrays required was plotted relative to the minimum detectable fold-change.
The authors are grateful to Marcelle Phaneuf for animal genotyping and project planning, Josée Guénette and DJ MacIntyre for conducting the inhalation exposures, Pat Goegan, Erica Blais, Lorraine Casavant, and Yunus Siddiqui for technical assistance, and Lynn Berndt for microarray advice. Errol Thomson was recipient of a scholarship from the Natural Sciences and Engineering Research Council of Canada. This work was supported by Health Canada (SEP 4320010; RV) and the Genomics R&D Initiative (RV, CY).
- Dockery DW, Pope CA III, Xu X, Spengler JD, Ware JH, Fay ME, Ferris BG Jr, Speizer FE: An association between air pollution and mortality in six U.S. cities. N Engl J Med 1993, 329: 1753–1759.PubMedView ArticleGoogle Scholar
- Burnett RT, Brook J, Dann T, Delocla C, Philips O, Cakmak S, Vincent R, Goldberg MS, Krewski D: Association between particulate- and gas-phase components of urban air pollution and daily mortality in eight Canadian cities. Inhal Toxicol 2000,12(Suppl 4):15–39.PubMedView ArticleGoogle Scholar
- Pope CA III, Burnett RT, Thurston GD, Thun MJ, Calle EE, Krewski D, Godleski JJ: Cardiovascular mortality and long-term exposure to particulate air pollution: epidemiological evidence of general pathophysiological pathways of disease. Circulation 2004, 109: 71–77.PubMedView ArticleGoogle Scholar
- Zanobetti A, Schwartz J, Gold D: Are there sensitive subgroups for the effects of airborne particles? Environ Health Perspect 2000, 108: 841–845.PubMed CentralPubMedView ArticleGoogle Scholar
- Goldberg MS, Burnett RT, Bailar JC III, Tamblyn R, Ernst P, Flegel K, Brook J, Bonvalot Y, Singh R, Valois MF, Vincent R: Identification of persons with cardiorespiratory conditions who are at risk of dying from the acute effects of ambient air particles. Environ Health Perspect 2001,109(Suppl 4):487–494.PubMed CentralPubMedView ArticleGoogle Scholar
- Sunyer J: Urban air pollution and chronic obstructive pulmonary disease: a review. Eur Respir J 2001, 17: 1024–1033.PubMedView ArticleGoogle Scholar
- Peel JL, Metzger KB, Klein M, Flanders WD, Mulholland JA, Tolbert PE: Ambient air pollution and cardiovascular emergency department visits in potentially sensitive groups. Am J Epidemiol 2007, 165: 625–633.PubMedView ArticleGoogle Scholar
- Donaldson K, Stone V, Seaton A, MacNee W: Ambient particle inhalation and the cardiovascular system: potential mechanisms. Environ Health Perspect 2001,109(Suppl 4):523–527.PubMed CentralPubMedView ArticleGoogle Scholar
- Vincent R, Bjarnason SG, Adamson IY, Hedgecock C, Kumarathasan P, Guenette J, Potvin M, Goegan P, Bouthillier L: Acute pulmonary toxicity of urban particulate matter and ozone. Am J Pathol 1997, 151: 1563–1570.PubMed CentralPubMedGoogle Scholar
- Adamson IY, Hedgecock C: Patterns of particle deposition and retention after instillation to mouse lung during acute injury and fibrotic repair. Exp Lung Res 1995, 21: 695–709.PubMedView ArticleGoogle Scholar
- Oberdorster G, Ferin J, Lehnert BE: Correlation between particle size, in vivo particle persistence, and lung injury. Environ Health Perspect 1994,102(Suppl 5):173–179.PubMed CentralPubMedView ArticleGoogle Scholar
- Adamson IY, Prieditis H: Silica deposition in the lung during epithelial injury potentiates fibrosis and increases particle translocation to lymph nodes. Exp Lung Res 1998, 24: 293–306.PubMedView ArticleGoogle Scholar
- Rahman I, Adcock IM: Oxidative stress and redox regulation of lung inflammation in COPD. Eur Respir J 2006, 28: 219–242.PubMedView ArticleGoogle Scholar
- Bouthillier L, Vincent R, Goegan P, Adamson IY, Bjarnason S, Stewart M, Guenette J, Potvin M, Kumarathasan P: Acute effects of inhaled urban particles and ozone: lung morphology, macrophage activity, and plasma endothelin-1. Am J Pathol 1998, 153: 1873–1884.PubMed CentralPubMedView ArticleGoogle Scholar
- Brook RD, Brook JR, Urch B, Vincent R, Rajagopalan S, Silverman F: Inhalation of fine particulate air pollution and ozone causes acute arterial vasoconstriction in healthy adults. Circulation 2002, 105: 1534–1536.PubMedView ArticleGoogle Scholar
- Samet JM, Silbajoris R, Huang T, Jaspers I: Transcription factor activation following exposure of an intact lung preparation to metallic particulate matter. Environ Health Perspect 2002, 110: 985–990.PubMed CentralPubMedView ArticleGoogle Scholar
- Sen B, Mahadevan B, DeMarini DM: Transcriptional responses to complex mixtures: a review. Mutat Res 2007, 636: 144–177.PubMedView ArticleGoogle Scholar
- Vincent R, Kumarathasan P, Goegan P, Bjarnason SG, Guenette J, Berube D, Adamson IY, Desjardins S, Burnett RT, Miller FJ, Battistini B: Inhalation toxicology of urban ambient particulate matter: acute cardiovascular effects in rats. Res Rep Health Eff Inst 2001, 5–54.Google Scholar
- Thomson E, Kumarathasan P, Goegan P, Aubin RA, Vincent R: Differential regulation of the lung endothelin system by urban particulate matter and ozone. Toxicol Sci 2005, 88: 103–113.PubMedView ArticleGoogle Scholar
- Thomson E, Kumarathasan P, Vincent R: Pulmonary expression of preproET-1 and preproET-3 mRNAs is altered reciprocally in rats after inhalation of air pollutants. Exp Biol Med (Maywood) 2006, 231: 979–984.Google Scholar
- Galatius-Jensen S, Wroblewski H, Emmeluth C, Bie P, Haunso S, Kastrup J: Plasma endothelin in congestive heart failure: a predictor of cardiac death? J Card Fail 1996, 2: 71–76.PubMedView ArticleGoogle Scholar
- Calderon-Garciduenas L, Vincent R, Mora-Tiscareno A, Franco-Lira M, Henriquez-Roldan C, Barragan-Mejia G, Garrido-Garcia L, Camacho-Reyes L, Valencia-Salazar G, Paredes R, Romero L, Osnaya H, Villarreal-Calderon R, Torres-Jardon R, Hazucha MJ, Reed W: Elevated plasma endothelin-1 and pulmonary arterial pressure in children exposed to air pollution. Environ Health Perspect 2007, 115: 1248–1253.PubMed CentralPubMedView ArticleGoogle Scholar
- Kooter I, Pennings J, Opperhuizen A, Cassee F: Gene expression pattern in spontaneously hypertensive rats exposed to urban particulate matter (EHC-93). Inhal Toxicol 2005, 17: 53–65.PubMedView ArticleGoogle Scholar
- Gunnison A, Chen LC: Effects of subchronic exposures to concentrated ambient particles (CAPs) in mice. VI. Gene expression in heart and lung tissue. Inhal Toxicol 2005, 17: 225–233.PubMedView ArticleGoogle Scholar
- Wise H, Balharry D, Reynolds LJ, Sexton K, Richards RJ: Conventional and toxicogenomic assessment of the acute pulmonary damage induced by the instillation of Cardiff PM10 into the rat lung. Sci Total Environ 2006, 360: 60–67.PubMedView ArticleGoogle Scholar
- Gold DR, Litonjua A, Schwartz J, Lovett E, Larson A, Nearing B, Allen G, Verrier M, Cherry R, Verrier R: Ambient pollution and heart rate variability. Circulation 2000, 101: 1267–1273.PubMedView ArticleGoogle Scholar
- Peters A, Dockery DW, Muller JE, Mittleman MA: Increased particulate air pollution and the triggering of myocardial infarction. Circulation 2001, 103: 2810–2815.PubMedView ArticleGoogle Scholar
- Thomson E, Goegan P, Kumarathasan P, Vincent R: Air pollutants increase gene expression of the vasoconstrictor endothelin-1 in the lungs. Biochim Biophys Acta 2004, 1689: 75–82.PubMedView ArticleGoogle Scholar
- Miyazaki Y, Araki K, Vesin C, Garcia I, Kapanci Y, Whitsett JA, Piguet PF, Vassalli P: Expression of a tumor necrosis factor-alpha transgene in murine lung causes lymphocytic and fibrosing alveolitis. A mouse model of progressive pulmonary fibrosis. J Clin Invest 1995, 96: 250–259.PubMed CentralPubMedView ArticleGoogle Scholar
- Fujita M, Shannon JM, Irvin CG, Fagan KA, Cool C, Augustin A, Mason RJ: Overexpression of tumor necrosis factor-alpha produces an increase in lung volumes and pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol 2001, 280: L39-L49.PubMedGoogle Scholar
- Fujita M, Mason RJ, Cool C, Shannon JM, Hara N, Fagan KA: Pulmonary hypertension in TNF-alpha-overexpressing mice is associated with decreased VEGF gene expression. J Appl Physiol 2002, 93: 2162–2170.PubMedView ArticleGoogle Scholar
- Lundblad LK, Thompson-Figueroa J, Leclair T, Sullivan MJ, Poynter ME, Irvin CG, Bates JH: Tumor necrosis factor-alpha overexpression in lung disease: a single cause behind a complex phenotype. Am J Respir Crit Care Med 2005, 171: 1363–1370.PubMed CentralPubMedView ArticleGoogle Scholar
- Kumarathasan P, Blais E, Goegan P, Yagminas A, Guenette J, Adamson IY, Crapo JD, Mason RJ, Vincent R: 90-day repeated inhalation exposure of surfactant Protein-C/tumor necrosis factor-alpha, (SP-C/TNF-alpha) transgenic mice to air pollutants. Int J Toxicol 2005, 24: 59–67.PubMedView ArticleGoogle Scholar
- Vincent R, Goegan P, Johnson G, Brook JR, Kumarathasan P, Bouthillier L, Burnett RT: Regulation of promoter-CAT stress genes in HepG2 cells by suspensions of particles from ambient air. Fundam Appl Toxicol 1997, 39: 18–32.PubMedView ArticleGoogle Scholar
- Takano H, Yanagisawa R, Ichinose T, Sadakane K, Inoue K, Yoshida S, Takeda K, Yoshino S, Yoshikawa T, Morita M: Lung expression of cytochrome P450 1A1 as a possible biomarker of exposure to diesel exhaust particles. Arch Toxicol 2002, 76: 146–151.PubMedView ArticleGoogle Scholar
- Morgan ET, Li-Masters T, Cheng PY: Mechanisms of cytochrome P450 regulation by inflammatory mediators. Toxicology 2002, 181–182: 207–210.PubMedView ArticleGoogle Scholar
- Guo L, Lobenhofer EK, Wang C, Shippy R, Harris SC, Zhang L, Mei N, Chen T, Herman D, Goodsaid FM, Hurban P, Phillips KL, Xu J, Deng X, Sun YA, Tong W, Dragan YP, Shi L: Rat toxicogenomic study reveals analytical consistency across microarray platforms. Nat Biotechnol 2006, 24: 1162–1169.PubMedView ArticleGoogle Scholar
- Kerr KF, Serikawa KA, Wei C, Peters MA, Bumgarner RE: What is the best reference RNA? And other questions regarding the design and analysis of two-color microarray experiments. OMICS 2007, 11: 152–165.PubMedView ArticleGoogle Scholar
- Adamson IY, Prieditis H, Hedgecock C, Vincent R: Zinc is the toxic factor in the lung response to an atmospheric particulate sample. Toxicol Appl Pharmacol 2000, 166: 111–119.PubMedView ArticleGoogle Scholar
- Adamson IY, Prieditis H, Vincent R: Pulmonary toxicity of an atmospheric particulate sample is due to the soluble fraction. Toxicol Appl Pharmacol 1999, 157: 43–50.PubMedView ArticleGoogle Scholar
- Reynolds LJ, Richards RJ: Can toxicogenomics provide information on the bioreactivity of diesel exhaust particles? Toxicology 2001, 165: 145–152.PubMedView ArticleGoogle Scholar
- Nadadur SS, Kodavanti UP: Altered gene expression profiles of rat lung in response to an emission particulate and its metal constituents. J Toxicol Environ Health A 2002, 65: 1333–1350.PubMedView ArticleGoogle Scholar
- Andre E, Stoeger T, Takenaka S, Bahnweg M, Ritter B, Karg E, Lentner B, Reinhard C, Schulz H, Wjst M: Inhalation of ultrafine carbon particles triggers biphasic pro-inflammatory response in the mouse lung. Eur Respir J 2006, 28: 275–285.PubMedView ArticleGoogle Scholar
- Yin D, Zhang Y, Stuart C, Miao J, Zhang Y, Li C, Zeng X, Hanley G, Moorman J, Yao Z, Woodruff M: Chronic restraint stress modulates expression of genes in murine spleen. J Neuroimmunol 2006, 177: 11–17.PubMedView ArticleGoogle Scholar
- Ghoshal K, Wang Y, Sheridan JF, Jacob ST: Metallothionein induction in response to restraint stress. Transcriptional control, adaptation to stress, and role of glucocorticoid. J Biol Chem 1998, 273: 27904–27910.PubMed CentralPubMedView ArticleGoogle Scholar
- Peng X, Wood CL, Blalock EM, Chen KC, Landfield PW, Stromberg AJ: Statistical implications of pooling RNA samples for microarray experiments. BMC Bioinformatics 2003, 4: 26.PubMed CentralPubMedView ArticleGoogle Scholar
- Kendziorski C, Irizarry RA, Chen KS, Haag JD, Gould MN: On the utility of pooling biological samples in microarray experiments. Proc Natl Acad Sci USA 2005, 102: 4252–4257.PubMed CentralPubMedView ArticleGoogle Scholar
- Ziccardi MH, Gardner IA, Denison MS: Application of the luciferase recombinant cell culture bioassay system for the analysis of polycyclic aromatic hydrocarbons. Environ Toxicol Chem 2002, 21: 2027–2033.PubMedView ArticleGoogle Scholar
- Hewett P, McCawley MA: A microcomputer spreadsheet technique for analyzing multimodal particle size distributions. Appl Occup Environ Hyg 1991, 6: 865–873.View ArticleGoogle Scholar
- Thomson E, Vincent R: Reagent volume and plate bias in real-time polymerase chain reaction. Anal Biochem 2005, 337: 347–350.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25: 402–408.PubMedView ArticleGoogle Scholar
- Montgomery DC: Design and Analysis of Experiments. 5th edition. John Wiley & Sons; 2001.Google Scholar
- Kerr MK: Design considerations for efficient and effective microarray studies. Biometrics 2003, 59: 822–828.PubMedView ArticleGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res 2002, 30: e15.PubMed CentralPubMedView ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2005.Google Scholar
- Wu H, Kerr MK, Cui X, Churchill GA: MAANOVA: A Software Package for the Analysis of Spotted cDNA Microarray Experiments. In The analysis of gene expression data: methods and software. Edited by: Parmigiani G, Garrett ES, Irizarry RA, Zeger S. Springer-Verlag; 2003:313–431.View ArticleGoogle Scholar
- Cui XG, Hwang JTG, Qiu J, Blades NJ, Churchill GA: Improved Statistical Tests for Differential Gene Expression by Shrinking Variance Components. Biostatistics 2005, 6: 59–75.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society 1995, 57: 289–300.Google Scholar
- Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 2003,4(5):P3.PubMedView ArticleGoogle Scholar
- Tibshirani R: A simple method for assessing sample sizes in microarray experiments. BMC Bioinformatics 2006, 7: 106.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.