An overlapping set of genes is regulated by both NFIB and the glucocorticoid receptor during lung maturation
© Lajoie et al.; licensee BioMed Central Ltd. 2014
Received: 2 December 2013
Accepted: 11 March 2014
Published: 25 March 2014
Lung maturation is a late fetal developmental event in both mice and humans. Because of this, lung immaturity is a serious problem in premature infants. Disruption of genes for either the glucocorticoid receptor (Nr3c1) or the NFIB transcription factors results in perinatal lethality due to lung immaturity. In both knockouts, the phenotype includes excess cell proliferation, failure of saccularization and reduced expression of markers of epithelial differentiation. This similarity suggests that the two genes may co-regulate a specific set of genes essential for lung maturation.
We analyzed the roles of these two transcription factors in regulating transcription using ChIP-seq data for NFIB, and RNA expression data and motif analysis for both. Our new ChIP-seq data for NFIB in lung at E16.5 shows that NFIB binds to a NFI motif. This motif is over-represented in the promoters of genes that are under-expressed in Nfib-KO mice at E18.5, suggesting an activator role for NFIB. Using available microarray data from Nr3c1-KO mice, we further identified 52 genes that are under-expressed in both Nfib and Nr3c1 knockouts, an overlap which is 13.1 times larger than what would be expected by chance. Finally, we looked for enrichment of 738 recently published transcription factor motifs in the promoters of these putative target genes and found that the NFIB and glucocorticoid receptor motifs were among the most enriched, suggesting that a subset of these genes may be directly activated by Nfib and Nr3c1.
Our data provide the first evidence for Nfib and Nr3c1 co-regulating genes related to lung maturation. They also establish that the in vivo DNA-binding specificity of NFIB is the same as previously seen in vitro, and highly similar to that of the other NFI-family members NFIA, NFIC and NFIX.
KeywordsLung development Nr3c1 Glucocorticoid receptor Nfib Regulation of transcription ChIP-seq analysis Expression analysis Motif analysis Transcription factor
Lung development is a complex developmental process initiated by budding of the lungs from the gut endodermal tube, multiple rounds of expansion and branching morphogenesis, and final maturation of the epithelial and endothelial components that comprise the airways, pulmonary circulation, and gas exchange surface [1, 2]. It is the final maturation of the lung epithelial cells that is frequently interrupted by premature birth, leading to both acute and chronic lung disease in premature infants [3, 4]. Here we demonstrate apparently related roles of the Nfib and glucocorticoid receptor (hereafter either Nr3c1 or GR) genes in lung maturation.
Previous studies showed that loss of Nfib resulted in perinatal lethality due to lung immaturity . The lungs of late fetal mice lacking Nfib showed reduced expression of Type I and Type II epithelial markers along with morphological immaturity exemplified by a failure of the formation of saccules, the precursor to the alveolar air exchange region. In addition, excess proliferation of both mesenchymal and epithelial cells is seen in Nfib null lungs. Surprisingly, while the phenotype is clearly related to the failure of epithelial cell maturation, loss of Nfib only in the mesenchymal cells of the lung yields a very similar phenotype , indicating that mesenchymal cells regulate late epithelial maturation through as yet unknown inductive mechanisms .
Prenatal administration of glucocorticoids has been shown to stimulate lung maturation in both mice and premature infants [8–10]. Conversely, deletion of Nr3c1, the gene encoding the glucocorticoid receptor, results in a phenotype remarkably similar to that of loss of Nfib, including excess cell proliferation, failure of saccularization and reduced expression of markers of epithelial differentiation . As with Nfib, loss of Nr3c1 only in the mesenchyme recapitulates much of this phenotype . The similarity in phenotype seen with the loss of either Nfib or Nr3c1, together with the shared cell-type expression requirement suggests that these two genes may co-regulate a specific set of genes essential for lung maturation. We therefore examined the lung genes regulated by Nfib and Nr3c1 and the specific binding targets of NFIB to determine how these genes may cooperate in the regulation of lung maturation.
ChIP-seq shows that NFIB binds to the known NFI motif in mouse fetal lung
Correlation of NFIB binding and expression of nearby genes
We studied the mechanism of transcriptional regulation by NFIB in fetal lung cells using our NFIB ChIP-seq data from E16.5 fetal lung cells and previously published gene expression data from E18.5 fetal lung cells in WT and Nfib-knockout mice .
We first sought for dysregulated genes in the Nfib-KO using a 2-fold expression change threshold and a maximal q-value of 0.05 for selection (see Methods).
We identified 631 genes, of which 412 are down-regulated and 219 are up-regulated. For convenience, we will refer to the down-regulated genes as “NFIB-activated”, and to the up-regulated as “NFIB-repressed”. Of course, we realize that the observed effect could be due to direct or indirect regulation of the gene in question (e.g., via NFIB interacting with another TF).
Number of NFIB-activated/repressed genes with an NFIB ChIP-seq peak near their TSS
If we extend the analysis to binding at up to 10 kbp and 100 kbp from the TSS, neither the NFIB-activated nor NFIB-repressed genes have a number of NFIB ChIP-seq peaks that differs significantly from the number expected by chance (p>0.05, two-tailed Fisher’s exact test).
The lack of evidence of a clear relationship between proximal NFIB binding and gene expression in Table 1 may be due to the fact that the expression data is from a later stage of fetal lung development than the ChIP-seq data (E18.5 vs. E16.5). It is quite possible that the set of genes bound by NFIB changes substantially between E16.5 and E18.5. Another confounding factor is that the gene expression data comes from embryonic lungs where Nfib is deleted from E10, but expression is not measured until E18.5, leaving ample time for compensatory changes in gene expression. In fact among the 631 genes identified as activated or repressed at day E18.5 in the Nfib-KO mouse, 28 are annotated as having “sequence-specific DNA binding transcription factor activity” in the Gene Ontology (GO) database . The changes in expression of these TFs will affect the expression of many genes, so many of the observed dysregulated genes may be indirect rather than direct targets of Nfib. Another possibility is that the majority of regulation by NFIB is via long-distance chromatin looping , but we consider this unlikely given the clear enrichment of NFIB binding events we observe in proximal promoter regions (Figure 1).
Promoters of genes activated by NFIB are enriched in NFIB motifs, but repressed ones are not
In the absence of ChIP-seq data in E18.5 mouse fetal lung, we turned to a motif-based analysis of the relationship between NFIB binding and gene expression. First, we tested for over-representation of putative NFIB binding sites (predicted using our new NFIB motif) in the promoters of NFIB-activated and NFIB-repressed genes (see Methods). We found a significant enrichment in the NFIB-activated genes (p a c t <0.0013), but not in the NFIB-repressed (p r e p >0.3). This suggest that many genes are activated by NFIB through direct interaction, but that repression generally results from indirect regulation.
Association of TF motifs with Nfib targets
p a c t
p r e p
As an additional control, we repeated the complete MAS analysis after replacing mouse promoter sequences by their ortholog from rat or human. For both species, NFI motifs were among the most strongly enriched within the set of NFIB-activated genes. In human, NFIX ranked second (p a c t <10-4), and NFIB ranked fourteenth (p a c t <10-3). In rat, NFIX ranked third (p a c t <10-3) and NFIB ranked tenth (p a c t <10-3). Overall, this strongly support the hypothesis that NFIB activates its targets during lung maturation through direct interactions near the promoter regions.
The presence of our novel, ChIP-derived NFIB motif in gene promoters shows less significant correlation with gene expression than the SELEX-based motifs for NFIA, NFIB and NFIX (Table 2). The new motif (NFIB_MEME) ranks twentieth according to MAS, despite being highly similar to the SELEX-based NFIB motif (Figure 2), from which it differs primarily in the preference for an ‘A’ in the right-most position. This may indicate that the new, ChIP-based motif is slightly less accurate than the SELEX-based motif, which, if true, could be due to numerous reasons. The accuracy of motifs derived from ChIP-seq experiments depends strongly on the number of sequences without the motif presented to the motif discovery algorithm. Such sequences can be due to imperfect antibody specificity or to indirect DNA-binding by the antibody via a protein complex or via chromatin loops bound jointly by the antibody and another DNA-binding protein . None of these issues are present in SELEX experiments, although they suffer from their own limitations. There is no guarantee that the DNA-binding specificity of the protein or DNA-binding domain assayed by SELEX is the same under the in vitro SELEX conditions as it is in in vivo. At any rate, our expression correlation results suggest that the existing SELEX-based NFI-family motifs are at least as accurate as our ChIP-derived motif for NFIB.
Other TFs may contribute to the Nfib-KO phenotype
For each motif in Table 2 we have included the expression fold change of the corresponding gene, if significant (p-value ≤0.05). The increased expression of Nfia in the Nfib-KO further suggests that there may be some compensatory mechanism at play between the two paralogs. Since the DNA-binding motif of NFIA is almost identical to that of NFIB, it is probable that both transcription factors bind to the same regulatory elements. Such an apparently compensatory change in one NFI family member upon loss of another was noted previously in Nfib-KO lungs and suggests some type of homeostatic regulation of total NFI levels . However, these data are not sufficient to indicate whether NFIA acts as an activator or a repressor, or whether the same genes are regulated by both NFIA and NFIB.
In addition to the NFI motifs, we note the large positive MAS of the Meis, SNAI2 and NR2F1 motifs, indicating that they are enriched in the NFIB-activated gene set (i.e. in genes that are under-expressed in the Nfib-KO). Because Snai2, Meis2 and Nr2f1 are over-expressed in the Nfib-KO, it is possible that some genes in our NFIB-activated set are repressed by these genes instead of being directly activated by NFIB. Repressor activity has been documented for each of these three factors [22–24]. For example Snai2, which represses transcription via the recruitment of histone deacetylases to target gene promoters , is known for its antiapoptotic activity and plays a role in epithelial-mesenchymal transition. While neither epithelial-mesenchymal transition nor altered apoptosis seem implicated in the phenotype of Nfib-/- lungs , other epigenetic changes mediated by Snai2 affecting cell proliferation and cell differentiation are clearly possible and will be investigated. Finally, we note the enrichment of the EBF1 motif, which regulates cell differentiation . However, according to the microarray data, the Ebf1 gene is not significantly dysregulated in the Nfib-KO.
Comparison of expression change for some TFs of interest in Nfib - and Nr3c1 -knockout lungs
Nfib and Nr3c1regulate an overlapping set of genes
Nfib-knockout mice show a phenotype very similar to that seen in glucocorticoid receptor (Nr3c1) knockout mice. To ask whether this may be due to a common set of dysregulated genes, we compared our Nfib-KO microarray data with available microarray data for Nr3c1-KO in fetal mice lung at 18.5 . Using the same selection threshold as for the Nfib-KO dataset (2-fold expression change and q-value ≤ 0.05), we identified 158 activated genes and 160 repressed genes by GR.
Regulatory targets of Nfib and Nr3c1 overlap significantly
Overlap fold enrichment
Overlap enrichment p-value
Promoters of genes that are activated by both NFIB and GR are enriched in NFIB and GR motifs
Association of TF motifs with Nfib and Nr3c1 common targets
FC N f i b
p a c t
p r e p
Among the motifs with a larger (but similar) positive MAS, we find the androgen receptor (which is almost identical to GR), the EBFI motif, motifs for two members of the T-box family of transcription factors and the ZNF410 motif. However, none of the mouse genes corresponding to these motifs shows a significant (p- value ≤0.05) expression change in both the Nr3c1 and Nfib knockouts. We also note the large negative MAS of Foxi1 and some estrogen-related receptors (Esrra, Esrrb). Estrogen controls many cellular processes such as growth and differentiation. While Esrrb is over-expressed in the Nr3c1-KO, we identified no such dysregulation in the Nfib-KO. We have no expression data for Foxi1 in Nfib-KO, and this gene is not dysregulated in the Nr3c1-KO. However, some other genes of the same Fox family are dysregulated in both knockouts, such as Foxp2 (over-expressed in both) and Foxn3 (under-expressed in Nr3c1-KO and over-expressed in Nfib-KO). Interestingly, it has been shown that loss of Foxp2 leads to defective postnatal lung alveolarization in mouse .
Similarly to what we did in a previous section, we repeated the MAS analysis using human and rat orthologous sequences. For human, NFIX ranked third (p a c t <10-4), and NFIB ranked nineteenth (p a c t <10-2). For rat, NFIX ranked sixth (p a c t <10-3), and NFIB ranked eighth (p a c t <10-2). This suggests that NFIB activates some NR3C1-activated genes through binding at their promoter sequences. However, NR3C1 did not show a significant enrichment for human, and it ranked only 51 st in rat (p a c t =0.044). These data suggest that in some instances mechanisms other than direct binding of promoter sequences by NR3C1 may mediate co-regulation by NR3C1 and NFIB. Consistent with this finding, previous studies have indicated that some functions of NR3C1 are mediated by mechanisms other than direct NR3C1 binding to DNA. For example, mice defective in DNA-binding by NR3C1 are viable while those deleted for NR3C1 die at birth . Thus it will be important to determine the fraction of co-regulated genes whose expression is regulated by direct binding of NR3C1 versus other indirect mechanisms of regulation.
Regulatory sub-network involving Nfib and Nr3c1
Firstly, as shown in Table 2 and discussed above, Nfib and Nfia transcripts are both repressed by NFIB, which is indicated by links 1 and 2 in Figure 8. For simplicity, and because both NFIA and NFIB are transcription factors and thus can act directly to affect transcription, we depict these as direct interactions, but they may well be indirect. The repression of Nfia by Nfib is further supported by the observation that Nfia is significantly under-expressed (p<10-5) while Nfib is significantly over-expressed (p=0.001) in the Nr3c1-KO (Table 3).
Secondly, both Nr3c1 and Nfib transcripts are significantly repressed by GR (p < 0.001 and p = 0.001, respectively, Table 3), which we indicate by links 3 and 4 in Figure 8. We once again depict these as direct interactions for simplicity.
would involve under-expression of Nr3c1 in the Nfib-KO, which is not observed. We indicate the direct cooperation hypothesis by link 5 in Figure 8.
Finally, since the DNA binding affinities of NFIA and NFIB are highly similar (see motif logos in Figure 5), we infer that NFIA may bind many of the same regulatory elements as NFIB. We therefore hypothesize that Nfia may regulate the genes in G52, which we indicate by link 6 Figure 8. However, it is not clear from the available data whether Nfia acts as a repressor or an activator of these genes.
GO analysis of putative common targets of Nfib and Nr3c1
To identify which biological processes could be activated by Nfib and Nr3c1, we tested for gene-annotation enrichment using GOrilla  and a ranked list approach (see Methods). We found significant enrichment for several general terms including “cell adhesion” (q-value = 0.03), “transport” (q-value = 0.0006) and “immune system process” (q-value = 0.004), but this required the inclusion of genes that are not in the set of 52 putative common targets (i.e. genes that appear under-expressed in both knockouts, but that do not meet our strict fold-change and q-value thresholds, see Methods). This suggests that both Nfib and Nr3c1 may regulate these processes in some manner during lung development. On the other hand, we found significant enrichment for three specific terms where Nfib and Nr3c1 are more likely to play a direct role since the associated genes are among the 52 common targets we have identified: First, “cellular defense response” (q-value = 0.03), with the Tnfrsf4 (tumor necrosis factor receptor superfamily, member 4) and Ncf1 (neutrophil cytosolic factor 1) genes; Second, “regulation of vascular endothelial growth factor signaling pathway” (q-value = 0.02), with the Myoc1c (myosin ic) and Xdh (xanthine dehydrogenase) genes; Third, “regulation of fibroblast proliferation” (q-value = 0.04), with Sphk1 (sphingosine kinase 1), Aqp1 (aquaporin 1) and Cdkn1a (cyclin-dependent kinase inhibitor 1a) genes. We are currently examining these genes in more detail to assess NFIB and GR binding to putative regulatory regions.
Lung maturation is a complex process dependent on cooperation between the mesenchyme, epithelium and endothelial cells to promote vascularization, epithelial differentiation, mesenchyme thinning and lung morphogenesis [1, 2]. Previous studies using transcriptional profiling have identified large numbers of genes whose expression levels change substantially during each of the characteristic stages of mouse lung development: 1) embryonic, 2) pseudoglandular, 3) canalicular, 4) saccular and 5) alveolar [31–33]. For example Mariani et al. demonstrated groups of extracellular matrix genes that exhibit stage-specific expression patterns in mouse lung [32, 33]. However the vast majority of these changes in gene expression represent stage-specific differentiation markers of epithelial or mesenchymal cells which define the cellular phenotype at each stage, but give little information on the regulatory mechanisms that control the differentiation process. Conversely, genetic studies have been instrumental in the discovery of gene regulatory pathways essential for early, and to a lesser extent, later stages of lung development. For example, reciprocal epithelial-mesenchymal inductions have been shown to be essential for both early and late stages of lung development [7, 34, 35]. In addition, while a number of signaling pathways including the FGF [36, 37], Shh [38, 39], BMP [40, 41], Wnt  and TGF [43, 44] pathways are known to signal within and between the epithelium and mesenchyme, the transcriptional networks that respond to and/or generate these pathways and mediate the maturation program remain largely unknown [2, 45].
Both the NFIB and GR transcription factors have been shown previously to regulate the transition from the canalicular to the saccular stage in mouse lung [5, 10, 46, 47]. Indeed mutations in either gene result in a very similar phenotype characterized by an excess of mesenchymal and epithelial cells at E18, severe reduction or failure of saccule formation, severe delay in type I and type II cell differentiation and resultant perinatal lethality with largely non-inflated lungs. More recently it was determined using conditional KO alleles that Nfib and Nr3c1 expression in the lung mesenchyme, not lung epithelium, is essential for normal lung maturation [6, 11]. The similarity in phenotype seen with loss of either Nfib or Nr3c1 initially suggested that these two transcription factors might be cooperating in a conserved pathway of lung maturation. Our determination of a significant overlap between the genes whose expression changes with the loss of either Nfib or Nr3c1 (Table 4 and Figure 6) in lung mesenchyme is consistent with this hypothesis. In addition, the motif analysis showed significant enrichment of NFI and GR binding motifs in genes whose expression decreased in both Nfib-KO and Nr3c1-KO lungs (Figure 7) suggesting cooperative activation of these genes by these factors.
Our computational analysis combining expression, binding and motif data provide the first evidence for Nfib and Nr3c1 co-regulating genes related to lung maturation. Although these data are consistent with the model of Nfib and Nr3c1 function shown in Figure 8, other models are possible. For example, it may be that the overlap in differentially expressed genes in the two mutants (Table 4) is a reflection of the overall phenotype (lung immaturity) rather than being causally dependent on direct co-regulation by Nfib and Nr3c1. To distinguish between such alternate hypotheses it will be necessary to simultaneous assess GR and NFIB binding to mouse lung target genes and assess the frequency of co-occupancy at potential target genes. While GR ChIP-seq data is available for human lung cancer cells , rat pheochromocytoma cells  and mouse adipocytes , because GR binding to DNA is highly dependent on chromatin context and thus shows high cell-type specificity  these data sets are of marginal use at best for determining GR targets in normal lung mesenchyme. In addition, while the ChIP-seq data set used here for NFIB was determined at E16.5, it will be necessary to assess both NFIB and GR binding at multiple stages of lung development to determine the temporal sequence during which NFIB and GR bind to target sites to regulate gene expression. However, even in the absence of such data our determination that both NFIB and GR binding motifs are among those most highly associated with genes whose expression is decreased by loss of either Nfib and Nr3c1 (Figure 8) provides a useful framework for future studies.
E16.5 wild type lungs were minced on ice in PBS, incubated in 1% formaldehyde at room temperature for 10 min., quenched with 0.125M glycine for 5 min. and prepared for chromatin isolation as described previously . Chromatin was sheared to ∼200–500 bp using a Branson Sonifier 250 sonicator and subjected to ChIP using a ChIP assay kit (Upstate Biotechnology) and NFIB antibody (Geneka Biotechnology). Immunoprecipitated chromatin was isolated, the crosslinks reversed, and the isolated DNA was used to prepare a sequencing library and sequenced at the Cornell University DNA Sequencing and Genotyping Laboratory, resulting in 8,717,818 unpaired reads of 42bp.
Reads were aligned to the mouse genome (mm9) using Bowtie2 software  with default parameters. 2,311,332 reads (26.51%) aligned at a unique location and 971,349 (11.14%) aligned to more than one location. Reads with a mapping quality score (MAPQ) below 30 were discarded. The remaining reads were then used to call peaks with the MACS1.4 software , using default parameter values except for the following: -g 1.87e9 (effective genome size) -m 5,30 (MFOLD parameter). The list of peak summits output by MACS is available in Additional file 4. De novo motif discovery was performed on the 100bp sequences centered around each peak summit with repeats masked (mm9, downloaded from UCSC), using the MEME software . DNA sequences were shuffled using the uShuffle software  with parameter k=3 to preserve trinucleotide frequencies.
For Nfib-KO, we used the mRNA expression profiling datasets published in . In this KO strain only exon 2 is deleted from the Nfib gene leaving the remaining exons present for detection in the microarray analysis. These datasets have been produced using Affymetrix arrays and are available at the NCBI Gene Expression Omnibus (GEO) data repository under Accession number GSE24465. The probe intensity signals were normalized using the GC-RMA algorithm.
For Nr3c1-KO, we used the datasets published in . These datasets have been produced using Codelink BioArrays and are available at the ArrayExpress data repository under Accession number E-MEXP-861. The signal of the probes with “Good quality” flag was extracted, log2-transformed and quantile normalized using R.
For both datasets, differential expression and p-values were calculated using the linear model implemented in the R Limma package (Smyth, 2004), and q-values were obtained from the p-values using the Benjamini-Hochberg method (they are referred to as adjusted p-values in the Limma output). For downstream analysis, we considered only the probe set with the most significant result for each gene. The complete results with all probes are available in Additional files 5 and 6. Note that differential expression is expressed as log2(KO/WT) in Additional files 5 and 6, but are simply expressed as ratio (KO/WT) in the main text.
In this paper we will refer to a gene as activated by TF X if its expression decreases at least 2-fold in the knockout of X, with a q-value ≤ 0.05. Similarly, we will say that a gene is repressed by X if its expression increases at least 2-fold in the KO of X, with a q-value ≤ 0.05. These two definitions define three sets of genes:
X-activated genes: K O/W T≤ 0.5 and q-value ≤ 0.05,
X-repressed genes: K O/W T≥ 2 and q-value ≤ 0.05, and
X-non-target genes: All other genes for which we have expression data.
Of course we are using the terms “activated” and “repressed” somewhat loosely here, since these sets of genes will include both direct and indirect regulatory targets of TF X.
We use q-values in the above definition to take multiple testing into account. However, when we refer to a particular gene that has been selected using another criteria (e.g. a motif analysis), we report unadjusted p-values of differential expression, and we say that this gene is dysregulated if its p-value ≤ 0.05, without enforcing any fold change threshold.
Motif enrichment in promoters of target genes
To determine how a given motif M is associated with the regulatory targets of TF X, we measure the enrichment of M in the proximal promoters of genes activated or repressed by X, relative to the promoters of non-target genes. (Note that the motif M can be any motif, not just that of the TF X).
First, we use motif M to scan the proximal promoter of each gene, defined as the region within 1000 bp of its TSS, and save the best score for each promoter. We then determine p a c t , the p-value of a Wilcoxon rank-sum test  with the null hypothesis that the promoter motif scores are no better in the activated genes compared to the non-target genes. Similarly, p r e p is the p-value when we test the scores of the repressed genes compared with the non-target genes. For analysis and plotting purposes, we combine these two p-values into a single motif association score (MAS), whose magnitude is log10(min(p a c t ,p r e p )), and whose sign is positive if the motif M is more significantly enriched in the activated genes, p a c t ≤p r e p , and negative otherwise.
In the current work, we first compute the MAS for Nfib targets using each of the 738 motifs in the compendium of SELEX-based motifs published by Jolma et al.. In a subsequent analysis, we compute the MAS for the common targets of Nfib and Nr3c1. In this case, the “activated” set of genes is simply the intersection of the Nfib-activated and Nr3c1-activated, likewise for the “repressed” set of genes. The set of “non-targets” is the intersection of the Nfib-non-targets and the Nr3c1-non-targets. We note that since p a c t and p r e p are not adjusted for multiple tests, MAS should be viewed as a score rather than as a true statistical confidence measure.
We extracted TSS coordinates from the UCSC browser KnownGene table. When a gene has multiple TSSs, we used the one corresponding to the shortest transcript.
Intersection between Nfib and Nr3c1targets
The size of the overlap z between two independent sets A and B, each sampled without replacement from a set C, follows a hypergeometric distribution with parameters m=|A|, n=|C|-|A| and k=|B|, where |X| denotes the number of elements in set X. If C is the set of all genes, and A and B are two sets of target genes, the probability that their intersection contains z or more genes (under the null model) can be obtained in the R programming language with a call to the hypergeometric distribution function phyper(z-1,m,n,k,lower.tail=FALSE). The expected size of the overlap is equal to m k/(n+m).
All genes were ranked according to their likeliness of being activated by both Nfib and Nr3c1, with the most likely at the top of the list. To achieve this, we used the maximum of the fold-change (KO/WT) observed in either the Nfib-KO or the Nr3c1-KO as a sorting key. This ranked list of genes was used as input to GOrilla with default parameters , and a q-value threshold of 0.05 was used to define significant results.
This work was supported in part by NIH/NHLBI grant HL080624 and NYSTEM grants C026429 and C026714 to RMG. ML and TLB were funded by National Institutes of Health grant R0-1 RR021692-05.
- Maeda Y, Davé V, Whitsett JA: Transcriptional control of lung morphogenesis. Physiol Rev. 2007, 87 (1): 219-244. 10.1152/physrev.00028.2006.PubMedView ArticleGoogle Scholar
- Morrisey EE, Hogan BL: Preparing for the first breath: genetic and cellular mechanisms in lung development. Dev Cell. 2010, 18 (1): 8-23. 10.1016/j.devcel.2009.12.010.PubMed CentralPubMedView ArticleGoogle Scholar
- Popova AP: Mechanisms of bronchopulmonary dysplasia. J Cell Commun Signal. 2013, 7 (2): 119-10.1007/s12079-013-0190-x.PubMed CentralPubMedView ArticleGoogle Scholar
- Roos AB, Berg T, Nord M: A relationship between epithelial maturation, bronchopulmonary dysplasia, and chronic obstructive pulmonary disease. Pulm Med. 2012, 1 (2): 119-Google Scholar
- Steele-Perkins G, Plachez C, Butz KG, Yang G, Bachurski CJ, Kinsman SL, Litwack ED, Richards LJ, Gronostajski RM: The transcription factor gene Nfib is essential for both lung maturation and brain development. Mol Cell Biol. 2005, 25 (2): 685-698. 10.1128/MCB.25.2.685-698.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Hsu Y-C, Osinski J, Campbell CE, Litwack ED, Wang D, Liu S, Bachurski CJ, Gronostajski RM: Mesenchymal nuclear factor IB regulates cell proliferation and epithelial differentiation during lung maturation. Dev Biol. 2011, 354 (2): 242-252. 10.1016/j.ydbio.2011.04.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Shannon JM, Hyatt BA: Epithelial-mesenchymal interactions in the developing lung. Annu Rev Physiol. 2004, 66: 625-645. 10.1146/annurev.physiol.66.032102.135749.PubMedView ArticleGoogle Scholar
- Seckl JR: Prenatal glucocorticoids and long-term programming. Eur J Endocrinol. 2004, 151 Suppl 3: 49-62.View ArticleGoogle Scholar
- Banks BA, Cnaan A, Morgan MA, Parer JT, Merrill JD, Ballard PL, Ballard RA: Multiple courses of antenatal corticosteroids and outcome of premature neonates. Am J Obstet Gynecol. 1999, 181 (3): 709-717. 10.1016/S0002-9378(99)70517-X.PubMedView ArticleGoogle Scholar
- Cole TJ, Blendy JA, Monaghan AP, Krieglstein K, Schmid W, Aguzzi A, Fantuzzi G, Hummler E, Unsicker K, Schütz G: Targeted disruption of the glucocorticoid receptor gene blocks adrenergic chromaffin cell development and severely retards lung maturation. Genes Dev. 1995, 9 (13): 1608-1621. 10.1101/gad.9.13.1608.PubMedView ArticleGoogle Scholar
- Habermehl D, Parkitna JR, Kaden S, Brügger B, Wieland F, Gröne H-J, Schütz G: Glucocorticoid activity during lung maturation is essential in mesenchymal and less in alveolar epithelial cells. Mol Endocrinol. 2011, 25 (8): 1280-1288. 10.1210/me.2009-0380.PubMedView ArticleGoogle Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- Smit A: RepeatMasker.1996-2010. http://www.repeatmasker.org,
- Jolma A, Yan J, Whitington T, Toivonen J, Nitta KR, Rastas P, Morgunova E, Enge M, Taipale M, Wei G, Palin K, Vaquerizas JM, Vincentelli R, Luscombe NM, Hughes TR, Lemaire P, Ukkonen E, Kivioja T, Taipale J: DNA-binding specificities of human transcription factors. Cell. 2013, 152 (1–2): 327-339.PubMedView ArticleGoogle Scholar
- Grant CE, Bailey TL, Noble WS: FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011, 27 (7): 1017-1018. 10.1093/bioinformatics/btr064.PubMed CentralPubMedView ArticleGoogle Scholar
- Landt SG, Marinov GK, Kundaje A, Kheradpour P, Pauli F, Batzoglou S, Bernstein BE, Bickel P, Brown JB, Cayting P, Chen Y, DeSalvo G, Epstein C, Fisher-Aylor KI, Euskirchen G, Gerstein M, Gertz J, Hartemink AJ, Hoffman MM, Iyer VR, Jung YL, Karmakar S, Kellis M, Kharchenko PV, Li Q, Liu T, Liu XS, Ma L, Milosavljevic A, Myers RM, et al: ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 2012, 22 (9): 1813-1831. 10.1101/gr.136184.111.PubMed CentralPubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Kadauke S, Blobel GA: Chromatin loops in gene regulation. Biochimica et Biophysica Acta (BBA)-Gene Regulatory Mechanisms. 2009, 1789 (1): 17-25. 10.1016/j.bbagrm.2008.07.002.View ArticleGoogle Scholar
- Farnham PJ: Insights from genomic profiling of transcription factors. Nat Rev Genet. 2009, 10 (9): 605-616. 10.1038/nrg2636.PubMed CentralPubMedView ArticleGoogle Scholar
- Mahony S, Auron PE, Benos PV: DNA familial binding profiles made easy: comparison of various motif alignment and clustering strategies. PLoS Comput Biol. 2007, 3 (3): 61-10.1371/journal.pcbi.0030061.View ArticleGoogle Scholar
- Gilbert DG: Phylodendron. 1990, http://iubio.bio.indiana.edu/treeapp/phylodendron-doc.html,Google Scholar
- Cobaleda C, Pérez-Caro M, Vicente-Dueñas C, Sánchez-García I: Function of the zinc-finger transcription factor SNAI2 in cancer and development. Annu Rev Genet. 2007, 41: 41-61. 10.1146/annurev.genet.41.110306.130146.PubMedView ArticleGoogle Scholar
- Capdevila J, Tsukui T, Esteban CR, Zappavigna V, Belmonte JCI: Control of vertebrate limb outgrowth by the proximal factorMeis2and distal antagonism of BMPs by gremlin. Mol Cell. 1999, 4 (5): 839-849. 10.1016/S1097-2765(00)80393-7.PubMedView ArticleGoogle Scholar
- Dai K, Hussain MM: NR2F1 disrupts synergistic activation of the MTTP gene transcription by HNF-4αand HNF-1α. J Lipid Res. 2012, 53 (5): 901-908. 10.1194/jlr.M025130.PubMed CentralPubMedView ArticleGoogle Scholar
- Pioli PD, Dahlem TJ, Weis JJ, Weis JH: Deletion of Snai2 and Snai3 results in impaired physical development compounded by lymphocyte deficiency. PloS one. 2013, 8 (7): 69216-10.1371/journal.pone.0069216.View ArticleGoogle Scholar
- Liberg D, Sigvardsson M, Åkerblad P: The EBF/Olf/Collier family of transcription factors: regulators of differentiation in cells originating from all three embryonal germ layers. Mol Cell Biol. 2002, 22 (24): 8389-8397. 10.1128/MCB.22.24.8389-8397.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Bird AD, Tan KH, Olsson PF, Zieba M, Flecknoe SJ, Liddicoat DR, Mollard R, Hooper SB, Cole TJ: Identification of glucocorticoid-regulated genes that control cell proliferation during murine respiratory development. J Physiol. 2007, 585 (1): 187-201. 10.1113/jphysiol.2007.136796.PubMed CentralPubMedView ArticleGoogle Scholar
- Shu W, Lu MM, Zhang Y, Tucker PW, Zhou D, Morrisey EE: Foxp2 and Foxp1 cooperatively regulate lung and esophagus development. Development. 2007, 134 (10): 1991-2000. 10.1242/dev.02846.PubMedView ArticleGoogle Scholar
- Reichardt HM, Kaestner KH, Tuckermann J, Kretz O, Wessely O, Bock R, Gass P, Schmid W, Herrlich P, Angel P, etal: Dna binding of the glucocorticoid receptor is not essential for survival. Cell. 1998, 93 (4): 531-541. 10.1016/S0092-8674(00)81183-6.PubMedView ArticleGoogle Scholar
- Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48-10.1186/1471-2105-10-48.PubMed CentralPubMedView ArticleGoogle Scholar
- Kho AT, Bhattacharya S, Mecham BH, Hong J, Kohane IS, Mariani TJ: Expression profiles of the mouse lung identify a molecular signature of time-to-birth. Am J Respir Cell Mol Biol. 2009, 40 (1): 47-10.1165/rcmb.2008-0048OC.PubMed CentralPubMedView 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 (5): 541-548. 10.1165/ajrcmb.26.5.2001-00080c.PubMedView ArticleGoogle Scholar
- Mariani TJ, Shapiro SD: Application of expression profiling to the developing lung: identification of putative regulatory networks controlling matrix production. Chest. 2002, 121 (3 Suppl): 42-44.View ArticleGoogle Scholar
- Chuang P-T, McMahon AP: Branching morphogenesis of the lung: new molecular insights into an old problem. Trends Cell Biol. 2003, 13 (2): 86-91. 10.1016/S0962-8924(02)00031-4.PubMedView ArticleGoogle Scholar
- Minoo P: Transcriptional regulation of lung development: emergence of specificity. Respir Res. 2000, 1 (2): 109-115. 10.1186/rr20.PubMed CentralPubMedView ArticleGoogle Scholar
- Weaver M, Dunn NR, Hogan B: Bmp4 and Fgf10 play opposing roles during lung bud morphogenesis. Development. 2000, 127 (12): 2695-2704.PubMedGoogle Scholar
- Perl A-KT, Hokuto I, Impagnatiello M-A, Christofori G, Whitsett JA: Temporal effects of sprouty on lung morphogenesis. Dev Biol. 2003, 258 (1): 154-168. 10.1016/S0012-1606(03)00106-4.PubMedView ArticleGoogle Scholar
- Chuang P-T, Kawcak T, McMahon AP: Feedback control of mammalian Hedgehog signaling by the Hedgehog-binding protein, Hip1, modulates Fgf signaling during branching morphogenesis of the lung. Genes Dev. 2003, 17 (3): 342-347. 10.1101/gad.1026303.PubMed CentralPubMedView ArticleGoogle Scholar
- Bellusci S, Furuta Y, Rush MG, Henderson R, Winnier G, Hogan B: Involvement of Sonic hedgehog (Shh) in mouse embryonic lung growth and morphogenesis. Development. 1997, 124 (1): 53-63.PubMedGoogle Scholar
- Weaver M, Yingling JM, Dunn NR, Bellusci S, Hogan B: Bmp signaling regulates proximal-distal differentiation of endoderm in mouse lung development. Development. 1999, 126 (18): 4005-4015.PubMedGoogle Scholar
- Hyatt BA, Shangguan X, Shannon JM: BMP4 modulates fibroblast growth factor-mediated induction of proximal and distal lung differentiation in mouse embryonic tracheal epithelium in mesenchyme-free culture. Dev Dyn. 2002, 225 (2): 153-165. 10.1002/dvdy.10145.PubMedView ArticleGoogle Scholar
- Li C, Xiao J, Hormi K, Borok Z, Minoo P: Wnt5aparticipates in distal lung morphogenesis. Dev Biol. 2002, 248 (1): 68-81. 10.1006/dbio.2002.0729.PubMedView ArticleGoogle Scholar
- Li C, Zhu N-L, Tan RC, Ballard PL, Derynck R, Minoo P: Transforming growth factor-βinhibits pulmonary surfactant protein B gene transcription through SMAD3 interactions with NKX2. 1 and HNF-3 transcription factors. J Biol Chem. 2002, 277 (41): 38399-38408. 10.1074/jbc.M203188200.PubMedView ArticleGoogle Scholar
- Shi W, Heisterkamp N, Groffen J, Zhao J, Warburton D, Kaartinen V: Tgf-β3-null mutation does not abrogate fetal lung maturation in vivo by glucocorticoids. Am J Physiol Lung Cell Mol Physiol. 1999, 277 (6): 1205-1213.Google Scholar
- Xu Y, Wang Y, Besnard V, Ikegami M, Wert SE, Heffner C, Murray SA, Donahue LR, Whitsett JA: Transcriptional programs controlling perinatal lung maturation. PloS one. 2012, 7 (8): 37046-10.1371/journal.pone.0037046.View ArticleGoogle Scholar
- 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. 2004, 30 (5): 613-619. 10.1165/rcmb.2003-0236OC.PubMedView ArticleGoogle Scholar
- Gründer A, Ebel TT, Mallo M, Schwarzkopf G, Shimizu T, Sippel AE, Schrewe H: Nuclear factor IB (Nfib) deficient mice have severe lung hypoplasia. Mech Dev. 2002, 112 (1): 69-77.PubMedView ArticleGoogle Scholar
- Reddy TE, Pauli F, Sprouse RO, Neff NF, Newberry KM, Garabedian MJ, Myers RM: Genomic determination of the glucocorticoid response reveals unexpected mechanisms of gene regulation. Genome Res. 2009, 19 (12): 2163-2171. 10.1101/gr.097022.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Polman JAE, Welten JE, Bosch DS, de Jonge RT, Balog J, van der Maarel SM, de Kloet ER, Datson NA: A genome-wide signature of glucocorticoid receptor binding in neuronal PC12 cells. BMC Neurosci. 2012, 13: 118-10.1186/1471-2202-13-118.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu C-Y, Mayba O, Lee JV, Tran J, Harris C, Speed TP, Wang J-C: Genome-wide analysis of glucocorticoid receptor binding regions in adipocytes reveal gene network involved in triglyceride homeostasis. PLoS One. 2010, 5 (12): 15188-10.1371/journal.pone.0015188.View ArticleGoogle Scholar
- John S, Sabo PJ, Thurman RE, Sung M-H, Biddie SC, Johnson TA, Hager GL, Stamatoyannopoulos JA: Chromatin accessibility pre-determines glucocorticoid receptor binding patterns. Nat Genet. 2011, 43 (3): 264-268. 10.1038/ng.759.PubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9 (4): 357-359. 10.1038/nmeth.1923.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nussbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9 (9): 137-10.1186/gb-2008-9-9-r137.View ArticleGoogle Scholar
- Jiang J, Chan Y-S, Loh Y-H, Cai J, Tong G-Q, Lim C-A, Robson P, Zhong S, Ng H-H: A core Klf circuitry regulates self-renewal of embryonic stem cells. Nat Cell Biol. 2008, 10 (3): 353-360. 10.1038/ncb1698.PubMedView ArticleGoogle Scholar
- Wolfe DA, Hollander M: Nonparametric Statistical Methods. 1973, New York: John WileyGoogle 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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.