Linkage analysis of anti-CCP levels as dichotomized and quantitative traits using GAW15 single-nucleotide polymorphism scan of NARAC families
© Yang et al; licensee BioMed Central Ltd. 2007
Published: 18 December 2007
Rheumatoid arthritis is a clinically and genetically heterogeneous disease. Anti-cyclic citrullinated (anti-CCP) antibodies have a high specificity for rheumatoid arthritis and levels correlate with disease severity. The focus of this study was to examine whether analyzing anti-CCP levels could increase the power of linkage analysis by identifying a more homogeneous subset of rheumatoid arthritis patients. We also wanted to compare linkage signals when analyzing anti-CCP levels as dichotomized (CCP_binary), categorical (CCP_cat), and continuous traits, with and without transformation (log_CCP and CCP_cont). Illumina single-nucleotide polymorphism scans of the North American Rheumatoid Arthritis Consortium families were analyzed for four chromosomes (6, 7, 11, 22) using nonparametric linkage (NPL) (rheumatoid arthritis and CCP_binary), regress (CCP_cat and Log_CCP), and deviates (CCP_cont) analysis options as implemented in Merlin. Similar linkage results were obtained from analyses of rheumatoid arthritis, CCP_binary, and CCP_cont. The only exception was that we observed improved linkage signals and a narrower region for CCP_binary as compared to a clinical diagnosis of rheumatoid arthritis alone on chromosome 7, a region which previously showed variation in linkage results with rheumatoid arthritis according to anti-CCP levels. Analyses of CCP_cat and Log_CCP had little power to detect linkage. Our data suggested that linkage analyses of anti-CCP levels may facilitate identification of rheumatoid arthritis genes but quantitative analyses did not further improve power. Our study also highlighted that quantitative trait linkage results are highly sensitive to phenotype transformation and analytic approaches.
Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease affecting about 1% of the population. A genetic component for RA has been well established, with the MHC region being the largest single contributing component. Other chromosome regions (11q, 10q, 14q, 6p, 6q, 16q, 12p, etc.) and candidate genes (PTPN22, CTLA4, PADI4) have been identified by whole-genome linkage scans and association studies [1–5]. Most recently, a high-density SNP analysis of 642 families affected with RA collected by the North American Rheumatoid Arthritis Consortium (NARAC), the largest single linkage study of RA, identified two new linkage regions, 11p and 2q . These findings reflect the genetic complexity of the disease and suggest that analysis of a more homogeneous RA phenotype might increase the power of linkage analysis. In addition, most previous studies have analyzed RA as a dichotomous trait, which could lead to a power loss if RA is a naturally quantitative trait .
Anti-cyclic citrullinated (anti-CCP) antibodies have a high specificity for RA  and the levels are correlated with disease severity [9, 10]. To examine whether the power of linkage analysis could be improved by analyzing a more homogeneous phenotype and by quantitative characterization of the trait, we performed linkage analysis of anti-CCP antibody levels for selected chromosome regions previously linked to RA using NARAC data.
Illumina SNP scans of the NARAC families were analyzed. We chose anti-CCP antibody levels as the phenotype of interest, and evaluated the effect of covariates including sex, age of onset, year of birth, ever/never smoking, and current smoking. Anti-CCP antibody levels were analyzed in three ways: dichotomized, categorical, and continuous. An antibody titer of 20 was used as a cut-off value to dichotomize anti-CCP levels into positive (>20) and negative (≤20). In addition, anti-CCP levels were characterized into multiple (four) categories (negative, 0–19.9; low, 20–49.9; medium, 50–99.9; high, ≥100). Anti-CCP levels were also analyzed as continuous measurements. Because the assay for anti-CCP antibody titer has an upper limit of 210, we recoded all measurements exceeding 210 to 210. A log transformation was applied to approximate normality of anti-CCP levels because the raw data were highly skewed.
Because the purpose of this study was to compare methods rather than search for a new locus for anti-CCP, we limited our linkage analyses to selected regions. Chromosomes were selected based on findings from a previous SNP scan of NARAC families for RA . Chromosome 6, which contains human lymphocyte antigen (HLA) locus, was chosen as a positive control region. Chromosome 22, which did not show evidence for linkage with RA, was selected as the negative control. Chromosomes 7 and 11, which showed suggestive and significant evidence for linkage with RA, respectively, were included in our linkage analyses as test regions. In particular, it has been shown that chromosome 7 might harbor a susceptibility locus that was more closely linked to anti-CCP positive disease .
SNPs on chromosomes 6, 7, 11, and 22 were analyzed for all four anti-CCP phenotypes (CCP_binary, CCP_cat, Log_CCP, and CCP_cont) as well as RA affection status. Linkage disequilibrium (LD) between markers was calculated and markers in LD defined by D' > 0.7 were removed using SNPLINK . CCP_binary and RA affection were analyzed by nonparametric (NPL) linkage analysis using Merlin. CCP_cat and Log_CCP were analyzed by regression analysis implemented in Merlin Regress, which uses trait-squared sums and differences to predict IBD sharing between sib pairs . To run Merlin Regress, it was necessary to specify some trait distribution parameters, such as mean, variance, and heritability in the general population. We did not use the sample mean and variance because the families were affected with RA and therefore had higher frequencies and levels of anti-CCP positives. Instead, we estimated the mean (0 for CCP_cat, 0.78 for Log_CCP) and variance (0.0088 for CCP_cat, 0.11 for Log_CCP) among individuals who were anti-CCP negative (≤20) to approximate the distribution in the general population. Heritability of 0.6 was estimated based on variance-component analysis using Merlin. Untransformed continuous anti-CCP levels were also analyzed using the deviates option implemented in Merlin, which makes no assumptions about the trait distribution. Again, we chose 0 to approximate the population mean for anti-CCP levels.
Distribution of anti-CCP levels among the 1499 individuals included in the linkage analyses
Anti-CCP levels by gender, age, year of birth, and smoking status
t test P
Age of onset
Year of birth
Never or former
Maximal LOD scores for RA and anti-CCP phenotypes in the four chromosome regions
In this study, we analyzed anti-CCP levels as dichotomized, categorical, and continuous traits in linkage analyses. On chromosome 7, we observed improved linkage signals and a narrower linkage peak for dichotomized anti-CCP levels as compared to a clinical diagnosis of RA alone, although neither analysis reached statistical significance. Chromosome 7 previously showed stronger evidence for linkage for RA with positive anti-CCP . However, for the three chromosomes we selected that demonstrated evidence for linkage with RA, analyzing anti-CCP levels as quantitative traits did not further improve power to detect linkage compared to RA. In addition, categorization and log transformation of the continuous phenotype resulted in significantly reduced power to detect linkage.
Mapping susceptibility genes is challenging for clinically and genetically heterogeneous diseases such as RA. Examining a more homogenous disease phenotype might increase the power of linkage analysis. Anti-CCP is highly specific for RA and its levels are correlated with disease severity [9, 10]. A previous study suggested genetic variability according to anti-CCP status on chromosomes 4, 5, 6, and 7 . Our data on chromosome 7 supports the hypothesis that anti-CCP status might represent a more genetically homogeneous phenotype of RA. In addition, although linkage signals on chromosome 11 were less significant for the CCP_binary phenotype than RA, linkage peak regions on both chromosomes 7 and 11 were narrower, suggesting that analyzing anti-CCP levels as a dichotomous phenotype might be helpful in narrowing candidate regions in fine mapping.
For naturally occurring continuous or polychotomous traits, dichotomization could lead to power loss for linkage analysis . In particular, if the underlying gene confers not only disease susceptibility but also disease severity, treating disease phenotypes as quantitative traits could provide additional information in linkage analysis. We hypothesized that performing quantitative linkage analyses might increase the power to detect linkage for anti-CCP levels. However, analyses of none of the three quantitative phenotypes provided improved linkage signals compared to NPL analysis of RA. In particular, regression-based quantitative linkage analyses of log-transformed or ordinal anti-CCP levels, with or without covariate adjustment, appeared to have reduced power to detect linkage compared to NPL analysis of dichotomized anti-CCP status. We also used the variance components and quantitative trait locus analysis options implemented in Merlin for CCP_cat and Log_CCP phenotypes and observed similar results. For this particular trait, transforming the phenotype to achieve normality reduced the signal due to the major locus. Alternate approaches to handle non-normality are needed for quantitative trait linkage analysis in selected samples. Consistent with another report , we found loss of power when treating CCP_cat as continuous. Novel approaches, such as the recently proposed proportional odds latent variable model by Feng et al. , need to be further explored. Finally, all these quantitative linkage analyses are very sensitive to the specification of population parameters, such as population mean and variance, chosen for the regress and deviates analyses. Replacement of the estimated population mean by the sample mean or default values set by the programs resulted in completely different LOD scores. These results further highlight the importance of precise estimates and appropriate analytical approaches in linkage analysis of quantitative traits.
Linkage analyses of anti-CCP levels may facilitate identification of RA genes by making the phenotype more homogeneous. However, quantitative analyses of anti-CCP levels did not improve power to detect linkage over analysis of RA. Furthermore, log transformation and categorization of anti-CCP levels resulted in reduced power to detect linkage.
This article has been published as part of BMC Proceedings Volume 1 Supplement 1, 2007: Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/1?issue=S1.
- Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Monteiro J, Kern M, Criswell LA, Albani S, Nelson JL, Clegg DO, Pope R, Schroeder HW, Bridges SL, Pisetsky DS, Ward R, Kastner DL, Wilder RL, Pincus T, Callahan LF, Flemming D, Wener MH, Gregersen PK: A genome-wide screen in multiplex rheumatoid arthritis families suggests genetic overlap with other autoimmune diseases. Am J Hum Genet. 2001, 68: 927-936. 10.1086/319518.View ArticlePubMed CentralPubMedGoogle Scholar
- Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Etzel C, Damle A, Xiao X, Chen D, Lum RF, Monteiro J, Kern M, Criswell LA, Albani S, Nelson JL, Clegg DO, Pope R, Schroeder HW, Bridges SL, Pisetsky DS, Ward R, Kastner DL, Wilder RL, Pincus T, Callahan LF, Flemming D, Wener MH, Gregersen PK, North American Rheumatoid Arthritis Consortium: Screening the genome for rheumatoid arthritis susceptibility genes. A replication study and combined analysis of 512 multicase families. Arthritis Rheum. 2003, 48: 906-916. 10.1002/art.10989.View ArticlePubMedGoogle Scholar
- Osorio Y, Fortéa J, Bukulmez H, Petit-Teixeira E, Michou L, Pierlot C, Cailleau-Moindrault S, Lemaire I, Lasbleiz S, Alibert O, Quillet P, Bardin T, Prum B, Olson JM, Cornélis F: Dense genome-wide linkage analysis of rheumatoid arthritis, including covariates. Arthritis Rheum. 2004, 50: 2757-2765. 10.1002/art.20458.View ArticleGoogle Scholar
- Tamiya G, Shinya M, Imanishi T, Ikuta T, Makino S, Okamoto K, Furugaki K, Matsumoto T, Mano S, Ando S, Nozaki Y, Yukawa W, Nakashige R, Yamaguchi D, Ishibashi H, Yonekura M, Nakami Y, Takayama S, Endo T, Saruwatari T, Yagura M, Yoshikawa Y, Fujimoto K, Oka A, Chiku S, Linsen SE, Giphart MJ, Kulski JK, Fukazawa T, Hashimoto H, Kimura M, Hoshina Y, Suzuki Y, Hotta T, Mochida J, Minezaki T, Komai K, Shiozawa S, Taniguchi A, Yamanaka H, Kamatani N, Gojobori T, Bahram S, Inoko H: Whole genome association study of rheumatoid arthritis using 27,039 microsatellites. Hum Mol Genet. 2005, 14: 2305-2321. 10.1093/hmg/ddi234.View ArticlePubMedGoogle Scholar
- Oliver JE, Worthington J, Silman AJ: Genetic epidemiology of rheumatoid arthritis. Curr Opin Rheumatol. 2006, 18: 141-146. 10.1097/01.bor.0000209425.84775.23.View ArticlePubMedGoogle Scholar
- Amos CI, Chen WV, Lee A, Li W, Kern M, Lundsten R, Batliwalla F, Wener M, Remmers E, Kastner DA, Criswell LA, Seldin MF, Gregersen PK: High-density SNP analysis of 642 Caucasian families with rheumatoid arthritis identifies two new linkage regions on 11p12 and 2q33. Genes Immun. 2006, 7: 277-286. 10.1038/sj.gene.6364295.View ArticlePubMedGoogle Scholar
- Corbett J, Gu CC, Rice JP, Reich T, Province MA, Rao DC: Power loss for linkage analysis due to the dichotomization of trichotomous phenotypes. Hum Hered. 2004, 57: 21-27. 10.1159/000077386.View ArticlePubMedGoogle Scholar
- Rantapaa-Dahlqvist S, de Jong BAW, Berglin W, Hallmans G, Wadell G, Stenlund H, Sundin U, van Venrooij WJ: Antibodies against cyclic citrullinated peptide and IgA rheumatoid factor predict the development of rheumatoid arthritis. Arthitis Rheum. 2003, 48: 2741-2749. 10.1002/art.11223.View ArticleGoogle Scholar
- Kastbom A, Strandberg G, Lindroos A, Skogh T: Anti-CCP antibody test predicts the disease course during 3 years in early rheumatoid arthritis (the Swedish TIRA project). Ann Rheum Dis. 2004, 63: 1085-1089. 10.1136/ard.2003.016808.View ArticlePubMed CentralPubMedGoogle Scholar
- del Val Del Amo N, Ibanez Bosch R, Fito Manteca C, Gutierrez Polo R, Loza Cortina E: Anti-cyclic citrullinated peptide antibody in rheumatoid arthritis: relation to disease aggressiveness. Clin Exp Rheumatol. 2006, 24: 281-286.PubMedGoogle Scholar
- Webb El, Sellick GS, Houlston RS: SNPLINK: multipoint linkage analysis of densely distributed DNP data incorporating automated linkage disequilibrium removal. Bioinformatics. 2005, 21: 3060-3061. 10.1093/bioinformatics/bti449.View ArticlePubMedGoogle Scholar
- Sham PC, Purcell S, Cherny SS, Abecasis GR: Powerful regression-based quantitative-trait linkage analysis of general pedigrees. Am J Hum Genet. 2002, 71: 238-253. 10.1086/341560.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang X, Ye Y, Zhang H: Family-based association tests for ordinal traits adjusting for covariates. Genet Epidemiol. 2006, 30: 728-736. 10.1002/gepi.20184.View ArticlePubMedGoogle Scholar
- Feng R, Leckman JF, Zhang H: Linkage analysis of ordinal traits for pedigree data. Proc Natl Acad Sci USA. 2004, 101: 16739-16744. 10.1073/pnas.0404623101.View ArticlePubMed CentralPubMedGoogle 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 cited.