Alterations in gene expression in T1α null lung: a model of deficient alveolar sac development
BMC Developmental Biology volume 6, Article number: 35 (2006)
Development of lung alveolar sacs of normal structure and size at late gestation is necessary for the gas exchange process that sustains respiration at birth. Mice lacking the lung differentiation gene T1α [T1α(-/-)] fail to form expanded alveolar sacs, resulting in respiratory failure at birth. Since little is known about the molecular pathways driving alveolar sacculation, we used expression microarrays to identify genes altered in the abnormal lungs and, by inference, may play roles in normal lung morphogenesis.
Altered expression of genes related to cell-cell interaction, such as ephrinA3, are observed in T1α(-/-) at E18.5. At term, FosB, Egr1, MPK-1 and Nur77, which can function as negative regulators of the cell-cycle, are down-regulated. This is consistent with the hyperproliferation of peripheral lung cells in term T1α (-/-) lungs reported earlier. Biochemical assays show that neither PCNA nor p21 are altered at E18.5. At term in contrast, PCNA is increased, and p21 is decreased.
This global analysis has identified a number of candidate genes that are significantly altered in lungs in which sacculation is abnormal. Many genes identified were not previously associated with lung development and may participate in formation of alveolar sacs prenatally.
Lung development starts in mice at embryonic day 9.5 (E9.5). By E16.5, airways have extensively grown and branched to form the bronchial tree. Between E16.5 and term (E20.5) lung cell proliferation is gradually reduced, and the distal lung undergoes significant morphogenetic changes to form the alveolar sacs. While a population of distal epithelial cells flattens, thins, and spreads to form type I cells, other distal epithelial cells remain cuboidal, acquire surfactant filled lamellar bodies and differentiate into type II cells. Differentiation of epithelial cells is accompanied by vascular remodeling and thinning of the mesenchyme, and results in enlargement of the diameter and surface area of the alveolar sacs. Overall this process is known as sacculation, and it is critical to increase the efficiency of fluid absorption and gas exchange processes at birth [1–3].
Very little is known about the molecular regulation of sacculation in normal animals. Abnormal sacculation has been reported in many genetically altered animals carrying null mutations, or transgenes that mis- or over-express growth factors, transcription factors, and other regulatory molecules. These molecular abnormalities result in formation of alveolar spaces that are either too small, as in glucocorticoid receptor (GR) , corticotropin releasing hormone (CRH) , and Sp3 knockout mice , and double p21(+/-)p57(+/-) and p21(-/-)p57(+/-) mice , or too large as in gp330 knockout mouse , the SP-C promoter-Bmp4 mouse , and SP-C promoter-GATA6 mouse . It is interesting that both extremes of alveolar sac size can result in death of the newborn shortly after birth due to respiratory failure. Collectively these observations suggest that formation of alveolar sacs of appropriate dimensions, surface area, and thickness is of fundamental importance in lung organogenesis and is critical for survival.
We have previously shown that mice carrying a null mutation of the T1α gene fail to form expanded alveolar sacs near term and die at birth due to an inability to inflate their lungs with the first few breaths . In normal late fetal and adult lungs, T1α protein is uniquely expressed in the apical membrane of type I alveolar epithelial cells, which form over 90% of the alveolar surface that is specialized for gas exchange [12–14]. In the absence of this protein the alveolar sacs still form but they are narrower than normal and do not properly expand at birth. This abnormality appears to be linked to deficient differentiation of type I cells. This was indicated by the presence of fewer attenuated type I cells and reduced expression of Aqp-5, another type I cell marker gene. Secreted surfactant and surfactant gene expression patterns indicate normal differentiation of type II cells in T1α (-/-) lungs.
Some insights into the process of alveolar sacculation in normal animals come from gene expression microarray data using lungs of normal mice at different developmental time points from embryonic day 9 through postnatal week 4 . This study shows marked changes in gene expression between fetal day 17 and newborn, a period that encompasses the process of sacculation. Among the genes altered are the transcription factors Pod1 and GATA6, the stress-related gene Cyr61, surfactant protein D, and caveolin-1α. However, this important survey was not designed specifically to study sacculation, which would require sampling at more frequent time points near term.
As an alternative approach for understanding the molecular regulation of sacculation, it is reasonable to compare lungs with normal alveolar sac formation to those with altered sac formation and infer from differences in gene expression a set of candidate genes involved in the process. Preliminary expression microarray analysis of the GR knockout lung  has explored the overall changes in gene expression that lead to alterations in alveolar sac formation. This approach has also been used to study the related process of alveologenesis, that involves formation of septae that 'subdivide' alveolar sacs into smaller units or true alveoli .
Here, we use a microarray approach to begin to understand the molecular mechanisms by which formation of alveolar sacs is altered in T1α null mice, animals that serve as a highly reproducible model of altered lung sacculation. We have identified genes that are preferentially altered in the T1α null lung prior to recognizable structural abnormalities, such as ephrinA3, and identified that a number of regulators of the cell cycle are altered at term. Our study also shows that p21, a negative cell cycle regulator, is decreased at term but not at E18.5 at both the mRNA and protein levels, a finding that is likely associated with the failure of sacculation. These genes are now candidates for future studies designed to understand the molecular mechanisms that regulate sacculation in the normal lung.
Altered alveolar sac morphology in the absence of T1α is evident on fetal day E18.5
T1α protein is expressed in the apical surface of lung epithelial cells throughout mouse development (Figure 1). To determine the timing of the initial phenotypic alterations in the absence of T1α, we analyzed the distal lung morphology of T1α(+/+) and T1α(-/-) lungs at E11.5, 15.5, 18.5, and 19.5. Histological analyses showed no differences between T1α(+/+) and T1α(-/-) lungs at E11.5 and E15.5 (Fig. 2A,B,E,F). Gross morphology of the lobes and size of the epithelial buds are similar in (+/+) and (-/-) lungs at E11.5. No significant difference in epithelial tube number and size is apparent at E15.5 (Fig. 2I,L), suggesting that branching morphogenesis is not altered in T1α(-/-) lungs. In normal lung development, thinning of the mesenchyme, epithelial type I cell flattening, and enlargement of the air spaces are observed between E16.5 and term, as the alveolar sacs form. In the absence of T1α, the abnormal phenotype is evident at E18.5 and becomes progressively more severe at E19.5 (Fig. 2C,D,G,H) and at term , as the alveolar sacs are very narrow and irregular and the mesenchyme appears thicker than normal (Fig. 2K,N).
Altered expression of PCNA and p21 genes
To determine whether morphological alterations observed in the absence of T1α at E18.5 are also accompanied by changes in cell proliferation, we compared expression of PCNA and p21 genes by northern and western blot analyses in (+/+), (+/-) and (-/-) lungs at E18.5 and term. Neither PCNA nor p21 protein level is altered at E18.5 in T1α(-/-) lungs compared to T1α(+/+) (Fig. 3A–C). However at term, p21 mRNA is reduced 3.6-fold in (-/-) lungs (p ≤ 0.05) (Fig. 3B) and p21 protein is reduced to 75% of the level in (+/+) lungs (p ≤ 0.05) (Fig. 3C). We have previously shown high numbers of PCNA-stained cells at term, and that PCNA protein is significantly up-regulated ~1.8 fold in (-/-) lungs . This trend is also observed in Fig. 3A (~1.7 fold, Student's t-test p = 0.059). The mRNA levels of other two members of the cdk inhibitor family, p27 and p57 at term are not altered (Fig. 3D). These data suggest that the increase in distal lung cell proliferation in the absence of T1α occurs after morphological alterations are evident, and is probably the consequence of earlier changes in gene expression triggered by the absence of T1α. In the microarray analysis neither p21 nor p27 are detected as present (detection p > 0.05) precluding the comparison to the northern blot analysis. PCNA and p57 are both detected as present in the microarrays (detection p ≤ 0.05) and their levels of expression do not change in the absence of T1α (-/-) compared to (+/+) at E18.5 or at term (data not shown), similarly to the data obtained by northern blots (Figure 3).
Changes in gene expression in T1α null mutant lungs
To identify genes potentially involved in the abnormal alveolar sac formation of the T1α null lungs we analyzed global gene expression profiles of T1α(-/-) vs. T1α(+/+) lungs, at E18.5 and at term, using oligonucleotide microarrays. The data were analyzed using Student's t-test and the more stringent Mann-Whitney analysis . We identified 239 genes that differ between (-/-) and (+/+) lungs at E18.5, with Student's t-test p ≤ 0.05, and median detection p value ≤ 0.05 at least in one genotype (Table S1, Additional file 3). Of those, 156 genes have a Mann-Whitney p value ≤ 0.05 (Table S3, Additional file 3). Genes in this group that are up (+) or down (-) regulated more than 1.5 fold with a Student's t-test p ≤ 0.05 and Mann-Whitney p ≤ 0.05 are shown in Table 1. At term, there are 50 genes that change significantly with a Student's t-test p ≤ 0.05 (Table S2, Additional file 3), Mann-Whitney p = 0.08 (Table S4, Additional file 3). The ability to detect genes that change significantly between T1α (-/-) and (+/+) at term was limited by the smaller sample size at that time point (n = 3). Known genes in this group changing more than 1.5 fold are shown in Table 2. The microarray quality control parameters analyzed, described in the Methods section, indicate that all hybridization experiments are comparable and of good quality.
We used EASE bioinformatic tool to identify over-represented biological processes, defined by the Gene Ontology database (GO), for 239 genes that have altered expression at E18.5 and for the 50 genes altered at term in the absence of T1α (data not shown). Genes involved in vasculogenesis, protein targeting, protein metabolism and mitosis are over-represented at E18.5 with an EASE score ≤ 0.05. Apoptosis, cell proliferation and regulation of transcription are over-represented at term with an EASE score ≤ 0.05. This analysis reflects the molecular events that occur in fetal and term lung in the absence of T1α during alveolar sac formation, and supports the finding that altered proliferation is a late gestation event that may be the result of earlier changes in gene expression in the absence of T1α. These data also correlate with the increased PCNA and reduced p21 protein expression observed at term. Changes in gene expression between T1α(-/-) and T1α(+/+) lungs either at E18.5 or at term are significant (p ≤ 0.05; not corrected for multiple comparisons) but moderate, considering the major morphological changes observed in the absence of T1α. We previously showed that the absence of T1α alters differentiation of type I epithelial cells but does not appear to affect other cell types of the distal lung at term, such as type II cells or endothelial cells . Also, a similar proportion of epithelial to mesenchymal cells in T1α(-/-) and (+/+) lungs is suggested by unchanged levels of expression of the epithelial markers cytokeratin 8 and 18 and the mesenchymal marker collagen 1α1 at E18.5 and term (Figure S1, Additional file 1). Thus, molecular changes in an epithelial cell population that accounts for ~10% of the cells of the lung are likely to provide minimal signals when analyzed in the context of total lung mRNAs or proteins. These small but statistically significant changes are likely to be biologically relevant. Likewise, the molecular alterations in the absence of T1α could be at the post-transcriptional level.
We evaluated by real time RT-PCR (QRT-PCR) changes in expression of selected genes identified by the microarray analyses (Figure 4). We confirmed the trend in expression observed in the microarrays for several genes, making them good candidates for further studies. Many of the genes have been previously described in the lung (see references for Tables 1 and 2, Additional file 2, or the UniGene database,), and are therefore potentially involved in the formation of functional alveolar sacs.
At E18.5 we identified genes that are involved in cell-cell interaction, cell motility, prostaglandin metabolism, and signal transduction (Table 1). One of the genes involved in cell-cell interaction and cell motility is ephrinA3, which is up regulated at the mRNA and protein level (Figures 4A,B and 5A). Another member of the ephrin family, ephrin B2, is down regulated (Table 1). Cytoplasmic dynein, a motor protein involved in microtubule-mediated forward cell migration , is up regulated in the absence of T1α (Figure 4A,B). Genes involved in cell adhesion are also altered, such as the extracellular matrix protein βig-h3, and VCAM1 (Table 1). Ptgs1 (prostaglandin synthase 1; Cox-1) mRNA and protein expression levels are reduced in the absence of T1α (Table 1, and Figure 5B). Overall these changes in gene expression suggest that a different reorganization of the cellular structure is taking place in (-/-) lungs compared to wild type at E18.5.
At term, expression of several genes involved in regulation of cell proliferation is altered (Table 2). Among these, the early response genes FosB, Egr1, MPK-1 and Nur77 are down regulated (Figure 4C,D).
Some genes tested by QRT-PCR showed a variable level of expression in different T1α(-/-) lungs, such as Sp1 and Rb (data not shown). This is likely due to differences in the penetrance of the phenotype  within independent samples isolated in different experiments. The down-regulation of the expression of Cyr61 mRNA observed in the microarray analysis was not validated by other methods and Cyr61 protein level is significantly up-regulated in the absence of T1α (data not shown). This finding underscores the fact that changes in mRNA levels do not always correlate with changes in protein expression.
We have used the T1α knock-out mouse lung as a model of failed alveolar sacculation to begin to discover what genes are involved in the regulation of this critical morphogenetic event before birth. We previously reported that at term, lungs of the T1α knockout mouse display a severe phenotype , with airspaces that are small and often tortuous. The small airspaces cannot be expanded under positive pressure, indicating that the peripheral lung surface area is underdeveloped. Despite this dramatic phenotype, the null lungs at term have many normal features, including a normal number of type II alveolar epithelial cells, normal expression levels of surfactant protein mRNAs (SP-A, -B, -C, and -D), of a bronchiolar cell marker (CCSP), and of selected growth factors including FGF7 and FGF10 . Of the genes tested, only aquaporin 5, another type I cell membrane protein, is reduced at both the mRNA and protein level. PCNA protein, a marker of post-proliferative cells, is increased at birth as was expected given the hyperproliferative phenotype of the term animals.
Several other genetically altered animals have shown a failure of sacculation resulting in lung histology that resembles the T1α deficient mice. The genes targeted in these animals therefore become candidates for key players in the process of sacculation. These include the transcription factors Sp3  and NF1B , hormones and hormone receptors CRH  and GR , and inhibitors of the cell cycle p21/p57 (double p21(+/-)p57(+/-) mice) . Mice deficient in these genes show failed sacculation at term, respiratory distress at birth, and usually neonatal death, while maintaining many normal features of the peripheral lung. The targeted genes that result in dysmorphic sacculation encode proteins of different functional categories, such as transcription factors or membrane receptors. It is noteworthy that null mutations of two important alveolar epithelial proteins, T1α and gp330/megalin [8, 11], expressed in the apical membranes of type I and type II cells respectively, disrupt sacculation. T1α loss inhibits formation of expanded saccules, while loss of gp330/megalin results in abnormally enlarged saccules with cystic or emphysematous-like features .
While the identity of the genes knocked-out in the above animals is obviously known, none of these models of altered sacculation has been explored in detail to define the molecular pathways that lead to their phenotype. T1α, Sp3, GR, CRH, and megalin, for example, are all expressed early in lung embryogenesis and in the lung thereafter [14, 22–25]. In those knock-outs there is considerable uncertainty about the timing of molecular events that, in the end, affect sacculation. One could envision that the accumulation of small molecular anomalies over time might block a major morphogenetic event. To address this question, we show here that the global gene expression patterns of T1α null lungs at gestational day 18.5 and term differ both from controls and from each other.
Comparison of microarray analyses of T1α-null and wild-type lungs on day 18.5 shows changes in expression of genes in several functional categories including transcription factors, cell-cell interaction, enzymes, channels, and others, that are reflected in overrepresentation of several biological processes including vasculogenesis, protein targeting, protein, metabolism, and mitosis. Among the cell-cell interaction genes ephrinA3 is significantly up-regulated in lungs showing narrow alveolar saccules. Ephrins are surface bound proteins that interact with eph tyrosine kinase receptors at sites of cell-cell interactions. The bidirectional signaling generated by this interaction regulates intercellular adhesion, cell shape and cell motility . Ephrins, semaphorins, netrins and slit protein families generate molecular cues to guide cell movement and tissue morphogenesis [ and references therein].
Semaphorins, netrins and slit proteins appear to participate in early lung branching morphogenesis. In contrast there is sparse information about the role of ephrins in the lung. Several ephrinAs and ephA receptors are expressed in adult lung suggesting that they can also participate in cell-cell interaction in this organ. In particular ephrinA3 is highly up-regulated in squamous lung cell carcinomas . The role of EphrinAs has been mostly studied in neurons and vascular formation [29, 30] but it is possible that increased ephrinA3 in T1α(-/-) lungs alters the motility and/or shape of distal lung cells resulting in abnormal alveolar sac morphogenesis, a possibility that will require further evaluation.
Several types of enzymes such as peptidases, ligases, and others are affected at E18.5. One enzyme of particular interest is Ptgs 1 (Cox-1) that is decreased in T1α knockout lungs. Ptgs1 catalyzes the production of PGE2 that is the major prostaglandin species produced in the lung. PGE2 is known to be a negative regulator of lung fibroblast proliferation in adult rodents. While many lung cells can produce PGE2, murine alveolar epithelial cells are immunoreactive for COX-1 protein and adult alveolar type II epithelial cells have been directly shown to produce PGE2. Although the pattern of expression of Ptgs1 in mouse fetal lung is not known , decreased epithelial-derived prostaglandins due to a near absence of flattened type I cells, might affect mesenchymal cells in term T1α(-/-) lungs [32, 33].
Biochemical data at term show that PCNA levels are higher and p21 levels are lower than normal in T1α(-/-) lungs. Microarray analyses of term T1α(-/-) lungs identify additional genes likely involved in the hyperproliferation phenotype of the term peripheral lung, such as FosB, Egr1, MKP-1 and Nur77. These immediate early genes are significantly down regulated. Changes in the expression of these genes were not observed at day 18.5. FosB belongs to the Ap-1 family of proteins regulated after mitogenic stimulation. c-jun/fosB dimers inhibit cell proliferation through binding to the Rb promoter . Egr1 is a transcription factor responsive to a variety of stress and mitogenic signals that inhibits proliferation by activation of p21 [35, 36]. Reduced expression of Egr1 in T1α (-/-) lungs may be involved in the down regulation of p21 transcription observed at term. The transcription factor Nur77 (NGFI-B) is an orphan nuclear receptor that modulates the retinoic acid signaling . MKP-1 is a member of the dual specificity protein phosphatases that inhibits G1 progression in response to growth factors. Decreases in MKP-1 have been associated with increased proliferation . MKP-1 transcription is up regulated by glucocorticoids through the GR receptor . The similar phenotype of GR null lungs and T1α null lungs suggests that MKP-1 could be a common gene altered in these two mutant mice that regulates distal lung differentiation. Although all these changes point to possible abnormalities in cell cycle regulation, exactly how this might work and what the relationship between excess cell division and failure of sacculation is, remain unclear. Based on studies of other cells, decreased p21 would be expected to result in increased cell proliferation in the affected cell population(s). However because the net effects on cell cycle regulation are likely complex and dependent on the context of the affected cells it will be important in future studies to determine which lung cell types show alterations in proliferation related genes. It is possible that proliferation is repressed in some cell types while others are hyperproliferative, and are PCNA positive. An important site of decreased p21 expression is likely the type II cell which, in the T1α knock-out mouse, continues to express PCNA inappropriately near term. Although the sites of p21 expression are not well mapped in the fetal mouse lung, p21 and pro-SP-B proteins colocalize in newborn baboon lung under some conditions , and p21 mRNA is expressed in peripheral lung of fetal mice as shown by in situ hybridization . If nascent type I cells are derived from post-mitotic type II cells, withdrawal from the cell cycle may be required before type I characteristics, including the distinctive morphology, can be expressed. This alone might affect sacculation. This possibility is consistent with the failed sacculation, neonatal respiratory death, and normal levels of surfactant proteins and CCSP in the double p21(-/-)p57(+/-) knockout mouse .
We have hypothesized that changes in gene expression might be few and of small magnitude based on previous northern analyses of null lungs at term . Given the striking phenotype of the term lung, the relatively small change in many of the affected genes identified by microarray analysis, often between 1.5 and 2 fold, is notable. Many of these changes are highly significant, however, and are included in our assembled data because of the low variance between samples. Although a >twofold difference has been conventional in considering genes as functionally important, a number of previous studies have also opted to use lower cut-off levels (1.2–1.5) [41–43] when the changes are highly significant. In the T1α null lung this issue may relate in part to the fact that the type I cell comprises only about 10% of total lung cells; thus even large fold changes in genes expressed by type I cell precursors could be masked by lack of change of the same gene in other cell populations. In this regard it would have been useful to determine the fold changes in Aqp-5 mRNA by microarray because it is expressed only in type I cells in the peripheral lung; however Aqp-5 mRNA is not represented on the arrays we used.
This global analysis of gene expression has identified a number of candidate genes that are significantly altered in lungs in which sacculation is diminished and has linked abnormal sacculation to alterations in the expression of cell cycle regulators. Some of these genes may be involved in the regulation of sacculation per se while others may be expressed as a consequence of failed sacculation. Our future studies will be directed towards further characterization of expression patterns genes identified, such as ephrinA3, Ptgs1 and the cell cycle regulators, as preliminary to test their functions in the sacculation process.
Generation of T1α null mutant [T1α(-/-)] mice was described previously . T1α heterozygous [T1α(+/-)] mice were bred to generate timed pregnancies. The presence of a vaginal plug at noon was considered day 0.5 of pregnancy (E0.5). Fetuses were collected at E11.5, E15.5, E18.5 and E19.5. Term animals were obtained within 10 min of normal delivery. Mice were euthanized following IACUC-Boston University Medical Center approved procedures. Genomic DNA purified from a small piece of tail tissue was used to genotype the animals by PCR. Briefly, the targeting vector and/or endogenous T1α were detected using the following conditions: we amplified ~450 bp of the T1α promoter, absent in the targeted allele, using genomic DNA, Taq polymerase (QIAGEN, Inc., Valencia, CA), solution Q, 200 nM dNTPs, 94°C, 1 min; 50°C, 1 min; 72°C 2 min; 30 cycles; forward oligonucleotide (F) TGAATGAATGTAGGAGCGAGCTG; reverse oligonucleotide (R) TTAGGCTCTGAACCCACAATGG. We also amplified ~750 bp of the targeting sequence containing part of the 5' arm of homology and of the neo sequence using the same genomic DNA, Taq polymerase (QIAGEN), 1 mM MgCl2, 600 nM dNTPs, 94°C, 1 min; 60°C, 1 min; 72°C, 2 min; 30 cycles; forward oligonucleotide CCTTCTTGACGAGTTCTTCTGAGG; reverse oligonucleotide GGTCTTATCAGTTAGCTTGTGG.
Histology and immunohistochemistry
E11.5 whole embryos or E15.5, 18.5, 19.5 and P1 dissected lungs were fixed in 4% paraformaldehyde/PBS pH 7.4 solution and embedded in paraffin as described previously . Sections of T1α(+/+) and (-/-) lungs (6 μm) were stained with H&E to evaluate morphological changes in the absence of T1α as described previously . T1α immunohistochemistry was performed as described by Kotton et al . Briefly, antigen retrieval was performed by heating sections to 90°C in a citric acid buffer (antigen retrieval solution; Vector Laboratories) for 20 minutes, slowly cooling to room temperature, and prior to quenching. Sections were blocked with 1% goat serum in PBS (60 minutes), incubated overnight (4°C) with a monoclonal hamster anti-mouse T1α antibody (Developmental Studies Hybridoma Bank, University of Iowa, Hybridoma #8.1.1,), and treated with HRP-labeled goat anti-hamster IgG (ICN, Costa Mesa, CA) for 30 minutes at room temperature. Signals were amplified using tyramide (TSA-Biotin System, NEN Life Science Products, Boston, MA) according to the manufacturer's protocol and exposed to diaminobenzidine for 6.5 minutes.
RNA isolation, labeling and microarray hybridization
Total RNA from E18.5 and term T1α(+/+) and T1α(-/-) lungs was isolated (TRIzol method, Invitrogen, Carlsbad, CA), re-purified using RNeasy Mini Kit (QIAGEN), and DNase treated on column with DNA-free (Ambion, Inc., Austin TX). First and second strand syntheses were performed as recommended by Affymetrix using 8 μg total RNA as template. Biotin-labeled cRNA was generated using the ENZO BioArray High Yield RNA transcript labeling kit (Affymetrix, Santa Clara, CA). Each biotin-labeled cRNA sample (15 μg) was hybridized (16 hours at 45°C and 60 rpm) in an Affymetrix GeneChip System to the Murine Genome U74 Av2 GeneChip (Affymetrix) that contains ~12,000 sequences of the UniGene database or ESTs . Confocal laser scanning was performed with a G2500A scanner (Agilent Technologies, Palo Alto, CA), to detect the streptavidin-labeled fluor.
Microarray data acquisition
Five E18.5 T1α(+/+) lungs, five E18.5 T1α(-/-) lungs, three term T1α(+/+) lungs and three term T1α(-/-) lungs were analyzed on microarrays. A single weighted mean expression level for each gene was derived using Microarray Suite 5.0 software (MAS 5.0, Affymetrix). Using a one-sided Wilcoxon signed-rank test, the MAS 5.0 software generated a detection p-value for each gene indicating whether or not the transcript was reliably detected. Data from each array was scaled in order to normalize the results for inter-array comparisons.
To filter out arrays of poor quality, several quality control parameters on each array were assessed. We reviewed the scanned image for significant artifacts, and the presence on the array of bacterial genes spiked into the hybridization mix (bioB and bioC). Background and noise measurements were below 100 in all arrays. The average 3' to 5' ratio of the signal for GAPDH was 1.01 ± 0.15 and for β-actin was 1.63 ± 0.09 in E18.5 arrays. In term arrays, the 3' to 5' ratio for GAPDH was 1.63 ± 0.82 and for β-actin was 1.72 ± 0.54. A ratio less than 3 is considered acceptable. An average of 45.4 ± 3.8% probe sets exceeded the probe pair threshold in all arrays (called 'present'). Correlation values were calculated among arrays in the same genotype and developmental group. The average correlation values were: E18.5 T1α(-/-) (5 arrays) R2 = 0.953 +/- 0.023, E18.5 T1α(+/+) (5 arrays) R2 = 0.939 +/- 0.026, term T1α(-/-) (3 arrays) R2 = 0.945 +/- 0.016, and term T1α(+/+) (3 arrays) R2 = 0.931 +/- 0.026.
Microarray data analysis
Prior to performing statistical analyses, we selected genes whose median detection p value was ≤ 0.05 (called 'present') in at least one of the 2 comparative groups (+/+ or -/-). To compare T1α(-/-) lung gene expression vs. T1α(+/+) at E18.5 (n = 5 each) or at term (n = 3 each), a Student's t-test was performed with a p-value less than 0.05 considered significant (Tables S1 and S2, Additional file 3). A Mann-Whitney test  was also performed given that expression levels of many genes did not follow a Gaussian distribution (Tables S3 and S4, Additional file 3). Due to the limited sample size, absolute correction for multiple comparisons was not performed. The magnitude of change in gene expression was calculated as the ratio of the average (arithmetic mean) signal in (-/-) lungs divided by the average signal in (+/+) at each time point and expressed as positive fold change (genes up regulated) or negative fold change (genes down regulated). Functional classification of differentially expressed genes was derived from the Gene Ontology annotation , and current literature. Microarray data have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database . GSE1363 is the GEO accession number for the series object that summarizes the experiments and provides links to all other relevant accessions .
Over-represented functional categories or "themes" altered in the absence of T1α at E18.5 and at term were identified using the EASE application . EASE analyzes a list of Affymetrix ID numbers of the genes under study (input list) and finds over-represented biological "themes" in the list of genes, compared to the total number of genes for each biological theme in the array. The significance of each category is determined by two statistical values that are used to sort the categories, i.e. the standard Fisher exact probability, and a conservative EASE score that identifies robust categories. We considered categories with an EASE score ≤ 0.05 as over-represented.
Real Time RT-PCR
Selected genes differentially expressed in the microarray analysis, not previously linked to alveolar sac formation, were evaluated by real time RT-PCR (QRT-PCR). E18.5 and term total RNA, treated with DNA-free DNase (Ambion), was reversed transcribed (RT) (1 μg in 25 μl) using AMV reverse transcriptase (Promega, Madison, WI). RT reactions diluted 1:32 were analyzed by QRT-PCR in an ABI 7000 instrument (Applied Biosystems, Foster City, CA). Primers were designed using PrimerExpress version 2.0 software (Applied Biosystems) or were obtained from Assays-on-Demand (Applied Biosystems). Reactions were performed in 50 μl using either TaqMan PCR universal master mix or SybrGreen PCR master mix (Applied Biosystems). For all primers and probes, optimal concentrations were found and dissociation curves were performed for primers used in SybrGreen reactions to confirm a single product. The relative concentration of RNA for each gene tested was obtained using calibration curves performed with wild type E18.5 or term total lung RNA (n = 3), respectively. Samples and calibration curves were reverse transcribed and amplified under equal conditions. E18.5 T1α (+/+), E18.5 T1α (-/-), term T1α (+/+), and term T1α (-/-) lungs were analyzed in duplicate (n = 3-6). Data were normalized to β-actin mRNA at E18.5 and at term; the level of expression of β-actin mRNA does not change between T1α(+/+) and T1α(-/-) samples at the respective time point. Oligonucleotides used in SybrGreen reactions were: Additional sex combs like (F) TCTTCCGCCATTTCATCTGA, (R) CACAACATGGACTTGGCCAA; Complement C1qa (F) GTGAAAGGCAATCCAGGCAA, (R) TGGTTGGTGAGGACCTTGTCA; Cyr61 Applied Biosystems Assays-on-Demand; Dual specificity phosphatase 1 (MKP-1) (F) ATGCTTGACACACCCACCAGT, (R) TCATACGCGCATGTCATCG; Dynein cytoplasmic (F) AAATCCACCCCAGCGTTTCTAC, (R) CCCCATTTGTTTTTGAGACGG; Early growth response 1 (F) TTTCAGCCTGAGTCCTTACCCA, (R) TGCATGGTCAGTTCCCAGAAC; Ephrin A3 (F) TGCCCCCTTCCTATGACAGAGAG, (R) GCCCAGAGAGATGTACAGTGCA; FosB (F) CGGCCTCTCCCTTTATCCTTT, (R) TGTCCCCGACAAATCCAGAA; Gas5 (F) TTCATTTGGCTGGCTTGCTT, (R) GGAATTCAAGGCCCCCAAT; Nuclear receptor subfamily 4A1/Nur77 (F) CCAGCCTGGTTTTGAGCTAAG, (R) CCATCTTCTGGAGTGAGGAGGT; OTS8/T1α (F) TGGCGTCTTGTTAGCCATTG, (R) GGCGAGACCTTCCAGAAATC; Prostaglandin-endoperoxide synthase 1 (F) CAAACCTACGTCTACGCCAA, (R) AGTCCATCTGTTCCCTCCAC; Rb1 (F) GAATTCTCAGAGATCTTGCTTCG, (R) AGTGGAGTGTCTGCAGCTTTAAG; Sp1 (F) CACCCAGTGCCAGAGACATATG, (R) CCAATTCTTCCTACCACCAAGC.
Total lung RNA (10 μg) was electrophoresed, blotted and hybridized as described previously . p21 probe was a gift of C. X. Deng, NIDDKD, NIH, Bethesda, MD). p27 and p57 probes were generated by RT-PCR (p27 forward oligonucleotide CGGGATATGGAAGAAGCGAGTC, reverse TTTACGTCTGGCGTCGAAGG; p57 forward oligonucleotide GCCAATGCGAACGACTTCTT, reverse AGAGTTCTTCCATCGTCCGCT). Data were normalized to β-actin .
E18.5 and term lung proteins were extracted using RIPA buffer . Protein extracts (10–20 μg) were electrophoresed in a 12% SDS/polyacrylamide gel, blotted on PVDF membranes, or nitrocellulose for p57, and analyzed under similar conditions using PCNA, p21, p27, and p57 antibodies (Santa Cruz Biotechnologies, Inc., Santa Cruz, CA) as described previously . EphrinA3 was analyzed in a 4% SDS/polyacrylamide gel, blotted on PVDF membranes and analyzed using ephrinA3 H-90 (1:1000) and cox-1 M-20 (1:500) antibodies (Santa Cruz Biothechnologies, Inc.). Histone H1 FL-219 (Santa Cruz), β-actin Clone AC-15 (Sigma, Saint Louis, MO) and tubulin-α (Sigma) were tested but discarded as potential controls for western blot loading because all three proteins show a higher expression in T1α(-/-) lungs compared to T1α(+/+) at term (data not shown). Therefore we used total amount of protein loaded as the normalization value, and Coomasie blue staining of the gel, where noted (n = 2-3, in triplicate).
Kauffman SL: Kinetics of pulmonary epithelial proliferation during prenatal growth of the mouse lung. Anat Rec. 1975, 183: 393-403. 10.1002/ar.1091830304.
O'Hare KH, Townes PL: Morphogenesis of albino rat lung: an autoradiographic analysis of the embryological origin of the type I and II pulmonary epithelial cells. J Morphol. 1970, 132: 69-87. 10.1002/jmor.1051320105.
Williams MC: Alveolar type I cells: molecular phenotype and development. Annu Rev Physiol. 2003, 65: 669-95. 10.1146/annurev.physiol.65.092101.142446.
Cole TJ, Solomon NM, Van Driel R, Monk JA, Bird D, Richardson SJ, Dilley RJ, Hooper SB: Altered Epithelial Cell Proportions in the Fetal Lung of Glucocorticoid Receptor Null Mice. Am J Respir Cell Mol Biol. 2003
Muglia LJ, Bae DS, Brown TT, Vogt SK, Alvarez JG, Sunday ME, Majzoub JA: Proliferation and differentiation defects during lung development in corticotropin-releasing hormone-deficient mice. Am J Respir Cell Mol Biol. 1999, 20: 181-8.
Bouwman P, Gollner H, Elsasser HP, Eckhoff G, Karis A, Grosveld F, Philipsen S, Suske G: Transcription factor Sp3 is essential for post-natal survival and late tooth development. Embo J. 2000, 19: 655-61. 10.1093/emboj/19.4.655.
Zhang P, Wong C, Liu D, Finegold M, Harper JW, Elledge SJ: p21(CIP1) and p57(KIP2) control muscle differentiation at the myogenin step. Genes Dev. 1999, 13: 213-24. 10.1101/gad.13.18.2388.
Willnow TE, Hilpert J, Armstrong SA, Rohlmann A, Hammer RE, Burns DK, Herz J: Defective forebrain development in mice lacking gp330/megalin. Proc Natl Acad Sci USA. 1996, 93: 8460-4. 10.1073/pnas.93.16.8460.
Bellusci S, Henderson R, Winnier G, Oikawa T, Hogan BL: Evidence from normal expression and targeted misexpression that bone morphogenetic protein (Bmp-4) plays a role in mouse embryonic lung morphogenesis. Development. 1996, 122: 1693-702.
Koutsourakis M, Keijzer R, Visser P, Post M, Tibboel D, Grosveld F: Branching and differentiation defects in pulmonary epithelium with elevated Gata6 expression. Mech Dev. 2001, 105: 105-14. 10.1016/S0925-4773(01)00386-0.
Ramirez MI, Millien G, Hinds A, Cao Y, Seldin DC, Williams MC: T1alpha, a lung type I cell differentiation gene, is required for normal lung cell proliferation and alveolus formation at birth. Dev Biol. 2003, 256: 61-72. 10.1016/S0012-1606(02)00098-2.
Ramirez MI, Cao YX, Williams MC: 1.3 kilobases of the lung type I cell T1alpha gene promoter mimics endogenous gene expression patterns during development but lacks sequences to enhance expression in perinatal and adult lung. Dev Dyn. 1999, 215: 319-31. 10.1002/(SICI)1097-0177(199908)215:4<319::AID-AJA4>3.0.CO;2-B.
Wetterwald A, Hoffstetter W, Cecchini MG, Lanske B, Wagner C, Fleisch H, Atkinson M: Characterization and cloning of the E11 antigen, a marker expressed by rat osteoblasts and osteocytes. Bone. 1996, 18: 125-32. 10.1016/8756-3282(95)00457-2.
Williams MC, Cao Y, Hinds A, Rishi AK, Wetterwald A: T1 alpha protein is developmentally regulated and expressed by alveolar type I cells, choroid plexus, and ciliary epithelia of adult rats. Am J Respir Cell Mol Biol. 1996, 14: 577-85.
Bonner AE, Lemon WJ, You M: Gene expression signatures identify novel regulatory pathways during murine lung development: implications for lung tumorigenesis. J Med Genet. 2003, 40: 408-17. 10.1136/jmg.40.6.408.
Kaplan F, MacRae T, Comber J, Gagnon S, Sladek R, Hudson TJ, Sweezey NB: Application of expression microarrays to the investigation of fetal lung development in a glucocorticoid receptor knockout mouse model. Chest. 2002, 121: 90S-10.1378/chest.121.3_suppl.90S.
Clerch LB, Baras AS, Massaro GD, Hoffman EP, Massaro D: DNA microarray analysis of neonatal mouse lung connects regulation of KDR with dexamethasone-induced inhibition of alveolar formation. Am J Physiol Lung Cell Mol Physiol. 2004, 286: L411-9. 10.1152/ajplung.00306.2003.
Troyanskaya OG, Garber ME, Brown PO, Botstein D, Altman RB: Nonparametric methods for identifying differentially expressed genes in microarray data. Bioinformatics. 2002, 18: 1454-61. 10.1093/bioinformatics/18.11.1454.
Dujardin DL, Barnhart LE, Stehman SA, Gomes ER, Gundersen GG, Vallee RB: A role for cytoplasmic dynein and LIS1 in directed cell movement. J Cell Biol. 2003, 163: 1205-11. 10.1083/jcb.200310097.
Grunder A, Ebel TT, Mallo M, Schwarzkopf G, Shimizu T, Sippel AE, Schrewe H: Nuclear factor I-B (Nfib) deficient mice have severe lung hypoplasia. Mech Dev. 2002, 112: 69-77. 10.1016/S0925-4773(01)00640-2.
Assemat E, Vinot S, Gofflot F, Linsel-Nitschke P, Illien F, Chatelet F, Verroust P, Louvet-Vallee S, Rinninger F, Kozyraki R: Expression and role of cubilin in the internalization of nutrients during the peri-implantation development of the rodent embryo. Biol Reprod. 2005, 72: 1079-86. 10.1095/biolreprod.104.036913.
Keegan CE, Herman JP, Karolyi IJ, O'Shea KS, Camper SA, Seasholtz AF: Differential expression of corticotropin-releasing hormone in developing mouse embryos and adult brain. Endocrinology. 1994, 134: 2547-55. 10.1210/en.134.6.2547.
Thompson A, Han VK, Yang K: Differential expression of 11beta-hydroxysteroid dehydrogenase types 1 and 2 mRNA and glucocorticoid receptor protein during mouse embryonic development. J Steroid Biochem Mol Biol. 2004, 88: 367-75. 10.1016/j.jsbmb.2003.12.014.
Zhao C, Meng A: Sp1-like transcription factors are regulators of embryonic development in vertebrates. Dev Growth Differ. 2005, 47: 201-11. 10.1111/j.1440-169X.2005.00797.x.
Kalo MS, Pasquale EB: Signal transfer by Eph receptors. Cell Tissue Res. 1999, 298: 1-9. 10.1007/PL00008807.
Hinck L: The versatile roles of "axon guidance" cues in tissue morphogenesis. Dev Cell. 2004, 7: 783-93. 10.1016/j.devcel.2004.11.002.
Hafner C, Schmitz G, Meyer S, Bataille F, Hau P, Langmann T, Dietmaier W, Landthaler M, Vogt T: Differential gene expression of Eph receptors and ephrins in benign human tissues and cancers. Clin Chem. 2004, 50: 490-9. 10.1373/clinchem.2003.026849.
Boyd AW, Lackmann M: Signals from Eph and ephrin proteins: a developmental tool kit. Sci STKE. 2001, 2001: RE20-
Nakamoto M: Eph receptors and ephrins. Int J Biochem Cell Biol. 2000, 32: 7-12. 10.1016/S1357-2725(99)00096-5.
Bauer AK, Dwyer-Nield LD, Malkinson AM: High cyclooxygenase 1 (COX-1) and cyclooxygenase 2 (COX-2) contents in mouse lung tumors. Carcinogenesis. 2000, 21: 543-50. 10.1093/carcin/21.4.543.
Pan T, Mason RJ, Westcott JY, Shannon JM: Rat alveolar type II cells inhibit lung fibroblast proliferation in vitro. Am J Respir Cell Mol Biol. 2001, 25: 353-61.
Panos RJ, Voelkel NF, Cott GR, Mason RJ, Westcott JY: Alterations in eicosanoid production by rat alveolar type II cells isolated after silica-induced lung injury. Am J Respir Cell Mol Biol. 1992, 6: 430-8.
Milde-Langosch K, Bamberger AM, Methner C, Rieck G, Loning T: Expression of cell cycle-regulatory proteins rb, p16/MTS1, p27/KIP1, p21/WAF1, cyclin D1 and cyclin E in breast cancer: correlations with expression of activating protein-1 family members. Int J Cancer. 2000, 87: 468-72. 10.1002/1097-0215(20000815)87:4<468::AID-IJC2>3.0.CO;2-R.
Ragione FD, Cucciolla V, Criniti V, Indaco S, Borriello A, Zappia V: p21Cip1 gene expression is modulated by Egr1: a novel regulatory mechanism involved in the resveratrol antiproliferative effect. J Biol Chem. 2003, 278: 23360-8. 10.1074/jbc.M300771200.
Krones-Herzig A, Adamson E, Mercola D: Early growth response 1 protein, an upstream gatekeeper of the p53 tumor suppressor, controls replicative senescence. Proc Natl Acad Sci USA. 2003, 100: 3233-8. 10.1073/pnas.2628034100.
Kang HJ, Song MR, Lee SK, Shin EC, Choi YH, Kim SJ, Lee JW, Lee MO: Retinoic acid and its receptors repress the expression and transactivation functions of Nur77: a possible mechanism for the inhibition of apoptosis by retinoic acid. Exp Cell Res. 2000, 256: 545-54. 10.1006/excr.2000.4832.
Li M, Zhou JY, Ge Y, Matherly LH, Wu GS: The phosphatase MKP1 is a transcriptional target of p53 involved in cell cycle regulation. J Biol Chem. 2003, 278: 41059-68. 10.1074/jbc.M307149200.
Kassel O, Sancono A, Kratzschmar J, Kreft B, Stassen M, Cato AC: Glucocorticoids inhibit MAP kinase via increased expression and decreased degradation of MKP-1. Embo J. 2001, 20: 7108-16. 10.1093/emboj/20.24.7108.
O'Reilly MA, Watkins RH, Staversky RJ, Maniscalco WM: Induced p21Cip1 in premature baboons with CLD: implications for alveolar hypoplasia. Am J Physiol Lung Cell Mol Physiol. 2003, 285: L964-71.
Jeong JW, Lee KY, Kwak I, White LD, Hilsenbeck SG, Lydon JP, DeMayo FJ: Identification of murine uterine genes regulated in a ligand-dependent manner by the progesterone receptor. Endocrinology. 2005, 146: 3490-505. 10.1210/en.2005-0016.
Lu J, Qian J, Izvolsky KI, Cardoso WV: Global analysis of genes differentially expressed in branching and non-branching regions of the mouse embryonic lung. Dev Biol. 2004, 273: 418-35. 10.1016/j.ydbio.2004.05.035.
McReynolds MR, Taylor-Garcia KM, Greer KA, Hoying JB, Brooks HL: Renal medullary gene expression in aquaporin-1 null mice. Am J Physiol Renal Physiol. 2005, 288: F315-21. 10.1152/ajprenal.00207.2004.
Kotton DN, Ma BY, Cardoso WV, Sanderson EA, Summer RS, Williams MC, Fine A: Bone marrow-derived cells as progenitors of lung alveolar epithelium. Development. 2001, 128: 5181-5188.
Website title. [http://www.uiowa.edu/~dshbwww/]
Website title. [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL81]
Website title. [http://www.ncbi.nlm.nih.gov/geo/]
Website title. [http://niaid.abcc.ncifcrf.gov/]
We thank Drs. Jerry Brody, Wellington Cardoso, Jining Lu and YuXia Cao for thoughtful comments about the manuscript, and Jennifer Beane, Gang Liu and Norman Gerry for their technical support in the analysis of the microarray data. This work was supported by NHLBI Program Project Grant HL47049, and Doris Duke Charitable Foundation Clinical Scientist Development Award (A.S.).
GM, has previously participated in the generation of the T1α null mutant mice. She obtained all the developmental lung samples for the studies and prepared samples for microarray hybridizations at the microarray core facility at Boston University. She also validated gene expression patterns by real time PCR; AH, performed histology and immunostaining analyses; JW, analyzed patterns of expression of cell cycle regulators by northern and western blot methods; AS participated in the design of the study, analyzed the microarray data, performed the statistical analyses, and provided his expertise to interpret the microarray data; MCW participated in the design of the study, in the interpretation of the data, and in the drafting of the manuscript; MIR participated in the design and coordination of the study, in the interpretation of the data, and in the drafting of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Figure S1. Ratios of epithelial vs. mesenchymal marker genes in T1α (+/+) and (-/-) lungs at E18.5 and term. (PDF 394 KB)
Additional File 3: Supplementary Tables. Table S1. Genes dysregulated in T1α (-/-) vs. (+/+) lungs at E18.5 (p ≤ 0.05). Table S2. Genes dysregulated in T1α (-/-) vs. (+/+) lungs at term (p ≤ 0.05). Table S3. Mann-Whitney p values for genes dysregulated in E18.5 T1α (-/-) lungs compared to wild type, with Student's t test p ≤ 0.05. Table S4. Mann-Whitney p values for genes dysregulated in term T1α (-/-) lungs compared to wild type, with Student's t test p ≤ 0.05. (PDF 99 KB)
About this article
Cite this article
Millien, G., Spira, A., Hinds, A. et al. Alterations in gene expression in T1α null lung: a model of deficient alveolar sac development. BMC Dev Biol 6, 35 (2006). https://doi.org/10.1186/1471-213X-6-35