Genetic variants associated with lung function: the long life family study
© Thyagarajan et al.; licensee BioMed Central Ltd. 2014
Received: 8 July 2014
Accepted: 16 October 2014
Published: 1 November 2014
Reduced forced expiratory volume in 1 second (FEV1) and the ratio of FEV1 to forced vital capacity (FVC) are strong predictors of mortality and lung function is higher among individuals with exceptional longevity. However, genetic factors associated with lung function in individuals with exceptional longevity have not been identified.
We conducted a genome wide association study (GWAS) to identify novel genetic variants associated with lung function in the Long Life Family Study (LLFS) (n = 3,899). Replication was performed using data from the CHARGE/SpiroMeta consortia. The association between SNPs and FEV1 and FEV1/FVC was analyzed using a linear mixed effects model adjusted for age, age2, sex, height, field center, ancestry principal components and kinship structure to adjust for family relationships separately for ever smokers and never smokers. In the linkage analysis, we used the residuals of the FEV1 and FEV1/FVC, adjusted for age, sex, height, ancestry principal components (PCs), smoking status, pack-years, and field center.
We identified nine SNPs in strong linkage disequilibrium in the CYP2U1 gene to be associated with FEV1 and a novel SNP (rs889574) associated with FEV1/FVC, none of which were replicated in the CHARGE/SpiroMeta consortia. Using linkage analysis, we identified a novel linkage peak in chromosome 2 at 219 cM for FEV1/FVC (LOD: 3.29) and confirmed a previously reported linkage peak in chromosome 6 at 28 cM (LOD: 3.33) for FEV1.
Future studies need to identify the rare genetic variants underlying the linkage peak in chromosome 6 for FEV1.
Pulmonary function, as determined by spirometry, indicates the physiological state of the airways and lung. Pulmonary function measures such as forced expiratory volume in one second (FEV1), forced vital capacity (FVC) and the ratio of the two measures (FEV1/FVC) are used to diagnose chronic obstructive pulmonary disease (COPD) and assist in the diagnosis of asthma, which are major causes of death and morbidity worldwide . Reduced (FEV1) and forced vital capacity (FVC), in healthy asymptomatic adults, are predictors of cardiovascular disease mortality and all cause mortality, independent of smoking history -. The Danish 1905 cohort study also showed that higher FEV1 and peak expiratory flow (PEF) in nonagenarians to be associated with lower mortality . Due to the strong and consistent association with mortality, pulmonary function has been viewed as a biomarker of aging itself . We have previously shown that the prevalence of self-reported COPD was 3 times lower in the Long Life Family Study (LLFS) as compared to the other similarly aged cohorts  suggesting that LLFS participants, who were selected for exceptional familial longevity, may have better lung function as compared to the general population. These findings are supported by a previous study that reported male offspring of long lived parents had higher lung function as compared to those with short lived parents . Several studies have shown that pulmonary function measures are heritable characteristics with estimates ranging from 38% for FEV1 to 37% for FEV1/FVC -. However, the 26 genetic loci associated with FEV1/FVC, FEV1 or both (23 loci associated with FEV1/FVC and 10 loci associated with FEV1) - in genome wide association studies explain only around 3% of the variance in FEV1/FVC and around 1.5% of the variance in FEV1. Furthermore, the mean age of Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE)/SpiroMeta consortia was 55 years; hence, it identified genes for lung function among older adults. Since findings from the LLFS  and a previous study  suggest that there may be genetic determinants of lung function among exceptionally long lived individuals and families, we conducted a genome wide association and linkage study among participants of the LLFS, a family based cohort of exceptional longevity, to identify novel genetic determinants of lung function in this unique sample.
Material and methods
Cohort description and study design
The LLFS study design has also been described in detail previously . Briefly, the LLFS is a family-based cohort study (n = 4,559) that enrolled long-lived probands and their siblings (n = 1,445), their offspring (n = 2,329) and spousal controls (n = 785) recruited from 3 U.S. field centers (Boston University Medical Center in Boston MA, Columbia College of Physicians and Surgeons in New York City NY, and the University of Pittsburgh in Pittsburgh PA) and the University of Southern Denmark to identify genetic determinants of longevity in these families. At the U.S. field centers, an initial recruitment brochure was mailed to all people in the Center for Medicare and Medicaid Services list of Medicare enrollees who were ≥89 years old on January 1, 2005, were not in end stage renal disease or hospice programs and lived in zip codes within 3 hours driving distance of one of the three U.S field centers. Mailings were conducted in collaboration with CMS and the NIA via an Intra-Agency Agreement. Study participants were also recruited from the local communities using mailed brochures, posters, web-based media and newspaper advertisements as well as community presentations at churches and senior centers. Additional mailing lists were obtained through local government agencies or purchased public domain lists from commercial vendors. The University of Southern Denmark used the Danish National Register of Persons to identify individuals who were ≥90 years during the study recruitment period without any restrictions on residence . Only families who had the proband, at least one living sibling, and one of their living offspring (minimum family size of 3) with a Family Longevity Selection Score (FLoSS) of 7 or higher that correlates well with later-observed longevity  and gave informed consent and were willing to participate in the interview and examination including donating a blood sample were eligible to participate in this study. This strategy led to the enrollment of families with the greatest potential utility for phenotypic and genetic studies of exceptional survival in families. All research was performed in compliance with the Helsinki Declaration. Written informed consent was obtained from all enrollees. In a few cases of cognitive impairment, family members were enrolled via proxy consent, provided that the participant was able to express assent at the time of the examination. This study was approved by the Institutional Review Boards at all the institutions in the United States of America and Denmark.
After excluding 15% of the participants due to presence of non-European ancestry (n = 6), low quality spirometry (defined as 2 or more acceptable spirometry maneuvers with reproducibility within 250 mL) (n = 295), self-reported pulmonary fibrosis (n = 11) obtained during an in-person interview, history of lung volume reduction surgery (n = 14), or missing genotypes (n = 344), a total of 3,889 participants were included in the present analysis.
Lung function measurements
The examinations were conducted in the home setting with portable equipment by centrally trained and certified research assistants using a standardized protocol. Lung function was measured with a portable spirometer (EasyOne Diagnostic, NDD Medical Technologies, Andover, MA) following American Thoracic Society guidelines . Calibration checks and the best 3 maneuvers were reviewed centrally by one investigator. Only spirometry tests with 2 or more acceptable maneuvers with reproducibility within 250 mL were selected for further data analyses.
Genotyping and imputation
The Human Omni chip 2.5 v1 (Illumina Inc., Ca), was used to genotype all the LLFS participants at the Center for Inherited Disease Research (CIDR). Ancestry principal components (PCs), to control for population structure, were produced with EIGENSTRAT  on 1,515 LLFS unrelated individuals using 120,093 tag SNPs, where in advance any SNPs with minor allele frequencies (MAF) <5%, Hardy Weinberg Equilibrium (HWE) p <1E-06, and with missing genotypes were excluded. Ancestry PCs produced from unrelated subjects were expanded, within EIGENSTRAT framework, to all members of LLFS. Genotype imputations were performed based on the cosmopolitan phased haplotypes of 1000 Human Genome (1000HG, version 2010–11 data freeze, 2012-03-04 haplotypes) using MACH and MINIMACH , and a total of 38,045,518 SNPs were imputed. When MAF ≥ 0.05 and r2 > 0.3 for imputed SNP filters were applied to the hybrid dataset for analysis, the number of SNPs for analysis is reduced to 6,522,421 (from a total of 38,245,546 SNPs), of which 1,204,935 SNPs were genotyped and 5,317,486 SNPs were imputed.
The statistical models used to test the association between the GWAS SNPs and lung function (FEV1 and FEV1/FVC) were identical to the models used by the CHARGE/SpiroMeta consortia , except that the LLFS study also included adjustment for kinship structure to facilitate replication of results in the CHARGE/SpiroMeta consortia. We employed a linear mixed effects model which adjusted for age, age2, sex, height, field center and ancestry PCs (PC1-20) in addition to the kinship matrix. The adjusted phenotypic residuals from these models (FEV1 and FEV1/FVC) were inverse normal transformed to normally distributed z-scores. These transformed residuals were then used as the phenotype for association testing under an additive genetic model, separately for ever smokers and never smokers. The associations between individual SNPs the FEV1 (milliliters) and FEV1/FVC (percent) were analyzed using a linear mixed effects model with kinship structure , to adjust for family relationships separately for ever smokers and never smokers. The effect estimates and standard errors for ever-smokers and never smokers were meta-analyzed using inverse-variance weighting.
We used the same criteria as the CHARGE/SpiroMeta consortia for assessing genetic associations . All SNPs that showed a borderline association (p < 5E-06) with FEV1 and FEV1/FVC in LLFS were evaluated for their association with FEV1 and FEV1/FVC in the CHARGE/SpiroMeta consortia  using statistical models identical to those described above (individual family-based studies within the CHARGE/SpiroMeta did adjust for kinship structure while other studies within the consortia did not adjust for kinship structure). Finally, genotypes from both LLFS and CHARGE/SpiroMeta were meta-analyzed using METAL to evaluate the overall association between individual SNPs and lung function in both studies. We also evaluated the replication of previously reported GWAS hits (from the CHARGE/SpiroMeta consortia) for FEV1 and FEV1/FVC in the LLFS population using the statistical models described above.
To calculate Identity by Descent (IBD) for the linkage analyses, the ZAPLO program was used to estimate haplotypes of SNPs in small regions (0.5 cM intervals) . The deCODE map was used to approximate the cM positions . We identified all SNPs in our GWAS scan with no Mendel inconsistencies and an average pedigree heterozygosity ≥0.1. Within each 0.5 cM interval we used the first five such SNPs to construct a haplotype and if there were fewer than 5, we took all such SNPs in the interval. With the resulting haplotypes, IBD at 1 cM intervals was estimated in the Loki program , which does chromosome-wide IBD estimation in intact pedigrees. These IBD estimates were then used in the SOLAR package  to conduct a variance-component linkage analysis. In the linkage analysis, we used the residuals of the FEV1 and FEV1/FVC phenotypes, adjusted for age, sex, height, ancestry PCs, smoking status (current, former, never), pack-years, and field center. The SOLAR package was also used to estimate heritability and empirical p of LOD.
Comparison of demographic and lung function variables * in LLFS and CHARGE/SpiroMeta consortia +
LLFS (n = 3889)
CHARGE/SpiroMeta (n = 48201)
Sex (% Male)
Smoking Status (% never smoker)
GWAS findings for FEV 1 in the LLFS (p <1E-07) and replication in the CHARGE/SpiroMeta consortia
Coded allele frequency
Coded allele frequency
Coded allele frequency
Chr 6: 144127026
GWAS findings for FEV 1 /FVC in the LLFS (p <1E-07) and replication in the CHARGE/SpiroMeta consortia
Coded allele frequency
Coded allele frequency
Coded allele frequency
This study identified a novel linkage peak in chromosome 2 for FEV1/FVC and confirmed a previously identified linkage peak in chromosome 6 for FEV1. This study also replicated some of the GWAS SNPs identified in previous studies. Though this study identified SNPs in the CYP2U1 and PHACTR2 genes that were associated with FEV1, these findings could not be replicated in independent datasets.
Two previous reports identified a linkage peak for FEV1 at the same locus (chromosome 6 at 28 cM (~280 kb)) as the present study ,. As compared to previous linkage studies we had improved marker density (average marker spacing = 0.5 cM) and improved accuracy by using haplotype-based IBD estimation. The availability of GWAS data also allowed us to evaluate whether GWAS SNPs under the linkage peak explained the linkage peak identified for FEV1. While we found a modest attenuation of the FEV1 linkage peak after adjustment for common GWAS SNPs under the linkage peak that were nominally associated with FEV1 (p < 0.001), they did not completely explain the linkage peak. In contrast, the novel linkage peak identified in chromosome 2 at 219 cM (~372 kb) for FEV1/FVC was almost completely explained, by adjustment for the common GWAS SNPs under the linkage peak that were nominally associated with FEV1/FVC (p < 0.001). These findings suggest that common variants alone are insufficient to explain some linkage peaks such as the linkage peak in chromosome 6 for FEV1. As shown in other diseases, the inability to identify association under linkage peaks could in part be attributable to the fact that only common variants are examined under the linkage peak whereas the linkage signal could be caused by multiple rare variants with higher penetrance . Hence, the contribution of multiple rare variants with high penetrance under the linkage peak towards FEV1 needs to be evaluated in future studies. The linkage peak identified in chromosome 6 at 28 cM (~280 kb) for FEV1 does not contain any known genes. However, there are several DNase I hypersensitivity sites and putative transcription factor binding sites (9H3K27Ac marks) that have been identified in cells from pulmonary epithelium and blood vessels derived from the pulmonary artery within this linkage peak (EnCode data) suggesting that regulatory elements in this region may play an important role in determining lung function. The linkage peak for FEV1/FVC in chromosome 2 at 219 cM (~373 kb) contains the genomic region that codes for the DIRC3 gene, a non coding RNA that is involved in the pathogenesis of familial renal cancers (EnCode data). Though DIRC3 is expressed in pulmonary tissue its role in determining lung function has not been evaluated. Furthermore, this region also contains DNase I hypersensitivity sites or putative transcription factor binding sites (9H3K27Ac marks) in cells derived from the pulmonary epithelium or vasculature. Thus the linkage peaks identified in this study on chromosomes 2 and 6 may indicate previously unidentified regulatory pathways that may influence longevity through their effect on lung function.
Though this study identified a few novel GWAS SNPs that were associated with FEV1 and FEV1/FVC, none of these findings could be replicated in the CHARGE/SpiroMeta consortia. However, only 38% of SNPs associated with FEV1 or FEV1/FVC in LLFS were available for replication in the CHARGE/SpiroMeta consortia. Thus, the remaining SNPs and insertion/deletions polymorphisms associated with FEV1 and FEV1/FVC in LLFS but not genotyped/imputed in the CHARGE/SpiroMeta consortia will need to be evaluated in future studies. Since the LLFS study participants were not randomly selected to represent the general population but were specifically selected for their family history of exceptional longevity, it is possible that there may be some unique genotypes associated with lung function that may not be replicated in studies that are more representative of the general population. One previous study showed that elderly male offspring (range: 65–89 years) with long lived parents (age at death of at least on parent >80 years) had FEV1 that was 330 mL larger than FEV1 for male offspring with short lived parents even after controlling for smoking . The findings of the present study along with previous findings  might indicate that long lived families have unique genetic variants that contribute to higher lung function among those with exceptional longevity. This study also replicated 12 of 31 GWAS hits for FEV1 and FEV1/FVC identified in previous meta-analysis (Additional file 2: Tables S5 and S6). The reasons for not confirming all previously identified variants in this study may include limited power, population specific genetic heterogeneity and differences in environmental exposures such as smoking. Genetic heterogeneity in contribution towards determination of lung function in elderly participants (average age among LLFS participants: 68.8 years) as compared to middle aged adults may also contribute to the differences in genetic association observed in the LLFS study as compared to previous studies. As shown in Additional file 2: Tables S3 and S4, several SNPs had highly significant p for heterogeneity between the LLFS and the CHARGE/SpiroMeta dataset supporting the idea that, at least for some loci, heterogeneity in genetic contribution may account for the lack of replication of certain loci.
The family-based cohort design of the LLFS with extensive genotype information and detailed lung function measurements makes this study a valuable resource to identify genetic determinants of lung function. In addition to confirming some of the previously identified GWAS SNPs and a previously identified linkage peak in chromosome 6 for FEV1, this study also identified a novel linkage peak in chromosome 2 for FEV1/FVC. Repeated measurements of lung function in this study population along with targeted resequencing under the observed linkage peaks in future studies may help clarify the role of genetic variants in determining preserved lung function among exceptionally long lived individuals.
BT participated in the design of the study, the data analysis, data interpretation and manuscript preparation. MW participated in the data analysis, data interpretation and manuscript preparation. RM, JS, LC and SB participated in the interpretation of the data and manuscript editing. RGB was involved in the Spirometry data collection and manuscript editing. Both the CHARGE and SpiroMeta consortia were involved in data replication in independent dataset and manuscript review. AN participated in the data analysis, data interpretation and manuscript review. All authors read and approved the final manuscript.
Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, SpiroMeta consortium.
The investigators thank LLFS participants and staff for their contributions.
- Hegewald MJCR, et al: Pulmonary function testing. Murray and Nadel’s Textbook of Respiratory Medicine. Edited by: Mason RJBV, Martin TR. 2010, Saunders Elsevier, Philadelphia, Pa, 5Google Scholar
- Ashley F, Kannel WB, Sorlie PD, Masson R: Pulmonary function: relation to aging, cigarette habit, and mortality. Ann Intern Med. 1975, 82: 739-745. 10.7326/0003-4819-82-6-739.PubMedView ArticleGoogle Scholar
- Beaty TH, Cohen BH, Newill CA, Menkes HA, Diamond EL, Chen CJ: Impaired pulmonary function as a risk factor for mortality. Am J Epidemiol. 1982, 116: 102-113.PubMedGoogle Scholar
- Chinn S, Gislason T, Aspelund T, Gudnason V: Optimum expression of adult lung function based on all-cause mortality: results from the Reykjavik study. Respir Med. 2007, 101: 601-609. 10.1016/j.rmed.2006.06.009.PubMedView ArticleGoogle Scholar
- Higgins MW, Keller JB: Predictors of mortality in the adult population of Tecumseh. Arch Environ Health. 1970, 21: 418-424. 10.1080/00039896.1970.10667260.PubMedView ArticleGoogle Scholar
- Hole DJ, Watt GC, Davey-Smith G, Hart CL, Gillis CR, Hawthorne VM: Impaired lung function and mortality risk in men and women: findings from the Renfrew and Paisley prospective population study. BMJ. 1996, 313: 711-715. 10.1136/bmj.313.7059.711. discussion 715–716PubMedPubMed CentralView ArticleGoogle Scholar
- Lange P, Nyboe J, Appleyard M, Jensen G, Schnohr P: Spirometric findings and mortality in never-smokers. J Clin Epidemiol. 1990, 43: 867-873. 10.1016/0895-4356(90)90070-6.PubMedView ArticleGoogle Scholar
- Schunemann HJ, Dorn J, Grant BJ, Winkelstein W, Trevisan M: Pulmonary function is a long-term predictor of mortality in the general population: 29-year follow-up of the Buffalo Health Study. Chest. 2000, 118: 656-664. 10.1378/chest.118.3.656.PubMedView ArticleGoogle Scholar
- Sin DD, Wu L, Man SF: The relationship between reduced lung function and cardiovascular mortality: a population-based study and a systematic review of the literature. Chest. 2005, 127: 1952-1959. 10.1378/chest.127.6.1952.PubMedView ArticleGoogle Scholar
- Jacobsen PK, Sigsgaard T, Pedersen OF, Christensen K, Miller MR: Lung function as a predictor of survival in very elderly people: the Danish 1905 cohort study. J Am Geriatr Soc. 2008, 56: 2150-2152. 10.1111/j.1532-5415.2008.01969.x.PubMedView ArticleGoogle Scholar
- Hollenberg M, Yang J, Haight TJ, Tager IB: Longitudinal changes in aerobic capacity: implications for concepts of aging. J Gerontol A Biol Sci Med Sci. 2006, 61: 851-858. 10.1093/gerona/61.8.851.PubMedView ArticleGoogle Scholar
- Newman AB, Glynn NW, Taylor CA, Sebastiani P, Perls TT, Mayeux R, Christensen K, Zmuda JM, Barral S, Lee JH, Simonsick EM, Walston JD, Yashin AI, Hadley E: Health and function of participants in the long life family study: a comparison with other cohorts. Aging (Albany NY). 2011, 3: 63-76.Google Scholar
- Jedrychowski W: Effects of smoking and longevity of parents on lung function in the apparently healty elderly. Arch Gerontol Geriatr. 1990, 10: 19-26. 10.1016/0167-4943(90)90040-D.PubMedView ArticleGoogle Scholar
- Givelber RJ, Couropmitree NN, Gottlieb DJ, Evans JC, Levy D, Myers RH, O’Connor GT: Segregation analysis of pulmonary function among families in the Framingham Study. Am J Respir Crit Care Med. 1998, 157: 1445-1451. 10.1164/ajrccm.157.5.9704021.PubMedView ArticleGoogle Scholar
- Wilk JB, Djousse L, Arnett DK, Rich SS, Province MA, Hunt SC, Crapo RO, Higgins M, Myers RH: Evidence for major genes influencing pulmonary function in the NHLBI family heart study. Genet Epidemiol. 2000, 19: 81-94. 10.1002/1098-2272(200007)19:1<81::AID-GEPI6>3.0.CO;2-8.PubMedView ArticleGoogle Scholar
- Zhou JJ, Cho MH, Castaldi PJ, Hersh CP, Silverman EK, Laird NM: Heritability of chronic obstructive pulmonary disease and related phenotypes in smokers. Am J Respir Crit Care Med. 2013, 188: 941-947. 10.1164/rccm.201302-0263OC.PubMedPubMed CentralView ArticleGoogle Scholar
- Hancock DB, Eijgelsheim M, Wilk JB, Gharib SA, Loehr LR, Marciante KD, Franceschini N, van Durme YM, Chen TH, Barr RG, Schabath MB, Couper DJ, Brusselle GG, Psaty BM, van Duijn CM, Rotter JI, Uitterlinden AG, Hofman A, Punjabi NM, Rivadeneira F, Morrison AC, Enright PL, North KE, Heckbert SR, Lumley T, Stricker BH, O'Connor GT, London SJ: Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet. 2010, 42: 45-52. 10.1038/ng.500.PubMedPubMed CentralView ArticleGoogle Scholar
- Soler Artigas M, Loth DW, Wain LV, Gharib SA, Obeidat M, Tang W, Zhai G, Zhao JH, Smith AV, Huffman JE, Albrecht E, Jackson CM, Evans DM, Cadby G, Fornage M, Manichaikul A, Lopez LM, Johnson T, Aldrich MC, Aspelund T, Barroso I, Campbell H, Cassano PA, Couper DJ, Eiriksdottir G, Franceschini N, Garcia M, Gieger C, Gislason GK, Grkovic I, et al: Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat Genet. 2011, 43: 1082-1090. 10.1038/ng.941.PubMedView ArticleGoogle Scholar
- Repapi E, Sayers I, Wain LV, Burton PR, Johnson T, Obeidat M, Zhao JH, Ramasamy A, Zhai G, Vitart V, Huffman JE, Igl W, Albrecht E, Deloukas P, Henderson J, Granell R, McArdle WL, Rudnicka AR, Barroso I, Loos RJ, Wareham NJ, Mustelin L, Rantanen T, Surakka I, Imboden M, Wichmann HE, Grkovic I, Jankovic S, Zgaga L, et al: Genome-wide association study identifies five loci associated with lung function. Nat Genet. 2010, 42: 36-44. 10.1038/ng.501.PubMedPubMed CentralView ArticleGoogle Scholar
- Pedersen CB, Gotzsche H, Moller JO, Mortensen PB: The Danish Civil Registration System. a cohort of eight million persons. Dan Med Bull. 2006, 53: 441-449.PubMedGoogle Scholar
- Sebastiani P, Hadley EC, Province M, Christensen K, Rossi W, Perls TT, Ash AS: A family longevity selection score: ranking sibships by their longevity, size, and availability for study. Am J Epidemiol. 2009, 170: 1555-1562. 10.1093/aje/kwp309.PubMedPubMed CentralView ArticleGoogle Scholar
- Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, Jensen R, Johnson DC, MacIntyre N, McKay R, Navajas D, Pedersen OF, Pellegrino R, Viegi G, Wanger J: Standardisation of spirometry. Eur Respir J. 2005, 26: 319-338. 10.1183/09031936.05.00034805.PubMedView ArticleGoogle Scholar
- Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.PubMedView ArticleGoogle Scholar
- Li Y, Willer C, Sanna S, Abecasis G: Genotype imputation. Annu Rev Genomics Hum Genet. 2009, 10: 387-406. 10.1146/annurev.genom.9.081307.164242.PubMedPubMed CentralView ArticleGoogle Scholar
- Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR: MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010, 34: 816-834. 10.1002/gepi.20533.PubMedPubMed CentralView ArticleGoogle Scholar
- Pinheiro JC, Bates D: Mixed-effects models in S and S-Plus. 2000, Springer, New YorkView ArticleGoogle Scholar
- Cortinas Abrahautes J, Burzykowski T: A version of the EM algorithm for proportional hazards models with random effects. Lecture Notes of the ICB Seminars. 2002, 15-20.Google Scholar
- O’Connell JR: Zero-recombinant haplotyping: applications to fine mapping using SNPs. Genet Epidemiol. 2000, 19 (Suppl 1): S64-S70. 10.1002/1098-2272(2000)19:1+<::AID-GEPI10>3.0.CO;2-G.PubMedView ArticleGoogle Scholar
- Kong A, Thorleifsson G, Gudbjartsson DF, Masson G, Sigurdsson A, Jonasdottir A, Walters GB, Jonasdottir A, Gylfason A, Kristinsson KT, Gudjonsson SA, Frigge ML, Helgason A, Thorsteinsdottir U, Stefansson K: Fine-scale recombination rate differences between sexes, populations and individuals. Nature. 2010, 467: 1099-1103. 10.1038/nature09525.PubMedView ArticleGoogle Scholar
- Heath SC: Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet. 1997, 61: 748-760. 10.1086/515506.PubMedPubMed CentralView ArticleGoogle Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMedPubMed CentralView ArticleGoogle Scholar
- Joost O, Wilk JB, Cupples LA, Harmon M, Shearman AM, Baldwin CT, O’Connor GT, Myers RH, Gottlieb DJ: Genetic loci influencing lung function: a genome-wide scan in the Framingham Study. Am J Respir Crit Care Med. 2002, 165: 795-799. 10.1164/ajrccm.165.6.2102057.PubMedView ArticleGoogle Scholar
- Wilk JB, DeStefano AL, Joost O, Myers RH, Cupples LA, Slater K, Atwood LD, Heard-Costa NL, Herbert A, O’Connor GT, Gottlieb DJ: Linkage and association with pulmonary function measures on chromosome 6q27 in the Framingham Heart Study. Hum Mol Genet. 2003, 12: 2745-2751. 10.1093/hmg/ddg311.PubMedView ArticleGoogle Scholar
- Eyre S, Ke X, Lawrence R, Bowes J, Panoutsopoulou K, Barton A, Thomson W, Worthington J, Zeggini E: Examining the overlap between genome-wide rare variant association signals and linkage peaks in rheumatoid arthritis. Arthritis Rheum. 2011, 63: 1522-1526. 10.1002/art.30315.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/4.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.