Published ahead of print on November 16, 2006, doi:10.1165/rcmb.2006-0314OC
© 2007 American Thoracic Society DOI: 10.1165/rcmb.2006-0314OC NPAS1 Regulates Branching Morphogenesis in Embryonic LungDepartment of Pediatrics, and Department of Pathology, Children's Hospital and Harvard Medical School, Boston, Massachusetts; and Department of Pathology, Duke University Medical Center, Durham, North Carolina Correspondence and requests for reprints should be addressed to Mary E. Sunday, M.D., Ph.D., Department of Pathology, Duke University Medical Center, Research Drive, Carl Building, Room 0043, Durham, NC 27710. E-mail: mary.sunday{at}duke.edu
Drosophila trachealess (Trl), master regulator of tracheogenesis, has no known functional mammalian homolog. We hypothesized that genes similar to trachealess regulate lung development. Quantitative (Q)RT-PCR and immunostaining were used to determine spatial and temporal patterns of npas1 gene expression in developing murine lung. Immunostaining for -smooth muscle actin demonstrated myofibroblasts, and protein gene product (PGP)9.5 identified neuroendocrine cells. Branching morphogenesis of embryonic lung buds was analyzed in the presence of antisense or sense oligodeoxynucleotides (ODN). Microarray analyses were performed to screen for changes in gene expression in antisense-treated lungs. QRT-PCR was used to validate the altered expression of key genes identified on the microarrays. We demonstrate that npas1 is expressed in murine embryonic lung. npas1 mRNA peaks early at Embryonic Day (E)10.5E11.5, then drops to low levels. Sequencing verifies the identity of npas1 transcripts in embryonic lung. NPAS1 immunostaining occurs in nuclei of parabronchial mesenchymal cells, especially at the tracheal bifurcation. Arnt, the murine homolog of Tango (the heterodimerization partner for Trl) is also expressed in developing lung but at constant levels. npas1- or arnt-antisense ODN inhibit lung branching morphogenesis, with altered myofibroblast development and increased pulmonary neuroendocrine cells. On microarrays, we identify > 50 known genes down-regulated by npas1-antisense, including multiple genes regulating cell migration and cell differentiation. QRT-PCR confirms significantly decreased expression of the neurogenic genes RBP-Jk and Tle, and three genes involved in muscle development: ig-h3, claudin-11, and myocardin. Npas1 can regulate myofibroblast distribution, branching morphogenesis, and neuroendocrine cell differentiation in murine embryonic lung.
Key Words: branching morphogenesis myofibroblasts smooth muscle actin cell migration neuroendocrine cells
Lung morphogenesis is the result of complex interactions among multiple transcription factors, growth factors, and signaling pathways (14). Tracheal development in Drosophila melanogaster provides a simple model for investigating analogous genes involved in respiratory system development (5). Trachealess (trl) is the master regulatory gene because null mutation of trl leads to complete absence of the fly respiratory system (6). Trl also regulates downstream branching morphogenesis of the developing tracheal system (6, 7). Trl is a transcription factor (TF) with a basic helix-loop-helix (bHLH) domain and two "PAS" domains (6). Several other PAS domain genes have been cloned, including periodicity (per), single-minded (sim), and aryl hydrocarbon nuclear translocator (arnt), all of which collectively led to the designation of this family of TFs as "PAS" (per-arnt-sim) domain proteins. We became particularly interested in npas1 because this gene is closely homologous to Drosophila trl. Amino acid alignment scores to the Drosophila Trl amino acid sequence on a pam250 matrix are for 393 for mouse Npas1, versus 335 for mouse Sim2 versus 210 for Arnt. We are interested in Arnt because Tango, its fly homolog, is required for Trl function via formation of heterodimers (8, 9), leading to regulation of the migration of the first tracheoblasts. This signaling initiates a developmental cascade critical for Drosophila tracheogenesis (9). Breathless (Drosophila fibroblast growth factor [FGF] receptor [FGFR]), branchless (Drosophila FGF), and sprouty, an FGF signaling antagonist, function downstream from Trl and also regulate tracheal branching (9). Notch, delta, and enhancer of split mutants have pleiotropic effects, including tracheal defects and increased neurogenesis (9). Similarly, pleiotropic signaling mechanisms with a role in both lung and tracheal development include FGF10 (2), the Nkx2.1 (1), and Notch-1 (10).
Functional information is now available for several PAS proteins. Murine NPAS1 (11, 12) is a PAS protein expressed in the central nervous system (CNS), where it regulates complex behaviors (1315). In previous studies in other laboratories, npas1 mRNA was not detected in either adult murine lung or whole embryos before Embryonic Day (E)13 using Northern blot analyses (12). Npas1-null mice did not have any apparent pulmonary phenotype, suggesting that either it has no role in lung development, or that the phenotype may occur only in certain genetic backgrounds, or that there is functional redundancy with other genes. However, there has been no known direct analysis of lung histopathology in npas1-null mice, nor any functional assay related to lung development in these mice. It is known that in vitro functional assays can be more sensitive than phenotypic analysis of knockout mice because short-term gene knock-down does not allow sufficient time for compensatory gene upregulation that can occur in knockout mice (16, 17). Similarly, in reports of npas2 and npas3-null mice, there has been no report of pulmonary abnormalities. Consistent with arnt being ubiquitously expressed, arnt-null mice are embryonic lethal before E14, but the hemizygous arnt-null mice do not have any obvious defects. Hypoxia-inducible factor (HIF)-1 None of the potential TRL homologs has been previously demonstrated in vivo to be associated with a mammalian respiratory phenotype. In the present study, we use a sensitive in vitro functional assay to test the hypothesis that NPAS1 can regulate lung development. This work was presented in part at the International Conferences of the American Thoracic Society in 2003 and 2004 (20, 21).
RNA Analyses Total RNA was isolated from developing lungs from timed-pregnant Swiss-Webster mice (Taconic Laboratories, Germantown, NY) using TRIzol reagent (Invitrogen, Carlsbad, CA). cDNA was synthesized with SuperScript-III First-strand Synthesis System (Invitrogen). For QRT-PCR, primer pairs were designed to span introns (see online supplement for primer sequences) using Primer Express-3.0 (Applied Biosystems, Foster City, CA). For real-time-PCR, cDNA from 200 ng RNA, 0.5 µM primer pair, and 0.5 µM probe (Integrated DNA Technologies, Coralville, IA) were combined in Taqman Universal PCR Master-Mix (Applied Biosystems) and run in ABI-PRISM 7,300 (Applied Biosystems) at 50°C for 2-min, 95°C for 10 min, then 45-cycles at 95°C for 15 s each and 60°C for 1 min. Sequences of primers and probes for QRT-PCR of npas1, sim2, arnt, and -actin are given in Table E1 in the online supplement. Sequences of primers and probes for QRT-PCR of claudin 11, myocardin, -ig-H, RBPsuh, Tle, npas1 and -actin are given in Table E2.
Antisense Oligodeoxynucleotide Treatment of Developing Lungs
Immunohistochemistry
RNA Microarray Analyses The bio-informatics pipeline for microarray data analysis and reporting was accomplished through integration with an analysis/workflow tool called GenePattern (http://www.broad.mit.edu/cancer/software/genepattern/). Normalized data for each array were exported to the Partners HPCGG Bioinformatics Server (http://portal.hpcgg.org), merged with updated gene annotation data for each probe set on the array, and downloaded for further data visualization and analysis.
QRT-PCR for NPAS1, Arnt, and Sim2 To determine temporal patterns and relative levels of gene expression, real-time quantitative RT-PCR (QRT-PCR) was performed for npas1 and arnt using mRNA from murine lung from E10E18 and postnatal day 1 (P1) (Figure 1a). Transcript levels were normalized for -actin and expressed relative to the unit level defined in E11.5 lung buds. Although arnt mRNA levels change minimally during development (Figure 1a), npas1 is expressed in embryonic lung during a specific developmental window. The level of npas1 mRNA is highest at E10.5E11.5, then markedly declines at E12.5 (Figure 1a and data not shown). As a specificity control, we also performed QRT-PCR for sim2, because sim2-null mice have normal lung development (22). Sim2 is expressed at high levels from E10.5E16.5, when mRNA levels drop markedly.
NPAS1 Immunostaining We then performed immunohistochemistry to localize NPAS1 protein in paraffin sections of embryonic lung (Figures 1b1c). E11.5 lung buds were cultured for 48 h with npas1 sense ODN. NPAS1 immunostaining is present in parabronchial mesenchymal cells and a few epithelial cells in the region of branchpoints (Figure 1b, large arrows). At higher magnification, NPAS1 is visualized in nuclei as strong, punctate immunostaining of nucleoli predominantly in mesenchymal cells immediately adjacent to branchpoints (Figure 1c, indicated by long thin arrows). NPAS1 immunostaining of another E11.5 lung bud from the same litter treated with npas1-AS ODN demonstrates abrogation of NPAS1 immunostaining in the parabronchial mesenchyme and the epithelium at the branchpoint (Figure 1d). A lower level of NPAS1 immunostaining remained in the bronchial ganglia (arrow), which is likely to represent relatively higher intracellular levels of NPAS1 in neurons, similar to neuropeptide levels observed using limiting dilution immunoperoxidase analysis (23). As an additional control, a section of the same npas1 sense-treated lung bud shown in Figures 1b and 1c was stained using polyclonal rabbit IgG at the same protein concentration as the affinity-purified rabbit anti-NPAS1. As shown in Figure 1e, this section is devoid of any immunostaining. We also performed immunostaining for NPAS1 at two later time points during fetal lung development. At E15.5, NPAS1 is strongly detected in nuclei of numerous mesenchymal cells, especially in the vicinity of branchpoints, in the cytoplasm of epithelial cells lining large airways (Figure 1f), and in nerve fibers and ganglia (data not shown). Lower levels of NPAS1 immunostaining are detected in cytoplasm of epithelial cells lining primitive alveoli and in nuclei of less than half of the airway epithelial cells. At E18.5, the only strong NPAS1 immunostaining is present in nerve fibers and ganglia, with weak to moderate staining of nuclei of epithelial cells lining the developing fetal airways (Figure 1g). These temporal and spatial patterns of NPAS1 antigen expression are consistent with a major role for this protein as a regulator of epithelialmesenchymal interactions during early lung embryogenesis.
NPAS1 and Branching Morphogenesis To verify that the ODN were taken up by potential target cells in the lung buds, we added Alex-488labeled npas1-antisense or -sense ODN to the cultures. After 24 h of incubation, both AS (Figures 2a and 2b) and sense (Figures 2c and 2d) ODN were concentrated predominantly in the mesenchyme (Figures 2b and 2d).
We then used unlabeled ODN for the functional studies. After 48 h, there is no effect of npas1-sense or arnt-sense ODN on branching morphogenesis (Figures 2e and 2g) compared with untreated control lung buds from the same litter (data not shown), whereas npas1-AS and arnt-AS ODN (Figures 2f and 2h) significantly inhibit branching morphogenesis compared with the corresponding sense controls. Compared with the corresponding sense controls, branching morphogenesis of E11.5 lung buds is decreased by 26% by npas1-AS (Figures 2f and 2i), and also reduced 32% by arnt-AS (Figures 2h and 2i). There is no significant difference in branching morphogenesis between lung buds treated with npas1-sense or arnt-sense and media-alone controls from the same litter (data not shown), in contrast to the abnormal morphology of lungs cultured with npas1- or arnt-antisense ODN (Figures 2f and 2h). These data are consistent with heterodimer formation by NPAS1 and Arnt to form a functional complex to regulate transcription, analogous to TRL and Tango in Drosophila (24). Thus, the effect of arnt-AS may be due to decreased NPAS1 function because Arnt is a ubiquitous dimerization partner for many PAS proteins (8, 9). Sim2-antisense did not alter branching morphogenesis compared with the sim2-sense control, whereas sim2-sense alone decreased branching by 29 ± 2% compared with media alone (P = 0.0003). However, sim2-AS did decrease sim2 mRNA levels by 24% compared with sim2-sense, indicating that sim2-AS does have a specific effect on sim2 mRNA levels, in spite of which there is no effect on branching morphogenesis. Thus, effects of antisense ODN are most valid when antisense groups are compared with the corresponding sense control to correct for nonspecific effects of the ODN (25). We used classical histopathologic methods to compare effects of npas1-AS and arnt-AS on cell proliferation, including both BrdU incorporation and immunostaining for PCNA. Extensive analyses of all sections from multiple experiments demonstrate that cell proliferation is not appreciably altered in npas1- or arnt-AStreated lungs compared with sense controls (data not shown). Similarly, there is no visible change in apoptosis by npas1-antisense treatment according to either in situ TUNEL analysis with the ApopTag kit or histologic detection of TUNEL-positive apoptotic bodies (data not shown).
Myofibroblasts in NPAS1-ASTreated Lung Buds
Neuroendocrine Cells in NPAS1-ASTreated Lung Buds Considering our earlier observations of altered lung branching morphogenesis together with increased numbers of neuroendocrine (NE) cells when Notch signaling is inhibited, we analyzed the expression of an NE cell marker in npas1-AStreated lung buds (10). We also chose to study NE cells because they are known to be the first epithelial cell to differentiate in ontogeny and can be detected after about one week of culture of E11.5 lung buds (26). After 6 d in culture, occasional NE cells can be visualized in npas1-sensetreated controls by immunostaining for the early neural/NE cytoplasmic marker PGP9.5 (Figures 4a and 4b). Other NE markers such as calcitonin gene-related peptide are not yet detectable (data not shown). These NE cells are distributed both proximally and distally in the primitive lung buds, extending from the mainstem bronchi down at least three airway generations. Lung buds treated with media alone without ODN for 6 d also had infrequent NE cells (data not shown). In contrast, lung buds treated with npas1-AS for 6 d (Figures 4c and 4d) have markedly increased numbers of NE cells, occurring both as single cells (thin arrows) and as clusters within the primitive airway epithelium (thick arrows). Occasionally linear arrays of NE cells are seen (double arrow in Figure 4c). In npas1-AStreated lung buds, the NE cells are distributed both proximally and distally, beginning in the mainstem bronchi, and throughout the epithelium for at least three airway generations. Quantification of relative numbers of NE cells in these lung buds normalized for the tissue area of each lung bud is shown in Figure 4e. In lung buds treated for 6 d with npas1-AS ODN, there is a highly significant increase in both the number of PGP+ foci (5.6-fold increased) and the total number of NE cells (9-fold increased) per in2 area of lung bud tissue (P < 0.0001). There is also a significant increase in the number of NE cells per focus (1.7-fold increased, P < 0.012).
mRNA Microarray Analyses To begin to investigate molecular mechanisms underlying the observed morphologic changes, we performed expression profiling of lung buds treated with npas1-AS versus -sense controls from the same litters. RNA was extracted from pooled lung samples from five separate experiments (three to six lung buds per experiment) after 48 h of culture. Microarray analysis was performed by probing five pairs of Affymetrix 430A-2.0 gene array chips (Affymetrix, Union City, CA) using cDNA prepared from npas1-AS versus npas1-sensetreated lung buds. Detailed protocols are given elsewhere (27). The cumulative data were analyzed by excluding all samples with an absent signal in either AS- or sense-treated mRNA samples. The remaining genes were sorted and selected for P < 0.05 and fold-change > 3-fold. We identified a total of 157 genes consistently down-regulated and 11 genes consistently up-regulated in all five experiments (listed in Table E3). Of these genes, 9 are duplicates and 62 are unidentified genes (24 RIKENs, 20 unknowns, 1 open reading frame, and 17 ESTs). Excluding these genes results in 93 genes down-regulated by 1/3 or less by npas1-AS and 5 genes up-regulated by 3-fold or more by npas1-AS. None of the up-regulated genes appear to be of functional significance. At least 13 of the genes down-regulated by npas1-AS have known functions in regulating cell shape, cell adhesion, cell migration, extracellular matrix formation, and/or smooth muscle cell differentiation (Table 1). At least 4 of these genes plus 12 additional genes are known to play a role in normal development and/or angiogenesis, including 5 transcription factors (Table 1) (also listed in Tables E3 and E4). Other down-regulated genes are components of signal transduction cascades, many of which are involved in G-protein signaling, which can regulate cell migration. Finally, some of the genes down-regulated by npas1-AS are likely to play a role in lung development, including sonic hedgehog (Shh) signaling molecules (hhip, Hoxb4, Tbx5, dachshund, IGF-I).
Some genes down-regulated by npas1-AS are known regulators of CNS development, including two neurogenic genes (RBPsuh, suppressor-of-hairless, and Tle, transducin-like enhancer-of-split), which function in Notch-related lateral inhibition of neural/NE differentiation (28). This gene family also plays a major downstream role in tracheal development in Drosophila (9). Decreased expression of neurogenic genes would provide a mechanistic basis for increased pulmonary NE cells in npas1-AStreated lungs.
To validate expression levels of several genes that would likely contribute to the altered morphogenesis that we observed in the lung buds, we performed QRT-PCR. The results of three pooled experiments are given in Figure 5. Npas1 mRNA is reduced
Gene expression for TGF- -induced 68 kD, also termed " -IG-H3" are reduced 30% (P < 0.005). BetaIG-H3 is likely to play a role in lung development, being expressed at high levels in distal lung during organogenesis (32). It is present in the extracellular matrix, cytoplasm, and nuclei of smooth muscle cells in human airways and pulmonary vasculature, especially in the septal tips of alveolar ducts and alveoli, suggesting that it may have a morphogenetic role (33). Claudin-11, also significantly decreased by npas1-AS with P < 0.005, is a transmembrane tight-junction protein that is homologous to Drosophila Megatrachea, which plays an important role in regulating morphogenesis of the Drosophila tracheal system (34). Claudin-11 deficiency could also lead to the distal airspace dilatation that we observe in npas1-AStreated lung buds. Claudin-3 is expressed in chick lung buds during development (35). Claudin-5 localizes to junctions between cardiomyocytes and its deficiency is associated with dilated cardiomyopathy (36), suggesting that decreased claudin-11 could be linked to the abnormal shape and lack of continuity between myofibroblasts along airways of npas1-AStreated lung buds. Finally, down-regulation of myocardin gene expression (P < 0.01) could be consistent with impaired development of smooth muscle along developing airways (37). Myocardin is a transcription factor that has been shown to be necessary and sufficient for smooth muscle cell differentiation (37, 38).
Conclusions Myofibroblast migration is recognized to be a major factor in normal lung morphogenesis, especially during alveolarization (39). Similarly in Drosophila TRL functions by regulating cell migration during formation of the primary tracheal tube and also during branching morphogenesis (6, 9). Previously we demonstrated that notch-1-AStreated lung buds have increased branching morphogenesis and increased numbers of NE cells (10). In Drosophila, Notch is implicated as a negative regulator of tracheal development downstream from TRL (9). The cumulative direct and circumstantial evidence suggests that npas1 may be one functional mammalian homolog of the Drosophila trl gene. Despite high sequence conservation, in the CNS several other PAS proteins have been found to be biologically essential and functionally nonredundant (8). However, functional redundancy is common in developmental/biological processes that are essential for survival. We propose that functional redundancy of multiple signaling pathways activated by npas1 could provide a selective advantage to preserve normal lung structure and function that are crucial for postnatal survival. Thus, our short-term in vitro studies are not incompatible with the apparent lack of a pulmonary phenotype in NPAS1 KO mice. In conclusion, the cumulative evidence indicates that NPAS1 can function as a proximal regulator of lung organogenesis, consistent with NPAS1 being a candidate mammalian homolog of TRL, the primary regulator of respiratory system development in Drosophila.
The authors thank Dr. Vance Morgan, Linette Mookanamparambil, and Deepti Anand in the Harvard Partners Center for Genetics and Genomics for carrying out the microarray analyses. The authors acknowledge the technical assistance of Fan Zhang.
* These authors contributed equally to this paper. This work was supported by NIH grant #2RO1 HL44984 (M.E.S.) and NRSA training grant (B.M.L.). This article has an online supplement, which is accessible from this issue's table of contents at www.atsjournals.org Originally Published in Press as DOI: 10.1165/rcmb.2006-0314OC on November 16, 2006 Conflict of Interest Statement: None of the authors has a financial relationship with a commercial entity that has an interest in the subject of this manuscript. Received in original form August 24, 2006 Accepted in final form October 2, 2006
This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||