Abstract
Only a fraction of all smokers develop chronic obstructive pulmonary disease (COPD), suggesting a large role for genetic susceptibility. The leptin receptor (LEPR) is present in human lung tissue and may play a role in COPD pathogenesis. The present study examined the association between genetic variants in the LEPR gene and lung function decline in COPD.
In total, 429 European Americans were randomly selected from the National Heart Lung and Blood Institute Lung Health Study. 36 single nucleotide polymorphisms (SNPs) in LEPR were genotyped using the IlluminaTM GoldenGate platform (Broad Institute, Cambridge, MA, USA). Mean annual decline in forced expiratory volume in 1 s % predicted over the 5-yr period was calculated using linear regression. Linear regression models were also used to adjust for potential confounders. In addition, in vivo expression of the receptor gene was assessed with immunohistochemistry on lungs from smoke-exposed inbred mice.
We identified significant associations (p<0.05) between lung function decline and 21 SNPs. Haplotype analyses confirmed several of these associations seen with individual markers. Immunohistochemistry results in inbred mice strains support a potential role of LEPR in COPD pathogenesis.
We identified genetic variants in the LEPR gene significantly associated with lung function decline in a population of smokers with COPD. Our results support a role for LEPR as a novel candidate gene for COPD.
Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of death in the USA and the fifth leading cause of death worldwide, and its prevalence is expected to increase in coming decades 1, 2. The overwhelming majority of COPD is caused by environmental exposures. In the USA, this exposure is primarily cigarette smoke (CS); however, only 15% of all smokers develop COPD. This suggests a large role for genetic susceptibility.
In addition to its role in obesity (appetite suppression), leptin has been shown to have multiple other functions, including increasing sympathetic nerve activity, maintaining reproductive function, immunity, angiogenesis, preserving normal respiratory function in the presence of obesity and cell proliferation of tracheal epithelial cells and lung growth 3. After adjustment for obesity, the leptin pathway has been associated with inflammatory markers as well as multiple inflammatory conditions, including cardiovascular disease 4, 5. Its potential role in the systemic inflammatory response in patients with COPD is evident from the correlation of leptin with other inflammatory markers 6–8. Leptin has been shown to directly stimulate phagocytic activity of macrophages and enhance endotoxin-induced production of tumour necrosis factor-α, interleukin (IL)-6, and IL-12, cytokines typically involved in COPD pathogenesis 9–12. Conversely, deficiency in leptin or its receptor may predispose sufferers to both immunodeficiency and infection 10. Serum leptin levels are elevated in exacerbations of COPD and have also been associated with bacterial pulmonary infections, which are associated with COPD exacerbations and decline in lung function 7, 13–15. Leptin exerts multiple effects through its leptin receptor (LEPR), located on human chromosome 1p31. The LEPR is produced in several alternatively spliced forms that share extracellular and transmembrane domains but have varying cytoplasmic residues 16, and has a wide tissue distribution, including lung tissue 17.
We tested for association between genetic variants in the LEPR gene and lung function decline in a subset of the multicentre National Heart, Lung and Blood Institute (NHLBI)-supported Lung Health Study (LHS) cohort. In addition, the potential role of LEPR in COPD pathogenesis was further evaluated in smoke-exposed AKR/J mice, which display marked airspace enlargement by histological and morphometric criteria 18. To determine whether the LEPR is differentially localised to resident lung cells under conditions that simulate human COPD, we performed immunohistochemistry on murine lung sections obtained from AKR/J mice exposed to 4 months of CS versus room air controls.
Identifying pathways and novel molecular targets that modify the clinical course of disease is fundamental to developing preventive strategies and novel therapies.
METHODS
Subjects
Subjects participating in the current study included 429 European Americans (EA) randomly selected from a group of 4,287 participants of the multi-centre NHLBI-supported LHS for whom DNA was available. Lung function was measured annually over 5 yrs and conducted according to American Thoracic Society guidelines using identical spirometers, software, procedures and reading centre personnel 19, 20. The quality of the spirometry testing conducted by the technicians was monitored centrally throughout the testing and comparison of baseline spirometry measures showed good reproducibility with very small mean short-term intra-individual variations in forced expiratory volume in 1 s (FEV1) 21. Lung function data from annual visit 1 to annual visit 5 was used for the current analyses and has been shown to have a good linear fit in previous LHS analyses 19, 22. Subjects with less than three annual lung function measurements were excluded from the analysis (n = 15).
Single nucleotide polymorphism selection and genotyping
Single nucleotide polymorphisms (SNPs) representing the LEPR gene were selected from Goldenpath (genome.ucsc.edu) and/or National Center for Biotechnology Information (NCBI; www.ncbi.nlm.nih.gov; Bethesda, MD, USA). Criteria included: 1) SNPs approximating inter-SNP distance as close to 5 kb as were available at the time of the dbSNP Build 124 (NCBI); 2) representation of SNPs in the promoter, coding and 3′-untranslated region; and 3) SNPs with acceptable design scores according to the Illumina Assay Design Tool for genotyping on the IlluminaTM GoldenGate platform (Broad Institute, Cambridge, MA, USA). Priority for selecting SNPs included: 1) regulatory and coding SNPs; 2) highly polymorphic SNPs, preferably ≥10% minor allele frequency; 3) validated SNPs; and 4) SNPs at intron/exon boundaries. A total of 36 LEPR SNPs spanning 228,294 bp on human chromosome 1p31 with an average inter-SNP distance of 6.53 kb (range 1.9–14 kb) are summarised in table 1⇓.
Statistical analysis
Mean annual decline in lung function (post-bronchodilator FEV1 % predicted) was calculated as a linear regression slope over the 5-yr study period. Linear regression models were used to adjust for potential confounders, including baseline characteristics, smoking history (pack-yrs), age, sex, FEV1 % pred and airway reactivity (AR). AR was calculated as a quantitative measure, using the two-point slope 22. In addition to adjusting for baseline characteristics, 1) change in body mass index (BMI), calculated as a linear regression slope and 2) smoking status at year 5, defined as “continuous smoker”, “intermittent smoker”, and “sustained quitter”, were also analysed.
Forward and backward selection was used to develop a parsimonious model. Residuals from the regression were included in a genetic additive model and the most common homozygote genotype for each SNP served as the reference category. The interaction term for smoking and genetic effect were tested using PLINK (Boston, MA, USA) 23, while testing for quantitative interaction. Each SNP locus was evaluated for Hardy–Weinberg equilibrium. All analyses were performed with StataSE, version 8.0 (Stata Corp, College Station, TX, USA) and PLINK 23.
Individual haplotypes were analysed using PLINK. Pairwise linkage disequilibrium (LD) based on the D′ statistic was measured using Haploview (Broad Institute) 24. LD blocks were defined using their default algorithm 25. Sliding windows of two to four adjacent SNPs were used to test for association. Haplotype estimates were computed using PHASE (Seattle, WA, USA) assuming no recombination 26.
Murine model
10-week-old AKR/J mice were exposed to CS 5 h·day−1, 5 days·week−1, for 6 months. The exposure was conducted by burning 2R4F reference cigarettes (University of Kentucky, Tobacco Research Institute, Lexington, KY, USA), using a TE-10 smoking machine (Teague Enterprises, Woodland, CA, USA). Each cigarette was puffed for 2 s, once every minute for a total of eight puffs, at a flow rate of 1.05 L·min−1, to provide a standard puff of 35 cm3. The smoke machine was adjusted to produce a mixture of sidestream smoke (89%) and mainstream smoke (11%). The smoke chamber was monitored daily for total suspended particles and carbon monoxide, with concentrations of 90 mg·m−3 and 350 ppm, respectively. Air-exposed control mice were housed in a filtered air environment.
Immunohistochemistry
5-μm paraformaldehyde-fixed, paraffin-embedded tissue sections from smoke-exposed and room air control AKR/J mice were deparaffinised and rehydrated in an ethanol series. Sections were blocked for nonspecific binding with 3% normal serum from chicken and incubated with the primary antibodies for 1 h at room temperature. Following incubation with the primary antibody overnight at 4°C, slides were washed with PBS Tween-20, incubated for 30 min at room temperature with an appropriate biotinylated secondary antibody and developed by using 3,3′-diaminobenzidine substrate and chromagen from Dako (Stockport, UK). Antigen retrieval was performed using citrate buffer for 30 min. Antibodies were used at the following concentrations: LEPR, ObR (goat polyclonal, I-17, 1:100; Santa Cruz Biotechnology, Santa Cruz, AZ, USA), donkey anti-goat secondary antibody (sc2020, 1:400, Santa Cruz Biotechnology).
RESULTS
Clinical characteristics of the subjects are presented in table 2⇓. There were no statistical differences in baseline characteristics between those included in the final analyses (n = 414) compared with those excluded because more than two data points in lung function were missing (n = 15, data not shown). Distributions of baseline characteristics and lung function measurements were also similar to those in the full LHS cohort from which this subset was selected (n = 5,887, data not shown). Age, baseline lung function, AR, smoking status at year 5, and change in BMI were independently associated with lung function decline, and were included in the final regression analyses. There was no significant interaction between smoking status at year 5 and LEPR polymorphisms.
Single-marker analyses
All 36 LEPR SNPs were in the Hardy–Weinberg equilibrium. All results of the two-point tests for association between LEPR markers and FEV1 % pred decline are presented in table 3⇓. We found evidence for significant associations between 21 SNPs in the LEPR gene and FEV1 % pred decline, spanning the length of the gene from intron 2 to the 5′ end. The minor alleles of most significant SNPs were associated with attenuation in lung function decline. Of particular interest, each G allele at the functional marker rs1137100, in exon 4, which creates a lysine to arginine amino acid change, was associated with a 0.33%·yr−1 attenuation in annual loss of FEV1 % pred.
Haplotype analyses
Considering the clustering of most genotyped SNPs within six LD blocks with high levels of LD also observed between SNPs across different blocks, a systematic sliding window approach was implemented, considering windows of two to four SNPs/window beginning with the first (5′) marker, and working across the gene, one marker at a time. Haplotype tests revealed multiple association signals in three specific regions that overlapped with the single-SNP results described above. 25 haplotypes showed stronger association than single SNP results (fig. 2⇓). One of the most compelling regions was at the 5′ end of LEPR (region 1) and the haplotype GCCT (rs7531867, rs1805096, rs1892535 and rs6691346), which revealed a strong association with lung function decline (p = 0.003) and includes the coding-synonymous SNP, rs1805096. Furthermore, one 3-SNP and two 4-SNP window haplotypes in region 2 spanning rs10443259 to rs10889562 (introns 2 through 5) were most strongly associated with lung function decline (p = 0.002), and included the functional marker rs1137100, as described previously.
Validation of LEPR as a candidate gene for COPD phenotypes in inbred mice
AKR/J mice, a well-characterised inbred strain, exposed to 4 months of CS not only exhibit airspace enlargement but also develop airway thickening and inflammation that is highly reminiscent of COPD 18. We examined whether CS exposure altered LEPR expression in the lungs of this strain. Using immunohistochemistry, we observed that LEPR expression was evident in macrophages, the airways and airspace compartments in room air-exposed AKR/J mice. However, upon smoke exposure, there was a reduction of staining in the airspace and airway wall (fig. 3⇓). These findings are consistent with CS-induced downregulation of LEPR expression in the epithelial compartments of CS-exposed mice.
DISCUSSION
We identified 21 SNPs in the LEPR gene that were significantly associated with lung function decline in an EA population with COPD and our haplotype analysis supported results from SNP analysis. The association signals observed across these SNPs likely represent two or three signals in the gene. Furthermore, the potential role of the LEPR in lung architecture and COPD phenotypes is supported by an AKR/J murine model showing decreased LEPR expression in airway wall and epithelium after smoke exposure. Our results identify LEPR as a novel candidate gene for COPD.
To our knowledge, no previous studies have investigated the role of genetic polymorphisms in the LEPR gene and COPD or its associated phenotypes. The phenotype that we studied was rate of decline of lung function rather than COPD per se. The rationale for this was that rate of decline in lung function is a more precise phenotype than a single cross-sectional measurement of lung function needed to define COPD, and may serve as a more genetically homogeneous phenotype. We identified 21 SNPs in the LEPR gene that were significantly associated with lung function decline in a EA population with COPD. One of the most compelling regions was at the 5′ end of LEPR (region 1), which includes part of the intracellular domain of the receptor 16 and may result in different signalling potency. Most variants identified in LEPR were associated with an attenuation of lung function decline, and notably, the G variant at the nonsynonymous SNP rs1137100, which codes for an amino acid substitution in the extracellular domain of the LEPR 27, was associated with a 0.33%·yr−1 attenuation in loss of FEV1 % pred. Airflow obstruction that occurs in COPD is caused by a mixture of small airway disease, parenchymal destruction and increased airways responsiveness that develops over decades 28. Thus, over time, these genetic variants may have a substantial impact on disease progression. Over a 40-yr period, a person homozygous for the G allele at rs1137100 will have a FEV1 % pred that is 26.4% higher compared with a homozygote for the major allele. Given a minor allele frequency of 0.26, as seen in our cohort, we might extrapolate 7% of the population to be homozygote for the minor allele (GG), 38% to be heterozygote (AG) and 55% to be homozygote for the major allele (AA). Even after adjusting for smoking status, this genetic variant could explain the difference between mild and moderate or moderate and very severe COPD under current Global Initiative for Chronic Obstructive Lung Disease (GOLD) criteria 29.
The dearth of well-characterised populations of COPD subjects with longitudinal lung function data comparable to the LHS adds to the uniqueness of our population, but limits the opportunity to replicate these results. However, finding multiple significant SNPs in a single gene reduces the likelihood of spurious results due to multiple testing, though it may also be due in part to the high degree of LD in this gene. Furthermore, when using a false discovery rate of 0.2, six SNPs would remain statistically significantly associated with lung function decline (p<0.006) in the present study. Previous studies have shown genetic variants in LEPR to be associated with markers of inflammation, including C-reactive protein (CRP) and fibrinogen levels, lending support to the hypothesis that the leptin pathway has a physiological influence on inflammatory traits. Specifically, the minor allele (T) for the rs1805096 locus was associated with lower levels of fibrinogen, CRP and IL-6 levels in a previous study of healthy EAs 30. Both lower CRP and IL-6 levels have been associated with attenuated lung function decline in the LHS cohort 31. Interestingly, the minor allele (T) in the coding SNP rs1805096 was also associated with lower rates of lung function decline in our study (0.295 %·yr−1; p = 0.007). Therefore, we hypothesise one potential mechanism of the disease-modifying effects of the LEPR gene may be mediated through inflammatory mechanisms. The notion that leptin may function as an immunomodulatory cytokine has become increasingly accepted. The association between leptin and muscle wasting and cachexia in COPD is well established 32; and the presence of leptin in induced sputum of patients with moderate COPD and its association with other inflammatory markers 33 suggests leptin may be involved in the local inflammatory response in COPD.
Inbred mice or guinea pigs subjected to chronic CS exposure have been shown to be invaluable models of CS-induced parenchymal lung disease 34, 35. Furthermore, the use of animal models of complex human diseases to parse candidate genes identified in broad genetic or genomic surveys is a standard approach to initial pathway validation, especially if additional populations are not available for replication and therapeutic targets are an ostensible goal 36, 37. We noted a reduction in staining in the airspace wall and airway epithelial compartment, but retained expression in macrophages in AKR/J smoke-exposed mice. This is consistent with findings of Bruno et al. 13 who found decreased expression of leptin and its receptor in smokers and subjects with mild to severe COPD as compared with healthy nonsmoking subjects. Accordingly, reduced leptin receptor expression (acquired or genetic) may be a critical predisposing factor to CS-induced lung disease.
In summary, the current study is the first to report an association between LEPR polymorphisms and COPD. Our results support the role of the leptin pathway, and particularly the LEPR in COPD and lung function decline. Specifically, individuals with the minor allele at the selected polymorphisms were less susceptible to loss of lung function and COPD progression. The SNP and haplotype tests point to three clusters of signal highlighting two potential loci considering the linkage disequilibrium between these signals: the functional SNP rs137100 and at the synonymous coding SNP rs1805096. While these two SNPs are not in LD with each other (D′ = 0.45, R2 = 0.125), haplotypes of significance that include each SNP overlap neighbouring LD blocks, making it somewhat difficult to tease apart these two signals. It is necessary to replicate our findings in other populations; however, our results identify LEPR as a novel target in COPD and lung function loss.
Support statement
This study was supported by the National Heart Lung and Blood Institute (grants HL076322, HL066583 and HL010342; Bethesda, MD, USA) and the National Institute on Ageing (AG21057; Bethesda). This study was also supported in part by the Intramural Research Program of the National Human Genome Research Institute and National Institutes of Health (Bethesda, MD, USA). Genotyping services were provided by the Johns Hopkins University (Baltimore, MD, USA) under US Federal Government contract number N01-HV-48195 from the NHLBI. K.C. Barnes was supported in part by the Mary Beryl Patch Turnbull Scholar Program (Baltimore, MD, USA).
Statement of interest
None declared.
Acknowledgments
We would like to thank: P. Chi, M. Campbell, A.L. Scott, A.E. Jedlicka and M.V. Mintz (Malaria Research Institute, Gene Array Core Facility, Johns Hopkins University, Baltimore, MD, USA); A.F. Scott, K. Doheny, R. Ashworth and C. Boehm (Center for Inherited Disease Research Institute, Johns Hopkins University); and H. Voelker and K. Farnell (LHS Data Coordinating Center, University of Minnesota, Minneapolis, MN, USA).
The principal investigators and senior staff of the clinical and coordinating centres, the National Heart Lung and Blood Institute (Bethesda, MD, USA) are as follows. M.D. Altose (Principal Investigator) and C.D. Deitz (Project Coordinator): Case Western Reserve University, Cleveland, OH, USA. M.S. Eichenhorn (Principal Investigator), K.J. Braden (Project Coordinator) and R.L. Jentons (Project Coordinator): Henry Ford Hospital, Detroit, MI, USA. R.A. Wise (Principal Investigator), C.S. Rand (Co-Principal Investigator) and K.A. Schiller (Project Coordinator): Johns Hopkins University School of Medicine, Baltimore, MD, USA. P.D. Scanlon (Principal Investigator), G.M. Caron (Project Coordinator), K.S. Mieras and L.C. Walters: Mayo Clinic, Rochester, MN, USA. A.S. Buist (Principal Investigator), L.R. Johnson (LHS Pulmonary Function Coordinator) and V.J. Bortz (Project Coordinator): Oregon Health Sciences University, Portland, OR, USA. W.C. Bailey (Principal Investigator) and L.B. Gerald (Project Coordinator): University of Alabama at Birmingham, Birmingham, AL, USA. D.P. Tashkin (Principal Investigator) and I.P. Zuniga (Project Coordinator): University of California, Los Angeles, CA, USA. N.R. Anthonisen (Principal Investigator, Steering Committee Chair), J. Manfreda (Co-Principal Investigator), R.P. Murray (Co-Principal Investigator), S.C. Rempel-Rossum (Project Coordinator): University of Manitoba, Winnipeg, MB, Canada. J.E. Connett (Principal Investigator), P.L. Enright, P.G. Lindgren, P. O'Hara (LHS Intervention Coordinator), M.A. Skeans and H.T. Voelker: University of Minnesota Coordinating Center, Minneapolis, MN, USA. R.M. Rogers (Principal Investigator) and M.E. Pusateri (Project Coordinator): University of Pittsburgh, Pittsburgh, PA, USA. R.E. Kanner (Principal Investigator) and G.M. Villegas (Project Coordinator): University of Utah, Salt Lake City, UT, USA. S.S. Hurd (Former Director, Division of Lung Diseases), J.P. Kiley (Former Project Officer and Director, Division of Lung Diseases), G. Weinmann (Former Project Officer and Director, Airway Biology and Disease Program, DLD), M.C. Wu (Division of Epidemiology and Clinical Applications): National Heart, Lung, and Blood Institute staff, Bethesda, MD, USA.
- Received August 5, 2008.
- Accepted December 19, 2008.
- © ERS Journals Ltd