Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

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



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.


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.


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.


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.


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 (FEV1) 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.

Evolving from the mouse phenome project, we previously assessed lung function in several inbred mouse strains [19,20,21]. The greatest differences were noted in C3H/HeJ as compared to JF1/MsJ mice. In male C3H/HeJ mice, total lung capacity (TLC) and specific TLC (TLC/body weight) = 1837 μl and 64 μl/g, respectively. In male JF1/MsJ mice, TLC = 904 μl and specific TLC = 40 μl/g. A similar difference was observed in female mice (C3H/HeJ: TLC = 1443 μl and specific TLC 65 μl/g and JF1/MsJ: TLC = 874 μl and specific TLC 53 μl/g). Contrasting C3H/HeJ and JF1/MsJ back- and inter-cross mice, lung function phenotyping was followed by quantitative trait loci (QTL) analysis identifying lung function QTLs (Lfnq1 – Lfnq7) [20]. Consistent with the early life impairment concept, JF1/MsJ mice were more susceptible to lung injury following exposure to carbon nanoparticle (CNP) and acrolein (a toxic component of tobacco smoke) compared to C3H/HeJ [22,23,24].

Examining the mouse lung function QTL regions Lfnq1 and Lfnq2 located on mouse chromosome (mCh) 5 [20], we previously demonstrated that superoxide dismutase 3 extracellular (Sod3), Kit oncogene (Kit) and secreted phosphoprotein 1 (Spp1) as determinants of lung function in mice [21, 25,26,27]. SOD3 variants were also associated with declined lung function in children [25]. When compared to strain-matched C57BL/6 J control mice, gene-targeted Sod3 (−/−) mice exhibited poor ventilation efficiency [21]; Kit mutant mice (Kit W-sh / W-sh) exhibited increased TLC, residual volume, lung compliance (CL), and alveolar size [27]; and gene-targeted Spp1 (−/−) mice exhibited decreased TLC and increased specific CL and alveolar size [26].

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 (40–73 Mbp) and Lfnq4 on mCh17 (10–39 Mbp).



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 H20) 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% OsO4 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 mm3) was performed. Tissue blocks were post-fixed with 1% OsO4 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 H2O), emptied in 0.1 ml steps, and volume and corresponding plateau pressures were determined after each step. CL was determined as the linear slope of the pressure-volume curve.

RNA and Protein extraction

Lung RNA was isolated (RNeasy #74104; Qiagen, Hilden, Germany), quantified (Thermofisher nanodrop 2100), and quality assessed (Agilent Bioanalyzer). For protein analysis, lung was homogenized using 50 mM TrisHCl, 150 nM NaCl, NP-40-1%; sodium desoxycholate – 0.5%W/V; SDS-0.1%, radio immunoprecipitation solution (RIPA buffer #RB4476 Biobasic; Ontario, Canada) in 1:100 ratio supplemented with protease inhibitors (# ab6562 Abcam, Cambridge, UK) and centrifuged (10,000×g, 15 min, 4 °C). The supernatant was aspirated. Total protein content was determined spectrophotometrically (620 nm Protein Assay Dye Reagent # 500–0006; BioRad, USA).

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., 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).


Lungs were fixed (20 cm H2O) in situ with 4% phosphate buffered formaldehyde, embedded into agar-agar blocks, and sectioned (2 μm) for SUR sampling. Deparaffinized sections were subjected to antigen retrieval by microwaving slides (3 X 5 min, 20% power) in 1 mM citrate solution (pH 6.0). Slides were rehydrated in phosphate-buffered saline, followed by immersion into a mixture of methanol and H2O2 (Cat. # 1.08597.1000; Merck, Darmstadt, Germany) (30 min) to inhibit endogenous peroxidases. Tissues were treated with 5% horse serum or 4% bovine serum albumin (BSA) in PBS (45 min, 20 °C) or goat serum in 1:10 ratio (20 min). Primary antibody omission was used as a negative staining control. Mouse WIF1 rabbit polyclonal primary antibody (Cat. #ab186845, Abcam), mouse FZD6 polyclonal goat antibody (Cat. # AF1526; R&D Systems, McKinley, MN), or follistatin (FST) polyclonal goat antibody (Cat. # AF669; R&D Systems) were used. Biotinylated horse anti-goat antibody (Cat. # BA9500; Vector Laboratories, Burlingame, CA) was used as secondary antibody. Vectastain elite ABC kit (Cat. # PK-6100; Vector Laboratories) was used with diaminobenzidine (DAB) reagent (Merck KGaA). Tissues were counterstained with Meyer’s hematoxylin (Cat. # 1.09249.1000; Merck KGaA).


Two-way ANOVA followed by all pairwise multiple-comparison 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 CLair, CLsaline, and specific CLsaline and ratio of CL-saline/CL-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.

Table 1 Morphometric evaluation of C3H/HeJ and JF1/MsJ mouse lung

Transcript analysis in the lung (E18, P28 and P70)

The number of increased and decreased transcripts (fold change ≥ |2|-fold and FDR < 0.10) in JF1/MsJ mice compared to C3H/HeJ was: (i) E18: 45 increased and 102 decreased; (ii) P28: 65 increased and 109 decreased; and (iii) P70: 46 increased and 127 decreased (Additional file 1: Table S1-S6). Using DAVID, decreased transcripts were found to be enriched in the gene ontology category “Extracellular region” (GO: 00005576) (Table 2) at E18 and P70. This was consistent with the observed strain specific differences in elastin to collagen ratio and lung architecture. Collagen type XXVIII, alpha 1 (COL28A1) transcript was decreased by −2.8 fold in JF1/MsJ lungs compared to C3H/HeJ at E18 and P70. At E18, P28, and P70, 12 transcripts increased and 22 decreased in JF1/MsJ compared to C3H/HeJ (Table 3). At stages E18 and P28, 16 transcripts increased and 28 decreased in JF1/MsJ compared to C3H/HeJ (Additional file 1: Table S7). At stages P28 and P70, 32 transcripts increased and 70 decreased lung transcripts in JF1/MsJ compared to C3H/HeJ (Additional file 1: Table S8).

Table 2 Extracellular region is an enriched gene ontology category (GO: 0005576) associated with the significantly decreased transcripts in JF1/MsJ vs C3H/HeJ mouse lung at embryonic day 18 (E18) and postnatal day 70 (P70)
Table 3 Lung transcripts (34 genes) consistently increased (12) or decreased (22) across embryonic stage 18 (E18), postnatal stage (P) 28 and P70 in JF1/MsJ (JF1) mice compared to C3H/HeJ (C3H)

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).

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 (CT) method (ΔΔCT) was used [ΔCT = CT (gene)-CT(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

Fig. 2

Localization of Wnt inhibitor factor 1 (WIF1) protein in the lungs of four weeks old female C3H/HeJ and JF1/MsJ mice (n = 5). Increased WIF1 immunostaining was detected in the b, e airway epithelium 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

Fig. 3

Localization of follistatin (FST) protein in the lungs of four weeks old female C3H/HeJ and JF1/MsJ mice (n = 5). Increased of FST 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

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

Expression quantitative trait loci (eQTL) analysis

To further assess candidate genes for Lfnq3 and Lfnq4 associated with TLC, we compared JF1/MsJ to C3H/HeJ lung transcript levels at any of the three investigated stages. Employing this strategy, we identified 9 eQTLs located on mCh15 (Lfnq3) and mCh17 (Lfnq4) (Table 4) not previously reported for lung function development in mice. The candidate for Lfnq3 was leucine rich repeat containing 6 (testis) (Lrrc6). Several candidates for Lfnq4 were identified that include radial spoke head 1 homolog (Chlamydomonas) (Rsph1), formyl peptide receptor 1 (Fpr1), ATP-binding cassette, sub-family G (WHITE), member 1 (Abcg1), 3 histocompatibility 2 genes including class II antigen E beta (H2-Eb1), D region locus 1 (H2-D1), Q region locus 4 (H2-Q4), surfactant associated 2 (Sfta2), and gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1).

Table 4 Candidate genes for total lung capacity (TLC) in mice identified through the expression quantitative trait loci (e-QTL) strategy comparing C3H/HeJ (C3H) and JF1/MsJ (JF1). [cut off for fold change ≥2.0 fold at least in one stage among embryonic stage 18, postnatal stage (P) 28 and P70; false discovery rate = 0.10]


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 CL-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 CL-Air/TLC and CL-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 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 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 FEV1 while 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 CHIL1-expressing 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 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 chitinase-like 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.


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.


Abcg1 :

ATP-binding cassette, sub-family G (WHITE), member 1


Beta actin


Bronchoalveolar lavage fluid

Chil :

Chitinase-like 1

Chil3 :

Chitinase-like 3

CL :



Chronic obstructive pulmonary disease


Embryonic day


Expression-quantitative trait loci

FEV1 :

Forced expiratory volume in 1 s

Fpr1 :

Formyl peptide receptor 1

Fst :


Fzd6 :

Frizzled homolog 6

Gabbr1 :

Gamma-aminobutyric acid (GABA) B receptor, 1


Gene ontology

H2-D1 :

Histocompatibility 2, D region locus 1

H2-Eb1 :

Histocompatibility 2, class II antigen E beta

H2-Q4 :

Histocompatibility 2, Q region locus 4


Institutional Animal Ethics Committee



Kit :

Kit oncogene


Lung function quantitative trait loci

Lrrc6 :

Leucine rich repeat containing 6 (testis)


Mouse chromosome


Developing lung characteristic subtranscriptome


Postnatal day


Quantitative real time polymerase chain reaction


Quantitative trait loci

Rsph1 :

Radial spoke head 1 homolog (Chlamydomonas)

Sfta2 :

Surfactant associated 2

Sod3 :

Superoxide dismutase 3, extracellular

Spp1 :

Secreted phosphoprotein 1


Systematic uniform random


Transmission electron microscope


Total lung capacity

Wif1 :

Wnt inhibitor factor 1


  1. 1.

    Krauss-Etschmann S, Bush A, Bellusci S, Brusselle GG, Dahlen SE, Dehmel S, Eickelberg O, Gibson G, Hylkema MN, Knaus P, Konigshoff M, Lloyd CM, Macciarini P, Mailleux A, Marsland BJ, Postma DS, Roberts G, Samakovlis C, Stocks J, Vandesompele J, Wjst M, Holloway J. Of flies, mice and men: a systematic approach to understanding the early life origins of chronic lung disease. Thorax. 2013;68:380–4.

  2. 2.

    Stocks J, Hislop A, Sonnappa S. Early lung development: lifelong effect on respiratory health and disease. Lancet Respir Med. 2013;1:728–42.

  3. 3.

    Stocks J, Sonnappa S. Early life influences on the development of chronic obstructive pulmonary disease. Ther Adv Respir Dis. 2013;7:161–73.

  4. 4.

    Martinez FD. The origins of asthma and chronic obstructive pulmonary disease in early life. Proc Am Thorac Soc. 2009;6:272–7.

  5. 5.

    Lange P, Celli B, Agusti A, Boje Jensen G, Divo M, Faner R, Guerra S, Marott JL, Martinez FD, Martinez-Camblor P, Meek P, Owen CA, Petersen H, Pinto-Plata V, Schnohr P, Sood A, Soriano JB, Tesfaigzi Y, Vestbo J. Lung-Function Trajectories Leading to Chronic Obstructive Pulmonary Disease. N Engl J Med. 2015;373:111–22.

  6. 6.

    Repapi E, Sayers I, Wain LV, Burton PR, Johnson T, Obeidat M, Zhao JH, Ramasamy A, Zhai G, Vitart V, Huffman JE, Igl W, Albrecht E, Deloukas P, Henderson J, Granell R, WL MA, Rudnicka AR, Wellcome Trust Case Control C, Barroso I, Loos RJ, Wareham NJ, Mustelin L, Rantanen T, Surakka I, Imboden M, Wichmann HE, Grkovic I, Jankovic S, Zgaga L, Hartikainen AL, Peltonen L, Gyllensten U, Johansson A, Zaboli G, Campbell H, Wild SH, Wilson JF, Glaser S, Homuth G, Volzke H, Mangino M, Soranzo N, Spector TD, Polasek O, Rudan I, Wright AF, Heliovaara M, Ripatti S, Pouta A, Naluai AT, Olin AC, Toren K, Cooper MN, James AL, Palmer LJ, Hingorani AD, Wannamethee SG, Whincup PH, Smith GD, Ebrahim S, TM MK, Pavord ID, AK ML, Morris AD, Porteous DJ, Cooper C, Dennison E, Shaheen S, Karrasch S, Schnabel E, Schulz H, Grallert H, Bouatia-Naji N, Delplanque J, Froguel P, Blakey JD, NRS T, Britton JR, Morris RW, Holloway JW, Lawlor DA, Hui J, Nyberg F, Jarvelin MR, Jackson C, Kahonen M, Kaprio J, Probst-Hensch NM, Koch B, Hayward C, Evans DM, Elliott P, Strachan DP, Hall IP, Tobin MD. Genome-wide association study identifies five loci associated with lung function. Nat Genet. 2010;42:36–44.

  7. 7.

    Yao TC, Du G, Han L, Sun Y, Hu D, Yang JJ, Mathias R, Roth LA, Rafaels N, Thompson EE, Loisel DA, Anderson R, Eng C, Arruabarrena Orbegozo M, Young M, Klocksieben JM, Anderson E, Shanovich K, Lester LA, Williams LK, Barnes KC, Burchard EG, Nicolae DL, Abney M, Ober C. Genome-wide association study of lung function phenotypes in a founder population. J Allergy Clin Immunol. 2014;133:248–55.e241-210.

  8. 8.

    Hancock DB, Eijgelsheim M, Wilk JB, Gharib SA, Loehr LR, Marciante KD, Franceschini N, van Durme YM, Chen TH, Barr RG, Schabath MB, Couper DJ, Brusselle GG, Psaty BM, van Duijn CM, Rotter JI, Uitterlinden AG, Hofman A, Punjabi NM, Rivadeneira F, Morrison AC, Enright PL, North KE, Heckbert SR, Lumley T, Stricker BH, O'Connor GT, London SJ. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet. 2010;42:45–52.

  9. 9.

    Soler Artigas M, Loth DW, Wain LV, Gharib SA, Obeidat M, Tang W, Zhai G, Zhao JH, Smith AV, Huffman JE, Albrecht E, Jackson CM, Evans DM, Cadby G, Fornage M, Manichaikul A, Lopez LM, Johnson T, Aldrich MC, Aspelund T, Barroso I, Campbell H, Cassano PA, Couper DJ, Eiriksdottir G, Franceschini N, Garcia M, Gieger C, Gislason GK, Grkovic I, Hammond CJ, Hancock DB, Harris TB, Ramasamy A, Heckbert SR, Heliovaara M, Homuth G, Hysi PG, James AL, Jankovic S, Joubert BR, Karrasch S, Klopp N, Koch B, Kritchevsky SB, Launer LJ, Liu Y, Loehr LR, Lohman K, Loos RJ, Lumley T, Al Balushi KA, Ang WQ, Barr RG, Beilby J, Blakey JD, Boban M, Boraska V, Brisman J, Britton JR, Brusselle GG, Cooper C, Curjuric I, Dahgam S, Deary IJ, Ebrahim S, Eijgelsheim M, Francks C, Gaysina D, Granell R, Gu X, Hankinson JL, Hardy R, Harris SE, Henderson J, Henry A, Hingorani AD, Hofman A, Holt PG, Hui J, Hunter ML, Imboden M, Jameson KA, Kerr SM, Kolcic I, Kronenberg F, Liu JZ, Marchini J, McKeever T, Morris AD, Olin AC, Porteous DJ, Postma DS, Rich SS, Ring SM, Rivadeneira F, Rochat T, Sayer AA, Sayers I, Sly PD, Smith GD, Sood A, Starr JM, Uitterlinden AG, Vonk JM, Wannamethee SG, Whincup PH, Wijmenga C, Williams OD, Wong A, Mangino M, Marciante KD, WL MA, Meibohm B, Morrison AC, North KE, Omenaas E, Palmer LJ, Pietilainen KH, Pin I, Pola Sbreve Ek O, Pouta A, Psaty BM, Hartikainen AL, Rantanen T, Ripatti S, Rotter JI, Rudan I, Rudnicka AR, Schulz H, Shin SY, Spector TD, Surakka I, Vitart V, Volzke H, Wareham NJ, Warrington NM, Wichmann HE, Wild SH, Wilk JB, Wjst M, Wright AF, Zgaga L, Zemunik T, Pennell CE, Nyberg F, Kuh D, Holloway JW, Boezen HM, Lawlor DA, Morris RW, Probst-Hensch N, International Lung Cancer C, Consortium G, Kaprio J, Wilson JF, Hayward C, Kahonen M, Heinrich J, Musk AW, Jarvis DL, Glaser S, Jarvelin MR, Ch Stricker BH, Elliott P, O'Connor GT, Strachan DP, London SJ, Hall IP, Gudnason V, Tobin MD. Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat Genet. 2011;43:1082–90.

  10. 10.

    Soler Artigas M, Wain LV, Repapi E, Obeidat M, Sayers I, Burton PR, Johnson T, Zhao JH, Albrecht E, Dominiczak AF, Kerr SM, Smith BH, Cadby G, Hui J, Palmer LJ, Hingorani AD, Wannamethee SG, Whincup PH, Ebrahim S, Smith GD, Barroso I, Loos RJ, Wareham NJ, Cooper C, Dennison E, Shaheen SO, Liu JZ, Marchini J, Medical Research Council National Survey of H, Development Respiratory Study T, Dahgam S, Naluai AT, Olin AC, Karrasch S, Heinrich J, Schulz H, TM MK, Pavord ID, Heliovaara M, Ripatti S, Surakka I, Blakey JD, Kahonen M, Britton JR, Nyberg F, Holloway JW, Lawlor DA, Morris RW, James AL, Jackson CM, Hall IP, Tobin MD, SpiroMeta C. Effect of five genetic variants associated with lung function on the risk of chronic obstructive lung disease, and their joint effects on lung function. Am J Respir Crit Care Med. 2011;184:786–95.

  11. 11.

    Obeidat M, Hao K, Bosse Y, Nickle DC, Nie Y, Postma DS, Laviolette M, Sandford AJ, Daley DD, Hogg JC, Elliott WM, Fishbane N, Timens W, Hysi PG, Kaprio J, Wilson JF, Hui J, Rawal R, Schulz H, Stubbe B, Hayward C, Polasek O, Jarvelin MR, Zhao JH, Jarvis D, Kahonen M, Franceschini N, North KE, Loth DW, Brusselle GG, Smith AV, Gudnason V, Bartz TM, Wilk JB, O'Connor GT, Cassano PA, Tang W, Wain LV, Soler Artigas M, Gharib SA, Strachan DP, Sin DD, Tobin MD, London SJ, Hall IP, Pare PD. Molecular mechanisms underlying variations in lung function: a systems genetics analysis. Lancet Respir Med. 2015;3:782–95.

  12. 12.

    Tang W, Kowgier M, Loth DW, Soler Artigas M, Joubert BR, Hodge E, Gharib SA, Smith AV, Ruczinski I, Gudnason V, Mathias RA, Harris TB, Hansel NN, Launer LJ, Barnes KC, Hansen JG, Albrecht E, Aldrich MC, Allerhand M, Barr RG, Brusselle GG, Couper DJ, Curjuric I, Davies G, Deary IJ, Dupuis J, Fall T, Foy M, Franceschini N, Gao W, Glaser S, Gu X, Hancock DB, Heinrich J, Hofman A, Imboden M, Ingelsson E, James A, Karrasch S, Koch B, Kritchevsky SB, Kumar A, Lahousse L, Li G, Lind L, Lindgren C, Liu Y, Lohman K, Lumley T, WL MA, Meibohm B, Morris AP, Morrison AC, Musk B, North KE, Palmer LJ, Probst-Hensch NM, Psaty BM, Rivadeneira F, Rotter JI, Schulz H, Smith LJ, Sood A, Starr JM, Strachan DP, Teumer A, Uitterlinden AG, Volzke H, Voorman A, Wain LV, Wells MT, Wilk JB, Williams OD, Heckbert SR, Stricker BH, London SJ, Fornage M, Tobin MD, O'Connor GT, Hall IP, Cassano PA. Large-scale genome-wide association studies and meta-analyses of longitudinal change in adult lung function. PLoS One. 2014;9:e100776.

  13. 13.

    Loth DW, Soler Artigas M, Gharib SA, Wain LV, Franceschini N, Koch B, Pottinger TD, Smith AV, Duan Q, Oldmeadow C, Lee MK, Strachan DP, James AL, Huffman JE, Vitart V, Ramasamy A, Wareham NJ, Kaprio J, Wang XQ, Trochet H, Kahonen M, Flexeder C, Albrecht E, Lopez LM, de Jong K, Thyagarajan B, Alves AC, Enroth S, Omenaas E, Joshi PK, Fall T, Vinuela A, Launer LJ, Loehr LR, Fornage M, Li G, Wilk JB, Tang W, Manichaikul A, Lahousse L, Harris TB, North KE, Rudnicka AR, Hui J, Gu X, Lumley T, Wright AF, Hastie ND, Campbell S, Kumar R, Pin I, Scott RA, Pietilainen KH, Surakka I, Liu Y, Holliday EG, Schulz H, Heinrich J, Davies G, Vonk JM, Wojczynski M, Pouta A, Johansson A, Wild SH, Ingelsson E, Rivadeneira F, Volzke H, Hysi PG, Eiriksdottir G, Morrison AC, Rotter JI, Gao W, Postma DS, White WB, Rich SS, Hofman A, Aspelund T, Couper D, Smith LJ, Psaty BM, Lohman K, Burchard EG, Uitterlinden AG, Garcia M, Joubert BR, WL MA, Musk AB, Hansel N, Heckbert SR, Zgaga L, van Meurs JB, Navarro P, Rudan I, Oh YM, Redline S, Jarvis DL, Zhao JH, Rantanen T, O'Connor GT, Ripatti S, Scott RJ, Karrasch S, Grallert H, Gaddis NC, Starr JM, Wijmenga C, Minster RL, Lederer DJ, Pekkanen J, Gyllensten U, Campbell H, Morris AP, Glaser S, Hammond CJ, Burkart KM, Beilby J, Kritchevsky SB, Gudnason V, Hancock DB, Williams OD, Polasek O, Zemunik T, Kolcic I, Petrini MF, Wjst M, Kim WJ, Porteous DJ, Scotland G, Smith BH, Viljanen A, Heliovaara M, Attia JR, Sayers I, Hampel R, Gieger C, Deary IJ, Boezen HM, Newman A, Jarvelin MR, Wilson JF, Lind L, Stricker BH, Teumer A, Spector TD, Melen E, Peters MJ, Lange LA, Barr RG, Bracke KR, Verhamme FM, Sung J, Hiemstra PS, Cassano PA, Sood A, Hayward C, Dupuis J, Hall IP, Brusselle GG, Tobin MD, London SJ. Genome-wide association analysis identifies six new loci associated with forced vital capacity. Nat Genet. 2014;46:669–77.

  14. 14.

    Young S, Sherrill DL, Arnott J, Diepeveen D, LeSouef PN, Landau LI. Parental factors affecting respiratory function during the first year of life. Pediatr Pulmonol. 2000;29:331–40.

  15. 15.

    Murray CS, Pipis SD, McArdle EC, Lowe LA, Custovic A, Woodcock A, National Asthma Campaign-Manchester A, Allergy Study G. Lung function at one month of age as a risk factor for infant respiratory symptoms in a high risk population. Thorax 2002;57: 388–392.

  16. 16.

    Mullane D, Turner SW, Cox DW, Goldblatt J, Landau LI, le Souef PN. Reduced infant lung function, active smoking, and wheeze in 18-year-old individuals. JAMA Pediatr. 2013;167:368–73.

  17. 17.

    Turner S, Fielding S, Mullane D, Cox DW, Goldblatt J, Landau L, le Souef P. A longitudinal study of lung function from 1 month to 18 years of age. Thorax. 2014;69:1015–20.

  18. 18.

    Skripak JM. Persistent effects of maternal smoking during pregnancy on lung function and asthma in adolescents. Pediatrics. 2014;134(Suppl 3):S146.

  19. 19.

    Reinhard C, Eder G, Fuchs H, Ziesenis A, Heyder J, Schulz H. Inbred strain variation in lung function. Mamm Genome. 2002;13:429–37.

  20. 20.

    Reinhard C, Meyer B, Fuchs H, Stoeger T, Eder G, Ruschendorf F, Heyder J, Nurnberg P, de Angelis MH, Schulz H. Genomewide linkage analysis identifies novel genetic Loci for lung function in mice. Am J Respir Crit Care Med. 2005;171:880–8.

  21. 21.

    Ganguly K, Stoeger T, Wesselkamper SC, Reinhard C, Sartor MA, Medvedovic M, Tomlinson CR, Bolle I, Mason JM, Leikauf GD, Schulz H. Candidate genes controlling pulmonary function in mice: transcript profiling and predicted protein structure. Physiol Genomics. 2007;31:410–21.

  22. 22.

    Ganguly K, Upadhyay S, Irmler M, Takenaka S, Pukelsheim K, Beckers J, De Angelis MH, Hamelmann E, Stoeger T, Schulz H. Impaired resolution of inflammatory response in the lungs of JF1/Msf mice following carbon nanoparticle instillation. Respir Res. 2011;12:94.

  23. 23.

    Leikauf GD, Concel VJ, Liu P, Bein K, Berndt A, Ganguly K, Jang AS, Brant KA, Dietsch M, Pope-Varsalona H, Dopico RA Jr, Di YP, Li Q, Vuga LJ, Medvedovic M, Kaminski N, You M, Prows DR. Haplotype association mapping of acute lung injury in mice implicates activin a receptor, type 1. Am J Respir Crit Care Med. 2011;183:1499–509.

  24. 24.

    Ganguly K, Upadhyay S, Irmler M, Takenaka S, Pukelsheim K, Beckers J, Hamelmann E, Schulz H, Stoeger T. Pathway focused protein profiling indicates differential function for IL-1B, −18 and VEGF during initiation and resolution of lung inflammation evoked by carbon nanoparticle exposure in mice. Part Fibre Toxicol. 2009;6:31.

  25. 25.

    Ganguly K, Depner M, Fattman C, Bein K, Oury TD, Wesselkamper SC, Borchers MT, Schreiber M, Gao F, von Mutius E, Kabesch M, Leikauf GD, Schulz H. Superoxide dismutase 3, extracellular (SOD3) variants and lung function. Physiol Genomics. 2009;37:260–7.

  26. 26.

    Ganguly K, Martin TM, Concel VJ, Upadhyay S, Bein K, Brant KA, George L, Mitra A, Thimraj TA, Fabisiak JP, Vuga LJ, Fattman C, Kaminski N, Schulz H, Leikauf GD. Secreted phosphoprotein 1 is a determinant of lung function development in mice. Am J Respir Cell Mol Biol. 2014;51:637–51.

  27. 27.

    Lindsey JY, Ganguly K, Brass DM, Li Z, Potts EN, Degan S, Chen H, Brockway B, Abraham SN, Berndt A, Stripp BR, Foster WM, Leikauf GD, Schulz H, Hollingsworth JW. c-Kit is essential for alveolar maintenance and protection from emphysema-like disease in mice. Am J Respir Crit Care Med. 2011;183:1644–52.

  28. 28.

    Hsia CC, Hyde DM, Ochs M, Weibel ER, Structure AEJTFoQAoL. An official research policy statement of the American Thoracic Society/European Respiratory Society: standards for quantitative assessment of lung structure. Am J Respir Crit Care Med. 2010;181:394–418.

  29. 29.

    Fehrenbach H, Zimmermann G, Starke E, Bratu VA, Conrad D, Yildirim AO, Fehrenbach A. Nitrogen dioxide induces apoptosis and proliferation but not emphysema in rat lungs. Thorax. 2007;62:438–46.

  30. 30.

    Madurga A, Mizikova I, Ruiz-Camp J, Vadasz I, Herold S, Mayer K, Fehrenbach H, Seeger W, Morty RE. Systemic hydrogen sulfide administration partially restores normal alveolarization in an experimental animal model of bronchopulmonary dysplasia. Am J Physiol Lung Cell Mol Physiol. 2014;306:L684–97.

  31. 31.

    Fehrenbach H, Ochs M. Studying lung ultrastructure. In: Uhlig S, Taylor A, editors. Methods in pulmonary research. Basel: Birkhäuser; 1998. p. 429–54.

  32. 32.

    Haidar A, Ryder TA, Mobberley MA, Wigglesworth JS. Two techniques for electron opaque staining of elastic fibres using tannic acid in fresh and formalin fixed tissue. J Clin Pathol. 1992;45:633–5.

  33. 33.

    Schulz H, Johner C, Eder G, Ziesenis A, Reitmeier P, Heyder J, Balling R. Respiratory mechanics in mice: strain and sex specific differences. Acta Physiol Scand. 2002;174:367–5.

  34. 34.

    R Development Core team. R: A language and environment for statistical computing. Vienna: the R Foundation for Statistical Computing; 2011.

  35. 35.

    Rainer J, Sanchez-Cabo F, Stocker G, Sturn A, Trajanoski Z. CARMA web: comprehensive R- and Bioconductor-based web service for microarray data analysis. Nucleic Acids Res. 2006;34:W498–503.

  36. 36.

    Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009;4(1):44–57.

  37. 37.

    Warburton D, El-Hashash A, Carraro G, Tiozzo C, Sala F, Rogers O, De Langhe S, Kemp PJ, Riccardi D, Torday J, Bellusci S, Shi W, Lubkin SR, Jesudason E. Lung organogenesis. In: Koopman P, editor. Organogenesis in development. 1st ed. Oxford: Academic; 2010. p. 73–158.

  38. 38.

    Schittny JC, Mund SI, Stampanoni M. Evidence and structural mechanism for late lung alveolarization. Am J Physiol Lung Cell Mol Physiol. 2008;294:L246–54.

  39. 39.

    Beauchemin KJ, Wells JM, Kho AT, Philip VM, Kamir D, Kohane IS, Graber JH, Bult CJ. Temporal dynamics of the developing lung transcriptome in three common inbred strains of laboratory mice reveals multiple stages of postnatal alveolar development. PeerJ. 2016;4:e2318.

  40. 40.

    Hagood JS, Ambalavanan N. Systems biology of lung development and regeneration: current knowledge and recommendations for future research. Wiley Interdiscip Rev Syst Biol Med. 2013;5:125–33.

  41. 41.

    Stabler CT, Morrisey EE. Developmental pathways in lung regeneration. Cell Tissue Res. 2017;367:677–85.

  42. 42.

    Xu B, Chen C, Chen H, Zheng SG, Bringas P Jr, Xu M, Zhou X, Chen D, Umans L, Zwijsen A, Shi W. Smad1 and its target gene Wif1 coordinate BMP and Wnt signaling activities to regulate fetal lung development. Development. 2011;138:925–35.

  43. 43.

    Fujiwara N, Doi T, Gosemann JH, Kutasy B, Friedmacher F, Puri P. Smad1 and WIF1 genes are downregulated during saccular stage of lung development in the nitrofen rat model. Pediatr Surg Int. 2012;28:189–93.

  44. 44.

    Jones AE, Price FD, Le Grand F, Soleimani VD, Dick SA, Megeney LA, Rudnicki MA. Wnt/beta-catenin controls follistatin signalling to regulate satellite cell myogenic potential. Skelet Muscle. 2015;5:14.

  45. 45.

    Matzuk MM, Lu N, Vogel H, Sellheyer K, Roop DR, Bradley A. Multiple defects and perinatal death in mice deficient in follistatin. Nature. 1995;374:360–3.

  46. 46.

    Zhao Y, Silbajoris R, Young SL. Identification and developmental expression of two activin receptors in baboon lung. Biochem Biophys Res Commun. 1996;229:50–7.

  47. 47.

    Aoki F, Kurabayashi M, Hasegawa Y, Kojima I. Attenuation of bleomycin-induced pulmonary fibrosis by follistatin. Am J Respir Crit Care Med. 2005;172:713–20.

  48. 48.

    Guo N, Hawkins C, Nathans J. Frizzled6 controls hair patterning in mice. Proc Natl Acad Sci U S A. 2004;101:9277–81.

  49. 49.

    Tang H, Fang Z, Sun Y, Li B, Shi Z, Chen J, Zhang T, Xiu Q. YKL-40 in asthmatic patients, and its correlations with exacerbation, eosinophils and immunoglobulin E. Eur Respir J. 2010;35:757–60.

  50. 50.

    Cunningham J, Basu K, Tavendale R, Palmer CN, Smith H, Mukhopadhyay S. The CHI3L1 rs4950928 polymorphism is associated with asthma-related hospital admissions in children and young adults. Ann Allergy Asthma Immunol. 2011;106:381–6.

  51. 51.

    Verlaan DJ, Ouimet M, Adoue V, Sirois-Gagnon D, Lariviere M, Ge B, Beaulieu P, Dias J, Lam KC, Koka V, Laprise C, Pastinen T, Sinnett D. Promoter polymorphisms in CHI3L1 are associated with asthma. J Allergy Clin Immunol. 2012;130:533–5.

  52. 52.

    Gomez JL, Crisafi GM, Holm CT, Meyers DA, Hawkins GA, Bleecker ER, Jarjour N, Severe Asthma Research Program I, Cohn L, Chupp GL. Genetic variation in chitinase 3-like 1 (CHI3L1) contributes to asthma severity and airway expression of YKL-40. J Allergy Clin Immunol. 2015;e10(136):51–8.

  53. 53.

    Li JM, Zhang HF, Shen XL, Xie H, Wu XD, Shen T, Wang Y. Association between CHI3L1 SNPs and susceptibility to childhood asthma. Zhongguo Dang Dai Er Ke Za Zhi. 2015;17:144–8.

  54. 54.

    Hector A, Kormann MS, Mack I, Latzin P, Casaulta C, Kieninger E, Zhou Z, Yildirim AO, Bohla A, Rieber N, Kappler M, Koller B, Eber E, Eickmeier O, Zielen S, Eickelberg O, Griese M, Mall MA, Hartl D. The chitinase-like protein YKL-40 modulates cystic fibrosis lung disease. PLoS One. 2011;6:e24399.

  55. 55.

    Santos CB, Davidson J, Covar RA, Spahn JD. The chitinase-like protein YKL-40 is not a useful biomarker for severe persistent asthma in children. Ann Allergy Asthma Immunol. 2014;113:263–6.

  56. 56.

    Konradsen JR, James A, Nordlund B, Reinius LE, Soderhall C, Melen E, Wheelock AM, Lodrup Carlsen KC, Lidegran M, Verhoek M, Boot RG, Dahlen B, Dahlen SE, Hedlin G. The chitinase-like protein YKL-40: a possible biomarker of inflammation and airway remodeling in severe pediatric asthma. J Allergy Clin Immunol. 2013;132:328–35.

  57. 57.

    Sakazaki Y, Hoshino T, Takei S, Sawada M, Oda H, Takenaka S, Imaoka H, Matsunaga K, Ota T, Abe Y, Miki I, Fujimoto K, Kawayama T, Kato S, Aizawa H. Overexpression of chitinase 3-like 1/YKL-40 in lung-specific IL-18-transgenic mice, smokers and COPD. PLoS One. 2011;6:e24177.

  58. 58.

    Letuve S, Kozhich A, Arouche N, Grandsaigne M, Reed J, Dombret MC, Kiener PA, Aubier M, Coyle AJ, Pretolani M. YKL-40 is elevated in patients with chronic obstructive pulmonary disease and activates alveolar macrophages. J Immunol. 2008;181:5167–73.

  59. 59.

    Otsuka K, Matsumoto H, Niimi A, Muro S, Ito I, Takeda T, Terada K, Yamaguchi M, Matsuoka H, Jinnai M, Oguma T, Nakaji H, Inoue H, Tajiri T, Iwata T, Chin K, Mishima M. Sputum YKL-40 levels and pathophysiology of asthma and chronic obstructive pulmonary disease. Respiration. 2012;83:507–19.

  60. 60.

    He CH, Lee CG, Dela Cruz CS, Lee CM, Zhou Y, Ahangari F, Ma B, Herzog EL, Rosenberg SA, Li Y, Nour AM, Parikh CR, Schmidt I, Modis Y, Cantley L, Elias JA. Chitinase 3-like 1 regulates cellular and tissue responses via IL-13 receptor α2. Cell Rep. 2013;4(4):830–41. Cell Rep. 2013 Nov 27;5(4):1156

  61. 61.

    Webb DC, McKenzie AN, Foster PS. Expression of the Ym2 lectin-binding protein is dependent on interleukin (IL)-4 and IL-13 signal transduction: identification of a novel allergy-associated protein. J Biol Chem. 2001;276:41969–76.

  62. 62.

    Cai Y, Kumar RK, Zhou J, Foster PS, Webb DC. Ym1/2 promotes Th2 cytokine expression by inhibiting 12/15(S)-lipoxygenase: identification of a novel pathway for regulating allergic inflammation. J Immunol. 2009;182:5393–9.

  63. 63.

    Kott E, Duquesnoy P, Copin B, Legendre M, Dastot-Le Moal F, Montantin G, Jeanson L, Tamalet A, Papon JF, Siffroi JP, Rives N, Mitchell V, de Blic J, Coste A, Clement A, Escalier D, Toure A, Escudier E, Amselem S. Loss-of-function mutations in LRRC6, a gene essential for proper axonemal assembly of inner and outer dynein arms, cause primary ciliary dyskinesia. Am J Hum Genet. 2012;91:958–64.

  64. 64.

    Raidt J, Wallmeier J, Hjeij R, Onnebrink JG, Pennekamp P, Loges NT, Olbrich H, Haffner K, Dougherty GW, Omran H, Werner C. Ciliary beat pattern and frequency in genetic variants of primary ciliary dyskinesia. Eur Respir J. 2014;44:1579–88.

  65. 65.

    Amsterdam A, Nissen RM, Sun Z, Swindell EC, Farrington S, Hopkins N. Identification of 315 genes essential for early zebrafish development. Proc Natl Acad Sci U S A. 2004;101:12792–7.

  66. 66.

    Sun Z, Amsterdam A, Pazour GJ, Cole DG, Miller MS, Hopkins N. A genetic screen in zebrafish identifies cilia genes as a principal cause of cystic kidney. Development. 2004;131:4085–93.

  67. 67.

    Kishimoto N, Cao Y, Park A, Sun Z. Cystic kidney gene seahorse regulates cilia-mediated processes and Wnt pathways. Dev Cell. 2008;14:954–61.

  68. 68.

    Zhao L, Yuan S, Cao Y, Kallakuri S, Li Y, Kishimoto N, DiBella L, Sun Z. Reptin/Ruvbl2 is a Lrrc6/Seahorse interactor essential for cilia motility. Proc Natl Acad Sci U S A. 2013;110:12697–702.

  69. 69.

    Inaba Y, Shinohara K, Botilde Y, Nabeshima R, Takaoka K, Ajima R, Lamri L, Takeda H, Saga Y, Nakamura T, Hamada H. Transport of the outer dynein arm complex to cilia requires a cytoplasmic protein Lrrc6. Genes Cells. 2016;21:728–39.

  70. 70.

    Kott E, Legendre M, Copin B, Papon JF, Dastot-Le Moal F, Montantin G, Duquesnoy P, Piterboth W, Amram D, Bassinet L, Beucher J, Beydon N, Deneuville E, Houdouin V, Journel H, Just J, Nathan N, Tamalet A, Collot N, Jeanson L, Le Gouez M, Vallette B, Vojtek AM, Epaud R, Coste A, Clement A, Housset B, Louis B, Escudier E, Amselem S. Loss-of-function mutations in RSPH1 cause primary ciliary dyskinesia with central-complex and radial-spoke defects. Am J Hum Genet. 2013;93:561–70.

  71. 71.

    Knowles MR, Ostrowski LE, Leigh MW, Sears PR, Davis SD, Wolf WE, Hazucha MJ, Carson JL, Olivier KN, Sagel SD, Rosenfeld M, Ferkol TW, Dell SD, Milla CE, Randell SH, Yin W, Sannuti A, Metjian HM, Noone PG, Noone PJ, Olson CA, Patrone MV, Dang H, Lee HS, Hurd TW, Gee HY, Otto EA, Halbritter J, Kohl S, Kircher M, Krischer J, Bamshad MJ, Nickerson DA, Hildebrandt F, Shendure J, Zariwala MA. Mutations in RSPH1 cause primary ciliary dyskinesia with a unique clinical and ciliary phenotype. Am J Respir Crit Care Med. 2014;189:707–17.

  72. 72.

    Walther A, Riehemann K, Gerke V. A novel ligand of the formyl peptide receptor: annexin I regulates neutrophil extravasation by interacting with the FPR. Mol Cell. 2000;5:831–40.

  73. 73.

    Grewal T, Enrich C. Annexins--modulators of EGF receptor signalling and trafficking. Cell Signal. 2009;21:847–58.

  74. 74.

    Cardini S, Dalli J, Fineschi S, Perretti M, Lungarella G, Lucattelli M. Genetic ablation of the fpr1 gene confers protection from smoking-induced lung emphysema in mice. Am J Respir Cell Mol Biol. 2012;47:332–9.

  75. 75.

    Stockley RA, Grant RA, Llewellyn-Jones CG, Hill SL, Burnett D. Neutrophil formyl-peptide receptors. Relationship to peptide-induced responses and emphysema. Am J Respir Crit Care Med. 1994;149:464–8.

  76. 76.

    Matheson M, Rynell AC, McClean M, Berend N. Cigarette smoking increases neutrophil formyl methionyl leucyl phenylalanine receptor numbers. Chest. 2003;123:1642–6.

  77. 77.

    Vacchelli E, Ma Y, Baracco EE, Sistigu A, Enot DP, Pietrocola F, Yang H, Adjemian S, Chaba K, Semeraro M, Signore M, De Ninno A, Lucarini V, Peschiaroli F, Businaro L, Gerardino A, Manic G, Ulas T, Gunther P, Schultze JL, Kepp O, Stoll G, Lefebvre C, Mulot C, Castoldi F, Rusakiewicz S, Ladoire S, Apetoh L, Bravo-San Pedro JM, Lucattelli M, Delarasse C, Boige V, Ducreux M, Delaloge S, Borg C, Andre F, Schiavoni G, Vitale I, Laurent-Puig P, Mattei F, Zitvogel L, Kroemer G. Chemotherapy-induced antitumor immunity requires formyl peptide receptor 1. Science. 2015;350:972–8.

  78. 78.

    Zhang Y, Ge C, Wang L, Liu X, Chen Y, Li M, Zhang M. Induction of DKK1 by ox-LDL negatively regulates intracellular lipid accumulation in macrophages. FEBS Lett. 2015;589(1):52–8. PMID: 25436422

  79. 79.

    Malur A, Huizar I, Wells G, Barna BP, Malur AG, Thomassen MJ. Lentivirus-ABCG1 instillation reduces lipid accumulation and improves lung compliance in GM-CSF knock-out mice. Biochem Biophys Res Commun. 2011;415:288–93.

  80. 80.

    Thomassen MJ, Barna BP, Malur AG, Bonfield TL, Farver CF, Malur A, Dalrymple H, Kavuru MS, Febbraio M. ABCG1 is deficient in alveolar macrophages of GM-CSF knockout mice and patients with pulmonary alveolar proteinosis. J Lipid Res. 2007;48:2762–8.

  81. 81.

    Baldan A, Tarr P, Vales CS, Frank J, Shimotake TK, Hawgood S, Edwards PA. Deletion of the transmembrane transporter ABCG1 results in progressive pulmonary lipidosis. J Biol Chem. 2006;281:29401–10.

  82. 82.

    Romero F, Shah D, Duong M, Penn RB, Fessler MB, Madenspacher J, Stafstrom W, Kavuru M, Lu B, Kallen CB, Walsh K, Summer R. A pneumocyte-macrophage paracrine lipid axis drives the lung toward fibrosis. Am J Respir Cell Mol Biol. 2015;53:74–86.

  83. 83.

    Baldan A, Gomes AV, Ping P, Edwards PA. Loss of ABCG1 results in chronic pulmonary inflammation. J Immunol. 2008;180:3560–8.

  84. 84.

    Baldan A, Gonen A, Choung C, Que X, Marquart TJ, Hernandez I, Bjorkhem I, Ford DA, Witztum JL, Tarling EJ. ABCG1 is required for pulmonary B-1 B cell and natural antibody homeostasis. J Immunol. 2014;193:5637–48.

  85. 85.

    Granada M, Wilk JB, Tuzova M, Strachan DP, Weidinger S, Albrecht E, Gieger C, Heinrich J, Himes BE, Hunninghake GM, Celedon JC, Weiss ST, Cruikshank WW, Farrer LA, Center DM, O'Connor GT. A genome-wide association study of plasma total IgE concentrations in the Framingham Heart Study. J Allergy Clin Immunol. 2012;129:840–5.

  86. 86.

    Mittal RA, Hammel M, Schwarz J, Heschl KM, Bretschneider N, Flemmer AW, Herber-Jonat S, Konigshoff M, Eickelberg O, Holzinger A. SFTA2--a novel secretory peptide highly expressed in the lung--is modulated by lipopolysaccharide but not hyperoxia. PLoS One. 2012;7:e40011.

  87. 87.

    Rausch F, Schicht M, Paulsen F, Ngueya I, Brauer L, Brandt W. “SP-G”, a putative new surfactant protein--tissue localization and 3D structure. PLoS One. 2012;7:e47789.

  88. 88.

    May-Simera HL, Kelley MW. Cilia, Wnt signaling, and the cytoskeleton. Cilia. 2012;1(1):7.

Download references


not applicable.


This study was supported by DST SERB: SB/SO/AS-026/2013; Department of Biotechnology, Government of India: BT/PR12987/INF/22/205/2015, VINNOVA (2016–01951) (K.G.); National Institutes of Health grants ES015675, HL077763, and HL085655 (G.D.L.); European Respiratory Society (ERS): ERS-LTRF-2015-3567 (SU); ERS-STRF-2014-7156 (AM); CSIR-SRF: [9/1045(0007) 2 K14-EMR-1], Fulbright Nehru Doctoral Research Fellowship (IIE grantee id 15,151,382 (TAT); Helmholtz Portfolio Theme ‘Metabolic Dysfunction and Common Disease’, The Helmholtz Alliance ‘Imaging and Curing Environmental Metabolic Diseases, ICEMED’ (J.B.); German Center of Lung Research (DZL) H.F., L.L.

Availability of data and materials

Microarray data has been submitted to the Genome Expression Omnibus (GEO) database at National Center for Biotechnology Information NCBI (GSE80078).

Author information

KG, MI, JB, HF, HS and GDL conceived the project and designed the experiments; LG, AM, TAT, MI, LL, DH, AMH, SU and KG performed the experiments; LG, TAT, AM, MI, SV,HF, SU, HS, GDL and KG analyzed the data; LG, MI, HF, HS, GDL and KG wrote the manuscript. All authors have read and approved the manuscript.

Correspondence to Koustav Ganguly.

Ethics declarations

Ethics approval and consent to participate

Human participants: human data or human tissue: not applicable.

Mice: The use of animals was in accordance with the German Law of Animal Protection and approved by the Bavarian Animal Research Authority and the Animal Research Authority of Schleswig-Holstein (reference number V312–72241.123-3. All procedures were also approved by IACUC of the University of Pittsburgh, PA, USA. Frozen and paraffin embedded tissues were procured to carry out experiments at SRM University, India according to the Institutional Animal Ethics Committee (IAEC) permission [79/IAEC/2013].

Consent for publication

not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1: Table S1.

Increased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at embryonic day 18. Table S2. Decreased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at embryonic stage 18. Table S3. Increased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at postnatal day 28. Table S4. Decreased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at postnatal day 28. Table S5. Increased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at postnatal day 70. Table S6. Decreased lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice at postnatal day 70. Table S7. Transcripts showing consistent pattern of expression lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice across E18/ P28 stages. Table S8. Transcripts showing consistent pattern of expression lung transcripts in JF1/MsJ (JF1) compared to C3H/HeJ (C3H) mice across P28/ P70 stages. Figure S1. Representative lung sections showing smaller alveoli in JF1/Msf (JF1) mice compared to C3H/HeJ (C3H) in both males and females. (DOCX 216 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

George, L., Mitra, A., Thimraj, T.A. et al. Transcriptomic analysis comparing mouse strains with extreme total lung capacities identifies novel candidate genes for pulmonary function. Respir Res 18, 152 (2017).

Download citation


  • Transcriptomics
  • Lung development
  • Chronic obstructive pulmonary disease
  • Asthma
  • WNT Signaling