<title>Introduction</title> <p>Rheumatoid arthritis is the most common systemic autoimmune disease affecting approximately 1% of the adult population worldwide, with prevalence varying from 0.2–0.3% in East Asians to 6% in Pima and Chippewa Indians ##REF##12110146##[1]##. The disease is characterized by inflammation of the synovial tissue and local articular damage ##REF##12748655##[2]##. Disability in this inflammatory polyarthritis primarily stems from progressive bone erosion and comorbidity with coronary artery disease, infection and lymphoma ##REF##16914706##[3]##,##REF##8605261##[4]##. As with many other autoimmune conditions, RA affects women more commonly than men.</p> <p>Although the etiology of RA is presently unknown, studies of RA heritability in two Northern European regions have demonstrated that an average of 60% of the disease variance can be attributed to genetic factors ##REF##10643697##[5]##. Through a combination of linkage and association studies, alleles segregating at the human leukocyte antigen (<italic>HLA</italic>) class II <italic>DRB1</italic> gene on chr 6p have consistently been shown to have strong RA-predisposing effects ##REF##147420##[6]##,##REF##7942852##[7]##. That said, studies suggest that <italic>HLA-DRB1</italic> accounts for at most 50% of the phenotypic variance due to genetic effects ##REF##2059076##[8]##; therefore, loci not linked to the HLA region may play a crucial role in RA susceptibility.</p> <p>Utilizing a variety of approaches such as positional mapping, candidate gene experiments and large-scale functional genetic association studies, several recent reports have yielded evidence for additional RA genes. The most robust, non-MHC, RA-associated marker is the R620W missense polymorphism in the <italic>PTPN22</italic> gene on chromosome 1p13, which has been repeatedly associated with RA in individuals of European ancestry ##REF##15208781##[9]##–##REF##16731003##[11]##. In addition, positional cloning work has suggested the peptidyl arginine deiminase gene cluster (including <italic>PADI4</italic>) underneath a linkage peak on chr 1p36 may harbor susceptibility variants ##REF##12833157##[12]##,##REF##17265154##[13]## while well-powered association studies identified RA-associated SNPs in <italic>STAT4</italic>
##REF##17804842##[14]## and the <italic>TNFAIP3</italic> region ##REF##17982456##[15]##,##REF##17982455##[16]##. A promoter polymorphism of the Fc receptor-like 3 gene, <italic>FCRL3</italic>, and a SNP within the RUNX1 binding site of <italic>SLC22A4</italic> have also been implicated in RA susceptibility ##REF##15838509##[17]##–##REF##14608356##[19]##, both with conflicting reports ##REF##17133581##[20]##,##REF##16821265##[21]##.</p> <p>Interestingly, some of these disease-associated polymorphisms appear to have heterogeneity in effect sizes across ethnic groups; for example, the disease-associated variants in <italic>PADI4</italic> and <italic>FCRL3</italic> have a strong effect in East Asians but little effect in whites of European descent ##REF##16380915##[10]##,##REF##17763442##[22]##. Similarly, the <italic>PTPN22</italic> W620 risk allele is virtually absent in East Asians and therefore plays no role in RA risk in these populations ##REF##16731003##[11]##. As RA is a major cause of disability and is correlated with increased mortality in severe cases, genetic studies promise to improve public health. Importantly, as predicted by careful meta-analyses of linkage studies ##REF##9707586##[23]##, some RA-susceptibility variants show pleiotropic effects across many autoimmune diseases [e.g. 11,14,24,25]. Consequently, further identification of RA genetic risk factors should aid in elucidating the underlying mechanisms of autoimmunity, in general, and may substantially impact drug discovery through the development of targeted diagnostics and therapeutics.</p> <p>Arguing that the power of linkage disequilibrium-based designs to map disease alleles is high compared to other approaches, Jorde ##REF##7825565##[26]##, Risch and Merikangas ##REF##8801636##[27]## and Long and colleagues ##UREF##0##[28]## helped motivate the recent wave of successful genome-wide disease association studies. Propelled by technological developments, this shift has recently transformed common, complex disease gene mapping resulting in a number of convincing susceptibility variants [e.g. 29–31]. We took a large-scale candidate SNP association approach, very similar to that used in our recent study of psoriasis ##REF##17236132##[32]##, to interrogate the genome for genetic variants that predispose individuals to RA. This genome-wide SNP panel (25,966 SNPs), which is primarily composed of missense (70%), acceptor/donor splice site and putative transcription-factor binding site SNPs, was applied to a multi-tiered, case-control association study of RA that incorporated replication of association effects as a key feature of the study design. By directly interrogating polymorphisms with higher likelihoods of producing biologically disruptive effects across multiple large sample sets, our aim was to maximize power to detect RA susceptibility genes.</p> <p>We previously reported the identity of the RA-associated <italic>PTPN22</italic> R620W variant which was discovered in the first step (quality control of all DNA samples) of our RA scan ##REF##15208781##[9]##,##REF##16175503##[33]##. Here, we report our finding of variants in the <italic>PHF19-TRAF1-C5</italic> region on chromosome 9q33.2 that show strong and consistent association across three independent RA case-control studies (1732 cases/2502 controls), paralleling and extending the results of a whole-genome association study ##UREF##1##[34]## and a candidate gene study ##REF##17880261##[35]##. Combining genetic information from <italic>HLA</italic>, <italic>PTPN22</italic> and <italic>TRAF1</italic> variants, we calculate the posterior probability of RA for every possible genotype combination. Results such as these may form the foundation for individualized prognosis and targeted medicine.</p>
<title>Materials and Methods</title> <title>Subjects and Samples</title> <p>All RA cases included in this study were white and met the 1987 American College of Rheumatology diagnostic criteria for RA ##REF##3358796##[69]##; informed written consent was obtained from every subject. Sample Set 1, which consisted of 475 RA cases and 475 individually-matched controls, was collected by Genomics Collaborative, Inc. All case samples were white North Americans of European descent who where rheumatoid factor (RF) positive. Control samples were healthy white individuals with no medical history of RA, also of European descent. A single control was matched to each case on the basis of sex, age (±5 years), and self-reported ethnic background. The 661 cases in Sample Set 2 were acquired from the North American Rheumatoid Arthritis Consortium (NARAC) (<ext-link ext-link-type="uri" xlink:href="http://www.naracdata.org/">http://www.naracdata.org/</ext-link>) and consisted of members from 661 white North American multiplex families ##REF##16175503##[33]##,##REF##11254450##[70]##,##REF##15022313##[71]##. Both RF-positive and RF-negative patients were included in this sample set. Controls for Sample Set 2 were selected from 20,000 healthy individuals enrolled in the New York Cancer Project ##REF##15136663##[72]##, a population-based prospective study of the genetic and environmental factors that cause disease (<ext-link ext-link-type="uri" xlink:href="http://www.amdec.org/">http://www.amdec.org/</ext-link>). Two control individuals were matched to a single, randomly chosen affected sibling from each NARAC family on the basis of sex, age (decade of birth), and self-reported ethnic background. Sample Set 3 was composed of 596 white RA patients from the Leiden University Medical Centre and 705 white controls from the same geographic region in The Netherlands ##REF##9189051##[73]##–##REF##7974338##[75]##. Both RF-positive and RF-negative patients were included in this sample set. ##TAB##0##Table 1## displays the clinical characteristics of all three sample sets and a detailed description of samples that overlap with published studies of this region ##UREF##1##[34]##,##REF##17880261##[35]##.</p> <title>Functional Genome-Wide Scan</title> <p>Our functional genome-wide scan included 25,966 gene-centric SNPs curated from dbSNP, the Applera Genome Initiative ##REF##16175505##[44]##,##REF##16237444##[76]## and the literature. SNPs were included if they appeared in more than one database and had a minor-allele frequency &gt;1%. Approximately seventy percent of the SNPs were annotated as missense polymorphisms. The majority of the remaining SNPs were either located within putative transcription-factor site motifs or within acceptor/donor splice site regions or were nonsense polymorphisms.</p> <title>Genotyping</title> <p>Allele-specific, real-time quantitative PCR ##REF##10673283##[77]## was used to amplify 3 ng of pooled DNAs and infer SNP allele frequencies as previously detailed ##REF##16175505##[44]##. Individual genotyping on SNPs was performed on 0.3 ng of DNA using a similar protocol. Blinded to case-control status, custom-made in-house software was used to call genotypes, followed by hand-curation. Individual genotyping accuracy has been estimated to be &gt;99.8% by comparison with an independent method. <italic>HLA-DRB1</italic> genotyping was performed using sequence-specific oligonucleotide probes as previously described ##REF##15208781##[9]##. Shared epitope (SE) status ##REF##2446635##[78]## was determined from the probe hybridization patterns. For this study, <italic>DRB1</italic> alleles positive for the SE include: 0101, 0102, 0401, 0404, 0405, 0408 and 1001.</p> <title>Fine-Scale Mapping SNP Selection</title> <p>To identify SNPs for inclusion in our fine-scale mapping effort of the 9q33.2 region, we first postulated two different disease models: 1) a model where the originally identified SNP is in linkage disequilibrium with one or more causative SNPs and 2) a model of allelic heterogeneity where several alleles at the locus independently predispose individuals to RA. To address both of these models, we first defined the region to be interrogated by calculating pairwise linkage disequilibrium (r<sup>2</sup>) values between the originally identified SNP 5′ of <italic>PHF19</italic>, rs1953126, and all HapMap-genotyped SNPs (<ext-link ext-link-type="uri" xlink:href="http://www.hapmap.org/">http://www.hapmap.org/</ext-link>) within 500 kb flanking either side for the CEPH samples (Utah residents with ancestry from northern and western Europe, or CEU individuals) ##REF##16255080##[36]##. With this information, we defined a broad region spanning 668 kb from <italic>MEGF9</italic>, 177 kb upstream of rs1953126 , to <italic>STOM</italic>, 491 kb downstream of rs1953126, for follow-up genotyping. SNPs within this region were partitioned into those in moderate to high LD (r<sup>2</sup>&gt;0.20) with rs1953126 to address the first model, and those in low LD (r<sup>2</sup>&lt;0.20) with rs1953126 to address the second model. The power-based SNP selection program Redigo ##REF##15297809##[79]## was then used on the low LD set of SNPs to identify a reduced number of SNPs (tagging SNP set) that retained high power to detect association. Those SNPs in moderate to high LD with the original SNP were reduced by selecting a subset of representative SNPs of any groups exhibiting extremely high inter-group LD (r<sup>2</sup>&gt;0.98). Further, any putative functional SNPs were automatically included in the fine-scale mapping effort if we were able to construct high-quality genotyping assays for them. The resulting set of 137 SNPs was genotyped in Sample Set 1 and the data analyzed. Additional removal of fine-scale mapping SNPs was performed for evaluation in subsequent sample sets on the basis of association results and refined LD patterns: a subset of 72 SNPs were selected for genotyping in Sample Set 2 and 42 SNPs were genotyped in Sample Set 3.</p> <title>Statistical Analyses</title> <p>The Cochran-Armitage trend test ##UREF##4##[80]## was used to calculate P-values for individual SNPs. A William's-corrected G-test ##UREF##2##[37]## was used to calculate P-values for genotypic association. P-values were corrected for multiple testing using the method of Dunn-Sidak ##UREF##2##[37]##. Odds ratios and confidence intervals were calculated according to standard procedures. Hardy-Weinberg equilibrium testing was accomplished through the exact test of Weir ##UREF##5##[81]##. P-values were combined across sample sets using the Fisher's combined P-value, or omnibus procedure ##UREF##6##[82]##. Likewise, Mantel-Haenszel common odds ratios ##REF##13655060##[83]## were calculated to combine data across sample sets. To avoid the small-count limitations of asymptotic-derived confidence intervals, a Monte Carlo simulation was written in XLISP-STAT to calculate 95% confidence intervals on the Mantel-Haenszel common odds ratios. We typically performed 20,000 iterations of the Monte Carlo for these confidence intervals. The standard measure of pairwise linkage disequilibrium (the r<sup>2</sup> statistic from estimated 2-site haplotypes) was used to characterize the genetic architecture of the region. The program LDMAX with an EM algorithm was used to perform the r<sup>2</sup> calculations ##REF##10842743##[84]##.</p> <title>Genetic Analyses</title> <title>Haplotype Analysis</title> <p>Haplotypes were estimated from unphased genotype data and evaluated for association with RA through the Haplo.Stats software package ##REF##11791212##[39]##. A sliding window of haplotype association was calculated using a window size of 5 SNPs. Global P-values (calculated across all haplotypes within a window) and haplotype-specific ORs and P-values were calculated. Additional haplotype analyses were performed using a combination of the Pseudo-Gibbs sampling algorithm in the program SNPAnalyzer (<ext-link ext-link-type="uri" xlink:href="http://snp.istech.info/snp/SNPAnalyzer.html">http://snp.istech.info/snp/SNPAnalyzer.html</ext-link>) ##REF##15980517##[40]## and the Haplo.Stats package.</p> <title>Genetic Background-Conditioned Analysis</title> <p>A panel of 749 SNPs previously selected to be informative for classifying individuals of European descent into northern and southern geographical groups was applied to case and control samples from the second sample set as described previously ##REF##17044734##[41]##. Applying this method, we placed 367 cases and 525 controls from Sample Set 2 into a northern European ancestry cluster. Each case or control individual had a greater than 0.95 probability of belonging to the northern European cluster. The remaining cases and controls from this study were binned into an “Other” category. A Breslow-Day analysis ##UREF##7##[85]## was applied to the stratified data to test for heterogeneity in ORs between the two groups for the 9q33.2-linked SNPs studied here. To test for association conditioned on these stratified data, we also calculated a Mantel-Haenszel P-value ##REF##13655060##[83]##.</p> <title>Subphenotype Analysis: Rheumatoid Factor</title> <p>Rheumatoid Factor (RF) levels were measured in cases as previously described ##REF##17195208##[86]##,##REF##8712877##[87]##. To test for heterogeneity of effect between RF-positive and RF-negative patients we used two different methods. In sample set 2, where case-control matching was part of the study design, we used the Breslow-Day ##UREF##7##[85]## test. Since individual matching was not incorporated into Sample Set 3, we used a Monte Carlo simulation to compare the effect size for RF-positive patients versus all controls to the effect size for RF-negative patients versus all controls. Similar to other tests of homogeneity of odds ratios, we constructed a test statistic measuring the departure between normalized odds ratios comparing two groups (see equation S1 in ##SUPPL##4##Text S1##) and ran a Monte Carlo simulation to account for correlated odds ratios in the null distribution. Monte Carlo P-values were calculated in the traditional manner.</p> <title>Logistic Regression</title> <p>Logistic regression models were performed to assess the relative importance of 27 SNPs chosen as distinct representatives of groups of SNPs with pairwise r<sup>2</sup> values &gt;.90. First, a logistic regression model for each unique pair of SNPs was performed. These pairwise models assumed a multiplicative effect on the risk of RA for each additional copy of an allele. P-values and odds ratios for the effect of each SNP, when controlling for each alternative SNP, were examined visually to determine if any SNP showed obvious patterns (attenuating the risk of each alternate SNP and retaining risk when adjusted for each alternate SNP). These types of patterns might be expected under a disease model of a single functional SNP. For models in which both SNPs remained strongly associated (P&lt;.01), additional models were performed to determine if adding a third SNP significantly improved the model. To examine multi-SNP relationships in a more automated fashion in a manner similar to that suggested by Cordell and Clayton ##REF##11719900##[88]##, both a forward as well as a backward stepwise logistic regression procedure was performed on each sample set individually as well as on the combined sample sets. The stepwise models were performed coding the genotypes with indicator variables and with a significance level of 0.05 for the two degree of freedom score test (for entry) or Wald test (for exit) on the effect of the SNP used as a threshold for entry or exit from the model. Models applied to the combined sample sets also forced sample set as a covariate in the model. The final model from each procedure was also applied to the other sample sets to assess consistency of the models across sample sets. The P-value from the likelihood ratio test of the global null hypothesis for each model is reported for single studies while for the combined study we report a P-value from a likelihood ratio test comparing the model containing the SNPs and the variable representing sample set to a model containing only the same sample set variable. All logistic regression models were performed using SAS version 9.</p> <title>Multi-Locus RA Risk Calculations</title> <p>Risk for RA given every possible 3-locus genotype combination at the <italic>HLA-DRB1</italic> shared epitope, the R620W SNP in <italic>PTPN22</italic>, and 3-SNP <italic>TRAF1</italic> diplotypes was calculated for sample set 1 using Bayes' theorem (see equations S2 and S3 in ##SUPPL##4##Text S1##) assuming conditional independence between loci (the commonly-used Naive Bayes model for predictive modeling) and a range of RA prevalence values (1%, 10% and 30%). Theoretical calculations (not shown) demonstrate that unless both sample sizes and epistatic effects are very large, probability estimates of the jointly-occurring genotypes have lower error rates assuming conditional independence between loci. Therefore, fully-factorizing the probability of multi-locus genotypes (using the conditional independence assumption) is warranted under a broad range of the parameter space. By estimating the posterior probability of RA for every possible multi-locus genotype combination, accurate individual-based prognosis is possible. Confidence intervals on the relative risk estimates were obtained through simulation. It is important to note that due to the selection of loci for inclusion in the model, some overfitting may be present.</p>
<title>Results</title> <title>Identification of the RA-Associated Chr 9q33.2 Region</title> <p>We are conducting three sequential case-control studies to identify SNPs associated with RA. In the first study, DNA samples from white North Americans with (N = 475 cases) and without (N = 475 controls) RA (sample set 1, see ##TAB##0##Table 1## for a breakdown of the clinical characteristics of each sample set) were genotyped for a set of 25,966 gene-centric SNPs utilizing disease-phenotype-based pooled DNA samples (pooled DNA samples were used to economically increase genotyping throughput while minimizing DNA consumption). The allele frequency of each SNP was determined in cases and controls as described in the <xref ref-type="sec" rid="s4">Methods</xref> and 1438 SNPs were significantly associated with RA using an allelic test (P&lt;0.05); 88 of these SNPs mapped to chr 6p21 between <italic>HLA-F</italic> and <italic>HLA-DPB1</italic> within the major histocompatibility complex (MHC). Of the 1350 non-MHC SNPs, 1306 were evaluated in a second independent white North American sample set (661 cases and 1322 controls) by use of a similar pooling strategy (44 SNPs were not genotyped due to insufficient primer quantities). Eighty-nine statistically compelling SNPs (P<sub>allelic</sub>&lt;0.05) with the same risk allele in these two sample sets were then individually genotyped in sample set 1 to verify the results from the pooled DNA phase of the experiment; 55 SNPs retained statistical significance (P<sub>allelic</sub>&lt;0.05) and 44 have been individually genotyped in sample set 2. Twenty-eight of these were significant (P<sub>allelic</sub>&lt;0.05) and are currently being evaluated in a third independent white Dutch sample set (596 cases and 705 controls).</p> <p>The most significant non-MHC SNP to emerge from a combined analysis of sample sets 1 and 2 after the <italic>PTPN22</italic> missense SNP, rs2476601 ##REF##15208781##[9]##, was rs1953126 an intergenic SNP located 1 kb upstream of the human homologue to the Drosophila polycomblike protein-encoding gene, <italic>PHF19</italic>, on chr 9q33.2 near two excellent candidate genes, <italic>TRAF1</italic> and <italic>C5</italic> (individual genotyping: Sample Set 1: OR = 1.30, 95% CI 1.08–1.58, trend P = 0.007; Sample Set 2: OR = 1.35, 95% CI 1.18–1.56, trend P = 1.69E-05). This SNP was genotyped in sample set 3 showing a nonsignificant trend towards association: OR = 1.16, 95% CI 0.99–1.36, trend P = 0.066 (##SUPPL##1##Table S1##). No significant deviations from Hardy-Weinberg equilibrium were observed for the genotypes of this SNP in the cases or controls in the three sample sets. The frequency of the minor allele was approximately 30.8% in white North American controls increasing to 37.3% in white North American cases and 34.9% in Dutch controls increasing to 38.3% in Dutch cases. A combined analysis across all three sample sets was highly significant (OR = 1.28, 95% CI 1.16–1.40, trend P<sub>comb</sub> = 1.45E-06).</p> <title>Chr 9q33.2 Fine-Mapping and LD Analyses</title> <p>To further explore the association signal in this region, we used patterns of LD from the CEU HapMap data (<ext-link ext-link-type="uri" xlink:href="http://www.hapmap.org">www.hapmap.org</ext-link>) ##REF##16255080##[36]## to define a broad 668 kb region, extending from <italic>MEGF9</italic> to <italic>STOM</italic> on chr 9q33.2, for follow-up individual genotyping. Postulating two different disease models, one where the originally identified SNP, rs1953126, is in LD with one or more causative SNPs and a second model of allelic heterogeneity where several alleles at a locus independently predispose individuals to disease, we selected a combination of 137 LD and tagging SNPs from this region for follow-up genotyping in Sample Set 1 (##FIG##0##Figure 1##) (A detailed description of SNP selection is outlined in the <xref ref-type="sec" rid="s4">Methods</xref>). Only four SNPs, all in the <italic>RAB14-GSN-STOM</italic> region, were mildly out of Hardy-Weinberg equilibrium (10<sup>−4</sup>&lt;P&lt;0.01) in the controls (##SUPPL##1##Table S1##). Including the original SNP, rs1953126, 38 of the 138 chr 9q33.2-region SNPs genotyped in Sample Set 1 were significant at the 0.01 level.</p> <p>To better understand these positive signals and select a subset of informative SNPs for genotyping in our other sample sets, we next investigated the LD architecture around rs1953126 by calculating pairwise r<sup>2</sup> values for all 138 SNPs genotyped in Sample Set 1. Evaluating cases and controls separately revealed very similar LD patterns across this region (##FIG##1##Figure 2A##). There were two primary haplotype blocks (LD Block 1 and LD Block 2) (here an LD block is defined as a region in which over 75% of all pairwise r<sup>2</sup> LD correlation values exceeded 0.3), with moderate LD between pairs of SNPs residing within each of the two blocks. LD Block 1, which contains the original SNP, rs1953126, and is approximately 70 kb, extends from rs10985070, an intronic SNP in the 5′ end of <italic>PHF19</italic>, across <italic>TRAF1</italic> into the <italic>TRAF1-C5</italic> intergenic region to rs2900180. Approximately 214 kb in length, LD Block 2 ranges from the middle of <italic>C5</italic> to the <italic>RAB14-GSN</italic> intergenic region. Given that haplotype block structures can have complex LD patterns within and between blocks and that we were focused on a single associated SNP in this region (rs1953126), we present a higher resolution plot shown in ##FIG##1##Figure 2B## where pairwise r<sup>2</sup> values were calculated for rs1953126 and each of the remaining 137 SNPs, revealing groups of highly correlated SNPs not readily visible in the LD heat-map.</p> <p>Integrating the Sample Set 1 association results with the LD measures, we found that the original SNP, rs1953126, was highly correlated (r<sup>2</sup>&gt;0.95) with 17 other SNPs (Group 1 in ##FIG##1##Figures 2Band 2C## in LD Block 1. As predicted, these 18 SNPs have similar association results increasing in frequency from approximately 30–31% in controls to 36–37% in cases (OR = 1.29–1.35, trend P∼0.002–0.009) (##SUPPL##1##Table S1##). Of interest was the observation that 20 non-Group 1 SNPs were associated with disease at equal or greater significance including 14 other SNPs from LD Block 1. Thirteen of these other LD Block 1 SNPs, which were highly correlated with one another (r<sup>2</sup>&gt;0.95) (Group 2 in ##FIG##1##Figures 2Band 2C## and reasonably correlated with the Group 1 SNPs (r<sup>2</sup> = 0.66–0.72), had minor allele frequencies of approximately 38% in controls increasing to 46% in cases (OR = 1.34–1.39, trend P≤0.002). The fourteenth significant SNP in LD Block 1, rs7021880, a <italic>TRAF1</italic> intronic SNP, was also highly significant (OR = 1.43, trend P = 3.12E-04) increasing in frequency from 27.1% in controls to 34.7% in cases. This SNP was in LD with both Group 1 (r<sup>2</sup> = 0.82–0.90) and Group 2 (r<sup>2</sup> = 0.59–0.64) SNPs (##FIG##1##Figure 2B##). The six other SNPs with P values &lt;0.01 lie upstream of LD Block 1 (n = 4) or downstream of LD Block 2 in <italic>GSN</italic> (n = 2) (##FIG##0##Figure 1##, ##SUPPL##1##Table S1##) and, with the exception of the <italic>PSMD5</italic> intronic SNP rs10760117, were not as significant as many of the LD Block 1 SNPs.</p> <p>Given the association results and the LD structure, we selected 72 of the 137 fine-scale mapping SNPs to genotype in Sample Set 2 (661 white North American RA patients and 1322 matched white North American controls) (##SUPPL##1##Table S1##). This subset of fine-scale mapping SNPs was chosen to reduce the genotyping load, while capturing the association signals and retaining full coverage of the genetic variation in this region. Two of these 72 SNPs, rs12683062 (in <italic>CEP110</italic>) in the cases and rs9409230 (a <italic>RAB14-GSN</italic> intergenic SNP) in the controls, were moderately out of Hardy-Weinberg equilibrium (P = 2.56E-04 and P = 0.003, respectively; ##SUPPL##1##Table S1##). Including the original SNP, rs1953126, 23 of these 72 SNPs were significant (trend P&lt;0.01) in Sample Set 2; however, the nine significant LD Block 1 SNPs in Sample Set 1 were the most significant, replicated SNPs in Sample Set 2 (##FIG##0##Figure 1##). Interestingly there were three SNPs in <italic>GSN</italic> (rs10985196, rs7046030 and rs12683459), all highly correlated with pairwise r<sup>2</sup> values &gt;0.90, which were highly significant (trend P&lt;10<sup>−6</sup>) in Sample Set 2 but only marginally significant in Sample Set 1 (trend P = 0.01–0.05). The difference between the two sample sets appears to be the result of disparate control allele frequencies – the case allele frequencies are nearly identical between the two sample sets (∼22%) but the control allele frequencies differ by 3% (18–19% in Sample Set 1 vs 15.5–16.5% in Sample Set 2) (##SUPPL##1##Table S1##).</p> <p>Forty-two SNPs were genotyped in Sample Set 3 (596 white Dutch RA patients and 705 white Dutch controls); none of these SNPs rejected HWE at the P&lt;0.01 significance level (##SUPPL##1##Table S1##). These 42 SNPs span over 600 kb and were selected to cover genetic variability, association patterns and gene boundaries. Four of the 42 SNPs, spanning 286 kb from <italic>TRAF1</italic> to <italic>RAB14</italic>, were significant at the 0.01 level (##FIG##0##Figure 1##). Of these four, two SNPs (rs4836834 and rs7021049) were members of Group 2 from LD Block 1, perfectly correlated (r<sup>2</sup> = 1) and both SNPs were highly significant in all three sample sets. The other two significant SNPs, rs1323472 and rs942152, were only moderately if at all significant in Sample Sets 1 and 2. The six Group 1 SNPs genotyped in Sample Set 3 were close to the 0.05 significance level, with the most significant of these being the synonymous P340P <italic>TRAF1</italic> SNP, rs2239657 and the <italic>TRAF1-C5</italic> intergenic SNP, rs2900180 (trend P = 0.052) (##SUPPL##1##Table S1##). The <italic>TRAF1</italic> intronic SNP, rs7021880, was not significant in this sample set (trend P = 0.102).</p> <p>In a combined analysis of the 43 SNPs genotyped in all three sample sets, including the original SNP, rs1953126, 25 SNPs, spanning a region of over 525 kb from rs7026635 within <italic>FBXW2</italic> to rs10818527 within <italic>GSN</italic>, were significantly associated with RA (trend P<sub>comb</sub>&lt;0.01) (##TAB##1##Table 2##). Several of these SNPs exhibited consistent and strong association across all three sample sets (##SUPPL##1##Table S1##). Using either a combined trend or genotypic P-value, the top-ranked five SNPs were: rs6478486, rs4836834, rs2239657, rs7021880 and rs7021049 (listed in order of position). All reside within or near <italic>TRAF1</italic> in LD Block 1, had common odds ratios of approximately 1.3 and were highly significant (trend P<sub>comb</sub>&lt;1.5E-07) (##TAB##1##Table 2##).</p> <title>Multiple Testing</title> <p>Since false-positive results can be problematic in any large-scale experiment in which modest nominal significance levels are used, we corrected the results from the combined analysis for multiple testing using the method of Dunn-Sidak ##UREF##2##[37]##. Seven SNPS, all within LD Block 1, survived a Dunn-Sidak correction for 25,966 SNPs at P&lt;0.01. The corrected trend P<sub>comb</sub> values for the five most significant SNPs were: 0.004 for rs6478486 and rs223957 (Group 1), 0.002 for rs4836834 and 0.001 for rs7021049 (Group 2), and 1.3E-04 for rs7021880.</p> <title>Haplotype Sliding Window</title> <p>Given that our fine-scale-mapping SNPs cluster into various groups based on their pairwise r<sup>2</sup> values and that under many models haplotypes can be more informative than single-markers ##REF##15368617##[38]##, we used the Haplo-Stats package ##REF##11791212##[39]## to run a 5-SNP sliding-window haplotype association analysis on the 43 SNPs genotyped in all three sample sets separately for each sample set and then combined the statistical evidence across all three sample sets. The combined analysis revealed a 29 kb-wide maximum peak of global association for haplotypes comprised of alleles segregating at rs6478486-rs4836834-rs2239657-rs7021880-rs7021049 in LD Block 1 (P<sub>comb</sub> = 4.15E-08) (##FIG##2##Figure 3##). This region ranges from 9 kb downstream of <italic>TRAF1</italic> in the <italic>PHF19-TRAF1</italic> intergenic region to intron 3 within <italic>TRAF1</italic>. Aside from this peak and a second highly significant peak in the <italic>TRAF1</italic> region (P<sub>comb</sub> = 5.45E-08; rs2239657-rs7021880-rs7021049-rs2900180-rs2269066), a second region of interest was centered over the <italic>RAB14-GSN</italic> region (P = 2.11E-06).</p> <p>Of these two regions, we view the disease association evidence to be stronger for the <italic>PHF19-TRAF1</italic> region for several reasons: First, combined analyses across all studies yielded the most significant results for both single markers and haplotypes in this region. Second, the association signal at this region shows a higher degree of consistency across the three studies. Indeed, Sample Set 3 haplotypes in the <italic>RAB14-GSN</italic> region show very little deviation from the null hypothesis (##FIG##2##Figure 3##). Finally, as discussed above, a subset of SNPs in the <italic>RAB14</italic>-<italic>GSN</italic> region (e.g. rs10985196, rs7046030, rs12683459) displayed substantial differences in control allele frequencies between the two North American groups (##SUPPL##1##Table S1##) drawing into question the validity of the association results for these SNPs.</p> <title>Haplotype Analyses of LD Block 1 Variants</title> <p>While both the single marker and sliding window haplotype analyses pointed to LD Block 1 as harboring RA-associated SNPs, these analyses did not identify a single SNP that was clearly the most significant across all three studies. The <italic>TRAF1</italic> intronic SNP, rs7021880, was the most significant SNP in Sample Sets 1 (trend P = 3.12E-04) and 2 (trend P = 5.09E-07) and in the combined analysis (trend P<sub>comb</sub> = 5.41E-09); however, this SNP was not significant in the Dutch sample set (trend P = 0.102) where the Group 2 SNPs, rs4836834 and rs7021049, were the most significant (trend P = 0.004 and 0.006, respectively) (##SUPPL##1##Tables S1## and ##SUPPL##2##S2##). Interestingly, these Group 2 SNPs ranked second in significance in Sample Set 1 and in the combined analysis while in Sample Set 2 they ranked third behind rs7021880 and the Group 1 SNPs.</p> <p>Given these results, we analyzed the haplotype structure of LD Block 1 using a subset of the nine SNPs from this region genotyped in all three studies. Taking into account the LD structure we picked rs2239657, the P340P <italic>TRAF1</italic> synonymous polymorphism to represent the six Group 1 SNPs; rs7021049, a <italic>TRAF1</italic> intronic SNP to represent the two Group 2 SNPs; and rs7021880 for these analyses. Haplotype frequencies for these three SNPs were estimated using the Haplo.Stats package ##REF##11791212##[39]## revealing the same four common haplotypes in each study (##TAB##2##Table 3##). Two of these haplotypes, AGT and GCG were strongly associated with disease (P<sub>comb</sub> =  3.08E-08 and 8.00E-09, respectively), with the former being protective – decreasing in frequency from ∼60.9% in North American controls to 53.8% in North American cases and 56.7% in Dutch controls to 51.2% in Dutch cases (OR<sub>common</sub> = 0.76, 95% CI 0.70–0.83); and the latter susceptible – 27.0% in North American controls increasing to 34.7% in North American cases and 33.2% in Dutch controls increasing to 36.0% in Dutch cases (OR<sub>common</sub> = 1.32, 95% CI 1.21–1.45). These haplotype P<sub>comb</sub>-values were not significantly different from those calculated for the individual SNPs (##TAB##1##Table 2##) suggesting there is no strong evidence for synergistic cis-acting effects between these variants.</p> <title>Dosage Effects</title> <p>To explore the effect of the number of copies of each haplotype at these three sites (rs2239657, rs7021880 and rs7021049) along with any dominant/recessive effects between haplotypes, we estimated diplotypes using the pseudo-Gibbs sampling algorithm from the program SNPAnalyzer ##REF##15980517##[40]##. Analyzing the diplotypes individually, two diplotype combinations achieved statistical significance (P&lt;0.01) when compared to all other diplotypes (##TAB##3##Table 4##). The AGT/AGT diplotype was strongly associated with protection against RA (OR<sub>Common</sub> = 0.68, 95%CI 0.59–0.78; P<sub>Comb</sub> = 5.35E-07), whereas the less frequent GCG/GCG diplotype was associated with predisposition (OR<sub>Common</sub> = 1.42, 95%CI 1.16–1.75; P<sub>Comb</sub> = 0.005).</p> <p>Assuming a disease prevalence of 1%, we calculated the relative risk of RA in those individuals carrying 2 copies of the protective AGT haplotype compared to those without the AGT haplotype (RR<sub>2 copies AGT</sub> = 0.77). This homozygous relative risk was substanftially reduced from the relative risk calculated for individuals carrying only one copy of the AGT haplotype (RR<sub>1 copy AGT</sub> = 1.06.). Similarly, we estimated the relative risks for the susceptible GCG haplotype (RR<sub>2 copies GCG</sub> = 1.38; RR<sub>1 copy GCG</sub> = 1.15).</p> <title>Genetic Background-Conditioned Results</title> <p>We used a collection of 749 SNPs informative for European substructure to stratify both the cases and controls in Sample Set 2 ##REF##17044734##[41]##. By partitioning cases and controls into similar genetic background groups (“Northern European” or “Other”), our aim was to interrogate the data for strata-specific effects – that is, whether or not association signals were specific to one of these genetic background groups – and avoid potential confounding by population stratification. Although two SNPs demonstrated moderately higher significance levels following stratification – rs16910233 in <italic>C5</italic> (P<sub>North</sub> = 0.019 compared to P<sub>Unstrat</sub> = 0.147) and rs12685539 in <italic>CEP110</italic> (P<sub>Other</sub> = 0.038 compared to P<sub>Unstrat</sub> = 0.115), a Breslow-Day test of effect heterogeneity comparing OR<sub>North</sub> and OR<sub>Other</sub> was not significant. Furthermore, a positional plot of Mantel-Haenszel P-values, testing for association given the genetic background stratification, was very similar to the unadjusted plot (##SUPPL##0##Figure S1##) suggesting that stratification of the case and control samples by SNPs informative for European substructure did not change the association patterns in Sample Set 2.</p> <title>Rheumatoid Factor (RF)</title> <p>Rheumatoid factor, a circulating antibody to immunoglobulin G, is a key serum analyte used in diagnosis of RA as well as an aid for the prognosis of RA-severity ##REF##12748655##[2]##. As the R620W missense polymorphism in <italic>PTPN22</italic> appears to have stronger susceptibility effects for RF-positive disease ##REF##15208781##[9]##–##REF##16731003##[11]## and since RF is clinically important, we investigated the role of RF status on the chr 9q33.2 association patterns for the three LD Block 1 SNPs, rs2239657, rs7021880 and rs7021049, testing for both strata-specific effects as well effect size differences between RF-positive and RF-negative disease.</p> <p>To explore the effect isolated to RF-positive patients compared to controls, we performed a strata-specific analysis for all sample sets using a genotypic test. The resulting combined P-values for the RF-positive stratum were highly significant (P<sub>rs2239657</sub> = 4.02E-05, P<sub>rs7021880</sub> = 7.10E-06, P<sub>rs7021049</sub> = 5.68E-06; ##TAB##4##Table 5##), which were slightly less significant when compared to the overall genotypic combined P-values (##TAB##1##Table 2##). A similar analysis of RF-negative disease in Sample Sets 2 and 3 yielded genotypic combined P-values of P<sub>rs2239657</sub> = 0.038, P<sub>rs7021880</sub> = 0.013 and P<sub>rs7021049</sub> = 0.082. Allelic odds ratios and 95% confidence intervals were also calculated for each individual sample set and the results did not demonstrate a clear pattern of strata-specific effects within a stratum or differential effects between the two strata (##TAB##4##Table 5##). A Breslow-Day test was performed on Sample Set 2 (individually matched cases and controls) to formalize the test of homogeneity of odds ratios, showing that none of the three SNPs exhibited significant differential effects (##TAB##4##Table 5##). Similarly, results for the analogous Monte Carlo-based test performed in Sample Set 3 (where cases and controls were not individually matched) also did not reveal significant heterogeneity between RF-positive and RF-negative effects.</p> <title>Logistic Regression</title> <p>To further dissect association signals from LD patterns, build predictive models and explore the relative effects of each SNP within the models constructed we used logistic regression. To accomplish this we first minimized the number of SNPs for these analyses by calculating pairwise r<sup>2</sup> values for the 43 SNPs genotyped in all three sample sets and divided the SNPs into distinct groups based on their LD structure. SNPs with pairwise r<sup>2</sup> values &gt;0.90 were grouped together resulting in 27 distinct groups (##TAB##5##Table 6##) and then the single most significant SNP from each group (P<sub>comb</sub> from ##TAB##1##Table 2##) was chosen for the logistic regression analyses.</p> <p>In the univariate analysis, the <italic>TRAF1</italic> intronic SNP rs7021049, which marks the Group 2 SNPs in LD Block 1, was the most significant SNP (P = 1.24E-06), followed by rs7021880 (1.39E-06) and then the Group 1 SNP, rs2239657 (P = 2.52E-06) (##TAB##5##Table 6##). In addition, 11 other SNPs were significant (P&lt;0.01). To assess whether other observed associations in the region were primarily a result of LD with the most significant SNP, we performed pairwise logistic regression on all 27 SNPs adjusting for rs7021049. One SNP retained modest statistical significance (P&lt;0.01): rs10985196 (Group 21), a <italic>GSN</italic> intronic SNP (P = 0.001). To test whether the combination of the rs7021049 and rs10985196 variants fully accounted for the association with RA, we repeated the logistic regression adjusting for both; none of the remaining groups of SNPs were significantly associated with RA. It should be noted, however, that analyses of each individual sample set suggested the evidence for association with rs10985196 (Group 21) was primarily driven by the data from sample set 2 (data not presented).</p> <p>To explore more complex models, we used both forwards and backwards stepwise logistic regression procedures separately on the same 27 SNPs in each individual sample set as well as in a combined analysis of all three sample sets. The final models generated from the stepwise procedures, however, were inconsistent across the sample sets (##SUPPL##2##Table S2##). In fact, seven distinct models were produced; the only instance where the same model was produced was for both the forwards and backwards models of Sample Set 2. Not surprisingly, the forward model for the combined samples, which included two SNPs, rs7021049 (the Group 2 <italic>TRAF1</italic> intronic SNP) and rs10985196 (the <italic>GSN</italic> intronic SNP), was consistent with the results of the pairwise logistic regression analysis presented above.</p> <title>Multi-Locus RA Risk Calculations</title> <p>Given that we have begun to witness the application of associated genetic variants to disease prognosis ##REF##15778708##[42]##,##REF##16936732##[43]## and thus far we have convincing evidence for three RA-predisposing loci in our studies (<italic>HLA-DRB1</italic>, <italic>PTPN22</italic> and the <italic>TRAF1</italic> region), we estimated the risk of RA given genotypes at these three loci under three different possible unconditional RA risk assumptions (i.e. RA disease prevalence values) using Bayes' theorem. In total, there were 18 multi-locus genotype combinations and RA risk was calculated for each combination using data from Sample Set 1 as described in the <xref ref-type="sec" rid="s4">Materials and Methods</xref>. Assuming a 1% RA prevalence, similar to that observed in the white North American general population, the results indicate that individuals with the protective genotype at all three loci (0SE for <italic>HLA-DRB1</italic>, CC genotype for <italic>PTPN22</italic> and the AGT/AGT <italic>TRAF1</italic> diplotype) have a substantially reduced predicted risk of RA (0.29% vs. 1%), whereas those individuals in the highest-risk category (<italic>HLA</italic>-2SE, TT or TC genotype at <italic>PTPN22</italic>, and the GCG/GCG <italic>TRAF1</italic> diplotype), have an estimated RA risk of 13.06% – representing more than a 45-fold increase in risk (##SUPPL##3##Table S3##). These data are presented as a 3-D plot in ##FIG##3##Figure 4## where the lowest risk value has been reset to 1 and the other values normalized accordingly. Approximately 19% of the general population will find themselves in the low-risk multi-locus genotype category and only 0.06% in the high risk group. In contrast, when the disease prevalence is increased to 30%, as might be observed in high-risk groups such as an early arthritis clinic, the range of risk drops to 7.88-fold, with the posterior probability of RA calculated to be 11% for the lowest-risk genotype combination and increasing to 86.4% in the highest risk category (##SUPPL##3##Table S3##).</p>
<title>Discussion</title> <p>We undertook a large-scale, multi-tiered association study of RA using a panel of putative functional SNPs that have been successfully applied to case-control studies of other disease phenotypes ##REF##17236132##[32]##, ##REF##16175505##[44]##–##REF##17461418##[47]##. The initial step of this large-scale RA association study, individual genotyping of 87 prioritized SNPs to evaluate DNA quality prior to constructing DNA pools for our scan, led to the identification of the <italic>PTPN22</italic> R620W SNP ##REF##15208781##[9]##. This SNP has been both widely replicated and associated with multiple autoimmune diseases ##REF##16731003##[11]##.</p> <p>The present study focuses on variants in the chr 9q33.2 region that were also convincingly correlated with RA status. In particular, three groups of SNPs, represented by rs2239657 (Group 1), rs7021049 (Group 2) and rs7021880 were highly significant and showed a localized effect to a 70 kb region extending from rs10985070, in intron 3 of <italic>PHF19</italic>, across <italic>TRAF1</italic> to rs2900180 in the <italic>TRAF1-C5</italic> intergenic region, but excluding the <italic>C5</italic> coding region (LD Block 1 in ##FIG##0##Figure 1##). Examination of the CEU HapMap data identified 16 additional SNPs that were highly correlated (r<sup>2</sup>&gt;0.95) with either the Group 1 or Group 2 SNPs genotyped in this study (##FIG##1##Figure 2C##) and all 16 fall within this 70 kb region (no such SNPs were found for rs7021880). Across sample sets, the evidence for association at these sites was stronger, maintaining statistical significance after correction for multiple testing, and more consistent than sites in neighboring regions. Additional analyses further buttressed the statistical support for these conclusions: (i) a haplotype sliding window analysis of all SNPs genotyped in the chr 9q33.2 region demonstrated strong statistical evidence for the <italic>TRAF1</italic>-region harboring RA risk variants (P<sub>comb</sub> = 4.15E-08) and (ii) haplotype analysis of SNPs within the 70 kb LD Block 1, identified a common protective haplotype (P<sub>comb</sub> = 3.08E-08) and a less frequent risk haplotype (P<sub>comb</sub> = 8.00E-09). The three representative SNPs (rs2239657, rs7021049 and rs7021880) were strongly associated with RF-positive disease and trended towards association in RF-negative disease although the small number of RF-negative patients in our study precludes drawing statistically meaningful conclusions about the role of these SNPs in this patient population.</p> <p>To tease apart association signals from LD patterns, we used logistic regression. The pairwise analyses of the combined datasets suggest there may be two independent statistical signals of association to RA at chr 9q33.2 – one in the <italic>TRAF1</italic> region represented by rs7021049 and one in the <italic>GSN</italic> region represented by rs10985196 (##TAB##5##Table 6##); however, analyses of the individual sample sets showed rs10985196 was independently associated with disease risk in Sample Set 2 only while rs7021049 showed consistent association across all three sample sets (data not presented). Consequently, additional samples are needed to determine whether the <italic>GSN</italic> region truly contains RA-predisposing effects.</p> <p>To explore more complex models and assess whether SNPs outside of LD Block 1 were incorporated into these models, we used both forwards and backwards stepwise logistic regression. The sets of SNPs included in the models chosen by the stepwise procedures were inconsistent indicating that the observed association in the region is not clearly explained by a single SNP or set of SNPs included in the tested models.</p> <p>Independently, Plenge and colleagues ##UREF##1##[34]##, using a whole genome association (WGA) study, and Kurreeman and coworkers ##REF##17880261##[35]##, using a candidate gene approach, have also shown this chr 9q33.2 region is associated with RA risk in whites of European descent. Although a partially overlapping subset of samples was used in all three studies (see ##TAB##0##Table 1## footnotes), each study employed unique experimental designs, analyses and presented different facets of the 9q33.2-RA association. Plenge and colleagues ##UREF##1##[34]## identified rs3761847, a <italic>TRAF1/C5</italic> intergenic SNP, as one of two non-MHC SNPs reaching genome-wide significance in their WGA study; not surprisingly, the other significant non-MHC SNP was in <italic>PTPN22</italic>. Their follow-up fine mapping of the chr 9q33.2 region with nine haplotype tag SNPs in four RA sample sets (2519 cases / 3627 controls) localized the region of interest to 100 kb extending from <italic>PHF19</italic> into <italic>C5</italic>. rs3761847, which is a Group 2 SNP in LD Block 1, remained the most significant SNP in their combined analysis (P = 4.00E-14) followed by rs2900180 (P = 8.00E-14), a member of the Group 1 SNPs in LD Block 1.</p> <p>Taking a candidate gene approach, Kurreeman and colleagues ##REF##17880261##[35]## studied 40 SNPs in a 300 kb region surrounding <italic>C5</italic> (from the 3′ UTR of <italic>PHF19</italic> to intron 25 of <italic>CEP110</italic>) in a staged-approach in four sample sets (2,719 cases / 1999 controls) and concluded that rs10818488, another <italic>TRAF1-C5</italic> intergenic SNP and member of the LD Block 1 Group 2 SNPs, was the SNP most significantly associated with a diagnosis of RA in their study. Association of a Group 1 SNP, rs2416806, was moderately significant in a combined analysis of three of their sample sets (P = 0.015). Neither the Plenge et al. nor Kurreeman et al. study included rs7021880.</p> <p>Our analyses, which included more SNPs and incorporated HapMap information for all SNPs highly correlated with SNPs genotyped in our study, permitted a comprehensive analysis of the genetic architecture of 9q33.2 region, allowing us to localize the RA-susceptibility effects to a 70 kb region (LD Block 1) that includes a portion of <italic>PHF19</italic>, all of <italic>TRAF1</italic> and the majority of the <italic>TRAF1-C5</italic> intergenic region, but excludes the <italic>C5</italic> coding region, narrowing the true region of interest. Our data, however, did not allow us to identify a single SNP or group of highly correlated SNPs (r<sup>2</sup>&gt;0.95) in this 70 kb region that unambiguously explained the association signal in all three independent sample sets. Other sample sets with different patterns of LD or functional studies will be required to resolve this issue.</p> <p>Interestingly, Potter and colleagues ##REF##17277003##[48]##, who studied 23 haplotype-tagging SNPs from the 6 <italic>TRAF</italic> genes, including three from <italic>TRAF1</italic>, in a UK case-control study (351 RA cases / 368 controls) failed to see association with both a Group1 (rs1468671) and a Group2 (rs4836834) SNP. Using an overlapping sample set to the Potter et al study, the recent Welcome Trust Case Control Consortium genome-wide association study of RA (1860 cases / 2930 controls) ##REF##17554300##[25]## also failed to identify RA-risk variants in this region. However, a more recent follow-up study from the same group of an independent RA sample set from the UK (3418 cases / 3337 controls) confirmed association with four LD Block 1 Group 2 SNPs although the effect size was less (meta analysis OR 1.08, 95%CI 1.03–1.14) ##UREF##3##[49]##.</p> <p>The original RA-associated, 9q33.2 SNP identified in our genome-wide scan, rs1953126, is located within LD Block 1, 1 kb upstream of the 5′ end of <italic>PHF19</italic>, the human homologue of the Drosophila polycomblike protein, <italic>PCL</italic>, gene. In Drosophila, the protein encoded by this gene is part of the 1MDa extra sex combs and enhancer of zeste [ESC-E(Z)] complex which is thought to mediate transcriptional repression by modulating the chromatin environment of many developmental regulatory genes such as homeobox genes. While the exact function of this gene in humans remains unknown, it encodes two nuclear proteins that appear to be upregulated in multiple cancers and preliminary evidence suggests that deregulation of these genes may play a role in tumor progression ##REF##15563832##[50]##.</p> <p>
<italic>TRAF1</italic> encodes a protein that is a member of the TNF receptor (TNFR) associated factor (TRAF) protein family that associates with, and mediates signal transduction from various receptors including a subset of the TNFR superfamily. There are six members of this family of adaptor proteins; however, TRAF1 is unique in that while it contains the hallmark carboxyl-terminal TRAF domain, it has a single zinc finger in the amino-terminal part and the N-terminal RING finger domain, required for NF-κB activation, is missing. TRAF1 appears to have both anti-apoptotic and anti-proliferative effects ##REF##9733516##[51]##,##REF##11672546##[52]##. In addition, this protein has been found to be elevated in malignancies of the B cell lineage ##REF##9885224##[53]##–##REF##12933571##[58]##. This observation is interesting given the risk of lymphoma, particularly diffuse large B cell lymphomas, appears to be increased in the subset of RA patients with very severe disease, independent of treatment ##REF##16508929##[59]##,##REF##18501101##[60]##. Although the precise mechanism of TRAF1 action in various signaling pathways has not been fully elucidated, it is clear that this molecule plays an important role in immune cell homeostasis making it an excellent candidate gene for RA. In fact, <italic>in vitro</italic> work suggests that TNFα-mediated synovial hyperplasia, a major pathophysiologic feature of RA, may be correlated with upregulation of TRAF molecules, particularly TRAF1 ##REF##12067756##[61]##. Given that TNF blockade has proved a highly effective therapy for RA ##REF##7934491##[62]##,##REF##9920948##[63]## and response to TNF-antagonists among RA patients is known to vary, investigation of whether the <italic>TRAF1</italic> variants identified in this study play a role in this differential response may be a fruitful pharmacogenetic avenue to pursue.</p> <p>
<italic>C5</italic> is also an excellent RA candidate gene and although our analyses allowed us to exclude the <italic>C5</italic> coding region, SNPs in LD Block 1 could differentially regulate the expression of this gene. <italic>C5</italic> encodes a zymogen that is involved in all three pathways of complement activation. Traditionally, the complement system has been viewed as a central part of the innate immune system in host defenses against invading pathogens and in clearance of potentially damaging cell debris; however, complement activation has also recently been implicated in the pathogenesis of many inflammatory and immunological diseases. Proteolytic cleavage of C5 results in C5a, one the most potent inflammatory peptides, and C5b, a component of the membrane attack complex (MAC) that can cause lysis of cells and bacteria. Genetic studies in various mouse models of RA, including collagen-induced arthritis (CIA) and the K/BXN T cell receptor transgenic mouse model of inflammatory arthritis have provided evidence that <italic>C5</italic> or a variant in strong LD, plays a role in disease ##REF##10698345##[64]##–##REF##11489951##[66]##. More striking is the observation that anti-C5 monoclonal antibody therapy can prevent and ameliorate disease in both mouse models ##REF##7568051##[67]##,##REF##11869678##[68]##.</p> <p>In summary, we have independently identified a region on chr 9q33.2 as a risk locus for RA. Although the evidence from the SNPs genotyped in our sample sets most strongly points towards <italic>TRAF1</italic> variants as being the most highly consistent with a disease model, the high LD that extends from the 5′ end of <italic>PHF19</italic> through <italic>TRAF1</italic> and into the <italic>TRAF1-C5</italic> intergenic region precludes conclusively determining causative genes or functional motifs through genetic means in these samples. Mapping studies in additional sample sets with a different LD architecture and/or functional studies will be required to resolve the molecular relevance of these findings.</p> <p>Aside from the possibility of developing targeted therapies with knowledge of predisposing variants underlying the onset of RA, the identification of RA susceptibility alleles may encourage earlier monitoring and provide an intervention avenue in advance of significant joint erosion. Our initial analysis of the three known genetic risk factors, <italic>HLA-SE</italic>, <italic>PTPN22</italic> and the chr 9q33.2 variants described here, suggests a &gt;45 fold difference in RA risk depending on an individual's genotype at these three loci. As additional markers are identified, the ability to accurately predict individuals at increased risk for developing RA, particularly within families with a history of RA, may prove useful. Finally, differential risk variants may prove to be drug response markers.</p>
<p><bold>¤a:</bold> Current address: Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America</p> <p><bold>¤b:</bold> Current address: Roche Diagnostics, Pleasanton, California, United States of America</p> <p>Conceived and designed the experiments: SS AB. Performed the experiments: MC JC. Analyzed the data: MC CR VG SS AB. Contributed reagents/materials/analysis tools: CR VG SS AvdH KA CA LC DK PG FK RT TH MS. Wrote the paper: CR SS AB.</p> <p>Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease affecting both joints and extra-articular tissues. Although some genetic risk factors for RA are well-established, most notably <italic>HLA-DRB1</italic> and <italic>PTPN22</italic>, these markers do not fully account for the observed heritability. To identify additional susceptibility loci, we carried out a multi-tiered, case-control association study, genotyping 25,966 putative functional SNPs in 475 white North American RA patients and 475 matched controls. Significant markers were genotyped in two additional, independent, white case-control sample sets (661 cases/1322 controls from North America and 596 cases/705 controls from The Netherlands) identifying a SNP, rs1953126, on chromosome 9q33.2 that was significantly associated with RA (OR<sub>common</sub> = 1.28, trend P<sub>comb</sub> = 1.45E-06). Through a comprehensive fine-scale-mapping SNP-selection procedure, 137 additional SNPs in a 668 kb region from <italic>MEGF9</italic> to <italic>STOM</italic> on 9q33.2 were chosen for follow-up genotyping in a staged-approach. Significant single marker results (P<sub>comb</sub>&lt;0.01) spanned a large 525 kb region from <italic>FBXW2</italic> to <italic>GSN</italic>. However, a variety of analyses identified SNPs in a 70 kb region extending from the third intron of <italic>PHF19</italic> across <italic>TRAF1</italic> into the <italic>TRAF1-C5</italic> intergenic region, but excluding the <italic>C5</italic> coding region, as the most interesting (trend P<sub>comb</sub>: 1.45E-06 → 5.41E-09). The observed association patterns for these SNPs had heightened statistical significance and a higher degree of consistency across sample sets. In addition, the allele frequencies for these SNPs displayed reduced variability between control groups when compared to other SNPs. Lastly, in combination with the other two known genetic risk factors, <italic>HLA-DRB1</italic> and <italic>PTPN22</italic>, the variants reported here generate more than a 45-fold RA-risk differential.</p> <title>Author Summary</title> <p>Rheumatoid arthritis (RA), a chronic autoimmune disorder affecting ∼1% of the population, is characterized by immune-cell–mediated destruction of the joint architecture. Gene–environment interactions are thought to underlie RA etiology. Variants within <italic>HLA-DRB1</italic> and the hematopoietic-specific phosphatase, <italic>PTPN22</italic>, are well established RA-susceptibility loci, and although other markers have been identified, they do not fully account for the disease heritability. To identify additional susceptibility alleles, we carried out a multi-tiered, case-control association study genotyping &gt;25,000 putative functional SNPs; here we report our finding of RA-associated variants in chromosome 9q33.2. A detailed genetic analysis of this region, incorporating HapMap information, localizes the RA-susceptibility effects to a 70 kb region that includes a portion of <italic>PHF19</italic>, all of <italic>TRAF1</italic>, and the majority of the <italic>TRAF1-C5</italic> intergenic region, but excludes the <italic>C5</italic> coding region. In addition to providing new insights into underlying mechanism(s) of disease and suggesting novel therapeutic targets, these data provide the underpinnings of a genetic signature that may predict individuals at increased risk for developing RA. Indeed, initial analyses of three known genetic risk factors, <italic>HLA</italic>, <italic>PTPN22</italic>, and the chromosome 9q33.2 variants described here, suggest a &gt;45-fold difference in RA risk depending on an individual's three-locus genotype.</p>
<title>Supporting Information</title>
<p>We are grateful to all the individuals with rheumatoid arthritis, collaborating physicians and to control individuals for their participation in this study; M. Schrodi for helpful comments on the manuscript; T.L. Bugawan and E. Trachtenberg for <italic>HLA</italic> typing; R. Lundsten, M. Kern, H. Khalili and A. Rodrigues-Brown for database and sample management of the NARAC collection; and J. Lemaire and S. Mahan for database and sample management at SeraCare. We would like to thank the Celera High Throughput and Computational Biology departments for their support; Y. Li for numerous invaluable discussions on mapping genetics and helpful comments on this manuscript; J. Martin and R. Venkatesh for useful discussions on machine learning methods; and T. White and J. Sninsky for constructive comments on the manuscript, support and scientific advice.</p>
<fig id="pgen-1000107-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.g001</object-id><label>Figure 1</label><caption><title>Case-control association results and linkage disequilibrium structure of the 9q33.2 region.</title><p>A physical map of the 668 kb surrounding the original associated SNP, rs1953126, with the location of all 138 markers genotyped in sample set 1 noted. The markers in red indicate the 43 SNPs genotyped in all three sample sets. The locations of LD Block 1 and LD Block 2 are indicated. Above the physical map, the trend P-values are displayed for the SNPs genotyped in each of the three sample sets. The red line indicates trend P = 0.01. The LD structure across the 668 kb region from <italic>MEGF9</italic> to <italic>STOM</italic>, based on pairwise D' values from the CEU HapMap, is displayed below the physical map.</p></caption></fig> <fig id="pgen-1000107-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.g002</object-id><label>Figure 2</label><caption><title>The LD architecture of the 9q33.2 region.</title><p>(A) Pairwise linkage disequilibrium values (r<sup>2</sup>) for all 138 SNPs genotyped in Sample Set 1. Cases and controls are shown separately. (B) Pairwise LD values between rs1953126 and each of the 137 other SNPs genotyped in Sample Set 1. Locations of the two main LD blocks are shown in bold. (C) SNPs within each of the two LD groups in Block 1. SNPs in black were genotyped in this study and are listed according to their position. SNPs in grey were not genotyped but highly correlated (r<sup>2</sup>&gt;0.93) with either Group 1 or Group 2 SNPs in the CEU HapMap data. Note that all of the SNPs in grey lie in LD Block 1 between rs10985070 and rs2900180.</p></caption></fig> <fig id="pgen-1000107-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.g003</object-id><label>Figure 3</label><caption><title>A five-SNP sliding window haplotype analysis of the 9q33.2 region.</title><p>Each sample set is shown separately with the combined analysis in bolded black. The approximate location of the <italic>PHF19</italic>, <italic>TRAF1</italic>, <italic>RAB14</italic> and <italic>GSN</italic> genes are listed above.</p></caption></fig> <fig id="pgen-1000107-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.g004</object-id><label>Figure 4</label><caption><title>Relative risk plotted as a function of the genetic load of three validated RA risk variants in <italic>HLA-DRB1</italic>, <italic>PTPN22</italic> and <italic>TRAF1</italic>.</title><p>Individuals are classified according to the number of copies of the <italic>HLA-DRB1</italic> shared epitope (0, 1 and 2) (SE-positive <italic>HLA-DRB1</italic> alleles found in this sample set were: 0101, 0102, 0401, 0404, 0405, 0408 and 1001), carriage of the W620 <italic>PTPN22</italic> missense SNP (TT + CT vs CC) and diplotypes at the TRAF1 SNPs, rs2239657, rs7021880 and rs7021049. The frequency of each combination of markers in cases and controls is highlighted in white (case/control).</p></caption></fig>
<table-wrap id="pgen-1000107-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t001</object-id><label>Table 1</label><caption><title>Demographic and clinical information.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Sample Set</td></tr><tr><td align="left" rowspan="1" colspan="1">Subphenotype</td><td align="left" rowspan="1" colspan="1">1<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">2<xref ref-type="table-fn" rid="nt102">b</xref>
</td><td align="left" rowspan="1" colspan="1">3<xref ref-type="table-fn" rid="nt103">c</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Genetic background</td><td align="left" rowspan="1" colspan="1">White (North American)</td><td align="left" rowspan="1" colspan="1">White (North American)</td><td align="left" rowspan="1" colspan="1">White (Dutch)</td></tr><tr><td align="left" rowspan="1" colspan="1">No. of cases</td><td align="left" rowspan="1" colspan="1">475</td><td align="left" rowspan="1" colspan="1">661</td><td align="left" rowspan="1" colspan="1">596</td></tr><tr><td align="left" rowspan="1" colspan="1">No. of controls</td><td align="left" rowspan="1" colspan="1">475</td><td align="left" rowspan="1" colspan="1">1322</td><td align="left" rowspan="1" colspan="1">705</td></tr><tr><td align="left" rowspan="1" colspan="1">Female:male</td><td align="left" rowspan="1" colspan="1">314:161</td><td align="left" rowspan="1" colspan="1">536:125</td><td align="left" rowspan="1" colspan="1">362:196<xref ref-type="table-fn" rid="nt104">d</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Average age of onset (years)</td><td align="left" rowspan="1" colspan="1">46.97±11.83</td><td align="left" rowspan="1" colspan="1">38.61±13.61</td><td align="left" rowspan="1" colspan="1">54.58±13.38<xref ref-type="table-fn" rid="nt105">e</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">% RF-positive</td><td align="left" rowspan="1" colspan="1">100%</td><td align="left" rowspan="1" colspan="1">82%</td><td align="left" rowspan="1" colspan="1">72%<xref ref-type="table-fn" rid="nt106">f</xref>
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000107-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t002</object-id><label>Table 2</label><caption><title>Combined analysis of 43 chr 9q33.2 SNPs genotyped in all three RA sample sets.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Combined Analysis</td></tr><tr><td align="left" rowspan="1" colspan="1">Marker</td><td align="left" rowspan="1" colspan="1">Gene</td><td align="left" rowspan="1" colspan="1">Type</td><td align="left" rowspan="1" colspan="1">Position &amp; Alleles<xref ref-type="table-fn" rid="nt107">a</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>common</sub> (95% CI)<xref ref-type="table-fn" rid="nt108">b</xref>
</td><td align="left" rowspan="1" colspan="1">Trend P<sub>comb</sub>
<xref ref-type="table-fn" rid="nt109">c</xref>
</td><td align="left" rowspan="1" colspan="1">Genotypic P<sub>comb</sub>
<xref ref-type="table-fn" rid="nt109">c</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">rs10760112</td><td align="left" rowspan="1" colspan="1">
<italic>MEGF9</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122507391T</td><td align="left" rowspan="1" colspan="1">1.17 (1.02–1.23)</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">0.136</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7026635</td><td align="left" rowspan="1" colspan="1">
<italic>FBXW2</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">G122589848A</td><td align="left" rowspan="1" colspan="1">1.24 (1.10–1.35)</td><td align="left" rowspan="1" colspan="1">0.001</td><td align="left" rowspan="1" colspan="1">0.012</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10760117</td><td align="left" rowspan="1" colspan="1">
<italic>PSMD5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T122626558G</td><td align="left" rowspan="1" colspan="1">1.26 (1.10–1.31)</td><td align="left" rowspan="1" colspan="1">2.79E-04</td><td align="left" rowspan="1" colspan="1">0.003</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10739575</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">G122645922A</td><td align="left" rowspan="1" colspan="1">1.16 (1.03–1.30)</td><td align="left" rowspan="1" colspan="1">0.081</td><td align="left" rowspan="1" colspan="1">0.349</td></tr><tr><td align="left" rowspan="1" colspan="1">rs933003</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">A122647650G</td><td align="left" rowspan="1" colspan="1">1.12 (0.79–1.40)</td><td align="left" rowspan="1" colspan="1">0.255</td><td align="left" rowspan="1" colspan="1">0.243</td></tr><tr><td align="left" rowspan="1" colspan="1">rs1837</td><td align="left" rowspan="1" colspan="1">
<italic>PHF19</italic>
</td><td align="left" rowspan="1" colspan="1">3′UTR</td><td align="left" rowspan="1" colspan="1">T122658050C</td><td align="left" rowspan="1" colspan="1">1.28 (1.12–1.36)</td><td align="left" rowspan="1" colspan="1">2.17E-04</td><td align="left" rowspan="1" colspan="1">0.002</td></tr><tr><td align="left" rowspan="1" colspan="1">rs1056567</td><td align="left" rowspan="1" colspan="1">
<italic>PHF19</italic>
</td><td align="left" rowspan="1" colspan="1">S181S</td><td align="left" rowspan="1" colspan="1">T122671866C</td><td align="left" rowspan="1" colspan="1">1.25 (1.12–1.35)</td><td align="left" rowspan="1" colspan="1">1.11E-04</td><td align="left" rowspan="1" colspan="1">0.002</td></tr><tr><td align="left" rowspan="1" colspan="1">rs1953126</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T122680321C</td><td align="left" rowspan="1" colspan="1">1.28 (1.16–1.40)</td><td align="left" rowspan="1" colspan="1">1.45E-06</td><td align="left" rowspan="1" colspan="1">4.24E-05</td></tr><tr><td align="left" rowspan="1" colspan="1">rs1609810</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">G122682172A</td><td align="left" rowspan="1" colspan="1">1.29 (1.19–1.42)</td><td align="left" rowspan="1" colspan="1">1.92E-07</td><td align="left" rowspan="1" colspan="1">5.24E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs881375</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T122692719C</td><td align="left" rowspan="1" colspan="1">1.27 (1.17–1.41)</td><td align="left" rowspan="1" colspan="1">4.69E-07</td><td align="left" rowspan="1" colspan="1">1.09E-05</td></tr><tr><td align="left" rowspan="1" colspan="1">rs6478486</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T122695150C</td><td align="left" rowspan="1" colspan="1">1.29 (1.19–1.42)</td><td align="left" rowspan="1" colspan="1">1.35E-07</td><td align="left" rowspan="1" colspan="1">3.75E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs4836834</td><td align="left" rowspan="1" colspan="1">
<italic>TRAF1</italic>
</td><td align="left" rowspan="1" colspan="1">3′UTR</td><td align="left" rowspan="1" colspan="1">T122705722A</td><td align="left" rowspan="1" colspan="1">1.32 (1.19–1.43)</td><td align="left" rowspan="1" colspan="1">8.13E-08</td><td align="left" rowspan="1" colspan="1">1.84E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2239657</td><td align="left" rowspan="1" colspan="1">
<italic>TRAF1</italic>
</td><td align="left" rowspan="1" colspan="1">P340P</td><td align="left" rowspan="1" colspan="1">G122711341A</td><td align="left" rowspan="1" colspan="1">1.29 (1.19–1.43)</td><td align="left" rowspan="1" colspan="1">1.49E-07</td><td align="left" rowspan="1" colspan="1">3.89E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7021880</td><td align="left" rowspan="1" colspan="1">
<italic>TRAF1</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122713711G</td><td align="left" rowspan="1" colspan="1">1.33(1.21–1.46)</td><td align="left" rowspan="1" colspan="1">5.41E-09</td><td align="left" rowspan="1" colspan="1">2.27E-07</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7021049</td><td align="left" rowspan="1" colspan="1">
<italic>TRAF1</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">G122723803T</td><td align="left" rowspan="1" colspan="1">1.32 (1.20–1.43)</td><td align="left" rowspan="1" colspan="1">4.09E-08</td><td align="left" rowspan="1" colspan="1">1.22E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2900180</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T122746203C</td><td align="left" rowspan="1" colspan="1">1.27 (1.18–1.41)</td><td align="left" rowspan="1" colspan="1">3.32E-07</td><td align="left" rowspan="1" colspan="1">7.62E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2269066</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T122776839C</td><td align="left" rowspan="1" colspan="1">1.29 (1.14–1.53)</td><td align="left" rowspan="1" colspan="1">1.68E-04</td><td align="left" rowspan="1" colspan="1">0.001</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2269067</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122776861G</td><td align="left" rowspan="1" colspan="1">1.27 (1.17–1.46)</td><td align="left" rowspan="1" colspan="1">1.71E-05</td><td align="left" rowspan="1" colspan="1">1.04E-04</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2159776</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122795981T</td><td align="left" rowspan="1" colspan="1">1.11 (0.99–1.19)</td><td align="left" rowspan="1" colspan="1">0.190</td><td align="left" rowspan="1" colspan="1">0.135</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7040033</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A122798865G</td><td align="left" rowspan="1" colspan="1">0.86 (0.80–0.96)</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">0.060</td></tr><tr><td align="left" rowspan="1" colspan="1">rs17611</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">I802V</td><td align="left" rowspan="1" colspan="1">A122809021G</td><td align="left" rowspan="1" colspan="1">0.84 (0.79–0.94)</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">0.040</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10985126</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">G385G</td><td align="left" rowspan="1" colspan="1">C122823755T</td><td align="left" rowspan="1" colspan="1">1.20 (1.11–1.39)</td><td align="left" rowspan="1" colspan="1">8.69E-04</td><td align="left" rowspan="1" colspan="1">0.001</td></tr><tr><td align="left" rowspan="1" colspan="1">rs2416811</td><td align="left" rowspan="1" colspan="1">
<italic>C5</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T122829455C</td><td align="left" rowspan="1" colspan="1">0.85 (0.79–0.95)</td><td align="left" rowspan="1" colspan="1">0.008</td><td align="left" rowspan="1" colspan="1">0.023</td></tr><tr><td align="left" rowspan="1" colspan="1">rs1323472</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">C122866156G</td><td align="left" rowspan="1" colspan="1">1.23 (1.12–1.34)</td><td align="left" rowspan="1" colspan="1">1.57E-04</td><td align="left" rowspan="1" colspan="1">7.06E-04</td></tr><tr><td align="left" rowspan="1" colspan="1">rs9657673</td><td align="left" rowspan="1" colspan="1">
<italic>CEP110</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T122900196C</td><td align="left" rowspan="1" colspan="1">0.86 (0.81–0.96)</td><td align="left" rowspan="1" colspan="1">0.019</td><td align="left" rowspan="1" colspan="1">0.052</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10081760</td><td align="left" rowspan="1" colspan="1">
<italic>CEP110</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A122924127G</td><td align="left" rowspan="1" colspan="1">1.15 (1.03–1.25)</td><td align="left" rowspan="1" colspan="1">0.049</td><td align="left" rowspan="1" colspan="1">0.066</td></tr><tr><td align="left" rowspan="1" colspan="1">rs12683062</td><td align="left" rowspan="1" colspan="1">
<italic>CEP110</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T122946625G</td><td align="left" rowspan="1" colspan="1">1.12 (1.00–1.33)</td><td align="left" rowspan="1" colspan="1">0.209</td><td align="left" rowspan="1" colspan="1">0.029</td></tr><tr><td align="left" rowspan="1" colspan="1">rs3747843</td><td align="left" rowspan="1" colspan="1">
<italic>CEP110</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A122954127G</td><td align="left" rowspan="1" colspan="1">1.13 (1.01–1.21)</td><td align="left" rowspan="1" colspan="1">0.108</td><td align="left" rowspan="1" colspan="1">0.304</td></tr><tr><td align="left" rowspan="1" colspan="1">rs3736855</td><td align="left" rowspan="1" colspan="1">
<italic>CEP110</italic>
</td><td align="left" rowspan="1" colspan="1">V1398V</td><td align="left" rowspan="1" colspan="1">A122956841T</td><td align="left" rowspan="1" colspan="1">0.87 (0.82–0.98)</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">0.191</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10760152</td><td align="left" rowspan="1" colspan="1">
<italic>RAB14</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A122987806C</td><td align="left" rowspan="1" colspan="1">1.15 (1.05–1.27)</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">0.024</td></tr><tr><td align="left" rowspan="1" colspan="1">rs942152</td><td align="left" rowspan="1" colspan="1">
<italic>RAB14</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122991506T</td><td align="left" rowspan="1" colspan="1">1.18 (1.11–1.32)</td><td align="left" rowspan="1" colspan="1">2.53E-04</td><td align="left" rowspan="1" colspan="1">0.002</td></tr><tr><td align="left" rowspan="1" colspan="1">rs9408928</td><td align="left" rowspan="1" colspan="1">
<italic>RAB14</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C122991738T</td><td align="left" rowspan="1" colspan="1">1.11 (0.93–1.38)</td><td align="left" rowspan="1" colspan="1">0.364</td><td align="left" rowspan="1" colspan="1">0.378</td></tr><tr><td align="left" rowspan="1" colspan="1">rs9409230</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T123007581A</td><td align="left" rowspan="1" colspan="1">1.14 (0.93–1.40)</td><td align="left" rowspan="1" colspan="1">0.499</td><td align="left" rowspan="1" colspan="1">0.217</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7030849</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">C123009655T</td><td align="left" rowspan="1" colspan="1">1.18 (1.08–1.29)</td><td align="left" rowspan="1" colspan="1">0.003</td><td align="left" rowspan="1" colspan="1">0.014</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10985196</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A123072865C</td><td align="left" rowspan="1" colspan="1">1.25 (1.18–1.46)</td><td align="left" rowspan="1" colspan="1">6.33E-07</td><td align="left" rowspan="1" colspan="1">4.12E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs306781</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C123082765T</td><td align="left" rowspan="1" colspan="1">0.68 (0.59–1.16)</td><td align="left" rowspan="1" colspan="1">0.119</td><td align="left" rowspan="1" colspan="1">0.284</td></tr><tr><td align="left" rowspan="1" colspan="1">rs7046030</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">C123087058T</td><td align="left" rowspan="1" colspan="1">1.26 (1.18–1.47)</td><td align="left" rowspan="1" colspan="1">2.05E-06</td><td align="left" rowspan="1" colspan="1">1.99E-05</td></tr><tr><td align="left" rowspan="1" colspan="1">rs12683459</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A123088119G</td><td align="left" rowspan="1" colspan="1">1.25 (1.18–1.47)</td><td align="left" rowspan="1" colspan="1">1.36E-06</td><td align="left" rowspan="1" colspan="1">9.79E-06</td></tr><tr><td align="left" rowspan="1" colspan="1">rs4837839</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T123111948C</td><td align="left" rowspan="1" colspan="1">0.85 (0.82–0.97)</td><td align="left" rowspan="1" colspan="1">0.021</td><td align="left" rowspan="1" colspan="1">0.076</td></tr><tr><td align="left" rowspan="1" colspan="1">rs306783</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T123112418C</td><td align="left" rowspan="1" colspan="1">1.11 (1.00–1.19)</td><td align="left" rowspan="1" colspan="1">0.198</td><td align="left" rowspan="1" colspan="1">0.405</td></tr><tr><td align="left" rowspan="1" colspan="1">rs306784</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T123112473G</td><td align="left" rowspan="1" colspan="1">1.15 (1.03–1.24)</td><td align="left" rowspan="1" colspan="1">0.049</td><td align="left" rowspan="1" colspan="1">0.131</td></tr><tr><td align="left" rowspan="1" colspan="1">rs10818527</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">A123115075G</td><td align="left" rowspan="1" colspan="1">1.21 (1.08–1.31)</td><td align="left" rowspan="1" colspan="1">0.001</td><td align="left" rowspan="1" colspan="1">0.004</td></tr><tr><td align="left" rowspan="1" colspan="1">rs12683989</td><td align="left" rowspan="1" colspan="1">
<italic>GSN</italic>
</td><td align="left" rowspan="1" colspan="1">intronic</td><td align="left" rowspan="1" colspan="1">T123125867C</td><td align="left" rowspan="1" colspan="1">1.17 (1.05–1.50)</td><td align="left" rowspan="1" colspan="1">0.016</td><td align="left" rowspan="1" colspan="1">0.010</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000107-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t003</object-id><label>Table 3</label><caption><title>Three SNP haplotypes for LD Block 1.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="4" align="left" rowspan="1">Sample Set 1</td><td colspan="4" align="left" rowspan="1">Sample Set 2</td><td colspan="4" align="left" rowspan="1">Sample Set 3</td><td colspan="2" align="left" rowspan="1">Combined</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="4" align="left" rowspan="1">Global P = 6.00E-04<xref ref-type="table-fn" rid="nt110">a</xref>
</td><td colspan="4" align="left" rowspan="1">Global P = 3.77E-05<xref ref-type="table-fn" rid="nt110">a</xref>
</td><td colspan="4" align="left" rowspan="1">Global P = 0.033<xref ref-type="table-fn" rid="nt110">a</xref>
</td><td colspan="2" align="left" rowspan="1">Global P<sub>comb</sub> = 1.81E-07<xref ref-type="table-fn" rid="nt111">b</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Haplotype<xref ref-type="table-fn" rid="nt112">c</xref>
</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">P<sub>comb</sub>
<xref ref-type="table-fn" rid="nt113">d</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>common</sub> (95% CI)<xref ref-type="table-fn" rid="nt114">e</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">AGT</td><td align="left" rowspan="1" colspan="1">507(0.539)</td><td align="left" rowspan="1" colspan="1">582(0.619)</td><td align="left" rowspan="1" colspan="1">5.08E-04</td><td align="left" rowspan="1" colspan="1">0.72</td><td align="left" rowspan="1" colspan="1">708(0.537)</td><td align="left" rowspan="1" colspan="1">1595(0.605)</td><td align="left" rowspan="1" colspan="1">4.01E-05</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">604(0.512)</td><td align="left" rowspan="1" colspan="1">794(0.567)</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.80</td><td align="left" rowspan="1" colspan="1">3.08E-08</td><td align="left" rowspan="1" colspan="1">0.76 (0.70–0.83)</td></tr><tr><td align="left" rowspan="1" colspan="1">GCG</td><td align="left" rowspan="1" colspan="1">326(0.347)</td><td align="left" rowspan="1" colspan="1">253(0.269)</td><td align="left" rowspan="1" colspan="1">2.13E-04</td><td align="left" rowspan="1" colspan="1">1.44</td><td align="left" rowspan="1" colspan="1">457(0.347)</td><td align="left" rowspan="1" colspan="1">714(0.271)</td><td align="left" rowspan="1" colspan="1">8.71E-07</td><td align="left" rowspan="1" colspan="1">1.43</td><td align="left" rowspan="1" colspan="1">425(0.360)</td><td align="left" rowspan="1" colspan="1">465(0.332)</td><td align="left" rowspan="1" colspan="1">0.133</td><td align="left" rowspan="1" colspan="1">1.13</td><td align="left" rowspan="1" colspan="1">8.00E-09</td><td align="left" rowspan="1" colspan="1">1.32 (1.21–1.45)</td></tr><tr><td align="left" rowspan="1" colspan="1">AGG</td><td align="left" rowspan="1" colspan="1">85(0.090)</td><td align="left" rowspan="1" colspan="1">71(0.075)</td><td align="left" rowspan="1" colspan="1">0.250</td><td align="left" rowspan="1" colspan="1">1.22</td><td align="left" rowspan="1" colspan="1">108(0.082)</td><td align="left" rowspan="1" colspan="1">232(0.088)</td><td align="left" rowspan="1" colspan="1">0.540</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">122(0.103)</td><td align="left" rowspan="1" colspan="1">120(0.086)</td><td align="left" rowspan="1" colspan="1">0.127</td><td align="left" rowspan="1" colspan="1">1.22</td><td align="left" rowspan="1" colspan="1">NC</td><td align="left" rowspan="1" colspan="1">1.09 (0.93–1.27)</td></tr><tr><td align="left" rowspan="1" colspan="1">GGG</td><td align="left" rowspan="1" colspan="1">22(0.023)</td><td align="left" rowspan="1" colspan="1">32(0.034)</td><td align="left" rowspan="1" colspan="1">0.168</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">41(0.031)</td><td align="left" rowspan="1" colspan="1">92(0.035)</td><td align="left" rowspan="1" colspan="1">0.539</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">25(0.021)</td><td align="left" rowspan="1" colspan="1">20(0.014)</td><td align="left" rowspan="1" colspan="1">0.135</td><td align="left" rowspan="1" colspan="1">1.49</td><td align="left" rowspan="1" colspan="1">NC</td><td align="left" rowspan="1" colspan="1">0.93 (0.70–1.21)</td></tr><tr><td align="left" rowspan="1" colspan="1">Other</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2(0.002)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4(0.003)</td><td align="left" rowspan="1" colspan="1">3(0.001)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">5(0.004)</td><td align="left" rowspan="1" colspan="1">1(0.001)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000107-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t004</object-id><label>Table 4</label><caption><title>Diplotype Analysis for the TRAF1-region SNPs rs2239657, rs7021880 and rs7021049.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="4" align="left" rowspan="1">Sample Set 1 Global P = 0.0069<xref ref-type="table-fn" rid="nt115">a</xref>
</td><td colspan="4" align="left" rowspan="1">Sample Set 2 Global P = 1.3E-04<xref ref-type="table-fn" rid="nt115">a</xref>
</td><td colspan="4" align="left" rowspan="1">Sample Set 3 Global P = 0.058<xref ref-type="table-fn" rid="nt115">a</xref>
</td><td colspan="2" align="left" rowspan="1">Combined Analysis Global P<sub>comb</sub> = 8.22E-06<xref ref-type="table-fn" rid="nt116">b</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">No. (Frequency) in</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Diplotype<xref ref-type="table-fn" rid="nt117">c</xref>
</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt118">d</xref>
</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt118">d</xref>
</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">Case</td><td align="left" rowspan="1" colspan="1">Control</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt118">d</xref>
</td><td align="left" rowspan="1" colspan="1">OR</td><td align="left" rowspan="1" colspan="1">P<sub>comb</sub>
<xref ref-type="table-fn" rid="nt119">e</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>common</sub>
<xref ref-type="table-fn" rid="nt120">f</xref> (95%CI)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">AGT/AGT</td><td align="left" rowspan="1" colspan="1">140(0.297)</td><td align="left" rowspan="1" colspan="1">180(0.383)</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">183(0.278)</td><td align="left" rowspan="1" colspan="1">482(0.366)</td><td align="left" rowspan="1" colspan="1">8.21E-05</td><td align="left" rowspan="1" colspan="1">0.67</td><td align="left" rowspan="1" colspan="1">153(0.259)</td><td align="left" rowspan="1" colspan="1">232(0.331)</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">5.35E-07</td><td align="left" rowspan="1" colspan="1">0.68(0.59–0.78)</td></tr><tr><td align="left" rowspan="1" colspan="1">AGT/Other</td><td align="left" rowspan="1" colspan="1">51(0.108)</td><td align="left" rowspan="1" colspan="1">64(0.136)</td><td align="left" rowspan="1" colspan="1">0.197</td><td align="left" rowspan="1" colspan="1">0.77</td><td align="left" rowspan="1" colspan="1">86(0.131)</td><td align="left" rowspan="1" colspan="1">204(0.155)</td><td align="left" rowspan="1" colspan="1">0.157</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">81(0.137)</td><td align="left" rowspan="1" colspan="1">73(0.104)</td><td align="left" rowspan="1" colspan="1">0.085</td><td align="left" rowspan="1" colspan="1">1.36</td><td align="left" rowspan="1" colspan="1">NC</td><td align="left" rowspan="1" colspan="1">0.94(0.78–1.13)</td></tr><tr><td align="left" rowspan="1" colspan="1">AGT/GCG</td><td align="left" rowspan="1" colspan="1">178(0.377)</td><td align="left" rowspan="1" colspan="1">158(0.336)</td><td align="left" rowspan="1" colspan="1">0.197</td><td align="left" rowspan="1" colspan="1">1.20</td><td align="left" rowspan="1" colspan="1">255(0.387)</td><td align="left" rowspan="1" colspan="1">426(0.324)</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">1.32</td><td align="left" rowspan="1" colspan="1">218(0.369)</td><td align="left" rowspan="1" colspan="1">257(0.367)</td><td align="left" rowspan="1" colspan="1">0.954</td><td align="left" rowspan="1" colspan="1">1.01</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">1.18(1.04–1.34)</td></tr><tr><td align="left" rowspan="1" colspan="1">GCG/GCG</td><td align="left" rowspan="1" colspan="1">51(0.108)</td><td align="left" rowspan="1" colspan="1">34(0.072)</td><td align="left" rowspan="1" colspan="1">0.068</td><td align="left" rowspan="1" colspan="1">1.55</td><td align="left" rowspan="1" colspan="1">76(0.115)</td><td align="left" rowspan="1" colspan="1">100(0.076)</td><td align="left" rowspan="1" colspan="1">0.004</td><td align="left" rowspan="1" colspan="1">1.59</td><td align="left" rowspan="1" colspan="1">77(0.130)</td><td align="left" rowspan="1" colspan="1">78(0.111)</td><td align="left" rowspan="1" colspan="1">0.304</td><td align="left" rowspan="1" colspan="1">1.19</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">1.42(1.16–1.75)</td></tr><tr><td align="left" rowspan="1" colspan="1">GCG/Other</td><td align="left" rowspan="1" colspan="1">46(0.098)</td><td align="left" rowspan="1" colspan="1">27(0.057)</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">1.77</td><td align="left" rowspan="1" colspan="1">50(0.076)</td><td align="left" rowspan="1" colspan="1">87(0.066)</td><td align="left" rowspan="1" colspan="1">0.452</td><td align="left" rowspan="1" colspan="1">1.16</td><td align="left" rowspan="1" colspan="1">54(0.091)</td><td align="left" rowspan="1" colspan="1">52(0.074)</td><td align="left" rowspan="1" colspan="1">0.309</td><td align="left" rowspan="1" colspan="1">1.25</td><td align="left" rowspan="1" colspan="1">0.086</td><td align="left" rowspan="1" colspan="1">1.32(1.04–1.66)</td></tr><tr><td align="left" rowspan="1" colspan="1">Other/Other</td><td align="left" rowspan="1" colspan="1">6(0.013)</td><td align="left" rowspan="1" colspan="1">7(0.015)</td><td align="left" rowspan="1" colspan="1">0.789</td><td align="left" rowspan="1" colspan="1">0.85</td><td align="left" rowspan="1" colspan="1">9(0.014)</td><td align="left" rowspan="1" colspan="1">18(0.014)</td><td align="left" rowspan="1" colspan="1">1.000</td><td align="left" rowspan="1" colspan="1">1.00</td><td align="left" rowspan="1" colspan="1">8(0.014)</td><td align="left" rowspan="1" colspan="1">8(0.011)</td><td align="left" rowspan="1" colspan="1">0.804</td><td align="left" rowspan="1" colspan="1">1.19</td><td align="left" rowspan="1" colspan="1">NC</td><td align="left" rowspan="1" colspan="1">1.01(0.56–1.72)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000107-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t005</object-id><label>Table 5</label><caption><title>Analysis of rs2239657, rs7021880 and rs7021049 stratified by the presence of rheumatoid factor.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="6" align="left" rowspan="1">rs2239657</td><td colspan="6" align="left" rowspan="1">rs7021880</td><td colspan="6" align="left" rowspan="1">rs7021049</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Genotypes</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Genotypes</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Genotypes</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">GG</td><td align="left" rowspan="1" colspan="1">GA</td><td align="left" rowspan="1" colspan="1">AA</td><td align="left" rowspan="1" colspan="1">MAF</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt121">a</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>Allelic</sub>
</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">CG</td><td align="left" rowspan="1" colspan="1">GG</td><td align="left" rowspan="1" colspan="1">MAF</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt121">a</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>Allelic</sub>
</td><td align="left" rowspan="1" colspan="1">GG</td><td align="left" rowspan="1" colspan="1">GT</td><td align="left" rowspan="1" colspan="1">TT</td><td align="left" rowspan="1" colspan="1">MAF</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt121">a</xref>
</td><td align="left" rowspan="1" colspan="1">OR<sub>Allelic</sub>
</td></tr></thead><tbody><tr><td colspan="19" align="left" rowspan="1">
<bold>Sample Set 1</bold>
<xref ref-type="table-fn" rid="nt122">b</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">RF-positive cases</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">224</td><td align="left" rowspan="1" colspan="1">184</td><td align="left" rowspan="1" colspan="1">0.370</td><td align="left" rowspan="1" colspan="1">0.008</td><td align="left" rowspan="1" colspan="1">1.35 (1.11–1.63)</td><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">225</td><td align="left" rowspan="1" colspan="1">195</td><td align="left" rowspan="1" colspan="1">0.347</td><td align="left" rowspan="1" colspan="1">0.001</td><td align="left" rowspan="1" colspan="1">1.43 (1.17–1.74)</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">229</td><td align="left" rowspan="1" colspan="1">140</td><td align="left" rowspan="1" colspan="1">0.461</td><td align="left" rowspan="1" colspan="1">0.002</td><td align="left" rowspan="1" colspan="1">1.39 (1.16–1.67)</td></tr><tr><td align="left" rowspan="1" colspan="1">Controls</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">195</td><td align="left" rowspan="1" colspan="1">229</td><td align="left" rowspan="1" colspan="1">0.304</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">187</td><td align="left" rowspan="1" colspan="1">249</td><td align="left" rowspan="1" colspan="1">0.271</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">222</td><td align="left" rowspan="1" colspan="1">180</td><td align="left" rowspan="1" colspan="1">0.381</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td colspan="19" align="left" rowspan="1">
<bold>Sample Set 2</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">RF-positive cases</td><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">268</td><td align="left" rowspan="1" colspan="1">206</td><td align="left" rowspan="1" colspan="1">0.373</td><td align="left" rowspan="1" colspan="1">5.60E-04</td><td align="left" rowspan="1" colspan="1">1.32 (1.14–1.55)</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">250</td><td align="left" rowspan="1" colspan="1">229</td><td align="left" rowspan="1" colspan="1">0.346</td><td align="left" rowspan="1" colspan="1">2.39E-04</td><td align="left" rowspan="1" colspan="1">1.27(1.08–1.49)</td><td align="left" rowspan="1" colspan="1">107</td><td align="left" rowspan="1" colspan="1">283</td><td align="left" rowspan="1" colspan="1">152</td><td align="left" rowspan="1" colspan="1">0.458</td><td align="left" rowspan="1" colspan="1">8.77E-04</td><td align="left" rowspan="1" colspan="1">1.30 (1.12–1.50)</td></tr><tr><td align="left" rowspan="1" colspan="1">Matched controls</td><td align="left" rowspan="1" colspan="1">106</td><td align="left" rowspan="1" colspan="1">457</td><td align="left" rowspan="1" colspan="1">520</td><td align="left" rowspan="1" colspan="1">0.309</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">87</td><td align="left" rowspan="1" colspan="1">425</td><td align="left" rowspan="1" colspan="1">571</td><td align="left" rowspan="1" colspan="1">0.277</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">175</td><td align="left" rowspan="1" colspan="1">505</td><td align="left" rowspan="1" colspan="1">403</td><td align="left" rowspan="1" colspan="1">0.395</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">RF-negative cases</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">57</td><td align="left" rowspan="1" colspan="1">41</td><td align="left" rowspan="1" colspan="1">0.406</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">1.63 (1.18–2.27)</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">56</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">0.368</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">1.74(1.24–2.44)</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">59</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">0.474</td><td align="left" rowspan="1" colspan="1">0.054</td><td align="left" rowspan="1" colspan="1">1.41 (1.02–1.93)</td></tr><tr><td align="left" rowspan="1" colspan="1">Matched controls</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">115</td><td align="left" rowspan="1" colspan="1">0.295</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">130</td><td align="left" rowspan="1" colspan="1">0.250</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">125</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">0.391</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Breslow-Day<xref ref-type="table-fn" rid="nt123">c</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.263</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.222</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.656</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td colspan="19" align="left" rowspan="1">
<bold>Sample Set 3</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">RF-positive cases</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">156</td><td align="left" rowspan="1" colspan="1">111</td><td align="left" rowspan="1" colspan="1">0.398</td><td align="left" rowspan="1" colspan="1">0.07</td><td align="left" rowspan="1" colspan="1">1.25 (1.03–1.52)</td><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">151</td><td align="left" rowspan="1" colspan="1">121</td><td align="left" rowspan="1" colspan="1">0.374</td><td align="left" rowspan="1" colspan="1">0.184</td><td align="left" rowspan="1" colspan="1">1.19(0.98–1.45)</td><td align="left" rowspan="1" colspan="1">73</td><td align="left" rowspan="1" colspan="1">164</td><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">0.494</td><td align="left" rowspan="1" colspan="1">0.019</td><td align="left" rowspan="1" colspan="1">1.28 (1.06–1.55)</td></tr><tr><td align="left" rowspan="1" colspan="1">RF-negative cases</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">0.364</td><td align="left" rowspan="1" colspan="1">0.483</td><td align="left" rowspan="1" colspan="1">1.09 (0.82–1.44)</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">0.357</td><td align="left" rowspan="1" colspan="1">0.312</td><td align="left" rowspan="1" colspan="1">1.11(0.83–1.47)</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">0.447</td><td align="left" rowspan="1" colspan="1">0.297</td><td align="left" rowspan="1" colspan="1">1.06 (0.81–1.39)</td></tr><tr><td align="left" rowspan="1" colspan="1">Controls</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">320</td><td align="left" rowspan="1" colspan="1">298</td><td align="left" rowspan="1" colspan="1">0.346</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">79</td><td align="left" rowspan="1" colspan="1">309</td><td align="left" rowspan="1" colspan="1">312</td><td align="left" rowspan="1" colspan="1">0.334</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">137</td><td align="left" rowspan="1" colspan="1">331</td><td align="left" rowspan="1" colspan="1">232</td><td align="left" rowspan="1" colspan="1">0.432</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Monte Carlo</bold>
<sup>d</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.218</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.645</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.116</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td colspan="19" align="left" rowspan="1">
<bold>Combined</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">RF-positive cases<xref ref-type="table-fn" rid="nt125">e</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4.02E-05</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">7.10E-06</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">5.68E-06</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">RF-negative cases <xref ref-type="table-fn" rid="nt126">f</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.082</td><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000107-t006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000107.t006</object-id><label>Table 6</label><caption><title>Pairwise logistic regression analysis of 27 chr 9q33.2 SNPs.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Group<xref ref-type="table-fn" rid="nt127">a</xref>
</td><td align="left" rowspan="1" colspan="1">Marker</td><td align="left" rowspan="1" colspan="1">r<sup>2</sup> with rs7021049<xref ref-type="table-fn" rid="nt128">b</xref>
</td><td align="left" rowspan="1" colspan="1">P<xref ref-type="table-fn" rid="nt129">c</xref>
</td><td align="left" rowspan="1" colspan="1">P adjusted for rs7021049</td><td align="left" rowspan="1" colspan="1">P adjusted for rs7021049 &amp; rs10985196</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">rs10760112</td><td align="left" rowspan="1" colspan="1">0.157</td><td align="left" rowspan="1" colspan="1">0.357</td><td align="left" rowspan="1" colspan="1">0.285</td><td align="left" rowspan="1" colspan="1">0.770</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">rs10760117</td><td align="left" rowspan="1" colspan="1">0.329</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.760</td><td align="left" rowspan="1" colspan="1">0.579</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">rs10739575</td><td align="left" rowspan="1" colspan="1">0.086</td><td align="left" rowspan="1" colspan="1">0.055</td><td align="left" rowspan="1" colspan="1">0.580</td><td align="left" rowspan="1" colspan="1">0.893</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">rs933003</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.757</td><td align="left" rowspan="1" colspan="1">0.420</td><td align="left" rowspan="1" colspan="1">0.448</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">rs1837, rs7026635</td><td align="left" rowspan="1" colspan="1">0.151</td><td align="left" rowspan="1" colspan="1">0.002</td><td align="left" rowspan="1" colspan="1">0.169</td><td align="left" rowspan="1" colspan="1">0.126</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">rs1056567</td><td align="left" rowspan="1" colspan="1">0.243</td><td align="left" rowspan="1" colspan="1">5.22E-04</td><td align="left" rowspan="1" colspan="1">0.200</td><td align="left" rowspan="1" colspan="1">0.208</td></tr><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">rs2239657, rs1953126, rs1609810, rs881375, rs6478486, rs2900180</td><td align="left" rowspan="1" colspan="1">0.685</td><td align="left" rowspan="1" colspan="1">2.52E-06</td><td align="left" rowspan="1" colspan="1">0.217</td><td align="left" rowspan="1" colspan="1">0.254</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">rs7021880</td><td align="left" rowspan="1" colspan="1">0.607</td><td align="left" rowspan="1" colspan="1">1.39E-06</td><td align="left" rowspan="1" colspan="1">0.104</td><td align="left" rowspan="1" colspan="1">0.072</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">rs7021049, rs4836834</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1.24E-06</td><td align="left" rowspan="1" colspan="1">–</td><td align="left" rowspan="1" colspan="1">–</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">rs2269066</td><td align="left" rowspan="1" colspan="1">0.114</td><td align="left" rowspan="1" colspan="1">0.002</td><td align="left" rowspan="1" colspan="1">0.115</td><td align="left" rowspan="1" colspan="1">0.094</td></tr><tr><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">rs2269067</td><td align="left" rowspan="1" colspan="1">0.261</td><td align="left" rowspan="1" colspan="1">7.64E-06</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">0.175</td></tr><tr><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">rs2159776</td><td align="left" rowspan="1" colspan="1">0.143</td><td align="left" rowspan="1" colspan="1">0.291</td><td align="left" rowspan="1" colspan="1">0.367</td><td align="left" rowspan="1" colspan="1">0.598</td></tr><tr><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">rs17611, rs7040033, rs2416811, rs9657673, rs3736855</td><td align="left" rowspan="1" colspan="1">0.328</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.716</td><td align="left" rowspan="1" colspan="1">0.450</td></tr><tr><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">rs10985126</td><td align="left" rowspan="1" colspan="1">0.206</td><td align="left" rowspan="1" colspan="1">1.86E-04</td><td align="left" rowspan="1" colspan="1">0.103</td><td align="left" rowspan="1" colspan="1">0.992</td></tr><tr><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">rs1323472, rs7030849</td><td align="left" rowspan="1" colspan="1">0.585</td><td align="left" rowspan="1" colspan="1">1.99E-04</td><td align="left" rowspan="1" colspan="1">0.935</td><td align="left" rowspan="1" colspan="1">0.415</td></tr><tr><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">rs12683062</td><td align="left" rowspan="1" colspan="1">0.113</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">0.696</td><td align="left" rowspan="1" colspan="1">0.327</td></tr><tr><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">rs3747843</td><td align="left" rowspan="1" colspan="1">0.337</td><td align="left" rowspan="1" colspan="1">0.112</td><td align="left" rowspan="1" colspan="1">0.123</td><td align="left" rowspan="1" colspan="1">0.059</td></tr><tr><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">rs10760152, rs10081760</td><td align="left" rowspan="1" colspan="1">0.297</td><td align="left" rowspan="1" colspan="1">0.007</td><td align="left" rowspan="1" colspan="1">0.933</td><td align="left" rowspan="1" colspan="1">0.790</td></tr><tr><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">rs942152</td><td align="left" rowspan="1" colspan="1">0.434</td><td align="left" rowspan="1" colspan="1">2.92E-05</td><td align="left" rowspan="1" colspan="1">0.161</td><td align="left" rowspan="1" colspan="1">0.919</td></tr><tr><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">rs9408928, rs9409230</td><td align="left" rowspan="1" colspan="1">0.063</td><td align="left" rowspan="1" colspan="1">0.270</td><td align="left" rowspan="1" colspan="1">0.955</td><td align="left" rowspan="1" colspan="1">0.307</td></tr><tr><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">rs10985196, rs7046030, rs12683459</td><td align="left" rowspan="1" colspan="1">0.089</td><td align="left" rowspan="1" colspan="1">6.17E-06</td><td align="left" rowspan="1" colspan="1">0.001</td><td align="left" rowspan="1" colspan="1">–</td></tr><tr><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">rs306781</td><td align="left" rowspan="1" colspan="1">0.015</td><td align="left" rowspan="1" colspan="1">0.905</td><td align="left" rowspan="1" colspan="1">0.661</td><td align="left" rowspan="1" colspan="1">0.147</td></tr><tr><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">rs4837839</td><td align="left" rowspan="1" colspan="1">0.079</td><td align="left" rowspan="1" colspan="1">0.171</td><td align="left" rowspan="1" colspan="1">0.988</td><td align="left" rowspan="1" colspan="1">0.667</td></tr><tr><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">rs306783</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.192</td><td align="left" rowspan="1" colspan="1">0.210</td><td align="left" rowspan="1" colspan="1">0.987</td></tr><tr><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">rs306784</td><td align="left" rowspan="1" colspan="1">0.009</td><td align="left" rowspan="1" colspan="1">0.054</td><td align="left" rowspan="1" colspan="1">0.144</td><td align="left" rowspan="1" colspan="1">0.876</td></tr><tr><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">rs10818527</td><td align="left" rowspan="1" colspan="1">0.02</td><td align="left" rowspan="1" colspan="1">0.007</td><td align="left" rowspan="1" colspan="1">0.044</td><td align="left" rowspan="1" colspan="1">0.368</td></tr><tr><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">rs12683989</td><td align="left" rowspan="1" colspan="1">0.019</td><td align="left" rowspan="1" colspan="1">0.009</td><td align="left" rowspan="1" colspan="1">0.054</td><td align="left" rowspan="1" colspan="1">0.573</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000107.s001"><label>Figure S1</label><caption><p>Genetic background-conditioned analysis. Allelic association for individual SNPs genotyped in Sample Set 2 (in blue) plotted as a function of position along with a Mantel-Haenszel P-value using stratifying information from an ancestry clustering procedure (in red).</p><p>(0.41 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000107.s002"><label>Table S1</label><caption><p>Minor allele frequencies and allele-based association of chr 9q33 SNPs with RA.</p><p>(0.10 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000107.s003"><label>Table S2</label><caption><p>Log likelihood ratio P-values for forward and backwards models using logistic regression.</p><p>(0.05 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000107.s004"><label>Table S3</label><caption><p>RA risk estimates for 3 loci – <italic>HLA-SE</italic>, <italic>PTPN22</italic>, and <italic>TRAF1</italic> – assuming a disease prevalence of 1%, 10% and 30%.</p><p>(0.11 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000107.s005"><label>Text S1</label><caption><p>Rheumatoid factor analysis and multilocus RA risk calculations.</p><p>(0.06 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>a</label><p>138 SNPs, including rs1953126, were genotyped in this sample set. Note, all 950 samples were genotyped for a single SNP, rs10818488, in the candidate gene study performed by Kurreeman et al ##REF##17880261##[35]##.</p></fn><fn id="nt102"><label>b</label><p>73 SNPs, including rs1953126, were genotyped in this sample set. Note, 475 of these patient samples were included in the initial whole genome association study performed by Plenge et al ##UREF##1##[34]##.</p></fn><fn id="nt103"><label>c</label><p>43 SNPs, including rs1953126, were genotyped in this sample set. Note, 436 patients and 94 controls samples were included in the candidate gene study performed by Kurreeman et al ##REF##17880261##[35]##.</p></fn><fn id="nt104"><label>d</label><p>Information on gender was available for 558 patients.</p></fn><fn id="nt105"><label>e</label><p>Information on age of onset was available for 306 patients.</p></fn><fn id="nt106"><label>f</label><p>Information on RF status was available for 440 patients.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt107"><label>a</label><p>Positions according to genomic contig NT_008470.18 (Entrez Nucleotide). The minor allele is listed first, followed by the position in National Center for Biotechnology Information Genome Build 36.2 and then the major allele.</p></fn><fn id="nt108"><label>b</label><p>Calculated for the minor allele using a Mantel-Haenszel common OR.</p></fn><fn id="nt109"><label>c</label><p>Calculated using Fisher's combined test.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt110"><label>a</label><p>The Haplo.Stats package was used to test for association between haplotypes and disease status.</p></fn><fn id="nt111"><label>b</label><p>Calculated with use of Fisher's combined test.</p></fn><fn id="nt112"><label>c</label><p>These haplotypes consist of SNPs: rs2239657, rs7021880 and rs7021049 respectively.</p></fn><fn id="nt113"><label>d</label><p>Calculated for haplotypes with the same effect (risk or protection) in all three sample sets, with use of Fisher's combined test.</p></fn><fn id="nt114"><label>e</label><p>Mantel-Haenszel common odds ratio with confidence intervals from Monte Carlo simulation.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt115"><label>a</label><p>Calculated using a Williams-corrected G test.</p></fn><fn id="nt116"><label>b</label><p>Calculated using Fisher's combined test.</p></fn><fn id="nt117"><label>c</label><p>Allele 1 rs2239657–allele 1 rs7021880–allele 1 rs7021049 / allele 2 rs2239657–allele 2 rs7021880–allele 2 rs7021049.</p></fn><fn id="nt118"><label>d</label><p>P-values calculated using Fisher's exact test.</p></fn><fn id="nt119"><label>e</label><p>Calculated for diplotypes with the same effect (risk or protection) in all three sample sets, with use of Fisher's combined test.</p></fn><fn id="nt120"><label>f</label><p>Mantel-Haenszel common odds ratio with confidence intervals from Monte Carlo simulation.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt121"><label>a</label><p>Genotypic P-values were calculated except where indicated.</p></fn><fn id="nt122"><label>b</label><p>All cases in this sample set were RF-positive.</p></fn><fn id="nt123"><label>c</label><p>Differential effects between RF-positive and RF-negative association were determined for sample set 2 using a Breslow-Day test (cases and controls were individually matched).</p></fn><fn id="nt124"><label>d</label><p>Differential effects between RF-positive and RF-negative association were determined for sample set 3 using a Monte Carlo simulation (cases and controls were not individually matched).</p></fn><fn id="nt125"><label>e</label><p>Includes all three sample sets.</p></fn><fn id="nt126"><label>f</label><p>Includes sample sets 2 and 3.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt127"><label>a</label><p>SNPs were grouped together if their pairwise r<sup>2</sup> values were &gt;0.90. The first SNP in each group was used for the analyses. With the exception of Groups 1 and 2, they are listed in order of appearance on the chromosome (for groups of SNPs we used the position of the first SNP). To avoid confusion, we retained the identity of the Group 1 and 2 SNPs assigned in ##FIG##1##Figure 2##.</p></fn><fn id="nt128"><label>b</label><p>Pairwise LD between rs7021049 and each of the 27 other SNPs as measured by r<sup>2</sup> in the cases and controls of the combined analysis of all three sample sets.</p></fn><fn id="nt129"><label>c</label><p>Univariate analysis using logistic regression.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>MC, CMR, VEG, SJS, JJC and ABB are employees of Celera and, as such, may own stock and/or stock options in Celera.</p></fn><fn fn-type="financial-disclosure"><p>Collection of the NARAC samples has been funded by a National Arthritis Foundation grant and the National Institutes of Health, acting through the National Institute of Arthritis and Musculoskeletal and Skin Diseases and the National Institute of Allergy and Infectious Diseases (contracts N01-AR-7-2232 and R01-AR44222). These studies were performed in part in the General Clinical Research Center, Moffitt Hospital, University of California, San Francisco, with funds provided by the National Center for Research Resources (5 M01 RR-00079, U.S. Public Health Service). Support was also provided by the National Institutes of Health through grant HG02275. The Department of Rheumatology, Leiden University Medical Center, received funding from the Dutch Arthritis Foundation.</p></fn></fn-group>
<graphic id="pgen-1000107-t001-1" xlink:href="pgen.1000107.t001"/> <graphic xlink:href="pgen.1000107.g001"/> <graphic xlink:href="pgen.1000107.g002"/> <graphic id="pgen-1000107-t002-2" xlink:href="pgen.1000107.t002"/> <graphic xlink:href="pgen.1000107.g003"/> <graphic id="pgen-1000107-t003-3" xlink:href="pgen.1000107.t003"/> <graphic id="pgen-1000107-t004-4" xlink:href="pgen.1000107.t004"/> <graphic id="pgen-1000107-t005-5" xlink:href="pgen.1000107.t005"/> <graphic id="pgen-1000107-t006-6" xlink:href="pgen.1000107.t006"/> <graphic xlink:href="pgen.1000107.g004"/>
<media xlink:href="pgen.1000107.s001.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000107.s002.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000107.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000107.s004.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000107.s005.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["28"], "element-citation": ["\n"], "surname": ["Long", "Grote", "Langley"], "given-names": ["AD", "MN", "CH"], "year": ["1997"], "article-title": ["Genetic analysis of complex diseases."], "source": ["Science"], "volume": ["275"], "fpage": ["1328"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["Plenge", "Seielstad", "Padyukov", "Lee", "Remmers"], "given-names": ["RM", "M", "L", "AT", "EF"], "year": ["2007"], "article-title": ["\n"], "italic": ["TRAF1-C5"], "source": ["N Eng J Med"], "volume": ["357"], "fpage": ["1199"], "lpage": ["1209"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Sokal", "Rohlf"], "given-names": ["RR", "FJ"], "year": ["1995"], "article-title": ["Biometry 3rd ed."], "publisher-loc": ["United States"], "publisher-name": ["W. H. Freeman and Company"]}, {"label": ["49"], "element-citation": ["\n"], "surname": ["Barton", "Thomson", "Ke", "Eyre", "Hinks"], "given-names": ["A", "W", "X", "S", "A"], "year": ["2008"], "article-title": ["Re-evaluation of putative rheumatoid arthritis susceptibility genes in the post-genome wide association study era and hypothesis of a key pathway underlying susceptibility."], "source": ["Hum Mol Genet"], "comment": ["Online pub April 22, 2008"]}, {"label": ["80"], "element-citation": ["\n"], "surname": ["Armitage"], "given-names": ["P"], "year": ["1955"], "article-title": ["Tests for linear trends in proportions and frequencies."], "source": ["Biometrics"], "volume": ["11"], "fpage": ["375"], "lpage": ["386"]}, {"label": ["81"], "element-citation": ["\n"], "surname": ["Weir"], "given-names": ["BS"], "year": ["1996"], "article-title": ["Genetic Data Analysis II: Methods for Discrete Population Genetic Data."], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer"]}, {"label": ["82"], "element-citation": ["\n"], "surname": ["Fisher"], "given-names": ["RA"], "year": ["1954"], "article-title": ["\n"], "italic": ["Statistical Methods for Research Workers."], "publisher-loc": ["Edinburgh"], "publisher-name": ["Oliver & Boyd"]}, {"label": ["85"], "element-citation": ["\n"], "surname": ["Breslow", "Day"], "given-names": ["NE", "NE"], "year": ["1980"], "article-title": ["Statistical methods in cancer research. Vol I. The analysis of case-control studies."], "publisher-loc": ["Lyon"], "publisher-name": ["International Agency for Research on Cancer"]}]
PLoS Genet. 2008 Jun 27; 4(6):e1000107