|
|
||||||||
TRANSLATIONAL PHYSIOLOGY
Departments of 1Anesthesiology and Critical Care Medicine, 2Pulmonary and Critical Care Medicine, and 3Pathology, The Johns Hopkins University, Baltimore, Maryland
Submitted 2 November 2005 ; accepted in final form 9 June 2006
| ABSTRACT |
|---|
|
|
|---|
adult respiratory distress syndrome; mechanical ventilation; canine; computed tomography; cross-species microarray
We believe that VALI does in fact have its origin in inflammatory and other cellular responses in lung tissue exposed to (and possibly, predisposed to injury from) mechanical stress and hypothesize that these responses will vary throughout the heterogeneous injured lung in relation to local mechanical events. To explore this hypothesis, we assessed regional cellular responses by genomic microarray analysis, correlating changes in gene expression with the specific regional mechanical stresses imposed on that local tissue, in a novel canine model of unilateral lung injury. With each animal serving as its own control, the microarray-based genomic approaches (cross-species hybridization to the Affymetrix Human U133A GeneChip) provide broad characterization of regional cellular responses which, combined with functional computed tomography (CT) imaging for the noninvasive measurement of regional mechanical stress, makes possible the investigation of in vivo local cellular responses in the lung associated with heterogeneous VALI. Our goal was to establish whether regional cellular responses to local conditions could be determined and provide a novel window into understanding the mechanism of development of VALI.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
All experimental procedures were approved by the Johns Hopkins University Animal Care and Use Committee. Three mongrel dogs (weight 21.3 ± 1.5 kg) were anesthetized with 25 mg/kg pentobarbital intravenously and instrumented with femoral arterial and venous catheters. Anesthesia was maintained with additional pentobarbital (5 mg/kg iv every hour and when indicated), and muscle relaxation was provided by pancuronium (3 mg bolus and 0.5 mg hourly iv). A 39 or 41 French double-lumen endobronchial tube (Mallinkrodt, St. Louis, MO) was placed via a tracheostomy, and position was confirmed by fiber-optic bronchoscopy. The animals were ventilated with a "dual piston" large animal ventilator (Harvard Apparatus, Holliston, MA) permitting independent control of tidal volume (Vt), inspired oxygen fraction (FIO2 1.0), positive end-expiratory pressure (PEEP), measurement of airway opening pressure (Paw), and end-tidal partial pressure of carbon dioxide (ETPCO2) for each lung. Oxygen saturation (SaO2) was continuously measured using a pulse oximeter applied to the tongue or ear, and ETPCO2, arterial blood pressure (Pa), Paw, and esophageal pressure (Pes) were continuously recorded. An infusion of lactated ringers (LR, 510 ml·kg1·h1) was given for maintenance fluid replacement. Rectal or pulmonary artery temperature was maintained at 36 ± 1°C with radiant heat lamps. At the conclusion of the study, the animals were killed by exsanguination after supplemental pentobarbital (10 mg/kg iv).
After instrumentation, the individual lung Vt were adjusted to provide an ETPCO2 of 3035 mmHg for each lung at a respiratory rate of 20 breaths/min. This procedure resulted in individual lung Vt (means ± SD) of 8.7 ± 0.3 ml/kg for the right lung and 6.1 ± 0.4 ml/kg for the left, approximating the normal distribution of lung volume. The left lung was then mildly injured by repeated lavage with warmed saline, 20 ml/kg repeated four times while switching the animal's position between right and left lateral decubitus, with gravity drainage. For comparison, typical lung lavage injury protocols use 4060 ml/kg (both lungs) for six to eight washes to achieve an arterial partial pressure of oxygen (PaO2) of
100 mmHg. Ventilation continued at baseline settings with the addition of 5 cmH2O PEEP to the control right lung for 5 h. At that time, the animal was killed by exsanguination, the chest was opened, and eight tissue samples were taken from four corresponding regions in both lungs (apex dependent and nondependent, base dependent and nondependent). The lung tissue samples were immersed in RNA later (Ambion, Austin, TX) and frozen for subsequent analysis. In some animals, additional tissue samples were taken and stored in 10% formalin for histology.
CT Imaging and Analysis
One animal was transported to a Toshiba Aquilion 16 multidetector CT scanner for imaging at 4 h postinjury. During the 15-min transport, both lungs were ventilated with a single transport ventilator, but independent lung ventilation was reestablished before imaging. Whole lung image sets gated to end-expiration and end-inspiration were obtained during steady-state ventilation, four 2-mm contiguous slices per breath, using the ECG-gating capability of the scanner and timing signals generated by a computer running LabView 6.1 (National Instruments, Austin, TX). Image density distributions were analyzed using Pulmonary Analysis Software Suite (27) developed at the University of Iowa Division of Physiologic Imaging.
Sample Preparation and Microarray Analysis
Expression data analysis. Each regional tissue sample was run on a separate Affymetrix U133A GeneChip. The signal intensity fluorescent images produced during Affymetrix GeneChip hybridizations were read using the Agilent Gene Array Scanner and converted into GeneChip probe result files (CEL) using MAS 5.0 software (Affymetrix, Santa Clara, CA). The analysis of the probe level data (available in the CEL files) was performed using the Bioconductor affy package (31) as described previously (22). Briefly, the MAS5 (30) module of this package was used for background correction, across array global normalization, and extraction of the probe level data. Poorly performing probe pairs revealed by probe level analysis were masked before converting CEL files into gene expression values (22, 23). This approach increased present call of a transcript (P < 0.04) on average by 45% compared with unadjusted probe set processing. The remaining absent calls (P > 0.06) for transcript abundance were assumed to be a result of undetectable message concentration [<1 pM (9)] rather than technical or detecting errors. Therefore, absent calls for each GeneChip were averaged and all absent calls for a given chip were replaced with the average value. This modified data set was used for further analysis.
Selecting significant gene expression changes. The gene expression profiles of control and injured lungs were paired by animal and lung region from which samples were collected. Signal intensity values for each pair were log transformed (base 2) and raw-wise normalized, thus rendering each profile pair an independent data module composed of normalized hybridization signals of corresponding control and injured samples. For apex/base comparisons, there were six modules in each group (dependent and nondependent samples for each apex and base location in 3 animals), and there were three modules in each group for dependent base/nondependent base comparisons. These modules were grouped and tested for similarities in gene expression either between base and apex regions (12 apex vs. 12 base) or dependent and nondependent areas of the lung base (6 nondependent vs. 6 dependent) using Significance Analysis of Microarray [SAM (53)]. To evaluate the significance of differences in gene expression between apex/base or nondependent/dependent regions, the control/injury modules were converted into injured/uninjured fold-change values (log base 2) and analyzed using SAM, comparing six apex vs. six base or three nondependent vs. three dependent, respectively. The injured/uninjured fold-change ratio for individual genes was computed by SAM as the antilog of the mean log-transformed individual dog-injured/uninjured ratios. All data are thus reported as the absolute (non-log-transformed) fold-change ratio. Genes with 70% change (±1.70-fold) in expression and false discovery rate (FDR) <2% were considered affected by lung injury based on the sample size analysis described below.
To determine the appropriate fold-change cutoff, we utilized a new concept for sample size determination in microarray experiments for class comparisons (11). The median (
2 = 0.169) of signal intensity log ratios in dog lung tissues was calculated from 10 arrays generated for uninjured lungs and was utilized as a combined variance following the recommendations of Yang and Speed (59). This derived variance value for inbred dogs is fairly typical for models employing homogeneous subject populations [e.g., inbred mice strains demonstrated
2 =
0.25 (11)]. We applied
2 = 0.169, 2% FDR (
= 0.02), and number of replicates n = 3 to the single-labeled microarray sample size-identifying formula and obtained a significance cut-off fold-change value
= ± 1.70 for changes in gene expression between dependent and nondependent base regions:
![]() | (1) |
/2 and z
are the 100
/2th and 100
th percentiles of the normal distribution, respectively;
is the distance between the class means that the study is being powered to detect; and
g2 is the combined gene expression variances (biological variance due to heterogeneity from individual to individual of gene expression within the phenotype and experimental error variance due to technical inaccuracies in the microarray measurement variances). Despite this theoretical approach to determine an appropriate fold-change cut-off, it should be recognized that the actual value used remains somewhat arbitrary and, because of the small number of animals in the study, data interpretation should be with caution. Semiquantitative Duplex RT-PCR
Total RNAs (0.5 µg) from each tissue sample were used as a starting material for the reverse transcriptase (RT) reaction. Primers for PBEF, THBS1, GADD45a, DEADH, GGH, and ribosomal protein S18 (housekeeping gene control) were designed based on published GenBank mRNA sequences (NM_005746.1, NM_003246.1, NM_001924.2, R0068, NM_003878, and NM_022551.2, respectively). One tenth of RT product was used in a typical 25-cycle PCR run. The PCR products were separated on 2% agarose gel electrophoresis and visualized by ethidium bromide staining. The image was captured on thermal paper (Fotodyne System, Hartland, WI) and scanned into Adobe Photoshop 7.0.1 using an HP Scanjet 7400c scanner, and band intensities were analyzed with AlphaEase FC Software (Alpha Innotech, San Leandro, CA). Results were expressed as a fold-change of injured tissues vs. control tissue.
| RESULTS |
|---|
|
|
|---|
Over the course of the 5-h observation period, the peak inspiratory pressure of the injured lungs approximately doubled, whereas that of the control lungs remained unchanged. SaO2 remained at 99100% at all times despite unilateral lung injury because of the use of a high FIO2 and the presence of the control lung. All animals remained normotensive throughout the study without intervention. Lung wet/dry weight ratios (w/d, means ± SD) were lower in control lungs (5.76 ± 1.01) compared with injured lungs (8.08 ± 1.86, P < 0.01), with w/d in injured lungs greater in base regions (9.14 ± 1.37) than apex regions (7.01 ± 1.70, P = 0.02). There was no significant difference in w/d between base and apex regions in the control lungs (6.15 ± 1.20 vs. 5.38 ± 0.66, P = 0.13).
Representative tissue histological sections show severe acute inflammation in the left (injured) lung base, with reduced inflammation in the left nondependent base (Fig. 1, A and B). In comparison, the right (control) lung sections revealed minimal inflammatory changes (Fig. 1, C and D). Apical sections of both control and injured lungs showed minimal inflammation.
|
|
A total of 24 samples from three animals were successfully hybridized to 24 Affymetrix HG-U133A microarray chips and are included in this analysis. Since samples were taken from corresponding regions of the injured and control lungs of each animal, results from each regional injury sample were normalized to the within-animal control, a unique feature of this unilateral injury model. Initial analysis of expression profiles generated using Affymetrix MicroArray Suite (MAS 5.0) showed that
14% (
3,000 of >22,000 probe sets analyzing 18,000 transcripts and variants) of the canine genome hybridized to corresponding human probe sets with an efficiency similar to human mRNA. Sequence comparison showed that this fraction of canine mRNAs represented canine genes that were highly homologous to their human counterparts. Overall transcript detection by the cross-species hybridization of canine mRNA to the HG-U133A chip was much less sensitive (14%) than hybridization of human mRNA (typically 48%). However, application of canine-specific masking techniques (22, 23) in which nonhomologous probe sets between species are excluded from the usual multiprobe set matching requirements increased transcript detection to >65%.
Comparing regions, there were 308 and 595 genes significantly changed [defined as FDR <2% and average regional injured/control fold-change (FC) >1.7 or <1.7] in the apex and base, respectively, of which 213 (69.2%) increased expression in the apex and 238 (40.0%) increased in the base; 367 genes were significantly differentially expressed between the apex and base (FDR <2% for apex/base difference and either apex and/or base injured/control FC >1.7 or <1.7). Cluster analysis of significantly changed genes revealed groupings of differential regional gene expression changes, along with the most represented gene ontologies (12) (Fig. 3).
|
|
|
|
|
GeneChip results for five different upregulated genes were validated using semiquantitative RT-PCR, showing good agreement in fold-change injured/control between samples taken from apex (n = 6) and base (n = 6) regions (Fig. 6). There were no statistically significant differences between microarray and RT-PCR measured fold-changes by paired t-test (P range 0.390.57).
|
| DISCUSSION |
|---|
|
|
|---|
Our goal was to establish whether regional cellular responses to local mechanical conditions could be determined as a first step towards using these regional events as a novel window into understanding the mechanism of development of VALI. We developed a unique canine model of unilateral lung injury and mechanical ventilation to allow the comparison of control and injured tissue samples from the same animal. We chose the canine saline lavage model of ALI because the insult is primarily mechanical and serves to amplify the effects of mechanical ventilation by removing surfactant and reducing stability of peripheral air spaces, although a secondary inflammatory response does ensue (48). The progressive increase in peak inspiratory pressure over the 5 h seen only in the lavage lung indicates the interaction between mechanical ventilation and the initial lavage insult that is characteristic of ventilator-associated lung injury. Previous experience with this model had demonstrated significant mechanical heterogeneity (14, 34) not unlike that observed in patients with ARDS (18, 20), with dependent, basal flooding and air space collapse, intermediate zones of apparent air space opening and closing, nondependent overdistension, and relatively preserved apical regions, behavior confirmed by our CT data from one animal. Furthermore, successful cross-species hybridization to human microarrays had been described (38, 40, 57) and canine-specific microarrays were under development (4, 6, 26, 50). Finally, the canine model is of sufficient size to allow lung isolation and independent treatment, permitting the use of each subject as its own control to reduce the baseline variation inherent in unbred subjects. The benefit of decreased intersubject variation from this control scheme should offset possible changes to the control lung from systemic "spillover" of inflammatory or other mediators from the injured lung (49, 51), although the possibility remains of false negative results from this effect.
Although it is known that there is regional mechanical heterogeneity within the injured lung (19), and we documented regional differences in tissue aeration, w/d ratio, and inflammation, there are potentially other factors that could influence these results. In addition to differences in regional ventilation, there are likely differences in regional perfusion within the injured lung. In the lavaged lung, dependent, poorly aerated regions reduce their perfusion via the homeostatic hypoxic pulmonary vasoconstriction response (15). Similarly, although postural repositioning was used in an attempt to uniformly distribute the saline lavage, the removal of surfactant and distribution of residual fluid may be nonuniform. Finally, the number and type of migratory inflammatory cells differ between regions. Thus the stimulus for the differential regional responses found may not be purely the local mechanical stresses but could reflect contributions of other heterogeneous aspects of this large animal injury model.
The extension of microarray technology to large animal models of disease has been limited by the lack of availability of species-specific gene chips, and thus cross-species hybridization has been increasingly utilized for this purpose (38, 40, 57). This approach is useful for gene products with an adequate degree of homology to their human counterparts, whereas failure to detect an expressed gene cannot distinguish between lack of stimulus effect vs. failure to hybridize because of species-specific sequence differences. The present call rate can be increased using masking techniques that either eliminate poorly functioning probe sets or relax the stringent criteria for detecting all the probe set sequences in Affymetrix-type arrays, since differences between species may be limited to only portions of the sequence (22, 23, 26). As our purpose was to demonstrate regional differences in gene expression in response to different mechanical stresses, false negatives are irrelevant with the acknowledgment that these results are but a subset of the gene expression changes that have occurred. Recent announcement of the availability of a commercial canine microarray from Affymetrix as well as the development of custom canine cDNA arrays (4, 6, 50) will facilitate more specific and exhaustive analysis in the future. In any event, the use of microarray data in this manner should be considered exploratory and any specific mechanistic inferences should be validated with more precise techniques.
Analysis of the microarray data for the injured vs. control lung showed significantly changed expression in 472 genes, of which 46% were increased. The functional groupings (ontologies) of these genes were similar to those we have previously described (22, 24), including angiogenesis, cell cycle arrest, inflammation, coagulation, and immune response. Confirming our initial hypothesis, there were large numbers of genes differentially regulated between the apex and base (Fig. 3) and between dependent and nondependent regions within the base (Table 1 and Figs. 4 and 5). These lung regions experience strikingly different mechanical environments, as evidenced by their CT density distributions (Fig. 2). Furthermore, imaging studies of ARDS survivors have shown that chronic lung damage, mostly bullous lesions, occurs predominantly in dependent and basal lung regions, whereas apical regions are relatively spared (18). The mechanical events related to artificial ventilation most commonly invoked as injurious to the lung are overdistension, airways repeatedly opening and closing with each breath, and to a lesser extent tissue collapse and/or flooding, although the data supporting these relationships in patients remain controversial (28). In our model, which used no PEEP in the injured lung, there was minimal evidence of overdistension based on the end-inspiratory criterion of >90% aeration by CT (54), and any overdistension present occurred toward the apex. Of course, it is also possible that differences in other regional properties, such as perfusion or migratory inflammatory cells, influenced the observed changes in gene expression.
Examination of the highly differentially regulated genes between apex and base regions reveals several genes commonly associated with ALI (vascular endothelial growth factor C, THBS1, tranforming growth factor-
) as well as novel genes not previously described in the ALI literature. For example, the gene encoding the proposed proinflammatory cytokine Pre-B cell colony enhancing factor (PBEF) was one of the most upregulated genes in the injured lung. PBEF was originally described for its role in the maturation of B cell precursors (42) and was subsequently found to be upregulated in amniotic membranes from patients undergoing premature labor (especially with amniotic infections) (37). Epithelial cells from amniotic membranes also upregulate PBEF expression when subjected to stretching in vitro (41), and human fetal membrane explants exposed to PBEF produce inflammatory cytokines and chemokines (43). Furthermore, PBEF was recently identified as a cell cycle regulator, delaying neutrophil apoptosis in experimental inflammation and clinical sepsis (32) and thus functions as a regulator of the inflammatory response. Thus PBEF represents a proinflammatory cytokine, with proposed roles in immune regulation and oxidative metabolism and known association to infected or stretched epithelial cells, which has now been identified as significantly upregulated in a primarily mechanical model of VALI. These observations led us to pursue in-depth investigation of its role in other ALI disease models and in patients with sepsis-related ALI, revealing a consistent upregulation of PBEF in murine and canine LPS-induced ALI, in cytokine and cyclic stretch lung endothelial cells, and increased PBEF protein in both bronchoalveolar lavage and serum of patients with ALI (60). Furthermore, analysis of two single nucleotide polymorphisms in patients with sepsis-associated ALI suggests that a susceptible haplotype (GC) in the PBEF promoter region is associated with a 7.7-fold higher risk of developing ALI (60), strongly supporting a role for PBEF both in the pathogenesis of ALI and as a biomarker for ALI susceptibility.
Among the genes differentially regulated within the lung base were several stress response proteins, notably the inducible 70-kDa heat shock protein (HSP) whose expression has been associated with protection from development of ALI (55). Three HSP family genes, along with related chaperonins CCT7 and CCT8 and early response genes EGR3 and IER5, were all upregulated in both dependent and nondependent base regions. The increased expression of all these genes was significantly greater in the nondependent regions (2.3- to 6.5-fold) compared with dependent regions (1.0- to 1.7-fold) (Fig. 5). In keeping with the reported protective and anti-inflammatory actions of these proteins (55, 58), expression of cytokines (TNF-
, IL-2, IL-6 receptor, and IL-7) and chemokines (IL-8, CCL, and CXCL families) was decreased in the nondependent regions (1.1- to 2.0-fold), whereas it increased in dependent regions (1.0- to 2.0-fold). Thus, in the dependent regions of the lung base, which are flooded and collapsed with minimal penetration of air throughout the respiratory cycle (Fig. 2), there is an enhanced inflammatory response. In contrast, the well-aerated and potentially overdistended nondependent regions experience a robust stress response along with reduced cytokine and chemokine expression. These disparate responses are consistent with the results of a recent study in rats in which bacterial endotoxin plus injurious ventilation revealed a similar inverse relationship between HSP70 levels and IL-6 and IL-1
mRNA expression (56). Studies examining the time course of these changes and further modulating mechanical conditions will help determine how these factors contribute to protection vs. injury in VALI.
To further illustrate the significance of these regional differences in gene expression, we compared our list of genes differentially regulated between dependent and nondependent regions of the lung base with a list of candidate genes identified from the literature (3, 10, 13, 22, 24, 29, 35, 60). Of the
226 ALI candidate genes identified, 87 (including closely related genes in the same family) met our criteria for significantly differentially changed between dependent and nondependent regions. Remarkably, expression of 19 of these genes were changed in opposite directions, mostly upregulated in the dependent regions while downregulated in the nondependent region of the same lung in the same animal. Thus, the particular pattern of gene expression in the injured lung is critically dependent on the sampling location. Whether these regional differences are specifically due to the different mechanical stresses, the population of cell types in these regions or some other aspect of the local physiology such as perfusion remains unknown. It should be apparent, however, that sampling approaches that average expression patterns across both apex and base, dependent and nondependent regions are likely to underestimate the impact of genes that exhibit regional-specific expression. Averaging out of regional differences may explain the increasing numbers of differentially regulated genes detected as the regions examined become more localized. Furthermore, insight into the earliest cellular responses constituting VALI, and the contribution of different regional mechanical phenomena to these responses (that may be modified by patient and ventilator management), must be sought in the context of the mechanical heterogeneity that characterizes human ALI.
In summary, we have demonstrated significant differences in gene expression between different regions of the lung in a canine model of ALI, regions that undergo very different mechanical stresses during mechanical ventilation. These findings have important implications for the design and interpretation of studies investigating the mechanisms of VALI. A hallmark of human ALI is the heterogeneity of tissue involvement, and management strategies that reduce heterogeneity, in particular extremes of mechanical behavior, are the focus of current clinical investigations. Correlating regional cellular responses with the local mechanical and biological milieu may help us to not only understand how different mechanical stresses cause or propagate existing lung injury but also identify therapeutic targets, associate specific biomarkers with different aspects of evolving injury, and find surrogate endpoints to better predict outcomes and guide therapy.
| GRANTS |
|---|
|
|
|---|
| 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 |
|---|
|
|
|---|
, LPS, and cyclic stretch. Physiol Genomics 19: 331342, 2004.This article has been cited by other articles:
![]() |
R. C. Geiger, C. D. Kaufman, A. P. Lam, G. R. S. Budinger, and D. A. Dean Tubulin Acetylation and Histone Deacetylase 6 Activity in the Lung under Cyclic Load Am. J. Respir. Cell Mol. Biol., January 1, 2009; 40(1): 76 - 82. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. M. Otto, K. Markstaller, O. Kajikawa, J. Karmrodt, R. S. Syring, B. Pfeiffer, V. P. Good, C. W. Frevert, and J. E. Baumgardner Spatial and temporal heterogeneity of ventilator-associated lung injury after surfactant depletion J Appl Physiol, May 1, 2008; 104(5): 1485 - 1494. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. R. Martin Interactions between Mechanical and Biological Processes in Acute Lung Injury Proceedings of the ATS, April 15, 2008; 5(3): 291 - 296. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. N. Grigoryev, M. Liu, H. T. Hassoun, C. Cheadle, K. C. Barnes, and H. Rabb The Local and Systemic Inflammatory Transcriptome after Acute Kidney Injury J. Am. Soc. Nephrol., March 1, 2008; 19(3): 547 - 558. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Frank and D. J. Erle Progress toward a systems biology approach to acute lung injury Am J Physiol Lung Cell Mol Physiol, August 1, 2007; 293(2): L290 - L291. [Full Text] [PDF] |
||||
![]() |
G. B. Allen, T. Leclair, M. Cloutier, J. Thompson-Figueroa, and J. H. T. Bates The response to recruitment worsens with progression of lung injury and fibrin accumulation in a mouse model of acid aspiration Am J Physiol Lung Cell Mol Physiol, June 1, 2007; 292(6): L1580 - L1589. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. He, B. Han, and M. Liu Long pentraxin 3 in pulmonary infection and acute lung injury Am J Physiol Lung Cell Mol Physiol, May 1, 2007; 292(5): L1039 - L1049. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. J. Meyer and J. G. N. Garcia Wading into the Genomic Pool to Unravel Acute Lung Injury Genetics Proceedings of the ATS, January 1, 2007; 4(1): 69 - 76. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. G. N. Garcia and L. Moreno Vinasco Genomic insights into acute inflammatory lung injury Am J Physiol Lung Cell Mol Physiol, December 1, 2006; 291(6): L1113 - L1117. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |