Smoking accelerates aging of the small airway epithelium
- Matthew S Walters1,
- Bishnu P De1,
- Jacqueline Salit1,
- Lauren J Buro-Auriemma1,
- Timothy Wilson1,
- Allison M Rogalski1,
- Lindsay Lief1,
- Neil R Hackett1,
- Michelle R Staudt1,
- Ann E Tilley1,
- Ben-Gary Harvey1,
- Robert J Kaner1,
- Jason G Mezey1,
- Beth Ashbridge1, 2,
- Malcolm A S Moore2 and
- Ronald G Crystal1Email author
© Walters et al.; licensee BioMed Central Ltd. 2014
Received: 3 December 2013
Accepted: 31 July 2014
Published: 24 September 2014
Aging involves multiple biologically complex processes characterized by a decline in cellular homeostasis over time leading to a loss and impairment of physiological integrity and function. Specific cellular hallmarks of aging include abnormal gene expression patterns, shortened telomeres and associated biological dysfunction. Like all organs, the lung demonstrates both physiological and structural changes with age that result in a progressive decrease in lung function in healthy individuals. Cigarette smoking accelerates lung function decline over time, suggesting smoking accelerates aging of the lung. Based on this data, we hypothesized that cigarette smoking accelerates the aging of the small airway epithelium, the cells that take the initial brunt of inhaled toxins from the cigarette smoke and one of the primary sites of pathology associated with cigarette smoking.
Using the sensitive molecular parameters of aging-related gene expression and telomere length, the aging process of the small airway epithelium was assessed in age matched healthy nonsmokers and healthy smokers with no physical manifestation of lung disease or abnormalities in lung function.
Analysis of a 73 gene aging signature demonstrated that smoking significantly dysregulates 18 aging-related genes in the small airway epithelium. In an independent cohort of male subjects, smoking significantly reduced telomere length in the small airway epithelium of smokers by 14% compared to nonsmokers.
These data provide biologic evidence that smoking accelerates aging of the small airway epithelium.
KeywordsAging Small airway epithelium Telomere Smoking
Aging is a complex natural process that affects all living organisms, and can be broadly defined as a progressive decline in homeostasis of biological processes that leads to a loss and impairment of physiological integrity and function with a subsequent increased risk of disease and death ,. The biologic mechanisms of aging are only partly understood, but include genomic instability, telomere attrition, cellular senescence, epigenetic alterations, loss of protein homeostasis (proteostasis), dysregulated nutrient sensing, mitochondrial dysfunction, altered intercellular communication and stem cell exhaustion ,. Similar to other organs, the lung demonstrates a progressive decrease in lung function in healthy individuals ,. Epidemiological studies have demonstrated that smoking results in accelerated lung function decline over time, suggesting smoking accelerates aging of the lung ,. In support of this concept, telomere length analysis of blood leukocytes has demonstrated that smoking is associated with shortened telomeres -.
Based on the data indicating that smoking accelerates the aging-related decline in lung function, and that smokers have shorter telomere lengths in blood leukocytes, we hypothesized that cigarette smoking accelerates aging of the small airway epithelium, the cell population directly exposed to the inhaled toxins from the cigarette smoke, and one of the primary sites in the lung to develop histologic abnormalities associated with cigarette smoking -. To assess the effect of cigarette smoking on small airway epithelium aging, we compared age matched healthy nonsmokers and healthy smokers with no clinical manifestation of lung disease using the sensitive biologic parameters of gene expression and telomere length. Analysis of small airway epithelium gene expression levels demonstrated that smoking significantly dysregulates expression of an aging gene signature . Further, telomere length analysis in an independent cohort of age matched male subjects revealed that smoking significantly reduced telomere length in the small airway epithelium. Overall, these data provide biologic evidence that before the manifestation of lung function decline, smoking accelerates aging of the small airway epithelium by dysregulation of age-related gene expression and enhanced telomere erosion.
Demographics of the study population and biologic samples for analysis of aging signature gene expression in the small airway epithelium 1
38.0 ± 12.4
38.8 ± 10.1
Smoking history (pack-yr)
17.8 ± 8.4
Urine nicotine (ng/ml)
1657 ± 1527
Urine cotinine (ng/ml)
1929 ± 1060
Pulmonary function parameters3
106.9 ± 12.1
111.1 ± 11.6
106.3 ± 12.2
109.4 ± 13.1
82.8 ± 4.9
81.0 ± 4.3
99.7 ± 17.3
97.6 ± 14.1
90.4 ± 9.6
86.0 ± 5.4
Body mass index (BMI)
26.4 ± 5.3
27.8 ± 4.8
% epithelial cells
98.9 ± 0.9
99.2 ± 0.8
% inflammatory cells
1.0 ± 0.8
0.8 ± 0.8
Differential cell count5
71.3 ± 5.3
62.9 ± 6.9
10.0 ± 4.5
12.3 ± 6.4
9.5 ± 6.6
9.6 ± 6.6
Undifferentiated columnar (%)
8.2 ± 4.9
14.5 ± 8.2
Demographics of the study population and biologic samples for telomere length analysis in the small airway epithelium 1
39.0 ± 12.9
43.9 ± 7.6
Smoking history (pack-yr)
22.6 ± 9.0
Urine nicotine (ng/ml)
1424 ± 1432
Urine cotinine (ng/ml)
1359 ± 993
Pulmonary function parameters3
103.1 ± 10.9
110.4 ± 8.5
105.4 ± 11.1
107.5 ± 8.2
83.9 ± 4.8
79.1 ± 3.7
102.0 ± 14.6
96.2 ± 12.7
93.3 ± 11.8
87.2 ± 7.7
Body mass index (BMI)
27.8 ± 4.5
26.0 ± 3.9
% epithelial cells
97.5 ± 1.4
99.1 ± 0.8
% inflammatory cells
2.5 ± 1.4
0.9 ± 0.8
Differential cell count5
62.3 ± 5.8
57.9 ± 7.1
11.0 ± 5.1
17.3 ± 6.0
5.8 ± 4.2
5.6 ± 4.5
Undifferentiated columnar (%)
18.4 ± 5.8
18.3 ± 8.6
Small airway epithelial sampling
Small airway epithelium (10th to 12th order) was collected by fiberoptic bronchoscopy by brushing as previously described ,,. Following withdrawal of the bronchoscope, the cells were dislodged from the brush by flicking the brush tip in 5 ml of ice-cold bronchial epithelium basal medium (BEBM, Lonza, Basel, Switzerland). An aliquot of all airway epithelial samples was used to quantify the total number of cells recovered, and to quantify the percentage of epithelial and inflammatory cells and the proportions of epithelial cell subtypes. Cells from a second aliquot were pelleted for either RNA or DNA extraction for subsequent transcriptome or telomere length analysis, respectively.
RNA and microarray processing
Pelleted small airway epithelial cells obtained by bronchoscopy were immediately processed for RNA extraction and microarray analysis as previously described ,,. Total RNA was extracted from the small airway epithelium of 29 nonsmokers and 29 smokers (Table 1) using the TRIzol method (Invitrogen) with subsequent clean-up using the RNeasy MinElute RNA purification kit (Qiagen, Valencia, CA) and stored in RNAsecure (Ambion, Austin, TX) at -80°C. RNA quality was assessed using a Bioanalyzer (Agilent Technologies, Santa Clara, CA) to determine RNA integrity (RIN value) with only samples displaying a RIN >6 included in this study and concentration quantified with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE). One μg of total RNA was prepared for microarray genome-wide transcriptome analysis using the 3 IVT Express Kit (Affymetrix, Santa Clara, CA) and samples were submitted to the Genomics Resources Core Facility at Weill Cornell Medical College for hybridization, washing and scanning of the Affymetrix HG U133 Plus 2.0 microarray GeneChips. All HG-U133 Plus 2.0 microarrays were processed according to Affymetrix protocols, hardware and software, including being processed by the Affymetrix fluidics station 450 and hybridization oven 640 and scanned with an Affymetrix Gene Array Scanner 3000 7G. The captured image data from the HG-U133 Plus 2.0 arrays were processed using the MAS5 algorithm in GeneSpring version 7.3 (Affymetrix Microarray Suite Version 5) and normalized per chip to the median expression value of each sample to generate P calls for filtering. Overall microarray quality was verified by the criteria: (i) 3’/5’ ratio for GAPDH ≤3; and (ii) scaling factor ≤10.0. CEL files were processed by Partek Genomics Suite Software version 6.6 (Partek, Inc., St Louis, MO) for quality control, identification of outliers and determination of expression level for all probesets, using the Robust Multi-chip Average (RMA) method with Partek default parameters.
Microarray data analysis
Genes containing U133 2.0 microarray probesets and an Affymetrix P call in ≥10% of samples were considered to be expressed in the small airway epithelium. To confirm the genome wide transcriptional changes in response to smoking the phenotypes were evaluated in Partek for sources of variation and an ANOVA was performed to assess the effect of smoking on gene expression corrected for all covariates, including age, gender, race and region (small airway epithelium from distal right vs distal left lower lobe of the lung). Smoking-dependent differences in expression identified by Students t test and a p < 0.05 (using a Benjamini-Hochberg correction to limit the false discovery rate) were considered to represent a significant difference in expression between groups. Global patterns of expression were examined using a euclidean dissimilarity and average linkage hierarchical cluster analysis. To assess the effect of smoking on aging of the small airway epithelium, we compared healthy smoker vs healthy nonsmoker small airway epithelium expression of a 73 gene aging signature identified by de Magalhaes et al.  by meta-analysis of age-related gene expression profiles using 27 datasets from mice, rats and humans representing a common signature of aging between different species and tissues. Using criteria described above to define small airway epithelial expression resulted in removal of 6 genes from the analysis with a final aging signature of 67 genes (Additional file 1: Table S1). The phenotypes were evaluated in Partek for sources of variation and an ANOVA was performed to assess the effect of smoking on gene expression corrected for all covariates, including age, gender, race and region (small airway epithelium from distal right vs distal left lower lobe of the lung). Smoking-dependent differences in expression identified by Students t test and a p < 0.05 (using a Benjamini-Hochberg correction to limit the false discovery rate) were considered to represent a significant difference in expression between groups. Global patterns of expression of the 67 aging-related genes were examined using principal component analysis (PCA). The raw data and full dataset are publically available at the Gene Expression Omnibus (GEO) site (http://www.ncbi.nlm.nih.gov/geo/), accession number GSE52237.
To analyze expression of aging genes in the small airway epithelium at the protein level immunohistochemistry analysis was performed. Commercially available biopsy sections of paraffin embedded normal human bronchus tissue (Cat No: HuFPT111, US Biomax, Inc., Rockville, MD) were purchased for 3 independent donors (Donor 1: 21 yr old female, Donor 2: 16 yr old female and Donor 3: 37 yr old male). The samples were first cleaned in xylene and rehydrated with graded ethanol followed by steaming the samples for 15 min in citrate buffer solution (Thermo Scientific), followed by cooling at 23°C for 20 min to recover antigens and enhance staining. Endogenous peroxidase activity was quenched using 0.3% H2O2 for 30 min, followed by incubation for 20 min with normal serum matched to the secondary antibody to reduce background staining. Samples were incubated overnight at 4°C with the following primary antibodies all purchased from Abcam (Cambridge, MA) at specific individual concentrations: rabbit polyclonal anti-C3 antibody (ab129945, 0.2 μg/ml), rabbit polyclonal anti-CX3CL1 antibody (ab85034, 0.2 μg/ml), rabbit monoclonal anti-CLU antibody (ab92548, 0.2 μg/ml), rabbit polyclonal anti-LAPTM5 antibody (ab108014, 0.2 μg/ml), rabbit monoclonal anti-MGST1 antibody (ab131059, 1 μg/ml) and rabbit monoclonal anti-SPP1 antibody (ab91655, 0.2 μg/ml). Isotype-matched IgG (Jackson ImmunoResearch Laboratories Inc., West Grove, PA) was used as a negative control. The Vectastain Elite ABC kit and AEC substrate kit (Dako North America Inc., Carpinteria, CA) were used to visualize antibody binding and slides were counterstained with Mayer’s hematoxylin (Polysciences Inc., Warrington, PA) and mounted using faramount mounting medium (Dako North America Inc.). Images were acquired using a Nikon Microphot microscope with a Plan 40 ± N.A. 0.70 objective lens and an Olympus DP70 CCD camera.
Analysis of telomere length
Pelleted small airway epithelial cells obtained by bronchoscopy or blood derived leukocytes were processed for DNA extraction using the Qiagen Puregene kit (Germantown, MD). DNA was extracted from the small airway epithelium of 21 male nonsmokers and 22 male smokers and matched blood derived leukocytes of a subset of 12 nonsmokers and 16 smokers (Table 2) and the telomere length for each sample determined by Southern analysis using the TeloTAGGG telomere length assay kit (Roche Applied Science, Indianapolis, IN). Briefly, 1 μg of purified DNA was digested with the restriction enzymes HinfI plus RsaI (20 U each) in a 20 μl reaction volume for 2 hr, 37°C. The reaction was stopped by adding gel loading buffer. Digested DNA was separated by electrophoresis on a 0.8% agarose gel in Tris-acetate ethylenediaminetetraacetic acid (EDTA) buffer at 5 V/cm for 4 h and transferred to a nylon membrane by Southern analysis using 3 M NaCl plus 0.3 M sodium citrate buffer pH 7.0 for 16 hr. The transferred DNA was fixed by UV-crosslinking (120 mJ) to the nylon membrane. The terminal restriction fragments (TRFs) were hybridized to a digoxigenin (DIG)-labeled probe specific for telomeric repeats for 1 hr, 42°C, followed by incubation with a DIG-specific antibody covalently coupled to alkaline phosphatase for 30 min, 25°C. Finally, the immobilized telomere probe on the membrane was detected by incubation with CDP-Star, a chemiluminescent substrate, followed by exposure to an X-ray film. The film was then scanned with a densitometer and the mean TRF length was calculated using Image quant TL software (GE Healthcare, Port Washington, NY). Samples were analyzed in an un-blinded fashion; however to minimize the effect of variability between experiments, equal numbers of nonsmoker and smoker samples were run on each gel (see Additional file 1: Figure S1 for representative gels).
Comparison of demographic parameters among groups was performed by 2-tailed Students t test. For the HG-U133 Plus 2.0 gene expression data, an ANOVA was performed to assess the effect of smoking significance on gene expression corrected for all covariates, including age, gender, ethnicity and region. Differences in expression were evaluated by Students t test and a p < 0.05, using a Benjamini-Hochberg correction to limit the false discovery rate, was considered to represent a significant difference in expression between groups. The difference in mean telomere (TRF) length between phenotypes was calculated by 2-tailed Students t test and a p value <0.05 was considered to represent a significant difference in length between groups. Correlations between mean telomere length and age as well as smoking history were done using both univariate and multivariate linear regression and calculated using commercial software (Stata 10.1, StataCorp, College Station, TX). A p < 0.05 was considered to represent a statistically significant correlation.
Effect of smoking on aging-related genes in the small airway epithelium
Aging genes differentially expressed in the small airway epithelium of healthy smokers vs healthy nonsmokers 1
Direction of expression with smoking
Direction of expression with aging4
Secreted phosphoprotein 1
Lysosomal protein transmembrane 5
Macrophage expressed 1
Microsomal glutathione S-transferase 1
Collagen, type III, alpha 1
Fc fragment of IgG, low affinity IIb, receptor (CD32)
Glutathione S-transferase alpha 1
Lipopolysaccharide-induced TNF factor
Growth hormone inducible transmembrane protein
Vesicle amine transport protein 1 homolog (T. californica)
Metallothionein 1 F
Proprotein convertase subtilisin/kexin type 6
Complement component 4A
Chemokine (C-X3-C motif) ligand 1
Complement component 3
Effect of smoking on telomere length in the small airway epithelium
Previous studies have shown that telomere length is shorter in peripheral blood leukocytes of smokers compared to nonsmokers -. Cell differential analysis of the small airway epithelial brushes of nonsmokers and smokers demonstrated epithelial cells were the dominant cell population with a small amount of inflammatory cells. However, there was a significantly lower level of inflammatory cells in the brushings of smokers compared to nonsmokers (0.9% vs 2.5%, p < 0.0001; Table 2) suggesting the reduced telomere length we observe in the small epithelium of smokers does not result from the presence of inflammatory cells and are specific for the epithelial cell populations. In support of this, we measured the telomere length in the peripheral blood leukocytes of a subset of age matched nonsmokers (n = 12/21) and smokers (n = 16/22) (40.8 ± 9.7 yrs nonsmokers and 46 ± 6.4 yrs smokers, p > 0.12). The results demonstrated no significant differences (p > 0.1) in the telomere length of peripheral blood leukocytes in smokers vs nonsmokers (7.04 ± 0.71 kb vs 7.42 ± 0.54 kb; Additional file 1: Figure S3). However, in the same samples we still observe a significant decrease (p < 0.002) in the telomere length of the small airway epithelium in smokers vs nonsmokers (6.49 ± 0.74 kb vs 7.95 ± 1.26 kb). Therefore, these data suggest that reduced telomere length observed in the small epithelium of smokers is specific for the epithelial cell populations. Overall, these data suggest that smoking prematurely ages the small airway epithelium with smokers displaying on average telomere lengths ~14% shorter than chronologically age-matched healthy nonsmokers.
During the natural aging process of healthy individuals, the lung undergoes a progressive decrease in lung function ,. Studies have demonstrated that cigarette smoking accelerates the rate of lung function decline, suggesting that smoking results in premature aging of the lung ,. The small airway epithelium consists of 4 major cell types that line the airway lumen ≥6th generations and take the initial brunt of inhaled toxins from the cigarette smoke ,,. The earliest abnormalities and pathology linked to smoking and development of COPD manifest in the small airway epithelium ,,, and biologic changes in the small airway epithelium in response to cigarette smoking can be observed long before the presence of pathology ,,. The focus of this study was to determine if cigarette smoking results in accelerated aging of the small airway epithelium at the molecular level before physical manifestation of lung disease or lung function decline. To address this question we used the sensitive molecular parameters of aging-related gene expression and telomere length to compare the aging process of the small airway epithelium in age matched healthy nonsmokers and healthy smokers. The data demonstrate that smoking significantly dysregulates expression of an aging gene signature in the small airway epithelium, with ~26% of the aging genes (18/67) up or down-regulated greater than 1.2-fold in response to smoking. Over 72% (13/18) of the 18 smoking dysregulated genes demonstrate the same direction of expression observed with aging, suggesting smoking is accelerating the aging process in the small airway epithelium. In support of this concept, telomere length analysis in the small airway epithelial cells of an independent cohort of age matched male subjects revealed smoking status had a significant impact on telomere length, with smokers displaying on average a ~14% decrease in telomere length compared to nonsmokers. One of the limitations of this study is the presence of multiple cell types and epithelial subsets in the brushings obtained during bronchoscopy which makes it difficult to determine the individual cell type responsible for the accelerated aging in response to smoking. However, the epithelial brushes we analyze are processed immediately upon collection thus reducing any potential artifacts from downstream processing and thus represent a true reflection of the in vivo environment. In addition, the differences in number of epithelial cell subsets (ciliated, secretory, basal and undifferentiated columnar) between nonsmokers and smokers may allow inference towards the cell populations responsible for the smoking dependent changes we observe at the gene expression level and in telomere length and suggest that specific cell populations are more susceptible to the aging effects of cigarette smoke. Currently the biological consequences of smoking-dependent changes in age related gene expression and telomere shortening in the small airway epithelium on lung function decline and/or clinical disease in smokers are unknown. In addition, it is currently unknown if the effects of aging gene dysregulation and telomere shortening are mutually exclusive with each contributing independently to development of specific lung disease (e.g. COPD vs lung cancer) or both processes are biologically connected and contribute equally to development of disease. Future studies dissecting the biological consequences of each molecular phenomenon will help clarify their role in development of specific lung diseases. However, based on the known function of the aging signature genes and telomeres we can postulate potential consequences of these changes.
Consequence of aging-related gene expression changes in the small airway epithelium
Characterization of the biological processes associated with the original 73 gene aging signature identified by de Magalhaes et al.  highlighted enrichment for genes involved in multiple cellular processes, including immune response, inflammation, mitochondrial function, metabolism, lysosomes, collagen, apoptosis, cell cycle and markers of cellular senescence. Pathway analysis of the 18 smoking dysregulated genes failed to identify enrichment of a single biological process. However, based on the known biological functions of a subset of these genes, we can hypothesize how they may play a role in regulating the response of the small airway epithelium to cigarette smoke exposure and their potential contribution to the development of long-term pathological changes associated with smoking.
Based on the hypothesis that age-related gene expression changes may have a protective role in mediating the harmful effects of cellular dysfunction over-time, smoking-induced expression of these genes may serve a similar role. For instance, increased expression of the anti-oxidant genes MGST1 and GSTA1 in smokers most likely protects against elevated levels of harmful oxidants -. In addition, elevated expression of the protein chaperone clusterin (CLU) may curtail the effects of protein mis-folding and aggregation in response to cigarette smoke ,.
Alternatively, dysregulated expression of some genes may play a key role in driving the pathological changes associated with cigarette smoking. SPP1, the most highly induced smoking-responsive aging gene, encodes osteopontin, a secreted cytokine-like molecule that mediates immune function and extracellular matrix (ECM) remodeling . A recent study demonstrated that cigarette smoke-dependent up-regulation of SPP1 expression mediates induction of emphysema in an experimental mouse model . In addition, asbestos exposure in a murine lung fibrosis model results in up-regulation of SPP1 in bronchiolar epithelial cells . Microarray analysis comparing lungs of wild-type and SPP1 knockout mice exposed to asbestos showed that SPP1 modulated expression of a number of genes involved in cell signaling, immune system/defense, cell cycle regulation and ECM remodeling . Therefore, elevated expression of SPP1 in the small airway epithelium of smokers may have significant long term pathological consequences and may play a key role in the process of airway remodeling characterized in smoking-induced respiratory disease such as COPD.
While a majority of the smoking dysregulated aging genes (13/18) displayed the same direction of smoking-induced expression as that observed with aging, the remaining 5 genes (COL3A1, GHITM, MT1F, PCSK6 and C3) demonstrated the opposite direction of expression observed with aging. At present, the long-term biological consequence of these alternate expression patterns is unknown. However, studies for COL3A1, which encodes for collagen 3A1, have shown that expression of this gene is positively regulated by SPP1 in a murine model of asbestos-induced lung fibrosis and that over-expression of the gene is linked to fibrosis and tissue remodeling . Therefore, while COL3A1 expression generally decreases with aging over time to possibly reduce age-related tissue remodeling, smoking-induced expression may contribute significantly to the increased incidence of airway remodeling in individuals with smoking-induced respiratory disease such as COPD ,,,,. Similarly to COL3A1, the putative tumor suppressor gene MT1F may contribute to smoking-induced pathological changes. MT1F expression is typically down-regulated in a majority of human colon tumor tissues and its enforced expression can result in increased apoptosis and inhibition of cell migration, invasion and in vivo tumorigenicity . Therefore, the general increase in expression of MT1F during aging may function as a protective mechanism to prevent development of cancer, whereas smoking-dependent down-regulation in the small airway epithelium may contribute to the increased risk of developing lung cancer in smokers . However, additional in vitro studies need to be performed to demonstrate a role of MT1F in development of lung cancer in addition to following subjects over time to confirm the correlation between MT1F expression and development of lung cancer and/or COPD.
Consequence of telomere length changes in the small airway epithelium
Based on the hypothesis that COPD is a disease of accelerated aging -, a number of different studies have analyzed the effect of smoking and COPD status on telomere length in blood leukocytes -. While some have demonstrated that smokers (independent of COPD status) have shorter blood leukocyte telomeres compared to nonsmokers ,,-, others have shown that COPD status, but not smoking is associated with shortened telomeres -. Recent findings from a large study group of 46,396 individuals from the Danish general population demonstrate that blood leukocyte short telomere length is associated with decreased lung function and increased risk of COPD . Furthermore, a meta-analysis of blood leukocyte telomere length and COPD from nine independent European studies demonstrated shorter telomere length in COPD subjects compared with healthy controls . Our study is the first to look at the effect of smoking on telomere length in the small airway epithelium and demonstrates that smoking accelerates telomere shortening in this cell population critical to lung function prior to the development of lung disease and lung function decline. In addition, our analysis of telomere length in the peripheral blood leukocytes in a subset of the nonsmokers and smokers in our study suggest the small airway epithelium is more sensitive to the telomere shortening effects of cigarette smoke exposure than peripheral blood leukocytes. In contrast to studies of blood leukocyte telomere length, we observed no negative correlation of telomere length with pack-yr. In our analysis we did not correct for confounders such as age and obesity as no significant differences are observed in these parameters between our nonsmoker and smoker groups. The presence of additional confounders may explain the lack of negative association of telomere length with pack-yr; alternatively it may result from the narrow range of pack-yr of our subjects which may prevent identification of an association. In addition, from a biological perspective a small amount of cigarette smoke exposure may result in significant effects on telomere length, therefore, masking any effect between individuals with medium to high pack-yr. The long term biological consequences of smoking-dependent telomere shortening in the small airway epithelium are unknown, however based on the function of telomeres, one may postulate a number of possible outcomes.
Telomeres are specialized regions at the ends of linear chromosomes comprised of TTAGGG repeats that play a key role in maintaining chromosome stability and are linked to aging . In germ cells, telomeres are maintained by the enzyme complex telomerase; however, lack of expression of telomerase in somatic tissues results in decreased telomere length following each round of cell division . Once telomere length shortens to a specific point, a DNA damage response is elicited that is largely p53-dependent, resulting in reduced replication potential of the cell and activation of senescence with subsequent removal of the cells by apoptosis -. The numbers of senescent cells are known to increase in tissues during the natural aging process either due to an increase in their generation or decrease in their removal . In contrast to normal cells, senescent cells are known to adopt an altered phenotype highlighted by loss of proliferative activity, persistence of active metabolism, altered secretome (enriched in pro-inflammatory cytokines and matrix metalloproteinases) and resistance to apoptosis which can lead to loss of tissue homeostasis ,,. Therefore, smoking-dependent accelerated telomere shortening may increase the rate at which cells of the small airway epithelium reach cell senescence. In support of this concept, in vivo studies have shown that alveolar epithelial cells from COPD smokers have increased numbers of senescent cells relative to healthy controls ,. Furthermore, during normal aging increased expression of the smoking-induced aging gene CLU is used as a marker for cell senescence . Therefore, our gene expression data supports the hypothesis that telomere shortening may increase the rate of senescence in the small airway epithelium. In addition, studies of lung endothelial cells from COPD smokers demonstrated they reached senescence earlier during in vitro culture relative to nonsmokers and in doing so, adopted an abnormal secretory phenotype with secretion of a diverse array of cytokines known to play a role in inflammation and airway remodeling . Senescent cells within the small airway may adopt a similar function and secrete mediators that further drive smoking-induced pathology and airway remodeling.
An alternative consequence of shorter telomere length in the small airway epithelium maybe a decline in resident stem/progenitor cell function similar to what is observed in other aging tissues ,. Basal cells are the stem/progenitor population of cells in the small airways that differentiate into the other specialized epithelial cell types of the airway during normal epithelial turnover and repair -. In vitro studies have demonstrated that continuous culture of these cells to a senescent state results in diminished differentiation capacity ,. Therefore, reduced telomere length in the airway basal cells of smokers may result in impaired differentiation capacity leading to disordered cell differentiation and airway epithelial remodeling . In support of this, the first manifestation of smoking-induced pathology in vivo of the small airway epithelium is basal cell hyperplasia which is characterized by increased proliferation and expansion of the basal cell population . However, at present we have no evidence suggesting telomere length in the basal cell population is reduced in response to cigarette smoke exposure. Future studies identifying the exact cell populations of the small airway epithelium that are more sensitive to smoking-dependent telomere shortening will help further our understanding of the biological consequences of this phenomenon.
To date, the evidence clearly shows that short telomere length is associated with incidence of specific lung diseases; however, it is unknown whether telomere shortening itself drives the pathological processes in the lung or whether this phenomenon itself is a bystander and just a consequence of other processes at work. A recent study using telomerase null mice with short telomeres demonstrated that short telomeres alone was not sufficient to induce spontaneous lung disease; however, when mice with short telomeres were experimentally exposed to cigarette smoke, they had a lower threshold to develop emphysema . These findings suggest that multiple factors contribute to the development of cigarette smoke-induced airway disease and that telomere length, at least in part, plays a partial active role in “sensitizing” the cells to further damage and destruction. Therefore, similar mechanisms may be conserved in humans, and individuals with short telomeres may be of increased risk for development of smoking-induced COPD.
In summary, these data provide evidence that smoking accelerates aging of the small airway epithelium at the molecular level. Future studies analyzing the function of smoking-induced aging-related genes and telomeres in small airway epithelium biology may help identify new therapeutic targets to treat smoking-induced lung disease.
MSW conceived of the study and its design, performed research, data analysis and interpretation, and manuscript writing. BPD conceived of the study and its design, performed research, data analysis and interpretation. JS performed data analysis and interpretation. LBA performed research. TW performed research. AR performed research. LL performed research. NRH performed data analysis and interpretation. MRS coordinated collection of clinical material. AET coordinated collection of clinical material and performed statistical analysis on data. BGH coordinated collection of clinical material. RJK coordinated collection of clinical material. JGM performed statistical analysis on data. BA performed research, data analysis and interpretation. MASM performed data analysis and interpretation. RGC conceived of the study and its design, performed data analysis and interpretation, manuscript writing, and gave final approval of manuscript. All authors read and approved the final manuscript.
We thank N. Mohamed and D.N. McCarthy for help in preparing this manuscript. These studies were supported, in part, by National Institutes of Health R01HL107882, P50 HL084936, UL1 TR000457 and UL1 RR024143. Lauren Buro-Auriemma and Lindsay Leif were supported, in part, by T32HL094284. Ann Tilley was supported, in part, by K23HL103837.
- Kirkwood TB: Understanding the odd science of aging. Cell. 2005, 120: 437-447. 10.1016/j.cell.2005.01.027.PubMedView ArticleGoogle Scholar
- Lopez-Otin C, Blasco MA, Partridge L, Serrano M, Kroemer G: The hallmarks of aging. Cell. 2013, 153: 1194-1217. 10.1016/j.cell.2013.05.039.PubMedPubMed CentralView ArticleGoogle Scholar
- Kerstjens HA, Rijcken B, Schouten JP, Postma DS: Decline of FEV1 by age and smoking status: facts, figures, and fallacies. Thorax. 1997, 52: 820-827. 10.1136/thx.52.9.820.PubMedPubMed CentralView ArticleGoogle Scholar
- Janssens JP, Pache JC, Nicod LP: Physiological changes in respiratory function associated with ageing. Eur Respir J. 1999, 13: 197-205. 10.1183/09031936.99.14614549.PubMedView ArticleGoogle Scholar
- Fletcher C, Peto R: The natural history of chronic airflow obstruction. Br Med J. 1977, 1: 1645-1648. 10.1136/bmj.1.6077.1645.PubMedPubMed CentralView ArticleGoogle Scholar
- Valdes AM, Andrew T, Gardner JP, Kimura M, Oelsner E, Cherkas LF, Aviv A, Spector TD: Obesity, cigarette smoking, and telomere length in women. Lancet. 2005, 366: 662-664. 10.1016/S0140-6736(05)66630-5.PubMedView ArticleGoogle Scholar
- Morla M, Busquets X, Pons J, Sauleda J, Macnee W, Agusti AG: Telomere shortening in smokers with and without COPD. Eur Respir J. 2006, 27: 525-528. 10.1183/09031936.06.00087005.PubMedView ArticleGoogle Scholar
- Rode L, Bojesen SE, Weischer M, Vestbo J, Nordestgaard BG: Short telomere length, lung function and chronic obstructive pulmonary disease in 46,396 individuals. Thorax. 2013, 68: 429-435. 10.1136/thoraxjnl-2012-202544.PubMedView ArticleGoogle Scholar
- Gansner JM, Rosas IO: Telomeres in lung disease. Transl Res. 2013, 162: 343-352. 10.1016/j.trsl.2013.04.001.PubMedView ArticleGoogle Scholar
- Mirabello L, Huang WY, Wong JY, Chatterjee N, Reding D, Crawford ED, De Vivo I, Hayes RB, Savage SA: The association between leukocyte telomere length and cigarette smoking, dietary and physical variables, and risk of prostate cancer. Aging Cell. 2009, 8: 405-413. 10.1111/j.1474-9726.2009.00485.x.PubMedPubMed CentralView ArticleGoogle Scholar
- McGrath M, Wong JY, Michaud D, Hunter DJ, De Vivo I: Telomere length, cigarette smoking, and bladder cancer risk in men and women. Cancer Epidemiol Biomarkers Prev. 2007, 16: 815-819. 10.1158/1055-9965.EPI-06-0961.PubMedView ArticleGoogle Scholar
- O’Donnell CJ, Demissie S, Kimura M, Levy D, Gardner JP, White C, D’Agostino RB, Wolf PA, Polak J, Cupples LA, Aviv A: Leukocyte telomere length and carotid artery intimal medial thickness: the Framingham Heart Study. Arterioscler Thromb Vasc Biol. 2008, 28: 1165-1171. 10.1161/ATVBAHA.107.154849.PubMedPubMed CentralView ArticleGoogle Scholar
- Aviv A, Chen W, Gardner JP, Kimura M, Brimacombe M, Cao X, Srinivasan SR, Berenson GS: Leukocyte telomere dynamics: longitudinal findings among young adults in the Bogalusa Heart Study. Am J Epidemiol. 2009, 169: 323-329. 10.1093/aje/kwn338.PubMedPubMed CentralView ArticleGoogle Scholar
- Albrecht E, Sillanpaa E, Karrasch S, Alves AC, Codd V, Hovatta I, Buxton JL, Nelson CP, Broer L, Hagg S, Mangino M, Willemsen G, Surakka I, Ferreira MA, Amin N, Oostra BA, Backmand HM, Peltonen M, Sarna S, Rantanen T, Sipila S, Korhonen T, Madden PA, Gieger C, Jorres RA, Heinrich J, Behr J, Huber RM, Peters A, Strauch K, et al: Telomere length in circulating leukocytes is associated with lung function and disease. Eur Respir J. 2014, 43: 983-992. 10.1183/09031936.00046213.PubMedView ArticleGoogle Scholar
- Houben JM, Mercken EM, Ketelslegers HB, Bast A, Wouters EF, Hageman GJ, Schols AM: Telomere shortening in chronic obstructive pulmonary disease. Respir Med. 2009, 103: 230-236. 10.1016/j.rmed.2008.09.003.PubMedView ArticleGoogle Scholar
- Lee J, Sandford AJ, Connett JE, Yan J, Mui T, Li Y, Daley D, Anthonisen NR, Brooks-Wilson A, Man SF, Sin DD: The relationship between telomere length and mortality in chronic obstructive pulmonary disease (COPD). PLoS One. 2012, 7: e35567-10.1371/journal.pone.0035567.PubMedPubMed CentralView ArticleGoogle Scholar
- Mui TS, Man JM, McElhaney JE, Sandford AJ, Coxson HO, Birmingham CL, Li Y, Man SF, Sin DD: Telomere length and chronic obstructive pulmonary disease: evidence of accelerated aging. J Am Geriatr Soc. 2009, 57: 2372-2374. 10.1111/j.1532-5415.2009.02589.x.PubMedView ArticleGoogle Scholar
- Savale L, Chaouat A, Bastuji-Garin S, Marcos E, Boyer L, Maitre B, Sarni M, Housset B, Weitzenblum E, Matrat M, Le Corvoisier P, Rideau D, Boczkowski J, Dubois-Rande JL, Chouaid C, Adnot S: Shortened telomeres in circulating leukocytes of patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2009, 179: 566-571. 10.1164/rccm.200809-1398OC.PubMedView ArticleGoogle Scholar
- Auerbach O, Forman JB, Gere JB, Kassouny DY, Muehsam GE, Petrick TG, Smolin HJ, Stout AP: Changes in the bronchial epithelium in relation to smoking and cancer of the lung; a report of progress. N Engl J Med. 1957, 256: 97-104. 10.1056/NEJM195701172560301.PubMedView ArticleGoogle Scholar
- Knight DA, Holgate ST: The airway epithelium: structural and functional properties in health and disease. Respirology. 2003, 8: 432-446. 10.1046/j.1440-1843.2003.00493.x.PubMedView ArticleGoogle Scholar
- Tam A, Wadsworth S, Dorscheid D, Man SF, Sin DD: The airway epithelium: more than just a structural barrier. Ther Adv Respir Dis. 2011, 5: 255-273. 10.1177/1753465810396539.PubMedView ArticleGoogle Scholar
- Beasley MB: Smoking-related small airway disease-a review and update. Adv Anat Pathol. 2010, 17: 270-276. 10.1097/PAP.0b013e3181e3bf97.PubMedView ArticleGoogle Scholar
- McDonough JE, Yuan R, Suzuki M, Seyednejad N, Elliott WM, Sanchez PG, Wright AC, Gefter WB, Litzky L, Coxson HO, Pare PD, Sin DD, Pierce RA, Woods JC, McWilliams AM, Mayo JR, Lam SC, Cooper JD, Hogg JC: Small-airway obstruction and emphysema in chronic obstructive pulmonary disease. N Engl J Med. 2011, 365: 1567-1575. 10.1056/NEJMoa1106955.PubMedPubMed CentralView ArticleGoogle Scholar
- Crystal RG, Randell SH, Engelhardt JF, Voynow J, Sunday ME: Airway epithelial cells: current concepts and challenges. Proc Am Thorac Soc. 2008, 5: 772-777. 10.1513/pats.200805-041HR.PubMedView ArticleGoogle Scholar
- Harvey BG, Heguy A, Leopold PL, Carolan BJ, Ferris B, Crystal RG: Modification of gene expression of the small airway epithelium in response to cigarette smoking. J Mol Med (Berl). 2007, 85: 39-53. 10.1007/s00109-006-0103-z.View ArticleGoogle Scholar
- de Magalhaes JP, Curado J, Church GM: Meta-analysis of age-related gene expression profiles identifies common signatures of aging. Bioinformatics. 2009, 25: 875-881. 10.1093/bioinformatics/btp073.PubMedPubMed CentralView ArticleGoogle Scholar
- Hackett NR, Butler MW, Shaykhiev R, Salit J, Omberg L, Rodriguez-Flores JL, Mezey JG, Strulovici-Barel Y, Wang G, Didon L, Crystal RG: RNA-Seq quantification of the human small airway epithelium transcriptome. BMC Genomics. 2012, 13: 82-10.1186/1471-2164-13-82.PubMedPubMed CentralView ArticleGoogle Scholar
- Buro-Auriemma LJ, Salit J, Hackett NR, Walters MS, Strulovici-Barel Y, Staudt MR, Fuller J, Stevenson CS, Hilton H, Ho MW, Crystal RG: Cigarette smoking induces small airway epithelial epigenetic changes with corresponding modulation of gene expression. Hum Mol Genet. 2013, 22: 4726-4732. 10.1093/hmg/ddt326.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrett EL, Richardson DS: Sex differences in telomeres and lifespan. Aging Cell. 2011, 10: 913-921. 10.1111/j.1474-9726.2011.00741.x.PubMedView ArticleGoogle Scholar
- Tilley AE, O°Connor TP, Hackett NR, Strulovici-Barel Y, Salit J, Amoroso N, Zhou XK, Raman T, Omberg L, Clark A, Mezey JG, Crystal RG: Biologic phenotyping of the human small airway epithelial response to cigarette smoking. PLoS One. 2011, 6: e22798-10.1371/journal.pone.0022798.PubMedPubMed CentralView ArticleGoogle Scholar
- Faux SP, Tai T, Thorne D, Xu Y, Breheny D, Gaca M: The role of oxidative stress in the biological responses of lung epithelial cells to cigarette smoke. Biomarkers. 2009, 14 (Suppl 1): 90-96. 10.1080/13547500902965047.PubMedView ArticleGoogle Scholar
- Whalen R, Boyer TD: Human glutathione S-transferases. Semin Liver Dis. 1998, 18: 345-358. 10.1055/s-2007-1007169.PubMedView ArticleGoogle Scholar
- Morgenstern R, Zhang J, Johansson K: Microsomal glutathione transferase 1: mechanism and functional roles. Drug Metab Rev. 2011, 43: 300-306. 10.3109/03602532.2011.558511.PubMedView ArticleGoogle Scholar
- Wyatt A, Yerbury J, Poon S, Dabbs R, Wilson M: Chapter 6: the chaperone action of clusterin and its putative role in quality control of extracellular protein folding. Adv Cancer Res. 2009, 104: 89-114. 10.1016/S0065-230X(09)04006-8.PubMedView ArticleGoogle Scholar
- Wei J, Rahman S, Ayaub EA, Dickhout JG, Ask K: Protein misfolding and endoplasmic reticulum stress in chronic lung disease. Chest. 2013, 143: 1098-1105. 10.1378/chest.12-2133.PubMedView ArticleGoogle Scholar
- O’Regan A: The role of osteopontin in lung disease. Cytokine Growth Factor Rev. 2003, 14: 479-488. 10.1016/S1359-6101(03)00055-8.PubMedView ArticleGoogle Scholar
- Shan M, Yuan X, Song LZ, Roberts L, Zarinkamar N, Seryshev A, Zhang Y, Hilsenbeck S, Chang SH, Dong C, Corry DB, Kheradmand F: Cigarette smoke induction of osteopontin (SPP1) mediates T(H)17 inflammation in human and experimental emphysema. Sci Transl Med. 2012, 4: 117-10.1126/scitranslmed.3003041. ra9Google Scholar
- Sabo-Attwood T, Ramos-Nino ME, Eugenia-Ariza M, Macpherson MB, Butnor KJ, Vacek PC, McGee SP, Clark JC, Steele C, Mossman BT: Osteopontin modulates inflammation, mucin production, and gene expression signatures after inhalation of asbestos in a murine model of fibrosis. Am J Pathol. 2011, 178: 1975-1985. 10.1016/j.ajpath.2011.01.048.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim V, Kelemen SE, Buel-Haija M, Gaughan JP, Sharafkaneh A, Evans CM, Dickey BF, Solomides CC, Rogers TJ, Criner GJ: Small airway mucous metaplasia and inflammation in chronic obstructive pulmonary disease. COPD. 2008, 5: 329-338. 10.1080/15412550802522445.PubMedView ArticleGoogle Scholar
- Tam A, Sin DD: Pathobiologic mechanisms of chronic obstructive pulmonary disease. Med Clin North Am. 2012, 96: 681-698. 10.1016/j.mcna.2012.04.012.PubMedView ArticleGoogle Scholar
- Yan DW, Fan JW, Yu ZH, Li MX, Wen YG, Li DW, Zhou CZ, Wang XL, Wang Q, Tang HM, Peng ZH: Downregulation of metallothionein 1 F, a putative oncosuppressor, by loss of heterozygosity in colon cancer tissue. Biochim Biophys Acta. 1822, 2012: 918-926.Google Scholar
- Punturieri A, Szabo E, Croxton TL, Shapiro SD, Dubinett SM: Lung cancer and chronic obstructive pulmonary disease: needs and opportunities for integrated research. J Natl Cancer Inst. 2009, 101: 554-559. 10.1093/jnci/djp023.PubMedPubMed CentralView ArticleGoogle Scholar
- Ito K, Barnes PJ: COPD as a disease of accelerated lung aging. Chest. 2009, 135: 173-180. 10.1378/chest.08-1419.PubMedView ArticleGoogle Scholar
- Lee J, Sandford A, Man P, Sin DD: Is the aging process accelerated in chronic obstructive pulmonary disease?. Curr Opin Pulm Med. 2011, 17: 90-97. 10.1097/MCP.0b013e328341cead.PubMedView ArticleGoogle Scholar
- Faner R, Rojas M, Macnee W, Agusti A: Abnormal lung aging in chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2012, 186: 306-313. 10.1164/rccm.201202-0282PP.PubMedView ArticleGoogle Scholar
- Chilosi M, Carloni A, Rossi A, Poletti V: Premature lung aging and cellular senescence in the pathogenesis of idiopathic pulmonary fibrosis and COPD/emphysema. Transl Res. 2013, 162: 156-173. 10.1016/j.trsl.2013.06.004.PubMedView ArticleGoogle Scholar
- Armanios M, Blackburn EH: The telomere syndromes. Nat Rev Genet. 2012, 13: 693-704. 10.1038/nrg3246.PubMedPubMed CentralView ArticleGoogle Scholar
- Jacobs JJ, de Lange T: p16INK4a as a second effector of the telomere damage pathway. Cell Cycle. 2005, 4: 1364-1368. 10.4161/cc.4.10.2104.PubMedView ArticleGoogle Scholar
- de Jesus BB, Blasco MA: Assessing cell and organ senescence biomarkers. Circ Res. 2012, 111: 97-109. 10.1161/CIRCRESAHA.111.247866.View ArticleGoogle Scholar
- Sahin E, DePinho RA: Axis of ageing: telomeres, p53 and mitochondria. Nat Rev Mol Cell Biol. 2012, 13: 397-404. 10.1038/nrm3352.PubMedPubMed CentralView ArticleGoogle Scholar
- Naylor RM, Baker DJ, van Deursen JM: Senescent cells: a novel therapeutic target for aging and age-related diseases. Clin Pharmacol Ther. 2013, 93: 105-116. 10.1038/clpt.2012.193.PubMedPubMed CentralView ArticleGoogle Scholar
- Campisi J, di d’Adda Fagagna F: Cellular senescence: when bad things happen to good cells. Nat Rev Mol Cell Biol. 2007, 8: 729-740. 10.1038/nrm2233.PubMedView ArticleGoogle Scholar
- Tsuji T, Aoshiba K, Nagai A: Alveolar cell senescence in patients with pulmonary emphysema. Am J Respir Crit Care Med. 2006, 174: 886-893. 10.1164/rccm.200509-1374OC.PubMedView ArticleGoogle Scholar
- Aoshiba K, Nagai A: Senescence hypothesis for the pathogenetic mechanism of chronic obstructive pulmonary disease. Proc Am Thorac Soc. 2009, 6: 596-601. 10.1513/pats.200904-017RM.PubMedView ArticleGoogle Scholar
- Trougakos IP, Gonos ES: Regulation of clusterin/apolipoprotein J, a functional homologue to the small heat shock proteins, by oxidative stress in ageing and age-related diseases. Free Radic Res. 2006, 40: 1324-1334. 10.1080/10715760600902310.PubMedView ArticleGoogle Scholar
- Amsellem V, Gary-Bobo G, Marcos E, Maitre B, Chaar V, Validire P, Stern JB, Noureddine H, Sapin E, Rideau D, Hue S, Le Corvoisier P, Le Gouvello S, Dubois-Rande JL, Boczkowski J, Adnot S: Telomere dysfunction causes sustained inflammation in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2011, 184: 1358-1366. 10.1164/rccm.201105-0802OC.PubMedView ArticleGoogle Scholar
- Hackett NR, Shaykhiev R, Walters MS, Wang R, Zwick RK, Ferris B, Witover B, Salit J, Crystal RG: The human airway epithelial basal cell transcriptome. PLoS One. 2011, 6: e18378-10.1371/journal.pone.0018378.PubMedPubMed CentralView ArticleGoogle Scholar
- Rock JR, Randell SH, Hogan BL: Airway basal stem cells: a perspective on their roles in epithelial homeostasis and remodeling. Dis Model Mech. 2010, 3: 545-556. 10.1242/dmm.006031.PubMedPubMed CentralView ArticleGoogle Scholar
- Xian W, McKeon F: Adult stem cells underlying lung regeneration. Cell Cycle. 2012, 11: 887-894. 10.4161/cc.11.5.19328.PubMedPubMed CentralView ArticleGoogle Scholar
- Gray TE, Guzman K, Davis CW, Abdullah LH, Nettesheim P: Mucociliary differentiation of serially passaged normal human tracheobronchial epithelial cells. Am J Respir Cell Mol Biol. 1996, 14: 104-112. 10.1165/ajrcmb.14.1.8534481.PubMedView ArticleGoogle Scholar
- Widdicombe JH, Sachs LA, Morrow JL, Finkbeiner WE: Expansion of cultures of human tracheal epithelium with maintenance of differentiated structure and function. Biotechniques. 2005, 39: 249-255. 10.2144/05392RR02.PubMedView ArticleGoogle Scholar
- Alder JK, Guo N, Kembou F, Parry EM, Anderson CJ, Gorgy AI, Walsh MF, Sussan T, Biswal S, Mitzner W, Tuder RM, Armanios M: Telomere length is a determinant of emphysema susceptibility. Am J Respir Crit Care Med. 2011, 184: 904-912. 10.1164/rccm.201103-0520OC.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.