Transcriptomic analysis comparing mouse strains with extreme total lung capacities identifies novel candidate genes for pulmonary function

Background Failure to attain peak lung function by early adulthood is a risk factor for chronic lung diseases. Previously, we reported that C3H/HeJ mice have about twice total lung capacity (TLC) compared to JF1/MsJ mice. We identified seven lung function quantitative trait loci (QTL: Lfnq1-Lfnq7) in backcross/intercross mice derived from these inbred strains. We further demonstrated, superoxide dismutase 3, extracellular (Sod3), Kit oncogene (Kit) and secreted phosphoprotein 1 (Spp1) located on these Lfnqs as lung function determinants. Emanating from the concept of early origin of lung disease, we sought to identify novel candidate genes for pulmonary function by investigating lung transcriptome in C3H/HeJ and JF1/MsJ mice at the completion of embryonic development, bulk alveolar formation and maturity. Methods Design-based stereological analysis was performed to study lung structure in C3H/HeJ and JF1/MsJ mice. Microarray was used for lung transcriptomic analysis [embryonic day 18, postnatal days 28, 70]. Quantitative real time polymerase chain reaction (qRT-PCR), western blot and immunohistochemical analysis were used to confirm selected differences. Results Stereological analysis revealed decreased alveolar number density, elastin to collagen ratio and increased mean alveolar volume in C3H/HeJ mice compared to JF1/MsJ. Gene ontology term “extracellular region” was enriched among the decreased JF1/MsJ transcripts. Candidate genes identified using the expression-QTL strategy include: ATP-binding cassette, sub-family G (WHITE), member 1 (Abcg1), formyl peptide receptor 1 (Fpr1), gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1); histocompatibility 2 genes: class II antigen E beta (H2-Eb1), D region locus 1 (H2-D1), and Q region locus 4 (H2-Q4); leucine rich repeat containing 6 (testis) (Lrrc6), radial spoke head 1 homolog (Rsph1), and surfactant associated 2 (Sfta2). Noteworthy genes selected as candidates for their consistent expression include: Wnt inhibitor factor 1 (Wif1), follistatin (Fst), chitinase-like 1 (Chil1), and Chil3. Conclusions Comparison of late embryonic, adolescent and adult lung transcript profiles between mouse strains with extreme TLCs lead to the identification of candidate genes for pulmonary function that has not been reported earlier. Further mechanistic investigations are warranted to elucidate their mode of action in determining lung function. Electronic supplementary material The online version of this article (doi:10.1186/s12931-017-0629-3) contains supplementary material, which is available to authorized users.


Background
The origin of many chronic lung diseases may be rooted early in life because of failure to achieve optimal peak lung function by early adulthood [1][2][3]. Lung dysfunction may develop in utero and during early life suggesting an underlying impairment of lung development, growth and maturation [4]. For example, persistent wheezing in the early years and lower lung function by age 6 years among children is associated with later onset of asthma. Similarly, following a cohort over 22 years, Lange et al. [5] reported that low forced expiratory volume in 1 s (FEV 1 ) before age 40 as a risk factor for chronic obstructive pulmonary disease (COPD). Diminished lung function and hindered lung development is likely due to multiple environmental (e.g., maternal smoking) and genetic factors [6][7][8][9][10][11][12][13][14][15][16][17][18]. Thus, accumulating evidence suggests that early life determinants of lung function are possible determinants of obstructive airways diseases later in life. Therefore, in this study we sought to identify genetic determinants of lung function development that could provide important insight into the predisposition of chronic lung disease.
In this study, we examined differences in lung architecture and lung transcripts in these highly divergent mouse strains. Transcriptomic analysis was performed at three strategic times covering completion of embryonic lung development, bulk alveolar formation and growth. Morphometric studies in the adult mouse lung were performed to identify any structural alterations arising from developmental differences between C3H/HeJ and JF1/ MsJ mice. Differential transcript expression was also used to identify genes associated with TLC QTLs located within Lfnq3 on mCh15 (40.3-72.5 Mbp) and Lfnq4 on mCh17 (10.9-39.0 Mbp) associated with TLC located within Lfnq3 on mCh15  and Lfnq4 on mCh17 .

Mice
All procedures were approved by the Bavarian Animal Research Authority, Animal Research Authority of Schleswig-Holstein, and University of Pittsburgh, PA. Frozen and paraffin embedded tissues were procured to carry out experiments at SRM University, India according to the Institutional Animal Ethics Committee (IAEC) permission. Mouse strains C3H/HeJ (#000659) and JF1/ MsJ (#003720) were purchased from Jackson laboratory (Bar Harbor, ME, USA) and housed in specific pathogen free conditions. Food and water were available ad libitum. Three developmental stages, namely, embryonic day (E) 18 (completion of embryonic lung development), postnatal day (P) 28 (completion of bulk alveolar formation) and P70 (completion of lung growth and maturity) were selected for the microarray studies. E18 embryos (days post coitum) were stage matched. Mice were anesthetized (4 mg/kg xylene and 188 mg/kg ketamine i.p.) and the posterior aorta was severed. To obtain tissue for mRNA and protein analysis (n = 5 mice/strain/stage, female), the diaphragm was punctured, and the chest cavity opened. Lungs were excised, frozen in liquid nitrogen, and stored (−80°C). Whole lung was used for RNA extraction in case of stage E18. The left lobe of the lung was used for RNA extraction and the right lobe was used for protein extraction for P28 and P70.

Design-based stereological analysis of lung structures and lung function measurements
Comparison of age and sex matched adult C3H/HeJ and JF1/MsJ mice was performed according to American Thoracic Society/European Respiratory Society guidelines [28]. Use of a cascade design allowed linking data retrieved at the light microscopic level using a PC-based software tool (newCAST, Visiopharm, Hersholm, Denmark) with the data obtained at the transmission electron microscope (TEM) [29]. Alveolar numbers were quantified as described previously [30]. Lungs were fixed (20 cm H 2 0) in situ by intra-tracheal instillation of a mixture of 1% glutaraldehyde (Serva, Germany), 1% paraformaldehyde in 0.1 M sodium cacodylate solution (Merck, Germany) [31]. After 20 min, the trachea was ligated and the lungs were fixed overnight (4°C). Lung volume was determined by fluid displacement. To obtain systematic uniform random (SUR) samples, lungs were embedded (2% agar-agar) and sectioned (2 mm). Tissues were randomly embedded into glycolmethacrylate (GMA, Technovit 7100, HeraeusKulzer, Germany) for light microscopy or epoxy resin (Araldite, Serva, Germany) for transmission electron microscopy. For light microscopy, samples were stained with 1% OsO 4 in 0.1 M sodium cacodylate solution (Merck, Germany), followed by half-saturated aqueous uranyl acetate (Agar Scientific, UK), acetone dehydrated GMA embedded, and sectioned (2 μm). For TEM analysis, SUR sampling of tissue blocks (1x1x1 mm 3 ) was performed. Tissue blocks were post-fixed with 1% OsO 4 in 0.1 M sodium cacodylate solution followed by 3% aqueous tannic acid (Mallinckrodt, USA) staining for elastin [32].

Lung function measurements
Lung functions were measured twice with air or twice with saline filled lung as described previously [33]. Saline enables exclusion of surface tension. Lungs were filled to TLC (30 cm H 2 O), emptied in 0.1 ml steps, and volume and corresponding plateau pressures were determined after each step. C L was determined as the linear slope of the pressure-volume curve.

Transcript and protein analyses
Transcriptome-wide analysis of steady state mRNA levels was analyzed by microarray (n = 5/strain/stage, GeneChip®-MouseExon 1.0 ST Array Cat. #902473, Affymetrix, Santa Clara, CA). Selected candidates were confirmed by quantitative real time -polymerase chain reaction (qRT-PCR) and western blot as described previously [21,25]. Total RNA was amplified (GeneChip®Whole Transcript (WT) Sense Target Labeling Assay) without rRNA depletion. Expression console (v.1.3.0.187, Affymetrix) was used for quality control and to obtain annotated normalized RMA data (setting: Gene Level -Extended: RMA-Sketch). The dataset was filtered for probesets with an annotation for gene symbols and only the probeset with the highest number of probes per gene was used in statistical analysis. Array data has been submitted to the Genome Expression Omnibus (GEO) database at National Center for Biotechnology Information NCBI (GSE80078). For qRT-PCR, Quantitech primers (Qiagen) for Wnt inhibitory factor 1 (WIF1 Cat. # QT01065848), frizzled homolog 6 (FZD6 Cat. # QT00109998), and beta actin (ACTB Cat. # QT00095242) were used. ACTB was used as the reference control. For western blot analysis primary antibody for chitinase-like 3 (CHI3L3, Cat. # ab192029, Abcam) and horse radish peroxidase conjugated goat anti-rabbit antibody (Cat. # ab97051, Abcam) was used as the secondary antibody. Clarity Western ECL substrate (Cat. # 170-50,605; Biorad) was used for immunodetection (Syngene G: Box). The blots were stripped using Restore western blot stripping buffer (Cat # 21059, Thermofisher, Waltham, MA) for reprobing with ACTB primary antibody (Cat. # ab8227, Abcam).

Statistics
Two-way ANOVA followed by all pairwise multiplecomparison procedures (SigmaStat 11.0 software; Systat Software Incorporation) were used to compare group means. P < 0.05 was considered as significant in all cases. For microarray data, statistical analyses were performed by utilizing the statistical programming environment R [34] implemented in CARMAweb [35]. Genewise testing for differential expression used t-test and Benjamini-Hochberg multiple testing correction [false discovery rate (FDR) ≤ 0.10]. Gene ontology (GO) term and pathway enrichment analyses were done with Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 (p < 0.05, FDR ≤ 0.05) [36].

Design-based stereological analysis of lung architecture
Previously, we noted that mice C3H/HeJ mice had increased TLC and TLC/Bw than age and sex-matched JF1/MsJ mice [20]. To assess these differences further, we employed design-based stereological analysis to examine the lung architecture of each strain. The light and TEM comparison of age-and sex-matched lungs of JF1/MsJ and C3H/HeJ mice revealed significantly different alveolar architecture at an age of 12-14 weeks ( Table 1, Additional file 1: Figure S1). Interestingly, C3H/HeJ mice had increased mean chord length and mean alveolar volume but decreased alveolar number density and elastin to collagen ratio compared to JF1/ MsJ mice. These structural differences are consistent with the functional measures of increased C Lair , C Lsaline, and specific C Lsaline and ratio of C L-saline /C L-air in C3H/ HeJ compared to JF1/MsJ mice. Increased compliance could result from decreased elasticity and surface tension of the larger alveoli of C3H/HeJ lung. Because alveolar number and size were complementary, the lungs of both mouse strains exhibit near equal total alveolar surface area. To achieve higher alveolar surface area for adequate gas exchange, the smaller lung volume of JF1/ MsJ mice would require more alveolar septa.

Evaluation of selected transcripts
Noteworthy transcripts altered in the later stages include several members of the wingless-type MMTV integration site family, member 1 (WNT1) protein family including increased Wnt inhibitor factor 1 (WIF1) and follistatin (FST) and decreased frizzled homolog 6 (FZD6) in JF1/MsJ compared to C3H/HeJ lungs. Decreased transcripts also included chitinase-like 1 (CHIL1), CHIL3, and hemolytic component (HC). qRT-PCR analysis of lung WIF1 was consistent with the microarray results ( Fig. 1a, b). Western blot analysis for CHIL3 in the total lung protein homogenate was also consistent with the microarray results ( Fig. 1c). Expression of the WNT pathway proteins (WIF1, FST and FZD6) were assessed by immunohistochemistry (Figs. 2, 3, 4). Increased immunoreactive WIF1 ( Fig. 2) and FST ( Fig. 3) were detected in the alveolar epithelial type I and II cells, bronchial epithelial cells, and macrophages of JF1/MsJ compared to C3H/HeJ mice. Decreased immunoreactive FZD6 was detected in the bronchial epithelial cells and endothelial cells of JF1/MsJ compared to C3H/HeJ mice (Fig. 4).

Discussion
Hindered lung development results in failure to attain peak lung function by early adulthood that, in turn, may increase susceptibility to lung disease. In humans, lung development begins at gestational week 3 and alveolarization continues until at least age 5 years [37]. In mice, lung development begins E9.5 and continues up to P28 and alveolar septal formation continues up to P36 [38]. Recently, Beauchemin et al. [39] described 4-stage postnatal alveolarization with episodic transcriptional activity of pulmonary vascular genes in mice. The processes controlling lung development and growth are believed to be recapitulated following lung injury as genetic subroutines for repair and remodeling processes [40,41]. Therefore, it is plausible that individuals having impaired lung development and growth that do not result in clinically significant symptom may have inefficient repair processes thereby predisposing them to subsequent chronic lung diseases. In this context, elucidating the genomics of lung function development through dissecting the genetics of lung development and growth is a promising approach to identify the predisposing factors. Thus, in this study, we sought to identify candidate genes for TLC by comparing the lung transcript profile of JF1/MsJ and C3H/HeJ strains at 3 critical times involving late embryonic, adolescent, and adult lung development. The lung architectural differences detected in JF1/MsJ and C3H/HeJ mice are suggestive of their difference in development, growth and maturation events. JF1/MsJ mice have decreased mean chord length and mean alveolar volume but increased alveolar number density compared to C3H/HeJ. These findings suggest that JF1/ MsJ mice with smaller alveoli (decreased mean chord length and alveolar volume) compensated by increased alveolar number density. Thus, the lungs of both strains exhibit near equal total alveolar surface area for gas exchange. Lung architectural differences detected between these mouse strains are consistent with their functional characteristics. Increased specific C L-saline/air ratio in JF1/ MsJ mice compared to C3H/HeJ suggests that surface tension contributes more than the elastic recoiling properties of JF1/MsJ lung in determining the lung compliance. This is in agreement with the smaller alveoli in JF1/MsJ mice. Increased elastin to collagen ratio and the decreased C L-Air /TLC and C L-Saline /TLC in JF1/MsJ mice indicates that elastic recoil due to tissue properties is higher in JF1/MsJ mice compared to C3H/HeJ.
In this study, we compared the lung transcript expression contrasting JF1/MsJ and C3H/HeJ mice at 3 stages, (I) embryonic day 18 (completion of embryonic lung development), (II) P28 (adolescent lung, bulk alveolar formation completed); and (III) P70 (mature lung). Noteworthy transcripts with consistent expression across the stages examined included those encoded by Wnt pathway and Values are mean fold difference of (JF1/MsJ)/(C3H/HeJ) for lung transcript levels (cut off for fold change ≥2 fold; false discovery rate 0.10) chitinase-like genes. Gene ontology cellular component "extracellular region" was enriched among decreased JF1/ MsJ transcripts. We employed an e-QTL strategy and consistent expression to identify novel candidate genes for lung function development. e-QTL genes previously not reported with TLC development include Lrrc6, Fpr1, Abcg1, Rsph1, H2-Eb1, H2-D1, H2-Q4, Sfta2, and Gabbr1.
WIF1 inhibits extracellular WNT signaling that plays a crucial role during the lung development. Loss of SMAD1 transcriptional activation of Wif1 has been associated with its decreased expression and increased Wnt/ β-catenin signalling, resulting in abnormal distal lung epithelial cell differentiation and branching morphogenesis [42]. Reduced SMAD1 and WIF1 expression in Table 3 Lung transcripts (34 genes) consistently increased (12)  nitrofen-induced hypoplastic lung resulted in developmental retardation during the saccular stage [43]. Our findings of increased lung WIF1 (>5-fold) transcript and protein levels in alveolar epithelial cells type I and II, as well as in bronchial epithelial cells during alveolar development in JF1/MsJ mice compared to C3H/HeJ strongly supports it to be a candidate for pulmonary function development. Consistent with WIF1 expression, FST is also increased (>4-fold) in JF1/MsJ lungs compared to C3H/HeJ. Canonical Wnt/β-catenin signalling controls FST signalling in satellite cell derived myoblasts [44]. FST-deficient mice exhibit severely retarded overall growth, breathing failure, and perinatal death with fluid filled lungs and poorly expanded alveolar spaces [45]. FST binds and inhibits activins, which are members of the transforming growth factor beta superfamily that plays a significant role in the lung developmental processes [46]. Therapeutic intervention with FST markedly reduced the number of infiltrating cells and ameliorated the destruction of lung architecture in bleomycin-treated rats [47].
FZD6 also functions as a negative regulator of the canonical Wnt/β-catenin signaling cascade. FZD6 controls macroscopic hair patterning in mice and studies show epithelial cells as the source of FZD6-dependent signaling [48]. Wnt signaling pathway was enriched for annotation in the murine Developing Lung Characteristic Subtranscriptome (mDLCS) involving C57BL/6 J, A/J and C3H/HeJ inbred strains. WIF1 and FZD6 transcripts also expressed significant strain dependent modulation during alveolar development [39].
The proteins encoded by Chil1 (aka Chi3i1 or YKL40) and Chil3 (aka Chi3l3 or Ym1) are similar to bacterial chitinases but lacks chitinase activity. JF1/MsJ mice with lower basal lung function have decreased lung CHIL1 (>2 fold) and CHIL3 (>15 fold) transcripts. In humans, the gene homolog for mouse Chil1 is CHIL1. Variations in CHIL1 have been associated with circulating and airway CHIL1 protein levels along with asthma prevalence, severity, hospital admissions, and lung function in children and young adults [49][50][51][52][53]. Variations in CHIL1 are also considered to modulate age-adjusted lung function in cystic fibrosis patients [54]. Elevated serum CHIL1 protein is associated with severe persistent asthma among adults but not in children [55]. However, increased CHIL1 levels are detected among children with severe, steroid resistant asthma [56]. Serum CHIL1 levels are negatively correlated with percent FEV 1 while A B C Fig. 1 Comparative expression of lung Wnt inhibitor factor 1 (Wif1) and frizzled homolog 6 (Fzd6) and chitinase 3 like 3 (CHIL3) in JF1/MsJ and C3H/HeJ mice. a Expression of lung Wif1 mRNA in JF1/Msf was increased compared with C3H/HeJ mice at postnatal (P) days P28 and P70. b Expression of lung Fzd6 mRNA in JF1/Msf was decreased compared to C3H/HeJ mice at P28 and P70. Quantitative real-time polymerase chain reaction was performed, and the comparative cycle number threshold (C T ) method (ΔΔC T ) was used [ΔC T = C T (gene)-C T (Actb)]. Data are presented as expression relative to time-matched C3H/HeJ level (means ± SE; n = 5 mice/strain/stage). *Significantly different from C3H/HeJ (ANOVA followed by all pairwise multiple-comparison procedures with Holm-Sidak method, P < 0.05). c Western blot analysis of CHIL3 from total lung homogenate of P28 C3H/HeJ and JF1/MsJ mice (n = 3 mice/strain; female). Lung transcript expression showed > 15 fold decreased Chil3 expression in Jf1/MsJ compared to C3H/HeJ. βactin (ACTB) was used as the control. Photograph is representative of at least three observations positively correlated with low attenuation area/total lung area percent in COPD patients [57]. Elevated CHIL1 protein has been detected in bronchoalveolar lavage fluid and sputum, with higher numbers of CHIL1expressing cells in bronchial biopsies of smokers with COPD [58,59]. CHIL1 can activate Wnt/β-catenin signaling in alveolar macrophage [60]. Amino acid sequence homology of CHIL3 to proteins associated with tissue remodeling together with their binding capacity with heparin sulfate and GlcN oligomers suggests their role in airway wall remodeling in the allergic lung [61]. CHIL3 protein is abundantly expressed in the allergic mouse lung and enhances Th2 cytokine production by inhibiting 12/15(S)-lipoxygenase which generates lipid metabolites that interferes with T cell proliferation [62]. However, a homologous sequence to Chil3 has not been found in the human genome.
The evaluation of eQTL located in Lfnq3 and Lfnq4 provided 9 candidate genes worthy of further analysis for their possible role in lung development and growth. Located in Lfnq3 region linked to TLC on mCh15, Lrrc6 is a gene expressed in the respiratory epithelium that is essential for proper cilia axonemal assembly of inner and outer dynein arms. Kott and colleagues identified an early frameshift in LRRC6 that is associated with primary ciliary dyskinesia (PCD) [63]. PCD is a group of autosomal-recessive developmental disorders resulting from cilia and sperm-flagella defects, which can lead to respiratory infections, situs inversus, and male infertility [64]. In zebrafish, lrrc6 (aka seahorse) mutants produce ventral body curvature and kidney cysts [65,66] and LRRC6 protein associates with disheveled, which constrains the canonical Wnt/β-catenin pathway and promotes the non-canonical Wnt pathway during gastrulation [67]. In gene-targeted mice lacking functional LRRC6, cilia microtubules remained normal but the outer dynein arms (ODAs), the structures essential for the ciliary beating, are absent. In the absence of LRRC6, ODA proteins that normally are assembled in the cytoplasm and transported to the ciliary axoneme, remain in the cytoplasm and are not transported to the ciliary axoneme [68,69].
Of the eight differentially expressed transcripts located within the Lfnq4 region in mCh17 that has been linked to TLC, Rsph1, Fpr1, Abcg1, and Sfta2 are most noteworthy. Like LRRC6 variants, splice site (275-2A > C) and non-sense (433 C > T) polymorphisms in RSPH1 are associated with PCD [64]. Although persons with these variants can have unexplained neonatal respiratory distress [70], the splice site variation in RSPH1 is associated with a lower prevalence of neonatal respiratory distress, later onset of daily wet cough, and better adult lung function than other forms of PCD [71].
Mitochondrial damage-associated molecular patterns (DAMPs) include formyl peptides that activate human neutrophils through FPR1. On human neutrophils, FPR1 interacts with annexin A1 (aka lipocortin I), which mediated anti-inflammatory activities of glucocorticoids and can regulate epidermal growth factor receptor (EGFR) localization and activity [72,73]. Gene-targeted Fpr1 −/− mice have a slight, but significant decrease in mean chord length compared to strain-matched (C57BL/6 J) control mice [74]. In addition, Fpr1 −/− mice have decreased inflammation and emphysema compared to control mice after cigarette smoke exposure. In humans, cigarette smoking increases the number of FPRs on the surface of neutrophils, which is greater in smokers with COPD than smokers without COPD [75,76]. In addition, a variant rs867228 (minor allele frequency = 0.21) producing a loss-of-function allele in FPR1 is associated with poor metastasis-free and overall survival in breast and colorectal cancer patients receiving adjuvant chemotherapy [77].
Patients with pulmonary alveolar proteinosis (PAP) display impaired surfactant clearance, foamy, lipid-filled alveolar macrophages, and increased cholesterol metabolites within the lung. Also regulated by Wnt/β-catenin signalling [78], ABCG1 lipid transporter is considered to be a key downstream target of colony stimulating factor 2 (aka GM-CSF) and is necessary for proper surfactant catabolism [79,80]. Mice lacking ABCG1 develop alveolar type II cell hypertrophy with lipid deposition, increased levels of surfactant [81], and develop more severe lung fibrosis after bleomycin [82]. Compared with strain-matched control mice, the lungs of ABCG1 null mice are inflamed with macrophage accumulation, lymphocytic infiltration, hemorrhage, eosinophilic crystals, and elevated levels of cytokines and cytokine receptors. BALF samples obtained from Abcg1 −/− mice are marked by increased foamy macrophages and leukocytes and the presence of multiple markers of inflammation including crystals of CHIL3 protein [83]. In mice, ABCG1 is also required for pulmonary B-1 B cell and natural antibody homeostasis [84]. Variations in ABCG1 (rs13050646) and HLA-A (rs3823343, rs2517725) have been associated to IgE dysregulation and atopy [85].
SFTA2 is a secretory protein predominantly expressed by alveolar type II cells and non-ciliated bronchiolar cells [86]. Similar to hydrophobic surfactant proteins B (SFTB) and SFTPC, SFTA2 protein may share particular physicochemical properties [87] and lung SFTA2 mRNA deceases in mice treated with lipopolysaccharide [86]. However, unlike SFTFB and SFTPC, SFTA2 does not co-localize to lamellar bodies, but is found to co-localize to the Golgi and clathrin vesicles like hydrophilic SFTPA1 and SFTPD [86]. Little is known about the disease phenotypes associated with SFTA2 variants, but human SFTA2 is located in a chromosomal locus associated with diffused panbronchiolitis.
Recent studies have implicated Wnt/β-catenin pathway to play an important role in alveolar development [39]. In an attempt to summarize the plausible effect of altered expression of the identified candidate genes based on the above reviewed information, our data suggest that Wnt/β-catenin signaling may be altered during postnatal lung development and growth of JF1/MsJ as compared to C3H/HeJ mice. This is supported by the Fig. 4 Localization of frizzled homolog 6 (FZD6) protein in the lungs of four weeks old female C3H/HeJ and JF1/MsJ mice (n = 5). Decreased FZD6 immunostaining was detected in the b, e airways and h alveolar type I and type II cells of JF1/MsJ lungs compared with those of c, f, i C3H/HeJ mice. a, d, g Sections incubated phosphate-buffered saline (PBS control) without primary antisera increased WIF1 transcript and protein levels in JF1/MsJ lungs. Primary ciliary proteins, like LRRC6 and RSPH, are required not only for regulation of Wnt/β-catenin signaling, but also act as downstream effectors of the Wnt/β-catenin pathway. [88] Alterations in a chitinaselike protein, CHIL1, an activin-binding protein, FST, and a lipid transporter, ABCG1, also implicates this pathway. However, elucidation of the precise mode of action of genetic variants in these genes in determining lung function development warrants functional and mechanistic investigations.

Conclusions
This study demonstrates that the divergent pulmonary function between C3H/HeJ and JF1/MsJ mice is associated with differences in transcription expression profiles in lung. We further dissected the altered genetic signature of divergent lung function development among C3H/HeJ and JF1/MsJ mice and identified several novel genes, especially those with roles in Wnt/β-catenin signalling, not previously associated to lung function. The study also provides a reference of transcripts altered between mouse strains with divergent TLC at the completion of embryonic lung development, bulk alveolar formation and lung growth. The generated data will provide a point base for future association, functional and mechanistic studies for pulmonary function.