|
|
||||||||
1Respiratory Disease Area, Novartis Institutes for Biomedical Research, Horsham, West Sussex, United Kingdom; 2Genome and Proteome Sciences, Novartis Institutes for Biomedical Research, Cambridge, Massachusetts; 3Friedrich Miescher Institute for Biomedical Research, Basel, Switzerland; 4Airway Disease, Department of Cardiothoracic Surgery, National Heart and Lung Institute, London, United Kingdom; and 5Department of Environmental Medicine, Division of Lung Biology and Disease, University of Rochester Medical Center, Rochester, New York
Submitted 19 March 2007 ; accepted in final form 13 August 2007
| ABSTRACT |
|---|
|
|
|---|
chronic obstructive pulmonary disease; metabolic pathways; oxidant stress
Because of these clinical limitations, animal models have played an important role in attempts to elucidate disease-related mechanisms. A variety of agents have been used to mimic COPD-like changes in the lungs of rodents, including elastase, sulphur dioxide, nitrogen dioxide, ozone, lipopolysaccharide, and coal dust (reviewed in Ref. 21). Although some of these models have been important for understanding of how some of these pathological changes may occur, how these models correlate with the changes that occur in response to cigarette smoking remains unclear. There have also been a number of descriptive studies investigating the acute and chronic effects of CS exposure in rodents (reviewed in Refs. 38 and 39, respectively). Most of these studies were largely focused on characterizing the inflammatory and/or pathological changes that resulted, whereas only a few have focused on global molecular changes (9). In addition, few studies have characterized the effects of acute and chronic exposures within the same study (7). The use of knockout and transgenic animals in smoking studies has provided important insights into potential molecular mechanisms that are linked to disease pathogenesis (10a, 19, 29). These studies are powerful but limited to demonstrating the protective or destructive effects attributed to the abrogation or gross overexpression of a single gene. Because COPD susceptibility is likely to be heterogeneous and polygenic, understanding how macro changes at the level of the gene expression translate into pathological and physiological changes in response to CS will be essential for understanding how the disease develops.
The aim of this study was to identify the molecular changes induced by CS inhalation that may drive the biological and pathological consequences leading to disease. To do this, Sprague-Dawley rats were exposed to CS for 1 day and up to 34 wk. In addition, we studied the effect of smoking cessation by exposing one group for 26 wk followed by 8 wk of no exposures. Using Gene Set Enrichment Analysis, global changes in gene expression were translated into alterations in sets of genes with common molecular functions. We then measured protein levels of some of the genes for which expression changed in response to CS. We also characterized inflammatory cell infiltration, oxidant-induced damage to cells in the lung, pathologies at every level of the airways, and lung function. Here, we demonstrate a clear distinction between the acute and chronic response to CS and provide a new hypothesis to explain why these changes occur, resulting in the pathologies linked to COPD.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Animal maintenance conditions. Male Sprague-Dawley rats (350–400 g; Charles River, Margate, Kent, United Kingdom) were housed in rooms maintained at constant temperature (21 ± 2°C) and humidity (55% ± 15%) with a 12-h light cycle and 15–20 air changes per hour. Two animals were housed per cage containing two nest packs filled with grade 6 sawdust (Datesand, Manchester, United Kingdom), nesting material (Enviro-Dri; Lillico), maxi fun tunnels, and Aspen chew blocks (Lillico) to provide environmental enrichment. Animals were allowed food, RM1 pellets (SDS), and water ad libitum.
CS exposure.
Animals were exposed to whole body smoke as previously described (32) with some minor modification. Briefly, animals were exposed to 30 ml of smoke every 60 s with fresh air being pumped in for the remaining time. The smoke was generated using 2R1 research cigarettes (University of Kentucky, Lexington, KY), and sham control animals were exposed to room air only for the same duration of time (
60 min per exposure period). Animals were killed with an overdose of terminal anesthetic (sodium pentobarbitone, 200 mg ip) followed by exsanguinations 2 h after the last smoke exposure period. The carboxyhemoglobin levels were
5% at the time of termination. After 3 wk of exposures, plasma cotinine levels varied from
450 ng/ml 1 h after exposure to
150 ng/ml 22 h after exposure. Animals were culled, and changes were examined after 1, 3, and 5 days and 2, 3, 4, 6, 8, 12, 16, 26, and 34 wk of CS exposures. In addition, one group was exposed for 26 wk followed by 8 wk with no exposures to look at the effects of smoking cessation (referred to in the following text as the smoking cessation group).
Microarray analysis.
For microarray analysis, animals that did not undergo unconscious lung function maneuvers or lavage were used. The treatment and exposure regimen of the animals was otherwise identical. At the prescribed end point, animals were culled as described above, and the lungs were removed and flash-frozen in liquid nitrogen. Lung tissue was homogenized, and RNA was extracted using TRIzol according to the manufacturer's instructions. Protein contamination was checked at 260 and 280 nm, and a ratio of <1.6 was accepted. RNA was hybridized against Affymetrix RG_U34A GeneChips. The raw data was analyzed using GeneSpring version 6. To derive differentially expressed genes, a Welch t-test (parametric, not assuming equal variances) was used. The statistical cut-off used was 95% confidence. The fold regulation cut-off used was a
2-fold up- or downregulation comparing smoke- vs. air-exposed controls. Gene lists for the 12 time points were annotated computationally using DAVID, NetAffx, MAPPfinder, GenMAPP, and by systematic searches of PubMed. A nonredundant list of 340 genes that significantly changed between time-matched CS-exposed and sham animals was generated by combining the gene lists from the 12 time points. Data complied with minimum information about microarray experiments (MIAME) standard and have been deposited in the GEO database (accession nos. GSM170276–GSM170484).
Gene Set Enrichment Analysis. An in-house implementation of the Gene Set Enrichment Analysis method was used to analyze the data as previously described (24, 35). As input, Gene Set Enrichment Analysis requires microarray data from two conditions to be compared (e.g., sham- vs. smoke-exposed animals at each time point). Genes with expression levels below 100 on >75% of the GeneChips were discarded. Differential expression for each probe set was estimated by its expression ratio between pairs of conditions. Differentially modulated gene sets were identified using Wilcoxon's rank sum test.
Preparation of bronchoalveolar lavage fluid and lung tissue. Rats not used for microarray analysis were culled as above, and lungs were lavaged with 3 x 4-ml aliquots of sterile PBS. All aliquots were combined for individual rats. Bronchoalveolar lavage (BAL) total cell counts and differentials were performed as previously described (32). The right lung lobes were tied off, snap-frozen in liquid nitrogen, and used for protein analysis. The trachea and left lung were inflated with 5 ml of 10% neutral buffered formalin (NBF) for histological analysis.
Immunohistochemistry for 4HNE. Immunohistochemistry for 4-hydroxy-2-nonenal (4HNE) was performed by microwave irradiation method. Briefly, lung sections were dewaxed in xylene, rehydrated, and endogenous peroxidase-inhibited with 0.5% hydrogen peroxide in methanol for 10 min. Sections were stained with anti-4HNE adduct (Calbiochem) overnight at 4°C. Immunodetection was performed using biotinylated rabbit anti-mouse Ig reagent (Dako Cytomation, Cambridgeshire, United Kingdom), SABC reagent (Dako Cytomation, Cambridgeshire, United Kingdom), and diaminobenzidine (Sigma, Dorset, United Kingdom). The nuclei were counterstained with Cole's hematoxylin solution. Tonsil was used as a positive control, and for negative controls the primary antibody was omitted from one section of each of the samples. Cells were identified by both positive staining and standard morphological criteria. Staining intensity of 4HNE was assessed in alveolar macrophages semiquantitatively in a blinded fashion and graded 0 = no staining, 1 = weak staining, 2 = moderate staining, and 3 = strong staining. A Mann-Whitney's nonparametric statistical test was performed to determine significance of data.
BAL fluid and tissue cytokine analysis. Cell-free BAL and tissue homogenate supernatants were used to measure cytokine levels by ELISA (R&D Systems, Abingdon, United Kingdom, except for GRO/CINC-1, Amersham). Tissue homogenate protein levels were measured using the Bio-Rad protein assay (Bio-Rad, Hertfordshire, United Kingdom), and cytokine values were normalized against protein levels for individual homogenate samples.
Determination of BAL fluid mucin concentrations. Samples were assayed using a mucin enzyme-linked lectin assay (ELLA) as previously described (3, 15).
Immunohistochemical staining of inflammatory cells. Macrophages and neutrophils were stained as previously described (22). Macrophages were stained for ED1, and neutrophils were identified with the substituted naphthol method for AS-D chloroacetate esterase. Cells were identified using an AxioHOME microscope by both positive staining and standard morphological criteria.
Immunohistochemical staining of mucin and quantification of goblet cell density. Both tracheal and primary bronchus mucin were detected by a two-stage immunoperoxidase method using Ulex europaeus agglutinin-1 (UEA-1; Sigma) validated by Jackson and colleagues (15) as previously described (32).
Staining for collagen and smooth muscle. Sections (3 µm) were incubated in 0.1% Sirius Red F3BA (VWR, Poole, United Kingdom) to identify collagen or with mouse anti-human smooth muscle actin, clone 1A4 (Dako Cytomation, Ely, United Kingdom; see supplementary METHODS online at the AJP-Lung Cellular and Molecular Physiology web site for further detail). The areas of collagen and smooth muscle actin staining were assessed under a Zeiss Axioplan 2 microscope (x10 magnification) with an Imaging Associates KS400 image analyzer (Imaging Associates, Bicester, United Kingdom). Values were determined by measuring the inside perimeter of the airway and dividing it by the area of pigmented pixels, which was then converted into micrometers to give an average thickness surrounding the airway. To avoid bias on slant cuts, where the airway wall is thicker, a constant pixel count away from the airway was always measured when calculating the area. This prevented scoring of pseudo-area of pigmented tissue owed to more staining on obliquely cut airways. Collagen and smooth muscle actin were quantified on areas of the upper, upper-middle, lower-middle, and lower segments of the lung. The mean thickness for each sample was determined, and then the group (smoke- or air-exposed) was averaged to give a mean and standard error of the mean (SE) for each time point. A Mann-Whitney's nonparametric statistical test was performed to determine significance of data.
Quantifying mean linear intercept. Mean linear intercept was measured manually using a method previously described (37) using a reticule with a Thurlbeck grid. See supplementary text for details.
Conscious lung function measurements by whole body plethysmography. Rats were placed in a plethysmograph chamber (Buxco, Wilmington, NC) and left to acclimatize for 5 min. Bias flow (air removed from chamber) was set at 1 l/min. The pressure wave forms generated by respiration of the rats were analyzed by Buxco software. PenH (area under the curve; PenH AUC), peak expiratory flow (PEF), and breathing frequency were calculated on a breath by breath basis as previously described (2, 6, 17) and averaged every 30 s. All measurements were calculated over a 5-min recording period.
Unconscious lung function measurements using forced maneuvers. Respiratory function was measured using methods similar to those described previously (1, 4). See supplementary text for details.
Statistical analysis. All data are presented as means ± SE. Unless otherwise stated, statistical significance of data was determined using a two-tailed t-test with each group being treated as an independent group and significance measured only against the time-matched, air-exposed controls. A result was considered significant if P < 0.05.
| RESULTS |
|---|
|
|
|---|
|
Pathway changes in response to smoke.
During the first week, metabolic functions, in particular genes involved in oxidative phosphorylation and glycolysis (Fig. 2 and Supplemental Fig. 2), were most significantly elevated by smoke exposure. Signal transduction gene sets were significantly reduced. Several metabolic pathways remained upregulated over the course of the study (e.g., peroxisome proliferator-activated receptor-
signaling and carbohydrate metabolism); however, by 3 wk, pathways related to immunity, defense, and proteolysis were most significantly increased (Fig. 3 and Supplemental Fig. 3). Both inflammation and stress response/detoxification genes were upregulated acutely. Genes including the rat GRO homologs, IL-1
, and nuclear factor-erythroid 2-related factor 2 (NF-E2 factor 2 or Nrf2) were increased from 3 days and stayed elevated over the course of the study. Many of these changes became more pronounced as the duration of exposures increased (Fig. 3). In particular, those genes associated with macrophage function (e.g., matrix metalloproteinase-12) and phagocytosis (e.g., CD36) progressively increased over the course of the model. This was also true for genes associated with the lymphocyte function such as Ly-6C antigen, IgG, and immunoglobulin receptor expression. By 34 wk, gene sets related to immunity and defense and B cell-mediated immunity were the most significantly increased (Fig. 3 and Supplemental Fig. 4). Major histocompatibility complex class II- and macrophage-mediated immunity pathways were also very significantly elevated. Although slightly lower than at earlier time points, inflammation and stress response/detoxification gene sets remained increased in the smoking cessation group (Fig. 3). In contrast, most the gene sets involved in metabolic processes returned to control levels (Fig. 2).
|
|
|
|
was also increased in response to CS, however, only in the lung homogenate supernatants (see Supplemental Fig. 5). The levels were elevated and greatest after 3 days of exposure, increasing from 173 ± 25 pg/mg tissue in sham animals to 375 ± 57 pg/mg tissue in the CS-exposed rats (P < 0.01). The levels remained significantly increased over the course of the study but returned to control levels after smoking cessation. In addition, IL-6, IFN-
, monocyte chemoattractant protein-1, TNF-
, and IL-4 were assayed for but were below the levels of detection. Levels of IL-13 were measured and, although variable, appeared to be increased at some but not all time points (data not shown). Effect of CS exposure on inflammatory cell infiltration. CS increased the number of neutrophils in the BAL fluid in a biphasic pattern similar to that of the neutrophil chemokines (Fig. 6A). After an initial peak at 5 days, the number of neutrophils increased again starting from 6 wk and reached their greatest levels at 34 wk, increasing from 0.5 ± 0.1 x 105 in the sham animals to 14.7 ± 2.3 x 105 neutrophils recovered in the BAL fluid of CS-exposed rats. A similar biphasic pattern of neutrophil infiltration was observed in the lung tissue (Fig. 6B). In both the BAL and tissue, the numbers of neutrophils decreased from their peak levels in the cessation group; however, neutrophil numbers remained elevated after smoking cessation in BAL fluid (P < 0.001).
|
The most prominent inflammatory changes occurred in the alveoli, where pigmented, some binucleated, macrophages began to accumulate by 3 days. These aggregations of macrophages began to fill alveolar spaces by 3 wk, initially in those spaces adjacent to the terminal bronchioles (Fig. 7B). At 6–8 wk, the aggregates progressively spread into surrounding spaces and began to incorporate neutrophils (Fig. 6C). Many of these cells showed signs of degeneration and cell death as determined by nuclear pyknosis and debris. These inflammatory lesions progressively worsened over the duration of the study (Fig. 7D). The interalveolar walls were multifocally distended with aggregates of lymphocytes and macrophages. In addition, by 8 wk, a moderate perivascular lymphocytic inflammation was present in all CS-exposed animals (data not shown).
|
Effect of CS exposure on the epithelium of the central airways. CS induced epithelial hypertrophy and hyperplasia (increased number of mitotic bodies) that peaked at 5 days but was subsequently attenuated thereafter (Fig. 8). These changes were also associated with a significant increase in the number of mucin-containing cells in both the trachea (Fig. 8, A–D) and the primary bronchus. The area of epithelium stained positive for mucin in both the trachea and the bronchus was increased from 3 days and progressively increased over the course of the study, reaching a plateau by 16 wk (Fig. 8E). These changes in mucus levels in the lung tissue correlated closely to the amount of mucin detected in the BAL fluid (Fig. 8F).
|
After 16 wk of CS exposure, collagen deposition around the airways (other than primary bronchus) increased from 26.9 ± 1.7 to 56.0 ± 2.5 µm (P < 0.0001; Fig. 9A). As time progressed, the amount of collagen around the airways in sham animals increased. Although there was still more collagen around the airways of CS-exposed animals, the increase was no longer significant. The amount of smooth muscle around the airways was not significantly changed compared with air-exposed controls after 16 or 34 wk (Fig. 9B).
|
20%. No difference was detected between the groups that received 26 wk of either air or CS exposures (data not shown).
|
| DISCUSSION |
|---|
|
|
|---|
The major features of the acute phase included an immediate increase in the expression of metabolic gene sets and stress response genes, inflammation, and changes to the epithelium. Increased expression of metabolic genes sets was a new finding that has not been previously reported in the response to smoke. We propose that these changes are due to increased energy demands for repairing smoke-mediated damage. In addition, these pathways are critical for producing important reducing equivalents that may be needed to protect cells from the oxidizing effects of smoke. There is evidence to support the concept that energy transfer pathways and oxygen detoxification pathways coevolved and are increased in concert in response to environmental stresses such as CS (33).
Stress response genes were also increased during the acute phase. Enzymes such as heme oxygenase-1 (HO-1), thioredoxin,
-glutamylcysteine synthetase (
-GCS), and manganese-containing SOD (MnSOD) and the transcription factor Nrf2 were immediately upregulated, consistent with previous reports (9). These genes appear to be critical in the response to smoke and have been implicated in the pathogenesis of COPD. HO-1 expression is upregulated in alveolar macrophages of smokers (20), and promoter polymorphisms for this gene have been identified in certain COPD patients (40). Furthermore, overexpression of HO-1 protects against lung injury in elastase and CS models (31, 36). HO-1 and other antioxidant enzymes are regulated by Nrf2 (5), which is essential for protection against CS-induced inflammation and emphysema (29).
The inflammation in this acute phase was primarily neutrophilic in nature. The observed decrease in BAL macrophage numbers is likely due to the activation of these cells. In addition, oxidant modification to matrix proteins by CS extract has also been proposed to enhance the adhesion of macrophages to matrix (16). Another basic feature of the acute response to CS was the epithelial thickening and goblet cell metaplasia, which peaked after 5 days of exposures.
Subsequent to this acute phase, components of the innate, host defense response significantly decreased but did not completely resolve. Conversely, tissue macrophage numbers began to increase. The mechanisms underlying this transition to the second, more chronic phase in the response to smoke is unknown; however, we propose two hypotheses to explain these observations that may occur independently or collectively. First, alveolar macrophages, which migrate into the lung to resolve inflammation, are damaged by oxidants and may become overloaded with particulate matter with continual exposure to CS. These changes can result in reduced migration and enhanced activation of macrophages (25, 26). This may trigger the chronic inflammatory response that is observed following 2 wk of exposures. This hypothesis is consistent with the observations made in this model, with increased 4HNE staining in alveolar macrophages that are particle-laden and accumulating in large aggregates in the alveoli and alveolar ducts by 2 wk.
Second, ultrafine particles may localize to the macrophage mitochondria inducing structural damage (18) thereby reducing the metabolic capacity needed for protection against CS exposure (33). This damage to the metabolic machinery leaves the cell no longer capable of adapting to meet the added demands of further CS exposures. This may account, at least in part, for the "transcriptional shift" that occurs in gene expression patterns after 2 wk of CS exposures, and we propose that these changes are responsible for initiating the chronic response to smoke. In summary, we hypothesize that eventually the system becomes overwhelmed with repeated oxidant and particulate challenges, resulting in metabolic insufficiency and subsequently a greater oxidant burden leading to greater inflammation.
The second chronic wave of inflammation, oxidant damage, and mucus production was progressive and greater than the initial acute response. As shown previously, the inflammation shifted from a primarily neutrophilic response to one that had greater macrophage involvement (12, 27). Lymphocytes were present in these aggregates, and markers of the adaptive immune response (e.g., B cell-mediated immunity pathway, Ly-6C antigen and IgG chains) were increased in the chronic phase; however, the presence of lymphocytes were not as apparent as in smoking models reported in other species such as the mouse (7).
Associated with this second wave of inflammation were the start of structural changes in the lung consistent with those observed in smokers and COPD patients including collagen deposition around the airways and macrophage aggregation. There has previously only been one other report that we know of demonstrating a change in airway collagen in rats in response to CS exposure (30). The increase in airway collagen was no longer significant at 34 wk due to the increased level of collagen in the air-exposed controls. After 34 wk of exposure, there was a significant increase in the mean linear intercept of the alveolar spaces in CS-exposed rats. The failure to see greater structural damage to the lung and subsequent changes to lung function is possibly related to a lack of any viral and/or bacterial infection in this model system.
In the smoking cessation group, the changes induced by CS exposure were resolving; however, many of the inflammatory and antioxidant genes remained significantly elevated and were associated with increased numbers of inflammatory cells as has been documented in COPD patients (13). Despite this, other phenotypic changes (e.g., mucus production) had returned to control levels.
One limitation to this study is that the microarray experiments were conducted using RNA generated from whole lung rather than single cell types. Therefore, changes in expression profiles may reflect changes in cellular constituency of the lung at each time point. Ideally, it would be useful to analyze the expression changes occurring in a more homogenous cell population using laser capture microdissection. However, this study is strengthened by the fact that samples from individual rats were used for each time point in contrast to previously published work (9) where pooled samples were used. The use of individual samples allowed us to differentiate CS- and air-exposed groups in a more significant and comprehensive fashion. In addition, for the first time in a rodent smoking model, we used both pathway- and gene-centric analysis of the microarray data in combination with pathological observations. This approach placed these changes in a broader biological context and led to the identification of changes in biochemical processes previously not associated with the response to CS inhalation or COPD.
In summary, this study has provided new insights into the molecular pathways activated in response to CS and possibly involved in the pathogenesis of COPD. After a transient acute response where metabolic pathways appeared most significantly changed, a transcriptional shift in the gene expression pattern after 2 wk appears to initiate the chronic phase of the response to smoke inhalation. This shift may be due to changes in metabolic capacity and/or macrophage function, similar to the situation found in human (10, 12). Chronic CS inhalation leads to greater activation of immunity and defense pathways and recapitulated many of the phenotypic changes observed in COPD patients. These include oxidative damage to macrophages, a chronic, low-grade inflammation consisting of neutrophils, macrophages, and lymphocytes, epithelial damage, mucus hypersecretion, airway fibrosis, and emphysema.
These data suggest that the rat model may provide insights into the molecular mechanisms underlying the response to chronic smoking that lead to the destructive lung pathologies observed in COPD patients. In addition, using "Omics" technologies to characterize both animal models and clinical samples may help identify the gaps and parallels in these systems, as well as uncover prospective biomarkers, both of which are vital for identifying and measuring the efficacy of drug candidates for COPD.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. Rangasamy, V. Misra, L. Zhen, C. G. Tankersley, R. M. Tuder, and S. Biswal Cigarette smoke-induced emphysema in A/J mice is associated with pulmonary oxidative stress, apoptosis of lung cells, and global alterations in gene expression Am J Physiol Lung Cell Mol Physiol, June 1, 2009; 296(6): L888 - L900. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Churg, S. Zhou, O. Preobrazhenska, H. Tai, R. Wang, and J. L. Wright Expression of Profibrotic Mediators in Small Airways versus Parenchyma after Cigarette Smoke Exposure Am. J. Respir. Cell Mol. Biol., March 1, 2009; 40(3): 268 - 276. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Wright, M. Cosio, and A. Churg Animal models of chronic obstructive pulmonary disease Am J Physiol Lung Cell Mol Physiol, July 1, 2008; 295(1): L1 - L15. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. F. Chung and I. M. Adcock Multifaceted mechanisms in COPD: inflammation, immunity, and tissue repair and destruction Eur. Respir. J., June 1, 2008; 31(6): 1334 - 1356. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Churg, M. Cosio, and J. L. Wright Mechanisms of cigarette smoke-induced COPD: insights from animal models Am J Physiol Lung Cell Mol Physiol, April 1, 2008; 294(4): L612 - L631. [Abstract] [Full Text] [PDF] |
||||
![]() |
I-M. Wang, S. Stepaniants, Y. Boie, J. R. Mortimer, B. Kennedy, M. Elliott, S. Hayashi, L. Loy, S. Coulter, S. Cervino, et al. Gene Expression Profiling in Patients with Chronic Obstructive Pulmonary Disease and Lung Cancer Am. J. Respir. Crit. Care Med., February 15, 2008; 177(4): 402 - 411. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |