Skip to main content

Comparative transcriptomics in human COPD reveals dysregulated genes uniquely expressed in ferrets



Chronic obstructive pulmonary disease (COPD) is a progressive lung disease with poor treatment options. However, most mouse models of COPD produce a primarily emphysematous disease not recapitulating clinically meaningful COPD features like chronic bronchitis.


Wild-type ferrets (Mustela putorius furo) were divided randomly into two groups: whole body cigarette smoke exposure and air controls. Ferrets were exposed to smoke from 1R6F research cigarettes, twice daily for six months. RNA-sequencing was performed on RNA isolated from lung tissue. Comparative transcriptomics analyses of COPD in ferrets, mice, and humans were done to find the uniquely expressed genes. Further, Real-time PCR was performed to confirmed RNA-Seq data on multiple selected genes.


RNA-sequence analysis identified 420 differentially expressed genes (DEGs) that were associated with the development of COPD in ferrets. By comparative analysis, we identified 25 DEGs that are uniquely expressed in ferrets and humans, but not mice. Among DEGs, a number were related to mucociliary clearance (NEK-6, HAS1, and KL), while others have been correlated with abnormal lung function (IL-18), inflammation (TREM1, CTSB), or oxidative stress (SRX1, AHRR). Multiple cellular pathways were aberrantly altered in the COPD ferret model, including pathways associated with COPD pathogenesis in humans. Validation of these selected unique DEGs using real-time PCR demonstrated > absolute 2-fold changes in mRNA versus air controls, consistent with RNA-seq analysis.


Cigarette smoke-induced COPD in ferrets modulates gene expression consistent with human COPD and suggests that the ferret model may be uniquely well suited for the study of aspects of the disease.


Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease diagnosed typically by irreversible airflow obstruction and persistent respiratory symptoms [1, 2]. COPD remains among the top five leading causes of death worldwide [3] and still lacks pharmacologic treatments that can mitigate COPD progression and persistent symptoms of the chronic bronchitis phenotype that include chronic cough and sputum production. Genetic studies can reveal novel mechanisms and suggest potential therapeutic strategies for investigation.

Cigarette smoking and chronic exposure to other inhaled pollutants are principal environmental risk factors driving COPD pathogenesis known to exert dysregulation of gene expression in the lung as well as extrapulmonary tissues [4,5,6]. However, transcriptomics findings in human studies cannot control for phenotypic variations that can be achieved using experimental systems, and animal studies have largely been limited to studies in mice [5, 6]. Mice are among the most utilized models to study COPD [7, 8], yet produce primarily very mild emphysematous disease and do not develop features pathognomonic for chronic bronchitis or mucus hypersecretion [9]. To remedy this, we developed a cigarette smoke exposure model in ferrets that recapitulates features of emphysema as well as COPD-associated chronic bronchitis including airway remodeling, mucous metaplasia, muco-obstructive lung disease and periodic infectious exacerbations [10,11,12,13].

In the current report, we conducted transcriptome profiling of ferrets with COPD and compared findings to prior results from human COPD lungs [14] and a mouse model of COPD [6]. Overall, we sought to identify whether our ferret model of COPD captures unique genetic signatures in comparison to mouse models that can help improve understanding of the molecular pathogenesis of human COPD and promote the development of new and effective therapies.

Materials and methods

Cigarette smoke exposure in ferrets

Age and sex-matched ferrets (Mustela putorious furo) were procured from Marshall BioResources and exposed to room air or cigarette smoke generated by an automated cigarette smoke generator (TSE Systems, Chesterfield, MO). Smoke exposure consisted of 1 h of smoke from 1R6F research cigarettes, twice daily for at least 6 months. Animals were continuously monitored for particulate matter (200 ug/L) and CO levels (~ 1% to ~ 3%). All animal experiments were approved by the Institutional Biosafety Committee of the University of Alabama at Birmingham (IACUC 20232).

RNA isolation

Following euthanasia, RNA was extracted from freshly frozen ferret lung samples using a RNeasy Mini Kit (Qiagen, Germantown, MD, USA) according to the manufacturer’s instructions. RNA integrity and concentration were examined prior to processing for RNA-seq.

Next-generation sequencing on Illumina platforms

mRNA sequencing was performed on an Illumina NextSeq 500 System as described by the manufacturer (Illumina Inc., San Diego, CA, USA). Briefly, the quality of the total RNA was assessed using the Agilent 2100 Bioanalyzer. RNA with an RNA Integrity Number (RIN) of 7.0 or above was used for sequencing library preparation. We used the Agilent SureSelect Strand Specific mRNA library kit as per the manufacturer’s instructions (Agilent, Santa Clara, CA, USA). Library construction began with two rounds of polyA selection using oligo dT containing magnetic beads. The resulting mRNA was randomly fragmented with cations and heat, which was followed by first strand synthesis using random primers with inclusion of Actinomycin D (2.4 ng/µL final concentration). Second strand cDNA production was done with standard techniques and the ends of the resulting cDNA were made blunt, A-tailed and adaptors ligated for indexing to allow for multiplexing during sequencing. The cDNA libraries were quantitated using qPCR in a Roche LightCycler 480 with the Kapa Biosystems kit for Illumina library quantitation (Kapa Biosystems, Woburn, MA, USA) prior to cluster generation. Cluster generation was performed according to the manufacturer’s recommendations for onboard clustering (Illumina). Paired end 75 bp sequencing runs were completed to allow for better alignment of the sequences to the reference genome.

RNA-seq analysis

RNA-seq data from 12 individual ferrets (six ferrets for air and smoke experimental conditions, respectively) were sequenced using an Illumina NextSeq 500 System (Additional file 1: Table S1). Paired end 75 bp reads were generated. Pre-alignment quality assessments of the raw fastq sequences were carried out using FastQC (version 0.11.7) (Andrews, 2010). The number of paired ends reads for the twelve samples range from 35 to 57 M (Additional file 2: Table S2, Additional file 3: Fig. S1). The raw fastq sequences were aligned to the Mustela putorius furo reference genome (GenBank assembly accession: GCA_000215625.1) [15]. The alignments were carried out using STAR (version 2.7.1a) [16]. (--outFilterType BySJout --outFilterMultimapNmax 20 --alignSJoverhangMin 8 --alignSJDBoverhangMin 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.04 --alignIntronMin 20 -alignIntronMax 1000000 --alignMatesGapMax 1000000; as per ENCODE STAR parameter options for RNA-seq). Post-alignment quality assessments were carried out with RSeQC (version 2.6.3) [17]. Samtools (version 0.0.19) [18] and IGV (version 2.6.2) were used for indexing and viewing the alignments respectively. Gene expression was quantified as gene level counts using the htseq-count function (version 0.12.3) [19]. The Ensembl gene annotations for Mustela putorius furo (genebuild-last-updated 2016-05) was used [20]. The htseq-count default parameters were used except for the strand parameter which was set to reverse to consider the strandedness of library. DESeq2 filtered the low counts; genes for which there are less than 3 samples with normalized counts greater than or equal to 4 were filtered out. For each sample, the read counts were normalized using the Median Ratio Normalization method. Batch effects were evaluated using the ComBat function in sva (version 3.33.1) [21]; principal component analysis plots with the normalized gene expression data were used. Differentially expressed genes (DEGs) were identified using DESeq2 (version 1.28) [22]. DESeq2 was configured by adding batch and sex to the design to control for these factor whilst DESeq2 tests for association due to the condition [22]. Genes were differentially expressed if the absolute log2Fold change was > 1, and the adjusted p-value was < 0.1. The transformed normalized gene expression data with batch removal and adjustments for sex, using the limma function removeBatchEffect (version 3.44.0) (Ritchie et al. 2015), was used for downstream analyses (such as clustering analysis and heatmap generation (Fig. 1). The complex heatmap package (version [23] was used to generate heatmaps. The volcano plots were generated with an R Bioconductor package EnhancedVolcano (version 1.12.0) [24]. Circos plots were drawn with Circos (version 0.699) [25]. Interactivenn was used to generate Venn diagrams [26].

Fig. 1
figure 1

Differential expression pattern of ferret lungs following cigarette smoke exposure. Age and sex-matched ferrets (Mustela putorious furo) were exposed to cigarette smoke or air control for 6 months, after which lung tissues were isolated and RNA was extracted. Total RNA was isolated from lung homogenate and RNA-Seq analysis was performed. A Principal component analysis of the RNA seq data showing the distinct grouping of air control and smoke-exposed ferret lung (experimental) samples. B Hierarchical clustering of differentially expressed genes (DEGs). Studies were performed in two groups of 6 synchronous ferrets, each with 3 males and 3 females per group

Functional enrichment analysis

To assess which biological functions were altered upon smoke exposure, over-representation enrichment analysis was performed using the WEB-based GEne SeT AnaLysis Toolkit (Webgestalt) [27]. The functional database selected was gene ontology (biological process); the genome was selected as the reference set. A list of gene symbols was supplied (those genes that mapped between the ferret and human [20]; the genes supplied were DEGs. The organism of interest selected was Human. The hypergeometric test was used for over-representation evaluation for the lists of statistically altered DEGs. The Benjamini and Hochberg method was used to calculate the adjusted p-values (q). Other functional databases used were the gene ontology for molecular function and the biological pathway (KEGG). The detection of over-represented gene ontologies in the differentially expressed genes was also performed using the clusterProfiler package (version 4.0) [28, 29].

Real-time PCR validation

For qRT-PCR analysis, RNA was isolated using miRNeasy Mini Kit as per manufacturer protocol (Qiagen, Cat No.: 217004, USA). 50 nanogram of RNA was used for further amplification by using TaqMan™ RNA-to-CT™ 1-Step Kit (Thermo Scientific). Analyses of mRNA levels were carried out using StepOnePlus Real-Time PCR System (Life Technologies, Carlsbad, CA, USA), TaqMan Universal PCR Mastermix (Applied Biosystems, Waltham, MA, USA) and gene-specific Prime PCR Probe Assays was designed and order from thermofisher (Assay ID: APEPYE2, APT2FKD, APU696A, APFVTYY, APGZMJW, APH6F4U, APKCAPR, APMF4AN, APNKXVK, APPRTFH and APRWKZF). GAPDH was used for reference housekeeping gene. Relative expressions were quantified with the 2−ΔΔCT method [30].

Data access

The experiments performed in this study are available in GEO under the accession number (GSE193749).


Evaluation of differentially expressed genes in smoke-exposed ferret lungs

The transcriptome profile of COPD ferrets has not been previously reported. To determine the gene expression signature of the ferret model, lung parenchyma was derived from cigarette smoke-exposed or air control ferrets and assessed using RNA-Seq. In the ferret transcriptomics data, the overall read mapping rate to the ferret reference genome was between ranged between 80 and 92% and did not significantly vary by smoking status. Principal component analysis of global transcriptomics demonstrated clear discrimination of smoke and control exposed ferrets (Fig. 1A). Hierarchical clustering of DEGs, shown by heat map after adjustment of transcript levels by batch and sex (Fig. 1B), demonstrated distinct expression patterns by exposure condition, with pronounced alterations of the transcriptome in response to smoke, in smoke exposed ferrets. Overall, a total of 420 DEGs were associated with smoke exposure in COPD ferrets (Fig. 2, Additional file 1: Table S1) Cytochrome P450 family 1 subfamily B member 1 (CYP1B1) was highly upregulated in the model (5.1-fold increase, Benjamini–Hochberg adjusted p-value = 2.41E−44), in line with previous observations in humans with COPD [31] and rodent exposure models [32], and thus served as an important quality control.

Fig. 2
figure 2

Volcano plot of differentially expressed genes (DEGs) in smoke-exposed ferrets. Differential expression volcano plot adjusted for batch and sex. Genes were differentially expressed if the absolute log2fold change > 1 and the adjusted p-value < 0.1; these are plotted as red dots. The known cigarette smoke exposure marker CYP1B1 and other genes that have biological or scientific basis in smoking or COPD are indicated by text

Gene ontology enrichment analysis

To support these pathway analyses, we further acquired an overview of Gene Ontology by evaluating the functional annotations of ferret DEGs based on GO terms. The top 15 associated functions in each component are presented in Fig. 3 and the complete annotation is provided in Additional file 4: Table S3. Disease relevant pathways that were altered in biological processes (GO: BP) prominently included neutrophil degranulation, neutrophil activation, and respiratory burst, whereas G protein-coupled receptor signaling pathway and cell adhesion were down regulated (Fig. 3A; Additional file 4: Table S3). Molecular functions (GO: MF) notable included receptor signaling pathways and oxidoreductase activity (Fig. 3B; Additional file 4: Table S3), whereas cellular component GO terms included extracellular matrix and NADPH oxidase (Fig. 3C; Additional file 4: Table S3).

Fig. 3
figure 3

Gene ontology (GO) annotation plots for differentially expressed genes (DEGs) identified in smoke-exposed ferrets. Bubble plots showing GO terms affected for A biological Process, B molecular function, and C cellular component categories, by DEG count and FDR in smoke exposed ferrets as compared to controls

Functional and gene ontology (GO) enrichment analysis of differentially expressed genes in smoke-exposed ferret lungs

DEGs associated with cigarette smoke exposure identified from the analysis of ferret lung were next subjected to functional enrichment analysis. Top ten upregulated and downregulated KEGG pathways are shown in Table 1 in rank order by p-value, noting many were not statistically significant by FDR.

Table 1 KEGG pathway analysis showing up-regulated (A) and down-regulated (B) pathways

Comparison of differentially expressed genes in smoke-exposed ferret lungs, mouse lungs, and human COPD

A unique property of the COPD ferret model is that it exhibits classical pathophysiological features of airway disease, including goblet cell hyperplasia, mucus hyperexpression, and delayed mucus clearance, as well phenotypical features of chronic bronchitis including cough [10,11,12]. We next sought to understand whether COPD ferrets can provide insights into the molecular basis of human COPD that otherwise cannot be examined in COPD mouse models. We next sought to evaluate the extent of overlap between the genetic signatures of DEGs in ferrets as compared to human and/or mice. To do that, overlap was assessed between the COPD ferret gene expression signature and DEGs in datasets from Obeidat et al. evaluating mice that were smoke-exposed for 24 weeks duration [6] and Bosse et al. that evaluated humans with severely affected COPD [14], a population most relevant to the duration and intensity of chronic smoke exposure experienced by ferrets [12]. Of the 420 DEGs identified in COPD ferrets, 266 DEGs also had analyzable gene expression data generated in human [14] and mouse [6] models enabling comparative transcriptomic analyses between ferret, human, and mouse lungs. As depicted in Fig. 4A Venn diagram and detailed in the Circos plot in Fig. 4B that illustrates commonalities between the ferret, human, and mouse genomes, 52 (18.5%) of these 266 DEGs in ferrets were common to all three species, whereas 77 (28.9%) were also differentially expressed in humans, and 90 (33.8%) were differentially expressed in mice (Top 25 listed in Table 2; Complete list in Additional file 5: Table S4). Interestingly, 25 DEGs were differentially expressed in human and ferret, but not mice (Table 3); among these genes uniquely modeled in ferrets were those previously determined to be relevant to COPD pathogenesis, with upregulation observed for IL-18, KL, SRXN1, TREM1 and downregulation observed for HAS1 and PNPLA3 [33,34,35,36]. Of note, the change in KL and PLPLA3 expression in ferret were discordant with changes in humans found by Bosse et al. Among the 25 most statistically significant DEGs common in ferret and human, three other genes also showed discordant regulation between the two species: GZMA, SNAP25, and EGFLAM (Table 3). These genes are important in airway goblet cells, although studies of its expression and localization in ferrets have not been performed.

Fig. 4
figure 4

A Venn diagram depicting overlap between differentially expressed genes (DEGs) in smoke-exposed ferret, human COPD, and smoke-exposed murine lungs. The number of DEGs in common between smoke-exposed ferrets, human COPD [14], and smoke-exposed mice [6] are shown. B Circos plot of smoking-related genes overlapping between smoke-exposed ferret, human COPD, and smoke-exposed murine lungs. The gene chromosomal positions are based on the human genome and shown inside the circles. The first pair of circles (light blue and dark blue) represents the dysregulated genes from the long term (24 weeks) smoke-exposed mouse [6]. The upregulated genes are plotted in the dark blue circle pointing outward whilst the downregulated genes are plotted in the light blue circle pointing inward. The second pair of circles (light red and dark red) show the ferret genes dysregulated due to smoke exposure. The upregulated genes are in the dark red circle and the downregulated genes are in the light red circle. The third pair of circles (light green and dark green) show the “human smoking” gene signature [6, 14]. The upregulated genes are in the dark green circle, whilst the downregulated genes are in the light green circle. The length of the line is proportional to the −log10 FDR values for differential expression in human and for the −log10 FDR values in murine and for the −log10 padj for the ferret. Genes highlighted are statistically significant from the respective studies

Table 2 The top 25 smoke related DEGs identified in the ferret, mouse and human lungs shown in order of statistically significance for the ferret study
Table 3 25 Smoke-related DEGs identified in the ferret and human lungs shown in order of statistically significance for the ferret study

RT-PCR validation analysis of key genes in COPD ferret model

To validate RNA-Seq results, we conducted RT-PCR using ferret tissues for selected DEGs (Fig. 5) that were uniquely modeled in ferrets, in addition to CYP1B1 as an exposure control, MMP12 as an important gene associated with COPD progression, and HAS1 as a disease of interest related to mucosal defense. Quantification of mRNA expression in smoke-exposed vs. air control ferrets by RT-PCR demonstrated expected changes with smoke exposure (Fig. 5), generally consistent with the RNA-seq analysis and providing validation for the utility of COPD ferrets in modeling a set of genes that can uniquely inform on the pathogenesis of COPD.

Fig. 5
figure 5figure 5

RT-PCR validation of differentially expressed genes (DEGs) unique to smoke-exposed ferret and human COPD lungs as compared to raw data from RNA-Seq. RT-PCR (right) and RNA-Seq (left) measurements for A CYP1B1; B MMP12; C IL-18; D KL; E TREM1; F TREM144; G SRXN1; H ATP6V0D2; I PNPLA3; and J HAS1 in smoke-exposed and air control ferret lungs. Mean ± SEM. RT PCR N = 6/condition, 3 males and 3 females per group, comparison by Mann–Whitney


Various studies using cigarette smoke induced animal models have demonstrated genes associated with smoking by transcriptomic analysis [4, 5]. However, comparative studies using smoke exposed mouse models of COPD have limitations, since mice do not exhibit all features of the complex human COPD sub-phenotypes. In this study, we assessed the transcriptome profile of lung parenchyma in a ferret model of cigarette smoke-induced COPD. Gene expression profiling has not been previously reported in COPD ferrets, which uniquely exhibit features of airway remodeling in comparison to rodent models of cigarette smoke exposure, while also exhibiting emphysema [10,11,12,13]. A total of 420 DEGs associated with smoke exposure were identified, many of which have been observed in mice, humans, or both. Furthermore, building upon a recent study by Obeidat et al. who compared the expression profiles of smoke exposed (both short term and long term) murine model with COPD human lungs [6], we identified multiple additional DEGs and associated pathways bearing pathophysiological relevance to COPD that were uniquely expressed in ferrets and humans, but not mice. To our knowledge, this is currently the only RNA-seq based transcriptomics analysis in the COPD ferret model, and due to the unique features of ferrets, provides opportunities for detection of novel pathways with human relevance.

A major goal of the analysis was to identify gene targets that might be uniquely modeled in ferrets, given the human-like phenotype it exhibits. Our analysis revealed 25 aberrantly regulated genes that were exclusive to humans and the ferret COPD model, but not mice. These included several DEGs that have previously been studied for their role in COPD and that may be implicated in mucus-related aspects of the disease. IL-18 has been correlated with pulmonary function in COPD patients [37] and, in an IL-18 transgenic mouse model, was identified as key driver regulating downstream responses including mucus metaplasia, airway fibrosis, and vascular remodeling [38]. KL-deficient mice have been observed to exhibit delayed mucociliary clearance [34], and downregulation of KL protein KLOTHO has been widely observed in COPD with cigarette exposure [39, 40]. That we observed increased expression of KL may indicate that smoke exposed ferrets are responding to delayed mucociliary clearance, as recently shown [34], The discordant expression of KL in ferrets compared to humans could be due to the possibility that higher expression of KL acts as a compensatory response to oxidative stress due to smoking, a finding reported by other studies [41], including in human smokers without COPD [42]. Other DEGs common to humans and ferret of interest included the antioxidant enzyme sulfiredoxin-1 (SRX1), which was shown to be upregulated in the lungs of nuclear factor erythroid-2-related factor 2 (Nrf+/+) mice upon acute exposure to cigarette smoke, conferring a protective effect against oxidative stress, although was decreased in lung tissue from patients with advanced COPD [43]. This may have been prominent since ferrets exhibit a mild form of disease, given their age and relative duration of cigarette smoke exposure. Upregulation of triggering receptor expressed on myeloid cells 1 (TREM1) has been observed in cigarette-smoke mice, exacerbating pyroptosis-mediated lung injury and inflammation [44]. These pathways could represent opportunities for further investigation in COPD pathogenesis using the ferret model.

Prior investigations evaluating gene expression profiles of smoked-exposed mouse [6] and human COPD [14] lungs enabled cross-species transcriptomics analysis for 266 of the 420 total DEGs detected in the COPD ferret model. There was substantial overlap in DEGs detected across all three species, with involvement seen in genes with widely studied roles in COPD such as CYP1B1 [31], matrix metalloproteinase 12 (MMP12) [45], and aryl hydrocarbon receptor repressor (AHRR) [46], underscoring the validity of the smoke exposure model in ferrets for evaluation of the lung transcriptome. Common DEGs between ferret, human, and mouse also reinforced the contribution of additional genes to COPD pathogenesis. Among other overlapping DEGs, TREM2 has been observed to promote macrophage survival in the lung [47], nimA-related protein kinase 6 (NEK6) has been implicated in COPD-associated ciliopathies [48], and cathepsin B (CTSB) induced emphysema and secretory cell hyperplasia upon intratracheal administration to hamster lungs [49]. Additional genes of interest that were substantially upregulated and have direct relevance to COPD and/or cigarette smoke exposure included SERPINA1 (2.17 log2 fold change (log2FC); padj = 5.13E−26), ATP6V0D2 (3.99 log2FC; padj = 3.42E−26), MMP-12 (2.69 log2FC ; padj = 6.20E−12), AHRR (1.74 log2FC; padj = 8.32E−05), and CTSB (1.83 log2FC; padj = 4.64E−18), among others [4, 50].

Downregulation of hyaluronic acid synthase 1 (HAS1) has been associated with abnormally low production of the high molecular weight form of HA (HMW-HA) in COPD [51]. While RNA-Seq revealed a trend, RT-PCR, conducted because of prior interest in this area, indicated downregulation of HAS1. Reduced HAS1, and thus lower HMW-HA, could also lead to abnormal HMW/LMW HA ratios, which would alter mucus clearance and also induce inflammation [51].

Multiple cellular pathways were aberrantly altered in the COPD ferret model, including pathways that are also associated with COPD pathogenesis in humans. Of the pathways were significant changes were observed (Table 1 and Additional file 4: Table S3), those that have been reported to be involved with COPD pathogenesis included activation of mitogen-activated protein kinase (MAPK), phosphatidylinositol 3 kinase/serine threonine protein kinase (PI3K-AKT), phagosome, cytokine-cytokine receptor interaction, extracellular matrix (ECM)-receptor interaction, chemokines, apoptotic signaling pathways, inhibition of transforming growth factor beta (TGF-ß), peroxisome proliferator-activated receptors (PPAR), and inflammatory mediator regulation of transient receptor potential (TRP) channel pathways [52].

Our study, while unique in the analysis of ferrets, is not without limitations. Ferrets are an outbred species, and thus genetic heterogeneity may have reduced the ability to detect some important DEGs, especially given the sample size encountered. We focused our analysis on parenchyma tissue, noting this was most similar to the specimens collected in humans and other species for which we compared, and this samples also include the small airways, where airway disease is prominent [12, 53]; nevertheless, the analysis of airway epithelia could be informative in the future, especially given pathological changes in the large airways of ferrets upon smoking [10,11,12]. Since ferrets are sexually dimorphic, it is not surprising that there were important differences by ferret sex, although this was controlled for in the differential gene analysis and included where possible in comparison to other data sets. Importantly, analysis of RNA has evolved throughout the study, providing a limitation in the inter-species comparisons. The RNA-Seq analysis and RT-PCR validation in this study was based on mRNA expression; validation by protein analysis will be useful in the future, noting some protein reagents in ferrets need further characterization. To assure specificity of our comparisons, we required at least one-fold change, in addition to statistical significance for the DEG to be included in the ferret study, which differed from those included in mice-human comparators, In a minority of cases for genes uniquely modeled in ferrets, the directional change in differential expression in ferrets and humans were discordant. Although this can be explained for KL due to the potential for evolving changes over time, further research will be needed for other genes where this was the case, such as GZMA, SNAP25, KL, PNPLA3 and EGFLAM.


In summary, our analysis provides a large RNA-seq based transcriptome analysis set in COPD and involves a ferret model that uniquely exhibits features of chronic bronchitis. Our data revealed 420 DEGs, including multiple genes established as abnormal among COPD patients, such as pathways involving neutrophil activation, immune regulation, and oxidative metabolism. We also identified 25 DEGs that are uniquely expressed in ferrets and humans, but not mice, some of which implicate abnormal mucus clearance, a finding uniquely modeled by ferrets. Selected genes were validated with real time PCR, and these were concordant with the RNA-Seq expression. Together, the transcriptomic profile and comparison with other available models of COPD identifies a unique set of genes that could be relevant to disease pathogenesis. This high-resolution view of the lung transcriptome associated with COPD may ultimately provide a valuable resource for the research community and warrants further mechanistic investigations.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.



Chronic obstructive pulmonary disease


Quantitative polymerase chain reaction


Differentially expressed genes


Kyoto Encyclopedia of Genes and Genomes


Cytochrome P450 family 1 subfamily B member 1


Gene Ontology


Biological process


Molecular functions






Triggering receptor expressed on myeloid cells 1


Nuclear factor erythroid-2-related factor 2


Matrix metalloproteinase 12


Aryl hydrocarbon receptor repressor


NimA-related protein kinase 6


Hyaluronic acid synthase 1


High molecular weight form of HA


Mitogen-activated protein kinase


Phosphatidylinositol 3 kinase/serine-threonine protein kinase


Extracellular matrix


Transforming growth factor beta


Peroxisome proliferator-activated receptors


Transient receptor potential


  1. Rabe KF, Watz H. Chronic obstructive pulmonary disease. Lancet. 2017;389:1931–40.

    Article  PubMed  Google Scholar 

  2. Chand HS, Vazquez-Guillamet R, Royer C, Rudolph K, Mishra N, Singh SP, Hussain SS, Barrett E, Callen S, Byrareddy SN, et al. Cigarette smoke and HIV synergistically affect lung pathology in cynomolgus macaques. J Clin Invest. 2018;128:5428–33.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Chronic obstructive pulmonary disease (COPD). World Health Organization. Accessed 29 Oct 2021.

  4. Morissette MC, Lamontagne M, Berube JC, Gaschler G, Williams A, Yauk C, Couture C, Laviolette M, Hogg JC, Timens W, et al. Impact of cigarette smoke on the human and mouse lungs: a gene-expression comparison study. PLoS ONE. 2014;9: e92498.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Yun JH, Morrow J, Owen CA, Qiu W, Glass K, Lao T, Jiang Z, Perrella MA, Silverman EK, Zhou X, Hersh CP. Transcriptomic analysis of lung tissue from cigarette smoke-induced emphysema murine models and human chronic obstructive pulmonary disease show shared and distinct pathways. Am J Respir Cell Mol Biol. 2017;57:47–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Obeidat M, Dvorkin-Gheva A, Li X, Bosse Y, Brandsma CA, Nickle DC, Hansbro PM, Faner R, Agusti A, Pare PD, et al. The overlap of lung tissue transcriptome of smoke exposed mice with human smoking and COPD. Sci Rep. 2018;8:11881.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Stevenson CS, Birrell MA. Moving towards a new generation of animal models for asthma and COPD with improved clinical relevance. Pharmacol Ther. 2011;130:93–105.

    Article  CAS  PubMed  Google Scholar 

  8. Vlahos R, Bozinovski S. Recent advances in pre-clinical mouse models of COPD. Clin Sci (Lond). 2014;126:253–65.

    Article  CAS  Google Scholar 

  9. Wright JL, Cosio M, Churg A. Animal models of chronic obstructive pulmonary disease. Am J Physiol Lung Cell Mol Physiol. 2008;295:L1-15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Stanford D, Kim H, Bodduluri S, LaFontaine J, Byzek SA, Schoeb TR, Harris ES, Nath HP, Bhatt SP, Raju SV, Rowe SM. Airway remodeling in ferrets with cigarette smoke induced COPD using microCT imaging. Am J Physiol Lung Cell Mol Physiol. 2020;319:L11–20.

    Article  CAS  PubMed Central  Google Scholar 

  11. Lin VY, Kaza N, Birket SE, Kim H, Edwards LJ, LaFontaine J, Liu L, Mazur M, Byzek SA, Hanes J, et al. Excess mucus viscosity and airway dehydration impact COPD airway clearance. Eur Respir J. 2020;55: 1900419.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Raju SV, Kim H, Byzek SA, Tang LP, Trombley JE, Jackson P, Rasmussen L, Wells JM, Libby EF, Dohm E, et al. A ferret model of COPD-related chronic bronchitis. JCI Insight. 2016;1: e87536.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Basil MC, Cardenas-Diaz FL, Kathiriya JJ, Morley MP, Carl J, Brumwell AN, Katzen J, Slovik KJ, Babu A, Zhou S, et al. Human distal airways contain a multipotent secretory cell that can regenerate alveoli. Nature. 2022;604:120–6.

    Article  CAS  PubMed  Google Scholar 

  14. Bosse Y, Postma DS, Sin DD, Lamontagne M, Couture C, Gaudreault N, Joubert P, Wong V, Elliott M, van den Berge M, et al. Molecular signature of smoking in human lung tissues. Cancer Res. 2012;72:3753–63.

    Article  CAS  PubMed  Google Scholar 

  15. Peng X, Alfoldi J, Gori K, Eisfeld AJ, Tyler SR, Tisoncik-Go J, Brawand D, Law GL, Skunca N, Hatta M, et al. The draft genome sequence of the ferret (Mustela putorius furo) facilitates study of human respiratory disease. Nat Biotechnol. 2014;32:1250–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    Article  CAS  PubMed  Google Scholar 

  17. Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28:2184–5.

    Article  CAS  PubMed  Google Scholar 

  18. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome Project Data Processing S: the sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.

    Article  CAS  PubMed  Google Scholar 

  20. Yates AD, Achuthan P, Akanni W, Allen J, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, Azov AG, Bennett R, et al. Ensembl 2020. Nucleic Acids Res. 2020;48:D682–8.

    CAS  PubMed  Google Scholar 

  21. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28:882–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32:2847–9.

    Article  CAS  PubMed  Google Scholar 

  24. Blighe KRS, Lewis M, EnhancedVolcano: publication-ready volcano plots with enhanced colouring and labeling. 2018.

  25. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Heberle H, Meirelles GV, da Silva FR, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015;16:169.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017;45:W130–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, Feng T, Zhou L, Tang W, Zhan L, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (N Y). 2021;2: 100141.

    CAS  Google Scholar 

  30. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

  31. Yang D, Yan Y, Hu F, Wang T. CYP1B1, VEGFA, BCL2, and CDKN1A affect the development of chronic obstructive pulmonary disease. Int J Chron Obstruct Pulmon Dis. 2020;15:167–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kim YS, Kokturk N, Kim JY, Lee SW, Lim J, Choi SJ, Oh W, Oh YM. Gene profiles in a smoke-induced COPD mouse lung model following treatment with mesenchymal stem cells. Mol Cells. 2016;39:728–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Nakanishi K, Nishida M, Taneike M, Yamamoto R, Moriyama T, Yamauchi-Takihara K. Serum Klotho levels contribute to the prevention of disease progression. Int J Gen Med. 2021;14:229–36.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Garth J, Easter M, Skylar Harris E, Sailland J, Kuenzi L, Chung S, Dennis JS, Baumlin N, Adewale AT, Rowe SM, et al. The effects of the anti-aging protein klotho on mucociliary clearance. Front Med (Lausanne). 2019;6:339.

    Article  Google Scholar 

  35. Papakonstantinou E, Roth M, Klagas I, Karakiulakis G, Tamm M, Stolz D. COPD exacerbations are associated with proinflammatory degradation of hyaluronic acid. Chest. 2015;148:1497–507.

    Article  PubMed  Google Scholar 

  36. Galdi F, Pedone C, McGee CA, George M, Rice AB, Hussain SS, Vijaykumar K, Boitet ER, Tearney GJ, McGrath JA, et al. Inhaled high molecular weight hyaluronan ameliorates respiratory failure in acute COPD exacerbation: a pilot study. Respir Res. 2021;22:30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Rovina N, Dima E, Gerassimou C, Kollintza A, Gratziou C, Roussos C. Interleukin-18 in induced sputum: association with lung function in chronic obstructive pulmonary disease. Respir Med. 2009;103:1056–62.

    Article  PubMed  Google Scholar 

  38. Nakajima T, Owen CA. Interleukin-18: the master regulator driving destructive and remodeling processes in the lungs of patients with chronic obstructive pulmonary disease? Am J Respir Crit Care Med. 2012;185:1137–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Krick S, Grabner A, Baumlin N, Yanucil C, Helton S, Grosche A, Sailland J, Geraghty P, Viera L, Russell DW, et al. Fibroblast growth factor 23 and Klotho contribute to airway inflammation. Eur Respir J. 2018;52:1800236.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Gao W, Yuan C, Zhang J, Li L, Yu L, Wiegman CH, Barnes PJ, Adcock IM, Huang M, Yao X. Klotho expression is reduced in COPD airway epithelial cells: effects on inflammation and oxidant injury. Clin Sci (Lond). 2015;129:1011–23.

    Article  CAS  Google Scholar 

  41. Verde Z, Gonzalez-Moro JM, Chicharro LM, Reinoso-Barbero L, Bandres F, Gomez-Gallego F, Santiago C. A paradox: alpha-Klotho levels and smoking intensity. Lung. 2017;195:53–7.

    Article  CAS  PubMed  Google Scholar 

  42. Nakanishi K, Nishida M, Harada M, Ohama T, Kawada N, Murakami M, Moriyama T, Yamauchi-Takihara K. Klotho-related molecules upregulated by smoking habit in apparently healthy men: a cross-sectional study. Sci Rep. 2015;5:14230.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Singh A, Ling G, Suhasini AN, Zhang P, Yamamoto M, Navas-Acien A, Cosgrove G, Tuder RM, Kensler TW, Watson WH, Biswal S. Nrf2-dependent sulfiredoxin-1 expression protects against cigarette smoke-induced oxidative stress in lungs. Free Radic Biol Med. 2009;46:376–86.

    Article  CAS  PubMed  Google Scholar 

  44. Wang L, Chen Q, Yu Q, Xiao J, Zhao H. TREM-1 aggravates chronic obstructive pulmonary disease development via activation NLRP3 inflammasome-mediated pyroptosis. Inflamm Res. 2021;70:971–80.

    Article  CAS  PubMed  Google Scholar 

  45. Churg A, Zhou S, Wright JL. Matrix metalloproteinases in COPD. Eur Respir J. 2012;39:197.

    Article  CAS  PubMed  Google Scholar 

  46. Kodal JB, Kobylecki CJ, Vedel-Krogh S, Nordestgaard Børge G, Bojesen SE. AHRR hypomethylation, lung function, lung function decline and respiratory symptoms. Eur Respir J. 2018;51:1701512.

    Article  PubMed  Google Scholar 

  47. Wu K, Byers DE, Jin X, Agapov E, Alexander-Brett J, Patel AC, Cella M, Gilfilan S, Colonna M, Kober DL, et al. TREM-2 promotes macrophage survival and lung disease after respiratory viral infection. J Exp Med. 2015;212:681–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Perotin J-M, Polette M, Deslée G, Dormoy V. CiliOPD: a ciliopathy-associated COPD endotype. Respir Res. 2021;22:74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Lesser M, Padilla ML, Cardozo C. Induction of emphysema in hamsters by intratracheal instillation of cathepsin B. Am Rev Respir Dis. 1992;145:661–8.

    Article  CAS  PubMed  Google Scholar 

  50. Kodal JB, Kobylecki CJ, Vedel-Krogh S, Nordestgaard BG, Bojesen SE. AHRR hypomethylation, lung function, lung function decline and respiratory symptoms. Eur Respir J. 2018;51:1701512.

    Article  PubMed  Google Scholar 

  51. Klagas I, Goulet S, Karakiulakis G, Zhong J, Baraket M, Black JL, Papakonstantinou E, Roth M. Decreased hyaluronan in airway smooth muscle cells from patients with asthma and COPD. Eur Respir J. 2009;34:616–28.

    Article  CAS  PubMed  Google Scholar 

  52. Hussain SS, George S, Singh S, Jayant R, Hu CA, Sopori M, Chand HS. A small molecule BH3-mimetic Suppresses cigarette smoke-induced mucous expression in airway epithelial cells. Sci Rep. 2018;8:13796.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Kaza N, Lin VY, Stanford D, Hussain SS, Libby EF, Kim H, Borgonovi M, Conrath K, Mutyam V, Byzek SA, et al. Evaluation of a novel CFTR potentiator in copd ferrets with acquired cftr dysfunction. Eur Respir J. 2021;60:2101581.

    Article  Google Scholar 

Download references


Not applicable.


These studies were supported by NIH National Heart, Lung, and Blood Institute (NHLBI) Grant R35 HL135816 (to S.M.R) and National Institute of Diabetes and Digestive and Kidney Diseases Grant P30 DK072482 (to S.M.R). Funded in part by a grant from the Alpha-1 Foundation 919960 (To S.S.H).

Author information

Authors and Affiliations



SSH, DDS, SMR conceived of and designed the research; SSH, YJKE, DS, performed the experiments; SSH, YJKE, EFL, MLM, DS, DDS, and SMR analyzed data and interpreted results; SSH, YJKE, EFL, and SMR prepared figures; SSH, YJKE, EFL, and SMR wrote the manuscript; YJKE, SVR and SMR supervised the project; all authors had an opportunity to edit the manuscript and approved of its submission. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Steven M. Rowe.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were approved by the Institutional Biosafety Committee of the University of Alabama at Birmingham (IACUC 20232).

Consent for publication

Not applicable.

Competing interests

The author(s) declare(s) that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1.

RNA-seq data for 12 samples (six biological replicates each for air and smoke experimental conditions) sequenced using an Illumina NextSeq 500 System. [GEO Accession Number GSE193749].

Additional file 2: Fig. S1.

The proportion of paired ends reads for the ferret lung samples, by individual ferrets.

Additional file 3: Table S2.

The number of paired ends reads for the ferret lung samples.

Additional file 4: Table S3.

Complete annotation GO TERMS in smoke exposed ferrets compared to controls.

Additional file 5: Table S4.

Complete list (52) of DEGs in common between ferret, mice, and human transcriptome.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hussain, S.S., Edwards, Y.J.K., Libby, E.F. et al. Comparative transcriptomics in human COPD reveals dysregulated genes uniquely expressed in ferrets. Respir Res 23, 277 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: