Genome-wide association analysis implicates the involvement of eight loci with response to tocilizumab for the treatment of rheumatoid arthritis.

Rheumatoid arthritis (RA) is an immune-mediated inflammatory disease affecting the joints. A heterogeneous response to available therapies demonstrates the need to identify those patients likely to benefit from a particular therapy. Our objective was to identify genetic factors associated with response to tocilizumab, a humanized monoclonal antibody targeting the interleukin (IL)-6 receptor, recently approved for treating RA. We report the first genome-wide association study on the response to tocilizumab in 1683 subjects with RA from six clinical studies. Putative associations were identified with eight loci, previously unrecognized as linked to the IL-6 pathway or associated with RA risk. This study suggests that it is unlikely that a major genetic determinant of response exists, and it illustrates the complexity of performing genome-wide association scans in clinical trials.

Factor V Leiden and thrombosis in patients with systemic lupus erythematosus: a meta-analysis.

The aim of this study was to perform a meta-analysis of the association between the factor V Leiden polymorphism (FVL) and thrombosis among patients with systemic lupus erythematosus (SLE) and/or antiphospholipid antibody (aPL) positivity. Included studies recruited patients based on SLE or aPL-positive status, confirmed subjects' SLE diagnosis as defined by the American College of Rheumatology, and documented thrombotic events. Excluded studies were non-English or considered only arterial thrombosis. Individual patient data, available from 5 studies, together with unpublished data from 1210 European-American SLE patients from the UCSF Lupus Genetics Collection genotyped for FVL, were further analyzed. Seventeen studies (n=2090 subjects) were included in the initial meta-analysis. Unadjusted odds ratios (OR) were calculated to assess association of FVL with thrombosis. The OR for association of thrombosis with FVL was 2.88 (95% confidence interval (CI) 1.98-4.20). In the secondary analysis with our individual patient dataset (n=1447 European-derived individuals), SLE subjects with the FVL polymorphism still had more than two times the odds of thrombosis compared to subjects without this polymorphism, even when adjusting for covariates such as gender, age and aPL status. SLE and/or aPL-positive patients with the FVL variant have more than two times the odds of thrombosis compared to those without this polymorphism.

Detailed genetic characterization of the interleukin-23 receptor in psoriasis.

Using a multi-tiered, case-control association design, scanning 25 215 gene-centric SNPs, we previously identified two psoriasis susceptibility genes: IL12B and IL23R. These results have recently been confirmed. To better characterize the IL23R psoriasis-association, we used a fine mapping strategy to identify 59 additional IL23R-linked SNPs, which were genotyped in our three independent, white North American sample sets (>2800 individuals in toto). A sliding window of haplotype association demonstrates colocalization of psoriasis susceptibility effects within the boundaries of IL23R across all sample sets, thereby decreasing the likelihood that neighboring genes, particularly IL12RB2, are driving the association at this region. Additional haplotype work identified two 5-SNP haplotypes with strong protective effects, consistent across our three sample sets (OR(common)=0.67; P(comb)=4.32E-07). Importantly, heterogeneity of effect was extremely low between sample sets for these haplotypes (P(Het)=0.961). Together, these protective haplotypes attain a frequency of 16% in controls, declining to 11% in cases. The characterization of association patterns within IL23R to specific predisposing/protective variants will play an important role in the elucidation of psoriasis etiology and other related phenotypes. Further, this work is essential to lay the foundation for the role of IL23R genetics in response to pharmaceutical therapy and dosage.

Strategies and technical challenges in allele level Class II typing in 2578 bone marrow transplantation donor-recipient pairs.

Human leukocyte antigen typing of 2578 donor-recipient pairs whose transplantation was facilitated by the National Marrow Donor Program allowed for an in-depth analysis of the accuracy of high-volume allele level testing data. The methods employed provided allele level typing at DRB1/3/5, DQA1, DQB1, DPA1, and DPB1 using sequence-specific oligonucleotide probe hybridization (SSOPH), polymerase chain reaction (PCR) restriction fragment length polymorphism analysis, sequence specific PCR, and direct sequence-based typing (SBT). Each typing was independently tested by two laboratories in Phase 1, and in subsequent phases targeted samples were typed in duplicate by SBT to monitor typing quality. Comparison with prior transplant center typing was also evaluated. SSOPH detected discrepancies ranged from 0.6% at DPB1 to 5.1% at DQB1 in Phase 1. The majority of discrepancies, 62%, resulted from human error such as sample handling, result interpretation, or clerical errors. Alleles that are frequently discrepant have been identified in this predominantly white population.

