Global gene expression patterns in the post-pneumonectomy lung of adult mice
© Paxson et al; licensee BioMed Central Ltd. 2009
Received: 23 June 2009
Accepted: 5 October 2009
Published: 5 October 2009
Adult mice have a remarkable capacity to regenerate functional alveoli following either lung resection or injury that exceeds the regenerative capacity observed in larger adult mammals. The molecular basis for this unique capability in mice is largely unknown. We examined the transcriptomic responses to single lung pneumonectomy in adult mice in order to elucidate prospective molecular signaling mechanisms used in this species during lung regeneration.
Unilateral left pneumonectomy or sham thoracotomy was performed under general anesthesia (n = 8 mice per group for each of the four time points). Total RNA was isolated from the remaining lung tissue at four time points post-surgery (6 hours, 1 day, 3 days, 7 days) and analyzed using microarray technology.
The observed transcriptomic patterns revealed mesenchymal cell signaling, including up-regulation of genes previously associated with activated fibroblasts (Tnfrsf12a, Tnc, Eln, Col3A1), as well as modulation of Igf1-mediated signaling. The data set also revealed early down-regulation of pro-inflammatory cytokine transcripts and up-regulation of genes involved in T cell development/function, but few similarities to transcriptomic patterns observed during embryonic or post-natal lung development. Immunohistochemical analysis suggests that early fibroblast but not myofibroblast proliferation is important during lung regeneration and may explain the preponderance of mesenchymal-associated genes that are over-expressed in this model. This again appears to differ from embryonic alveologenesis.
These data suggest that modulation of mesenchymal cell transcriptome patterns and proliferation of S100A4 positive mesenchymal cells, as well as modulation of pro-inflammatory transcriptome patterns, are important during post-pneumonectomy lung regeneration in adult mice.
Pulmonary emphysema is an example of a chronic disease with parenchymal destruction, where repair is relatively ineffectual . To provide effective therapies for treating this disease, a better understanding of the cellular and molecular processes that govern the phenomenon of lung regeneration, and in particular alveolar regeneration, is crucial. An important approach is the analysis of tissues from animal species that retain a high degree of regenerative capacity. Adult mice are capable of regenerating functional alveoli following either lung resection or injury to a greater degree than the regenerative capacities observed in larger adult mammals [2, 3]. In healthy adult mice (or rats) unilateral pneumonectomy evokes compensatory lung regeneration from the remaining lung lobes , in part through neoalveolarization within the existing parenchyma . This regenerative process is characterized by restoration of lung volume, surface area, morphometry, DNA and protein content within 14 days, as demonstrated by our lab as well as others [6, 7]. Despite a plethora of macrophysiologic and morphometric studies on lung regeneration in rodents and larger animals , the cellular and molecular mechanisms that regulate this process are not well understood. Previous studies using gene expression have focused on specific pathways  rather than global transcriptomic approaches. For example, using a gene array (588 genes) designed for analysis of transcription factors, Landenberg et al demonstrated a 2-fold or greater up-regulation of six genes, including early-growth response gene-1 (Egr-1), Nurr77, tristetraprolin, the inhibitor of kB-alpha (IkB-alpha), Klf-4 (GKLF) and LRG-21, all within two hours of pneumonectomy in mice . The authors concluded that expression of early transcription factors (i.e. early immediate genes) activated by mechanical stress trigger a cascade of growth signals that promote lung regeneration. Likewise, repeated over-inflation of the murine lung soon after pneumonectomy (30 min) was associated with over-expression of the proto-oncogenes c-fos and junB , underscoring the ability of pneumonectomy-induced mechanical stretch to evoke transcription factors. Indeed, lung stretch by mechanical ventilation without pneumonectomy induces similar early immediate gene transcription .
While past studies have revealed early immediate genes that participate in the activation of lung regeneration, the majority of the regenerative process takes place over a prolonged period (7-14 days). Gene expression patterns corresponding to important biologic processes such as cellular proliferation, matrix formation, angiogenesis, and progenitor cell differentiation have not been fully characterized. It is also unclear from past studies why processes such as matrix formation and angiogenesis occur during the remodeling process, but are not associated with fibrosis in this context.
The objective of this study was to measure the effects of pneumonectomy (vs. sham surgery) on gene transcriptome patterns that are robustly expressed (fold change ≥1.5, or ≤ -1.5) at multiple time points during lung regeneration. Therefore, we analyzed the transcriptome (over 39,000 genes) from mouse lung tissues following unilateral pneumonectomy using Affymetrix GeneChip microarray technology. Samples were taken at four time points (6 hours, 1 day, 3 days and 7 days) post-pneumonectomy spanning the period during which the bulk (>80%) of lung regeneration occurs, as measured by changes in vital capacity . Following analysis of the transcriptomic patterns, an immunohistochemical study of the regenerating parenchyma using the fibroblast markers S100A4 and αSMA was also performed at two points during lung regeneration (3 days and 7 days) to further elucidate the role of fibroblasts in this process.
Animals used for microarray analysis and q RT-PCR
Mice used in this study were adult (10-12 week) female C57BL/6 (20-25 g) obtained from Jackson Laboratories. All experiments were performed in accordance with NIH guidelines, as dictated by Institutional Animal Care and Use Committee at Tufts University. For each of the four time points, the mice were divided into two groups: (1) pneumonectomy (PNY) and (2) sham operated (SHAM - thoracotomy without lung resection), with eight animals in each group. Mice were anesthetized by intraperitoneal injection of ketamine (50-75 mg/kg) and xylazine (5 mg/kg), and then received 2 ml of warmed normal saline and 100 mg/kg sodium ampicillin subcutaneously. Orotracheal intubation was performed under direct visualization using a 20-gauge catheter (BD Insyte catheter; Becton, Dickinson and Co, Franklin Lakes NJ) over a flexible stylet. Mice were secured in supine position, and mechanically ventilated (AUT6110, Buxco Electronics, Wilmington, NC) at 200 tidal breaths of 0.3 ml of room air per minute, at positive end-expiratory pressure of 3 cm H2O during surgery and recovery.
After achieving adequate anesthetic depth (absence of response to toe-pinch) the left thoracic wall was clipped and disinfected. The skin, chest wall and pleura were incised at the 5th intercostal space, and the left lung was gently lifted through a ~5-7 mm incision and ligated at the hilum with 4-0 silk (Sofsilk, Synture Norwalk Ct). The lungs were then inflated to 30 cmH20 airway pressure, and the chest wall closed during this inflation with a single interrupted suture. The skin was closed with 5-0 PDS in a simple interrupted pattern. Mice were extubated at the onset of vigorous spontaneous breathing. The mice recovered from surgery in a warmed cage, and post-operative pain was managed with buprenorphine subcutaneously (0.05 mg/kg) as soon as mice showed conscious motor control, and every 12 hours thereafter as needed (<3 days). Chow, nutrient gel (on the cage floor), and water were provided ad libitum. Sham pneumonectomy animals underwent an identical procedure, except that after the thoracotomy, the chest was left open for 5 minutes to simulate the conditions of the pneumonectomy group without removal of the left lung, then closed as described.
Tissue preparation and RNA isolation
The mice were anesthetized as above at 6 hours, 1 day, 3 days and 7 days after surgery (PNY or SHAM) then euthanized by cervical dislocation. The pulmonary vasculature was perfused with ice cold Hanks balanced salt solution, the trachea cannulated, and the lungs removed en bloc. RNA preservation was achieved by flooding the lung intratracheally with RNAlater solution (Qiagen #76104), followed by storage of lung tissue samples in RNAlater at -80°C.
RNA isolation and microarray analysis
Equal amounts of lung tissue were pooled from eight animals in each group (PNY or SHAM at each of 4 time points) to minimize biological variability . Total RNA was prepared from the dissected lung tissue using the Qiagen RNAeasy mini kit (Qiagen #74104) according to the manufacturer's directions. The total RNA samples from the primary purification were purified a second time on Qiagen RNAeasy columns according to the manufacturer's instructions. Total RNA concentrations, A260/A280 and A260/A230 ratios were determined using a NanoDrop ND1000 spectrophotometer. All microarray analysis was performed as described in the Affymetrix GeneChip Expression Analysis Technical Manual using Affymetrix Mouse Genome 430 2.0 GeneChips and the One-Cycle cDNA Synthesis and HT IVT Labeling kits (Affymetrix Inc.). The complete microarray dataset is available (accession number GSE15999 at: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE15999).
Microarray data analysis on PNY versus SHAM animals (time-independent)
The initial goal in the analysis was to identify genes differentially regulated in the comparison between pneumonectomized and sham-operated animals. By treating all time points as replicates within their respective groups two datasets were created (PNY and SHAM), each with an n = 4. The corresponding Affymetrix CEL files were background corrected, summarized and quantile-normalized using the RMA library in BioConductor http://www.bioconductor.org, yielding one expression value per probe set for each of the 8 arrays . Based on the 'Rank Products' algorithm proposed by Brietling, et. al , the RankProd library was employed to find differentially expressed genes. This algorithm works by performing comprehensive pair-wise comparisons to calculate a rank statistic RPg, defined as the probability of seeing the observed, pair-wise expression patterns for any given gene g. As a vehicle for measuring statistical significance a non-parametric P-value is also calculated, using 1000 permutations to determine how often the calculated RPg statistic would occur by chance alone. Finally, the RankProd library compares average expression between the two groups to derive a fold-change value. Genes with a reported P-value < 0.001 and a fold-change ≥ 1.5 or ≤ -1.5 were selected for further investigation.
Microarray data analysis to identify temporal changes in lung regeneration (time dependent)
In a second analysis, the focus was shifted to temporal changes in gene expression during lung regeneration as opposed to overall transcriptomic patterns. With only one array per experimental condition at each time point, derivation of statistical measures and the subsequent search for truly differentially expressed genes can be challenging. However, the S-Score algorithm described by Zhang et al. (2002) and Kerns et al. (2003) provides a method for determining statistical significance when biological replicates are not available by applying pair-wise comparisons to probe-level data [16–19]. On average, the Affymetrix 3' IVT platform contains 22 probes for every transcript represented on the array. Using this information directly, the S-Score algorithm has shown good sensitivity when compared to many other existing analysis methods without sacrificing specificity (including RMA, dChip and MAS5), and can produce accurate results when no biological replicates are present [18, 19]. This is particularly applicable and appropriate to our individual time point datasets in which we have only one paired array set for each time point. Using the S-Score algorithm, the relative change in probe pair intensity is calculated to convert the probe pair signal differences into multiple measurements with equalized errors. The relative changes for each probe pair are then summed to form the S-score, which represents a single measure of the significance of change for the gene in question. By definition, S-score is related to P-value by an exponential relation, and a value of 3 corresponds to a P-value of 0.003 [16, 19, 19]. Genes with an S-score ≥ 3.0 or ≤ - 3.0 (P ≤ 0.003) were selected for further analysis.
Ingenuity Pathways Analysis
For selected genes (genes with a P-value < 0.001 and a fold-change ≥ 1.5 or ≤ -1.5 for the time-independent analysis; genes with an S-score ≥ 3.0 or ≤ - 3.0 (P ≤ 0.003) for the time-dependent analysis), Ingenuity Pathway Analysis (IPA) version 2.0 (Ingenuity® Systems Inc, Redwood City, CA; http://www.ingenuity.com) was used to search for biological functions and interrelationships between significantly modulated genes in PNY versus SHAM mice. IPA provides a large manually curated database containing over 200,000 full text articles and information about thousands of human, mouse and rat genes  with which experimental data sets can be statistically compared. Genes from the dataset were overlaid onto a global molecular network developed from information contained within the IPA database, and networks of genes in the dataset were then algorithmically generated based on their connectivity (both direct and indirect relationships). Each network displays the type of relationship between two gene products, including genes that are not significantly altered in the user's microarray data set. The networks are ranked depending on the number of significantly expressed genes they contain, based on a P-value that indicates the likelihood of the genes in a network being found together due to chance. A score of 2 indicates a 1 in 100 chance that the focus genes of interest were linked in the network by chance rather than a direct biological relationship. Therefore, scores of 2 or higher have at least a 99% confidence level of not being generated by random chance alone .
Quantitative reverse transcription PCR validation
Validation of the microarray data using quantitative rt-PCR
SABiosciences Catalog #
1.4 (1.0 - 1.9)
4.0 (3.0 - 4.9)
1.1 (-1.6 - 2.0)
-4.0 (-8.0 - -4.0)
2.8 (1.8 - 4.3)
5.0 (2.6 - 9.8)
1.2 (-1.1 - 1.6)
1.4 (1.1 - 2.2)
Animals used for immunohistochemical study
Mice used for the immunohistochemistry study were also adult (10-12 week) female C57BL/6 (20-25 g) obtained from Jackson Laboratories. For each time point (3 days and 7 days), the mice were divided into two groups: (1) pneumonectomy (PNY) and (2) sham operated (SHAM - thoracotomy without lung resection), with three animals in each group. Unilateral pneumonectomy or sham thoracotomy were performed as described above, and after recovery, all mice were fed BrdU in drinking water (0.8 mg/ml) between days 0-3 (day 3 time point) or 4-7 (day 7 time point) before euthanasia.
The mice were anesthetized as above at 3 days and 7 days after surgery (PNY or SHAM) then euthanized by cervical dislocation. Following median sternotomy, the pulmonary vasculature was perfused with ice cold Hanks balanced salt solution, the trachea cannulated, and the lungs removed en bloc. Tissue fixation was achieved with intratracheal 10% buffered formalin at 25 cmH20 overnight. The trachea was then ligated, and the lung was embedded in paraffin. Immunofluorescent staining (IF) was performed on 5 μm paraffin sections. Primary antibodies included the monoclonal mouse antibody anti-BrdU (Santa Cruz, dilution 1:100), the monoclonal rabbit antibody anti-S100A4 (AbCam, dilution 1:100), and the monoclonal mouse antibody anti-αSMA (Santa Cruz, dilution 1:100). Tissue sections were deparaffinized and hydrated using standard methods, and antigen retrieval was performed using a citrate buffer (pH 6.0) and microwave heating (5 mins at high, 15 mins at 40% power). Tissues were washed (TBS with 0.1% Tween) three times before a 20 minute protein block (Dako, Carpinteria, CA), and then exposed to the primary antibodies (15-18 hours at 4 degrees Celsius). Detection of the primary antibodies was achieved using donkey anti-mouse Alexa Fluor 594 (red) for BrdU and donkey anti-mouse or donkey anti-rabbit Alexa Fluor 488 (green) for αSMA and S100A4 respectively, both at 1:200 (30 mins at 37 degrees Celsius). The appropriate isotype control assays were also performed; non-specific staining was not observed.
To examine the proliferation of S100A4 positive parenchymal cells during lung regeneration, 20 randomly selected high power fields (400×) were photographed digitally for each sample (3 PNY and 3 SHAM animals at each time point). Cells were counted (averaging 50-100 nucleated cells/HPF) and divided into four categories (nucleated cells (DAPI), S100A4 positive, BrdU positive, and double positive (S100A4+BrdU), and the mean percentage of S100A4 cells/nucleated cells, BrdU cells/nucleated cells and S100A4+BrdU/nucleated cells were obtained. A two-way ANOVA and independent t-tests were performed to test for significance (P < 0.05) between PNY and SHAM, and between 3 day and 7 day time points.
Microarray analysis and validation
Microarray data was collected from pooled lung samples (n = 8) at four time points (6 hr, 1 day, 3 day and 7 day) during post-pneumonectomy lung regeneration from both PNY and SHAM animals. As mentioned in the methods section, the data obtained from these microarrays were analyzed in two different ways. First, data from each time point was combined in a non-parametric replicate analysis generating a time-independent data set, PNY vs. SHAM, with four replicates. In this time-independent analysis to identify consistently regulated genes, 179 genes were identified as differentially expressed between PNY and SHAM (P < 0.001) with fold changes of ≥ 1.5 or ≤ -1.5 (Table S1, additional file 1). Second, global gene expression patterns in the lung were analyzed independently at each of the four time points following pneumonectomy (6 hours, 1 day, 3 day and 7 day). In this time-dependent analysis, 346, 472, 556 and 733 genes were differentially expressed between PNY and SHAM at 6 hour, 1 day, 3 day and 7 day post-pneumonectomy respectively (complete data not shown), with an S-score of ≥ 3 or ≤ -3 (equivalent to P < 0.003). In addition, validation of the microarray data was performed using quantitative rt-PCR. PCR was performed using 8 genes that showed modulated expression across several different areas of interest at the 1 day time point, as well as in the time-independent analysis. The expression patterns observed using qRT-PCR are similar to patterns observed by microarray (see Table 1).
Time-independent transcriptomic patterns during lung regeneration
A complete list of genes with significantly (P < 0.001) altered expression (fold change ≥ 1.5, or ≤ -1.5) in the time-independent analysis of post-pneumonectomy lung regeneration is compiled in Table S1, additional file 1. This table is organized into biological functions that are prevalent during lung regeneration including cell cycle/cell division, DNA synthesis or repair, cell proliferation, extracellular matrix, cytoskeleton, inflammatory, fibrotic, and immune responses, and protein phosphorylation, and miscellaneous biological functions. From this table, several important transcriptomic patterns emerge. First is the significant up-regulation of many cell cycle, cell division, and DNA synthesis-related genes. Many genes involved in cell proliferation are also differentially expressed, including members of Igf1 signaling (up-regulation of Igf1, Igfbp2, Cyr61 and Pappa2 and down-regulation of Igfbp3), as well as Ctgf, Hbegf and Tnfrsf12a. Components of the extracellular matrix including Tnc, Eln, Fbn1, Col3A1, Col5a2 and Vcan are up-regulated, as well as two members of the Adamts metalloproteinase family (Adamts2, Adamts9). Interestingly, genes relating to goblet cell hyperplasia and mucous production are also up-regulated (Clca3, Agr2, Slc26a4), with differential expression of other genes associated with inflammatory, fibrotic or immune responses (up-regulation of Reg3g, Ear11, Retnla, Nappa, and Ccr9; down-regulation of Arg1, CD5L and Alox15).
Time-dependent transcriptomic patterns during lung regeneration
Proliferation of S100A4 positive parenchymal cells following pneumonectomy
The murine response to pneumonectomy represents a unique biological phenomenon among adult mammals. By understanding the molecular signaling that defines alveolar regeneration in this species, novel targets to promote lung regeneration in higher mammals may be revealed. Transcriptomic patterns and cell proliferation data during post-pneumonectomy lung regeneration suggest that mesenchymal cells appear to be critical for rebuilding the infrastructure of the lung in this setting. Their proliferation corresponds to an increase in expression of extracellular matrix structural genes and proteases that are essential for tissue remodeling. While it may seem obvious that mesenchymal cells rebuild the post-pneumonectomy lung, the specific pattern of transcription and associated gene networks have not been systematically evaluated in this model. Patterns of growth factor expression suggest that these cells influence the growth and migration of epithelial and endothelial cells in a paracrine fashion. Furthermore there is an increase in expression of "anti-inflammatory" cytokines which may be important in preventing an over-exuberant innate immune response, and for prevention of fibrosis. Which of these processes are specifically carried out by mesenchymal cells warrants further investigation. However, our initial analysis shows that S100A4+ fibroblast-like cells, but not αSMA+ myofibroblasts, are proliferating cells that contribute to the transcriptomic pattern observed here.
Mesenchymal cell proliferation and immunomodulation during lung regeneration
The transcriptomic patterns suggest that mesenchymal cells are activated during lung regeneration, resulting in increased expression of common extracellular matrix components, in concert with transcriptional control of inflammatory and immunomodulatory pathways. The end result of this well orchestrated process is a non-fibrogenic form of 'wound healing'. The role of fibroblasts and myofibroblasts in alveolar regeneration was further investigated using immunohistochemical analysis of two markers, S100A4 and αSMA that have been previously associated with fibroblasts and myofibroblasts respectively [22–24]. This analysis confirmed that S100A4-positive fibroblastic cells proliferate in the lung parenchyma during lung regeneration (3 days), but are later down-regulated (7 days). Although S100A4 was not up-regulated at the transcriptional level during lung regeneration, this is consistent with regulation via post-transcriptional modulation as seen in other systems . In contrast, αSMA staining was only associated with vessels and not proliferating parenchymal cells, suggesting that myofibroblasts are not highly active in the process of alveolar regeneration. These data suggest that while fibroblast proliferation is important during lung regeneration, subsequent down-regulation of fibroblast proliferation and minimization of myofibroblastic differentiation may be equally important to avoid parenchymal fibrosis.
Persistent proliferation of S100A4-positive fibroblasts has been associated with fibrotic diseases such as murine bleomycin-induced fibrosis, in which S100A4-positive cell numbers peak at 2-3 weeks and are still above baseline at 4 weeks . Therefore, it is crucial to better understand how fibroblast proliferation can participate in lung regeneration without myofibroblastic differentiation and/or subsequent fibrosis. One possibility is that a combination of reduced Tgfβ signaling activity and reduced pro-inflammatory Th-2 cytokine production controls this response. For example, atrial natriuretic peptide (ANP, encoded by Nappa) is up-regulated during post-pneumonectomy lung regeneration, but down-regulated in OVA-induced asthma, a process associated with a fibrotic response . ANP is anti-fibrotic, and acts through inhibition of Tgfβ-induced fibroblast transformation . The combined effects of different growth factors can also modulate the actions of Tgfβ on cell proliferation and collagen production in fibroblast trans-differentiation. For example, Tgfβ generally induces myofibroblast differentiation, with concurrent increased production of αSMA and type I and III collagen . However, the presence of either Egf or Igf can influence fibroblasts to undergo cell proliferation and DNA synthesis (Igf), or differentiation and αSMA production (Egf) . The actions of Igf1 depend largely on its binding to extracellular proteins such as Igfbp2 and Igfbp3 [29–31]. While Igf1 and Igfbp2 are up-regulated during lung regeneration, Igfbp3 is down-regulated. Increased Igfbp3 expression is associated with emphysema  and senescent fibroblasts, where it results in Igf1 sequestration and reduced cell proliferation [33, 34]. Conversely, Igfbp2 binds to extracellular matrix or fibroblasts in the presence of Igf1 or Igf2 and increases their local bioavailability [35, 36]. Since Igf1 signaling is highly invoked during post-pneumonectomy lung regeneration, the actions of Igf1 might contribute to a proliferative rather than fibrotic fibroblast phenotype.
Although fibroblast proliferation is a common feature of lung regeneration and pulmonary fibrotic diseases, the early (6 hour and 1 day) down-regulation of cytokines and chemokines involved in inflammatory cell migration (Figure 2B), and up-regulation of genes involved in T cell development and function (Figure 3B) is unique to lung regeneration. Arginase 1 (Arg1) expressed by Th2-induced fibroproliferative M2 macrophages  is also up-regulated during OVA-induced asthma , but down-regulated during post-pneumonectomy lung regeneration. These data suggest an important role for reduced inflammatory response, immunomodulation, and M1 rather than M2 macrophages in successful lung repair/regeneration. These patterns may represent important mechanisms for balancing cell proliferation and extracellular matrix remodeling without incurring fibrotic scarring resulting from the action of exuberant pro-inflammatory cascades. In contrast to the early anti-inflammatory transcriptomic pattern observed during lung regeneration in this study, another recent study describes up-regulation of select inflammatory cytokines (Hmgb1 and INF-y) during post-pneumonectomy lung regeneration, although most cytokines studied were unaffected by the procedure .
Does adult alveolar regeneration recapitulate lung alveolar development?
Both the observed transcriptomic patterns and the results of the immunohistochemical analysis suggest that adult alveolar regeneration is not identical to lung alveolar development. Myofibroblasts have been previously described as essential effector cells in post-natal secondary septation and alveolarization, with characteristic α-SMA synthesis and localization to developing alveolar septa [39, 40]. In contrast, little αSMA expression is present in the parenchyma during post-pneumonectomy lung regeneration, and immunohistochemical localization of αSMA suggests that activated myofibroblasts are absent from alveolar structures. Igf signaling and mRNA expression patterns are up-regulated during alveolar development [41, 42], a trend which is also apparent in lung regeneration transcriptomic patterns. However, other common signaling pathways such as Fgf [43–45], Pdgf signaling  and HoxA5  have been implicated in alveolar development, but are not significantly present in transcriptomic patterns during lung regeneration. These differences could represent alterations between the ontogeny of alveolar development and post-pneumonectomy lung regeneration, or may reflect the broad experimental design of our study, which may have reduced our ability to detect subtle and/or temporally restricted changes.
Analysis of transcriptomic patterns at four time points during post-pneumonectomy lung regeneration reveals significant insight into the regeneration of normal lung tissue after partial pneumonectomy. By including both time-independent and time-dependent analyzes, the data provides insight into important checks and balances that are activated to facilitate growth of functional lung tissue without fibrogenesis. Analysis of the data set using IPA revealed two themes that are important in the process of lung regeneration. The first is the transcriptomic patterns consistent with activation of mesenchymal cells, and the second is transcriptomic patterns consistent with anti-inflammatory immunomodulatory activity. The presence of proliferating mesenchymal cells in the alveolar parenchyma was also demonstrated immunohistochemically. Although proliferating S100A4+ cells (fibroblasts) have been previously associated with fibrotic scarring, this data demonstrates that S100A4+ cells can actively participate in non-fibrogenic tissue regeneration in the lung . In lung regeneration, it may be the influence of immune modulation and modulation of the inflammatory response that is responsible for balancing cell proliferation and extracellular re-modeling against fibrosis.
Taken together with previous reports examining the mechanism of lung regeneration, we can speculate that from the initial stimulus of increased mechanical stress and hypoxia , post-pneumonectomy lung regeneration occurs through a combination of early immediate gene expression [9, 10], up-regulation of genes important in cell cycle regulation and cell proliferation, and the careful orchestration of fibroblast proliferation, extracellular matrix deposition and immunomodulation to prevent excessive fibrosis. Our study provides a unique description of transcriptomic patterns and mesenchymal cell proliferation during post-pneumonectomy lung regeneration in the adult mouse. However, better understanding of biologically important mechanisms using transgenic models, lineage tagging, and transplant models will be important to further understand the process of lung regeneration.
Our thanks to Alisha Gruntman for performing the pneumonectomy and sham surgeries. Funding was provided through grants awarded to Dr. Edward P. Ingenito (HL 072780-02 and HL 090145-02). We would also like to thank Dr. Tom Mariani (University of Rochester) for his insight and discussions. All microarray experiments and analysis were carried out the Tufts Center for Neuroscience Research funded by Tufts University and NINDS (P30 NS047243).
- Snider GL, Kleinerman JL, Thurlbeck WM, Bengali ZH: The definition of emphysema. Report of a National Heart, Lung, and Blood Institute, Division of Lung Diseases workshop. Am Rev Respir Dis. 1985, 132: 182-185.Google Scholar
- Hsia CC: Signals and mechanisms of compensatory lung growth. J Appl Physiol. 2004, 97: 1992-1998. 10.1152/japplphysiol.00530.2004.PubMedView ArticleGoogle Scholar
- Laros CD, Westermann CJ: Dilatation, compensatory growth, or both after pneumonectomy during childhood and adolescence. A thirty-year follow-up study. J Thorac Cardiovasc Surg. 1987, 93: 570-576.PubMedGoogle Scholar
- Brody JS, Burki R, Kaplan N: Deoxyribonucleic acid synthesis in lung cells during compensatory lung growth after pneumonectomy. Am Rev Respir Dis. 1978, 117: 307-316.PubMedGoogle Scholar
- Fehrenbach H, Voswinckel R, Michl V, Mehling T, Fehrenbach A, Seeger W, Nyengaard JR: Neoalveolarization contributes to compensatory lung growth following pneumonectomy in mice. Eur Respir J. 2008, 31: 515-522. 10.1183/09031936.00109407.PubMedView ArticleGoogle Scholar
- Brown LM, Rannels SR, Rannels DE: Implications of post-pneumonectomy compensatory lung growth in pulmonary physiology and disease. Respir Res. 2001, 2: 340-347. 10.1186/rr84.PubMedPubMed CentralView ArticleGoogle Scholar
- Voswinckel R, Motejl V, Fehrenbach A, Wegmann M, Mehling T, Fehrenbach H, Seeger W: Characterisation of post-pneumonectomy lung growth in adult mice. Eur Respir J. 2004, 24: 524-532. 10.1183/09031936.04.10004904.PubMedView ArticleGoogle Scholar
- Landesberg LJ, Crystal RG: Molecular Response to Pneumonectomy. Lung Development and Regeneration. Edited by: Massaro D, Massaro GD, Chambon P. 2004, New York: Marcel Dekker, 455-View ArticleGoogle Scholar
- Landesberg LJ, Ramalingam R, Lee K, Rosengart TK, Crystal RG: Upregulation of transcription factors in lung in the early phase of postpneumonectomy lung growth. Am J Physiol Lung Cell Mol Physiol. 2001, 281: L1138-49.PubMedGoogle Scholar
- Gilbert KA, Rannels DE: Increased lung inflation induces gene expression after pneumonectomy. Am J Physiol. 1998, 275: L21-9.PubMedGoogle Scholar
- Tremblay L, Valenza F, Ribeiro SP, Li J, Slutsky AS: Injurious ventilatory strategies increase cytokines and c-fos m-RNA expression in an isolated rat lung model. J Clin Invest. 1997, 99: 944-952. 10.1172/JCI119259.PubMedPubMed CentralView ArticleGoogle Scholar
- Nolen-Walston RD, Kim CF, Mazan MR, Ingenito EP, Gruntman AM, Tsai L, Boston R, Woolfenden AE, Jacks T, Hoffman AM: Cellular kinetics and modeling of bronchioalveolar stem cell response during lung regeneration. Am J Physiol Lung Cell Mol Physiol. 2008, 294: L1158-65. 10.1152/ajplung.00298.2007.PubMedPubMed CentralView 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. 10.1073/pnas.0500607102.PubMedPubMed CentralView ArticleGoogle Scholar
- Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.PubMedPubMed CentralView ArticleGoogle Scholar
- Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573: 83-92. 10.1016/j.febslet.2004.07.055.PubMedView ArticleGoogle Scholar
- Zhang L, Wang L, Ravindranathan A, Miles MF: A new algorithm for analysis of oligonucleotide arrays: application to expression profiling in mouse brain regions. J Mol Biol. 2002, 317: 225-235. 10.1006/jmbi.2001.5350.PubMedView ArticleGoogle Scholar
- Kerns RT, Zhang L, Miles MF: Application of the S-score algorithm for analysis of oligonucleotide microarrays. Methods. 2003, 31: 274-281. 10.1016/S1046-2023(03)00156-7.PubMedView ArticleGoogle Scholar
- Kennedy RE, Kerns RT, Kong X, Archer KJ, Miles MF: SScore: an R package for detecting differential gene expression without gene expression summaries. Bioinformatics. 2006, 22: 1272-1274. 10.1093/bioinformatics/btl108.PubMedView ArticleGoogle Scholar
- Kennedy RE, Archer KJ, Miles MF: Empirical validation of the S-Score algorithm in the analysis of gene expression data. BMC Bioinformatics. 2006, 7: 154-10.1186/1471-2105-7-154.PubMedPubMed CentralView ArticleGoogle Scholar
- Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, Tschoeke SK, Miller-Graziano C, Moldawer LL, Mindrinos MN, Davis RW, Tompkins RG, Lowry SF, Inflamm and Host Response to Injury Large Scale Collab. Res. Program: A network-based analysis of systemic inflammation in humans. Nature. 2005, 437: 1032-1037. 10.1038/nature03985.PubMedView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaarteenaho-Wiik R, Paakko P, Sormunen R: Ultrastructural features of lung fibroblast differentiation into myofibroblasts. Ultrastruct Pathol. 2009, 33: 6-15. 10.1080/01913120802608430.PubMedView ArticleGoogle Scholar
- Lawson WE, Polosukhin VV, Zoia O, Stathopoulos GT, Han W, Plieth D, Loyd JE, Neilson EG, Blackwell TS: Characterization of fibroblast-specific protein 1 in pulmonary fibrosis. Am J Respir Crit Care Med. 2005, 171: 899-907. 10.1164/rccm.200311-1535OC.PubMedView ArticleGoogle Scholar
- Strutz F, Okada H, Lo CW, Danoff T, Carone RL, Tomaszewski JE, Neilson EG: Identification and characterization of a fibroblast marker: FSP1. J Cell Biol. 1995, 130: 393-405. 10.1083/jcb.130.2.393.PubMedView ArticleGoogle Scholar
- Ambartsumian N, Klingelhofer J, Grigorian M, Karlstrom O, Sidenius N, Georgiev G, Lukanidin E: Tissue-specific posttranscriptional downregulation of expression of the S100A4(mts1) gene in transgenic animals. Invasion Metastasis. 1998, 18: 96-104. 10.1159/000024502.PubMedView ArticleGoogle Scholar
- Di Valentin E, Crahay C, Garbacki N, Hennuy B, Gueders M, Noel A, Foidart JM, Grooten J, Colige A, Piette J, Cataldo D: New asthma biomarkers: lessons from murine models of acute and chronic asthma. Am J Physiol Lung Cell Mol Physiol. 2009, 296: L185-97. 10.1152/ajplung.90367.2008.PubMedView ArticleGoogle Scholar
- Mohapatra SS: Role of natriuretic peptide signaling in modulating asthma and inflammation. Can J Physiol Pharmacol. 2007, 85: 754-759. 10.1139/Y07-066.PubMedView ArticleGoogle Scholar
- Grotendorst GR, Rahmanie H, Duncan MR: Combinatorial signaling pathways determine fibroblast proliferation and myofibroblast differentiation. FASEB J. 2004, 18: 469-479. 10.1096/fj.03-0699com.PubMedView ArticleGoogle Scholar
- Alami N, Page V, Yu Q, Jerome L, Paterson J, Shiry L, Leyland-Jones B: Recombinant human insulin-like growth factor-binding protein 3 inhibits tumor growth and targets the Akt pathway in lung and colon cancer models. Growth Horm IGF Res. 2008, 18: 487-496. 10.1016/j.ghir.2008.04.002.PubMedView ArticleGoogle Scholar
- Grimberg A, Cohen P: Role of insulin-like growth factors and their binding proteins in growth control and carcinogenesis. J Cell Physiol. 2000, 183: 1-9. 10.1002/(SICI)1097-4652(200004)183:1<1::AID-JCP1>3.0.CO;2-J.PubMedPubMed CentralView ArticleGoogle Scholar
- Grimberg A, Coleman CM, Shi Z, Burns TF, MacLachlan TK, Wang W, El-Deiry WS: Insulin-like growth factor factor binding protein-2 is a novel mediator of p53 inhibition of insulin-like growth factor signaling. Cancer Biol Ther. 2006, 5: 1408-1414.PubMedPubMed CentralView ArticleGoogle Scholar
- Muller KC, Welker L, Paasch K, Feindt B, Erpenbeck VJ, Hohlfeld JM, Krug N, Nakashima M, Branscheid D, Magnussen H, Jorres RA, Holz O: Lung fibroblasts from patients with emphysema show markers of senescence in vitro. Respir Res. 2006, 7: 32-37. 10.1186/1465-9921-7-32.PubMedPubMed CentralView ArticleGoogle Scholar
- Grigoriev VG, Moerman EJ, Goldstein S: Overexpression of insulin-like growth factor binding protein-3 by senescent human fibroblasts: attenuation of the mitogenic response to IGF-I. Exp Cell Res. 1995, 219: 315-321. 10.1006/excr.1995.1234.PubMedView ArticleGoogle Scholar
- Izumi K, Kurosaka D, Iwata T, Oguchi Y, Tanaka Y, Mashima Y, Tsubota K: Involvement of insulin-like growth factor-I and insulin-like growth factor binding protein-3 in corneal fibroblasts during corneal wound healing. Invest Ophthalmol Vis Sci. 2006, 47: 591-598. 10.1167/iovs.05-0097.PubMedView ArticleGoogle Scholar
- Conover CA, Khosla S: Role of extracellular matrix in insulin-like growth factor (IGF) binding protein-2 regulation of IGF-II action in normal human osteoblasts. Growth Horm IGF Res. 2003, 13: 328-335. 10.1016/S1096-6374(03)00092-3.PubMedView ArticleGoogle Scholar
- Arai T, Busby W, Clemmons DR: Binding of insulin-like growth factor (IGF) I or II to IGF-binding protein-2 enables it to bind to heparin and extracellular matrix. Endocrinology. 1996, 137: 4571-4575. 10.1210/en.137.11.4571.PubMedGoogle Scholar
- Strieter RM: What differentiates normal lung repair and fibrosis? Inflammation, resolution of repair, and fibrosis. Proc Am Thorac Soc. 2008, 5: 305-310. 10.1513/pats.200710-160DR.PubMedPubMed CentralView ArticleGoogle Scholar
- Tajima A, Kohno M, Watanabe M, Izumi Y, Tasaka S, Maruyama I, Miyasho T, Kobayashi K: Occult injury in the residual lung after pneumonectomy in mice. Interact Cardiovasc Thorac Surg. 2008, 7: 1114-1120. 10.1510/icvts.2007.170456.PubMedView ArticleGoogle Scholar
- Burri PH: Fetal and postnatal development of the lung. Annu Rev Physiol. 1984, 46: 617-628. 10.1146/annurev.ph.46.030184.003153.PubMedView ArticleGoogle Scholar
- Burri PH: Structural aspects of postnatal lung development - alveolar formation and growth. Biol Neonate. 2006, 89: 313-322. 10.1159/000092868.PubMedView ArticleGoogle Scholar
- Liu H, Chang L, Rong Z, Zhu H, Zhang Q, Chen H, Li W: Association of insulin-like growth factors with lung development in neonatal rats. J Huazhong Univ Sci Technolog Med Sci. 2004, 24: 162-165.PubMedView ArticleGoogle Scholar
- Silva D, Venihaki M, Guo WH, Lopez MF: Igf2 deficiency results in delayed lung development at the end of gestation. Endocrinology. 2006, 147: 5584-5591. 10.1210/en.2006-0498.PubMedView ArticleGoogle Scholar
- Mariani TJ, Reed JJ, Shapiro SD: Expression profiling of the developing mouse lung: insights into the establishment of the extracellular matrix. Am J Respir Cell Mol Biol. 2002, 26: 541-548.PubMedView ArticleGoogle Scholar
- Chailley-Heu B, Boucherat O, Barlier-Mur AM, Bourbon JR: FGF-18 is upregulated in the postnatal rat lung and enhances elastogenesis in myofibroblasts. Am J Physiol Lung Cell Mol Physiol. 2005, 288: L43-51. 10.1152/ajplung.00096.2004.PubMedView ArticleGoogle Scholar
- Powell PP, Wang CC, Horinouchi H, Shepherd K, Jacobson M, Lipson M, Jones R: Differential expression of fibroblast growth factor receptors 1 to 4 and ligand genes in late fetal and early postnatal rat lung. Am J Respir Cell Mol Biol. 1998, 19: 563-572.PubMedView ArticleGoogle Scholar
- Lindahl P, Karlsson L, Hellstrom M, Gebre-Medhin S, Willetts K, Heath JK, Betsholtz C: Alveogenesis failure in PDGF-A-deficient mice is coupled to lack of distal spreading of alveolar smooth muscle cell progenitors during lung development. Development. 1997, 124: 3943-3953.PubMedGoogle Scholar
- Mandeville I, Aubin J, LeBlanc M, Lalancette-Hebert M, Janelle MF, Tremblay GM, Jeannotte L: Impact of the loss of Hoxa5 function on lung alveogenesis. Am J Pathol. 2006, 169: 1312-1327. 10.2353/ajpath.2006.051333.PubMedPubMed CentralView ArticleGoogle Scholar
- Strutz F: The great escape--myofibroblasts in fibrosis and the immune system. Nephrol Dial Transplant. 2008, 23: 2477-2479. 10.1093/ndt/gfn171.PubMedView 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.