Variants in the 5q31 cytokine gene cluster are associated with psoriasis.

A multitiered genetic association study of 25 215 single-nucleotide polymorphisms (SNPs) in three case-control sample sets (1446 patients and 1432 controls) identified three IL13-linked SNPs (rs1800925, rs20541 and rs848) associated with psoriasis. Although the susceptibility effects at these SNPs were modest (joint allelic odds ratios (ORs): 0.76 to 0.78; P(comb): 1.3E-03 to 2.50E-04), the association patterns were consistent across the sample sets, with the minor alleles being protective. Haplotype analyses identified one common, susceptible haplotype CCG (joint allelic OR=1.27; P(comb)=1.88E-04) and a less common, protective haplotype TTT (joint allelic OR=0.74; P(comb)=7.05E-04). In combination with the other known genetic risk factors, HLA-C, IL12B and IL23R, the variants reported here generate an 11-fold psoriasis-risk differential. Residing in the 5q31 cytokine gene cluster, IL13 encodes an important T-cell-derived cytokine that regulates cell-mediated immunity. These results provide the foundation for additional studies required to fully dissect the associations within this cytokine-rich genomic region, as polymorphisms in closely linked candidate genes, such as IRF1, IL5 or IL4, may be driving these results through linkage disequilibrium.

Risk of progression from undifferentiated arthritis to rheumatoid arthritis: the effect of the PTPN22 1858T-allele in anti-citrullinated peptide antibody positive patients.

OBJECTIVES: Anti-citrullinated peptide antibodies (ACPA) and the C1858T missense single-nucleotide polymorphism (SNP) in the PTPN22 gene are both associated with the development of rheumatoid arthritis (RA). We investigated whether the combination of these two biomarkers yielded better test characteristics to predict progression from undifferentiated arthritis (UA) to RA compared with ACPA alone. METHODS: A total of 394 individuals with UA from a Dutch population-based inception cohort were included in this study. At baseline, ACPA were measured and the PTPN22 C1858T and HLA-DRB1 genotypes determined. Progression to RA was monitored at 1 yr after entry into the cohort. RESULTS: A priori, UA patients had a 35% (95% CI 30-40%) risk of developing RA, which increased to 66% (95% CI 57-75%) in patients who were ACPA-positive. There was an additional, although non-significant (P = 0.34), increase in RA risk to 76% (95% CI 57-90%) when patients were positive for both ACPA and the PTPN22 1858T-allele. The area under the receiver operator characteristic curve increased from 0.68 for ACPA-status alone to 0.70 for the combination of ACPA-status and the PTPN22 C1858T polymorphism. In logistic regression analysis, ACPA predicted RA-development independent of PTPN22, while the PTPN22 polymorphism had no independent effect. In HLA-DRB1 shared epitope positive, ACPA-positive UA patients, ACPA-levels were significantly increased in PTPN22 1858T allele carriers compared with non-1858T carriers. CONCLUSIONS: In this Dutch cohort of UA-patients, the PTPN22 1858T allele does not markedly improve individual decision-making to predict RA-development over ACPA alone, but it is associated with higher ACPA-levels.

The 620W allele is the PTPN22 genetic variant conferring susceptibility to RA in a Dutch population.

OBJECTIVES: A missense SNP, C1858T, in PTPN22 has been identified as a genetic risk factor for rheumatoid arthritis (RA). Subsequent work has suggested that other variants in this gene, in particular a haplotype marked by the minor allele of rs3789604, are associated with RA in white North Americans independent of C1858T. We tested this hypothesis in an independent white Dutch study. METHODS: A total of 667 RA patients and 286 controls were genotyped for 13 PTPN22 single nucleotide polymorphisms (SNPs) by allele-specific kinetic polymerase chain reaction. rs3789604 was genotyped in an additional 410 RA and 270 UA patients participating in the Leiden early arthritis inception cohort. We conducted single-marker and haplotype association tests. RESULTS: The sole haplotype strongly associated with RA in our Dutch population carries the PTPN22 1858T allele. A second haplotype identical at all other SNPs tested except 1858 was not associated with disease. No significant association of the haplotype tagged by the 3' PTPN22 SNP, rs3789604, with RA susceptibility (P = 0.134) was observed in our sample set. CONCLUSION: We conclude that C1858T is the sole PTPN22 variant predisposing to RA in our white Dutch sample set.

Finnish case-control and family studies support PTPN22 R620W polymorphism as a risk factor in rheumatoid arthritis, but suggest only minimal or no effect in juvenile idiopathic arthritis.

Several studies have identified the PTPN22 allelic variant 1858 C/T that encodes the R620W amino-acid change as a putative susceptibility factor in autoimmune diseases. The current study was undertaken to examine a large cohort of Finnish rheumatoid arthritis (RA) and juvenile idiopathic arthritis (JIA) subjects using both population control and, importantly, family-based association methods. The latter is particularly important when, as is the case for the 1858 C/T polymorphism, the frequency of the variant allele (T) differs in both major ancestral populations and in subpopulations. The analysis of rheumatoid factor-positive 1030 RA probands from Finland provides strong support for association of this variant in both population studies (allele specific odds ratio (OR)=1.47, 95% confidence interval (CI)=1.27-1.70, P=3 x 10(-7)) and in family studies (P<10(-6)). In contrast, no allelic association was seen with JIA (230 probands) and only weak evidence for a genotypic effect of 1858T homozygotes was observed in this population.

The PTPN22 R620W polymorphism associates with RF positive rheumatoid arthritis in a dose-dependent manner but not with HLA-SE status.

We have recently described the association between rheumatoid arthritis and a coding single-nucleotide polymorphism in the intracellular protein tyrosine phosphatase, PTPN22. The disease-associated polymorphism, 1858 C/T (rs2476601), encodes an amino-acid change (R620W) in one of four SH3 domain binding sites in the PTPN22 molecule. We have now extended our initial studies to address three questions: (1) Is the association with rheumatoid arthritis limited to rheumatoid factor (RF) positive disease? (2) Does homozygosity for PTPN22 R620W substantially increase disease susceptibility? (3) Is there an interaction between PTPN22 and the rheumatoid arthritis (RA)-associated HLA-DRB1 shared epitope alleles? A total of 1413 Caucasian rheumatoid arthritis patients and 1401 Caucasian controls were genotyped. The results support the view that PTPN22 was strongly and preferentially associated with RF positive disease (OR=1.75, 95% CI 1.46-2.10, P=1.3 x 10(-9)). The PTPN22 risk allele was not significantly associated with RF negative disease (OR=1.19, 95% CI 0.92-1.53, P=0.18), although a very weak association cannot be completely excluded. There was a strong dose effect on disease risk; two copies of the PTPN22 R620W allele more than doubles the risk for RF positive RA (OR=4.57, 95% CI 2.35-8.89). There was no evidence of a genetic association between PTPN22 and HLA susceptibility alleles.

Expression levels for many genes in human peripheral blood cells are highly sensitive to ex vivo incubation.

Monitoring of gene and protein expression in peripheral blood cells has significant potential for improving the diagnosis and therapy of many human diseases. As genomic-scale microarray and proteomic technologies are applied to peripheral blood, it is important to consider the variables that may affect interpretation of data. Here we report experiments performed to identify genes that are particularly sensitive to ex vivo handling prior to RNA extraction for gene expression microarrays or quantitative real-time RT-PCR assays. We examined Affymetrix gene expression in samples from eight normal individuals where blood was processed for RNA either immediately after blood draw or the next day following overnight incubation. These studies identified hundreds of genes that are sensitive to ex vivo handling of blood, and suggest that this is an important variable to consider when designing and interpreting human PBMC experiments.

HLA-DR2 dose effect on susceptibility to multiple sclerosis and influence on disease course.

Models of disease susceptibility in multiple sclerosis (MS) often assume a dominant action for the HLA-DRB1*1501 allele and its associated haplotype (DRB1*1501-DQB1*0602 or DR2). A robust and phenotypically well-characterized MS data set was used to explore this model in more detail. A dose effect of HLA-DR2 haplotypes on MS susceptibility was revealed. This observation suggests that, in addition to the role of HLA-DR2 in MS, two copies of a susceptibility haplotype further increase disease risk. Second, we report that DR2 haplotypes modify disease expression. There is a paucity of benign MS and an increase of severe MS in individuals homozygous for DR2. Concepts of the molecular mechanisms that underlie linkage and association of the human leukocyte antigen (HLA) region to MS need to be revised to accommodate these data.

Two new HLA-B alleles, B*4422 and B*4704, identified in a study of families with autoimmunity.

HLA-B typing of approximately 1 262 individuals from a study of 372 simplex families with multiple sclerosis has led to the identification of two new alleles (HLA-B*4422 and HLA-B*4704). Sequencing confirmed that both of these new alleles represent novel combinations of previously described sequence motifs, reinforcing the notion that recombination and/or gene conversion-like events play an important role in generating HLA allelic diversity. The identification of these new alleles brings the total number of HLA-B alleles to 465.

Effect of HLA class II gene disparity on clinical outcome in unrelated donor hematopoietic cell transplantation for chronic myeloid leukemia: the US National Marrow Donor Program Experience.

The clinical importance of HLA class II gene disparity in unrelated stem cell transplantation is not entirely known. The impact was evaluated of matching donors and recipients for HLA-DR, HLA-DQ, and HLA-DP genes on clinical outcome after stem cell transplantation for chronic myeloid leukemia (CML) performed between 1988 and 1997. HLA-DRB1, -DQA1, -DQB1, -DPA1, and -DPB1 alleles were identified in 831 transplant pairs using a combination of sequence-specific oligonucleotide probes, sequence-specific priming, and sequencing methods. Among the 831 pairs, 696 (84%) were HLA-A and -B serologically matched; of these, 565 (81%) were also matched for HLA-DRB1. HLA-DRB1 matching correlated with significantly improved survival (relative risk [RR], 1.29 [95% confidence interval (CI), 1.02-1.64; P =.04]) independently of HLA-DQA1 or HLA-DQB1 (RR, 1.01 [95% CI, 0.81-1.26; P =.94]) and HLA-DPA1 or HLA-DPB1 (RR, 1.11 [95% CI, 0.84-1.48; P =.46]). Single-locus HLA-DQ or HLA-DP disparity was not associated with significantly poorer survival. For patients who underwent transplantation in the first chronic phase (CP) from HLA-A, B matched donors, the presence of DRB1 allele mismatching was independently associated with increased incidence of grades III-IV acute graft-versus-host disease (GVHD). No significant associations of class II allele mismatching with risk for delayed engraftment or chronic GVHD disease were detected. This study clearly demonstrates the importance of precise matching of HLA-DRB1 alleles for successful transplantation. Furthermore, a good-risk population of patients whose transplantations were performed in the first CP of disease from HLA-A, B, DRB1 matched unrelated donors can be shown to have superior survival.

Genetic variability and linkage disequilibrium within the HLA-DP region: analysis of 15 different populations.

In order to understand the forces governing the evolution of the genetic diversity in the HLA-DP molecule, polymerase chain reaction (PCR)-based methods were used to characterize genetic variation at the DPA1 and DPB1 loci encoding this heterodimer on 2,807 chromosomes from 15 different populations including individuals of African, Asian, Amerindian, Indian and European origin. These ethnically diverse samples represent a variety of population substructures and include small, isolated populations as well as larger, presumably admixed populations. Ten DPA1 and 39 DPB1 alleles were identified and observed on 87 distinct DP haplotypes, 34 of which were found to be in significant positive linkage disequilibrium in at least one population. Some haplotypes were found in all ethnic groups while others were confined to a single ethnic group or population. Strong positive global linkage disequilibrium (Wn) between DPA1 and DPB1 was present in all 15 populations. The African populations displayed the lowest values of Wn whereas the Amerindian populations displayed near absolute disequilibrium. Analysis of the distribution of haplotypes using the normalized deviate of the Ewens-Watterson homozygosity statistic, F, suggests that DP haplotypes encoding the functional heterodimer are subject to much lower degrees of balancing selection than other loci within the HLA region. Finally, neighbor joining tree analyses demonstrate the power of haplotype diversity for inferring the relationships between the different populations.

Three new DP alleles identified in a study of 800 unrelated bone marrow donor-recipient pairs.

HLA-DP genotyping of 800 unrelated donor-recipient pairs in phase 5 of a retrospective analysis of unrelated bone marrow transplantation, sponsored by the National Marrow Donor Program (NMDP), has identified two new DPB1 alleles (DPB1*8701 and DBP1*8801) and one new DPA1 (DPA1*0108) allele. Sequencing confirmed that all three of these new alleles represent novel combinations of previously described sequence motifs, reinforcing the notion that "gene conversion-like" events play an important role in generating HLA allelic diversity. The identification of these new alleles brings the total number of DPA1 alleles to 20 and the total number of DPB1 alleles to 94.

Evolution of Pacific/Asian populations inferred from HLA class II allele frequency distributions.

The allele frequency distributions for the HLA class II loci, DRB1, DQB1 and DPB1, in eight Pacific/Asian populations: Hawaiian, Samoan, Malay, Papua New Guinea (PNG) Highlands, and two Indonesian and PNG Lowland groups, were determined using high-resolution polymerase chain reaction/sequence-specific oligonucleotide probe (PCR/SSOP) typing methods. The allele frequency distributions for the HLA-DRB1 locus were determined for a third Indonesian population as well as for an additional Filipino population. DRB1 alleles in the DR2 serogroup (or allelic lineage) are very common in this region; in some populations, more than 50% of the alleles belong to this serogroup. The DRB1*1502 allele is frequent in nine of the ten populations studied, reaching a frequency of 0.48 in one Indonesian population and among Filipinos. Extensive DR-DQ haplotype diversity was detected in these populations. Seven different DR2-DQB1 haplotypes were observed in the Indonesian and PNG Lowland populations, eight in the PNG Highlands and ten in Malays and Filipinos. The DRB1*0410 allele, commonly observed in Australia, is observed in the PNG Highlands at a low frequency (f=0.03) and is absent in the other populations. Two additional DRB1 alleles commonly observed in Australia, DRB1*0405 and *1407, are also observed in the PNG Highlands at high frequencies (f=0.132 and 0.126), while they are rare in the PNG Lowlands (f=0.039 and 0.013). These alleles are generally rare or absent in the other populations. The DPB1*0501 allele, common in Chinese and Japanese populations, is most frequent in the Samoan, Hawaiian, Indonesian, and Malay populations, and the *0401 allele is the most frequent DPB1 allele in the PNG Lowlands. Both of these alleles have the same very high frequency (f=0.34) in the PNG Highlands. Analyses of homozygosity (the Ewens-Watterson F statistic) in these and other populations indicate that, while most allele frequency distributions are consistent with balancing selection, values of F for the Indonesian and Javan populations may reflect positive directional selection. Phylogenetic trees constructed using the allele frequencies at the DRB1 locus of the populations reported here, as well as those for additional Pacific, Asian, and Australian populations, indicate that the PNG Highland population is more closely related to Australian populations than to PNG Lowland populations, while the PNG Lowlands are more closely related to other Melanesian populations.

The extent of HLA class II allele level disparity in unrelated bone marrow transplantation: analysis of 1259 National Marrow Donor Program donor-recipient pairs.

A comprehensive analysis of the HLA-D region loci, DRB1, DRB3, DRB5, DQA1, DQB1, DPA1 and DPB1, was performed to determine allelic diversity and underlying HLA disparity in 1259 bone marrow recipients and their unrelated donors transplanted through the National Marrow Donor Program. Although 43.0% of DRB1 alleles known to exist at the beginning of the study were found in this predominantly Caucasian transplant population, a few alleles predominated at each locus. In recipients, 67.1% of DRB1 alleles identified were one or two of six common DRB1 alleles. Only 118 (9.4%) donor-recipient pairs were matched for all alleles of DRB1, DQA1, DQB1, DPA1 and DPB1. While 79.4% of the pairs were matched for DRB1, only 13.2% were matched for DPB1 alleles. Almost 66% of pairs differed by more than one allele mismatch and 59.0% differed at more than one HLA-D locus. DQB1 was matched in 85.9% of DRB1-matched pairs. In contrast, only 13.9% of the pairs matched for DRB1, DQA1 and DQB1 were also matched for DPA1 and DPB1. This database, highlighting the underlying HLA disparity within the pairs, forms the foundation of an ongoing study to establish the relationship between HLA matching and successful outcome in unrelated allogeneic stem cell transplant.

T cell receptor BV gene rearrangements in the spinal cords and cerebrospinal fluid of patients with amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal disorder whose etiology and pathogenesis remain unknown. Recent studies, however, have demonstrated the presence of inflammatory infiltrates within ALS spinal cord and suggested the possibility of an immune-mediated process in motor neuron degeneration. We have analyzed the diversity of T-cells in the spinal cord in ALS. Reverse transcriptase polymerase chain reaction (RT-PCR) with variable (V) region sequence specific oligonucleotide primers was used to amplify T-cell receptor (TCR)BV transcripts from spinal cords obtained at autopsy from patients with ALS, patients who died without inflammatory disease of the central nervous system, brains from patients with ALS, and brains from patients who died with inflammatory CNS disease. Sequencing was then performed on the amplified transcripts. An overall increase in the level of TCRBV 2 transcripts was detected in ALS specimens when compared to controls. This result was independent of the HLA genotype of the individual. Furthermore, enrichment of TCRBV2-positive T cells could be demonstrated in cerebrospinal fluid derived from patients with ALS, using PCR analysis and a T cell stimulation assay with toxic shock syndrome toxin-1 (TSST-1), a Vbeta2-specific superantigen. Our results suggest that an immunological process involving the specific expansion of Vbeta2 TCR-positive T-cells may be important in the pathogenesis of ALS.