Psoriasis Patients Are Enriched for Genetic Variants That Protect against HIV-1 Disease

An important paradigm in evolutionary genetics is that of a delicate balance between genetic variants that favorably boost host control of infection but which may unfavorably increase susceptibility to autoimmune disease. Here, we investigated whether patients with psoriasis, a common immune-mediated disease of the skin, are enriched for genetic variants that limit the ability of HIV-1 virus to replicate after infection. We analyzed the HLA class I and class II alleles of 1,727 Caucasian psoriasis cases and 3,581 controls and found that psoriasis patients are significantly more likely than controls to have gene variants that are protective against HIV-1 disease. This includes several HLA class I alleles associated with HIV-1 control; amino acid residues at HLA-B positions 67, 70, and 97 that mediate HIV-1 peptide binding; and the deletion polymorphism rs67384697 associated with high surface expression of HLA-C. We also found that the compound genotype KIR3DS1 plus HLA-B Bw4-80I, which respectively encode a natural killer cell activating receptor and its putative ligand, significantly increased psoriasis susceptibility. This compound genotype has also been associated with delay of progression to AIDS. Together, our results suggest that genetic variants that contribute to anti-viral immunity may predispose to the development of psoriasis.

We were intrigued by the observation that several of the most highly significant SNPs from psoriasis GWAS were identical to the top SNPs from GWAS of HIV-1 virologic control, a clinical phenotype whereby certain HIV-1 infected individuals, termed ''HIV-1 controllers,'' are able to maintain low levels of plasma HIV-1 RNA in the absence of antiretroviral therapy and who generally do not develop clinical symptoms [12]. For example, rs2395029 within the MHC gene HCP5 and a proxy for HLA-B*57, was identified in a psoriasis GWAS as the SNP with the largest odds ratio, OR = 4.1, p = 2.13610 226 [8]. This same SNP has been shown to be the first or second most significant SNP in three GWAS of HIV-1 control [13][14][15]. The two most significant SNPs identified in multiple psoriasis GWAS are rs10484554 and rs12191877 near HLA-C (r 2 = 1 with each other in Europeans) [8][9][10]. rs10484554 and rs12191877 were found to be associated with HIV-1 control [Supplementary Materials in [13,15]] and both are in moderate linkage disequilibrium (r 2 = 0.33) with rs9264942, another top SNP for HIV-1 control [14][15][16].
The relationship between psoriasis and HIV-1 is also interesting because of the clinical observation that HIV-1 infection can exacerbate existing psoriasis or trigger new-onset psoriasis [17]. As HIV-1 infection progresses and CD4+ T cell counts decrease, psoriasis can worsen [18,19]. This has puzzled dermatologists and infectious disease clinicians because it has been convincingly established that psoriasis is an immune disorder that is mediated through activation of T cells. Several explanations for this ''psoriasis HIV-1 paradox'' have been proposed, including HIV-1 induced destruction of regulatory CD4+ T cells [20], an increase in number of memory CD8+ T cells late in disease [21], HIV-1 proteins acting as superantigens [22], or co-stimulation through traditional antigenic presentation [20].
Due to these genetic and clinical observations, we pursued a more in-depth analysis of the HLA region in psoriasis to determine whether patients with psoriasis are enriched for the major genetic determinants of HIV-1 control. The psoriasis data generated in this study were compared to the largest GWAS for HIV-1 control performed to date, involving 516 cases and 1,196 controls and for which detailed HLA allele information was available [15].

Accurate Imputation of HLA Alleles
We imputed to four-digit resolution the HLA class I alleles (-A, -B, -C) and HLA class II alleles (-DQA1, -DQB1, -DRB1) of 1,727 psoriasis cases and 3,581 healthy controls which were obtained from 3 separate case-control cohorts of European ancestry (Table  S1). Imputation was performed using the software HLA*IMP, which has been shown to have an accuracy of at least 96% for class I loci and 92% for Class II loci [23]. To further validate the accuracy of our imputation, we compared the imputed HLA alleles to empirically obtained HLA class I alleles for a subset of our samples (n = 98). The concordance was 566/581 alleles (97.4%), indicating that the imputation was of high accuracy. A sensitivity analysis examining the imputation accuracy of low frequency HLA alleles (allele frequency between 1% and 5%) demonstrated similar high accuracy (177/179 alleles = 98.9%). Only HLA alleles with a minor allele frequency greater than 1% in the control group were used for subsequent analyses.

HLA Association Testing Identifies a Similar Genetic Architecture between Psoriasis and HIV-1 Control
We tested all imputed HLA alleles for association with psoriasis using logistic regression, adjusting for gender, ancestry, and cohort. The top ten HLA associations for psoriasis are shown in Table 1 (Full four-digit and two-digit results in Tables S2 and S3, respectively). Overall, we observed a striking pattern in which the HLA alleles which are enriched in psoriasis patients are also enriched in HIV-1 controllers, and the HLA alleles which have decreased frequency in psoriasis patients are also decreased in HIV-1 controllers. We found that psoriasis patients are highly enriched for HLA-B*57:01 (12.5% in cases vs 3.9% in controls, p = 5.50610 242 , OR = 3.61), which in multiple studies has been shown to be the most significant predictor of both HIV-1 control and delayed progression time to AIDS [14,15,[24][25][26][27]. Psoriasis patients also display a significant enrichment of the HIV-1 control allele B*13:02, whereas they display a relative paucity of B*07:02, B*40:01, and C*04:01 which are associated with lack of virologic control [15]. The HLA allele B*35, almost always seen with C*04:01, and the most significant HLA allele associated with rapid progression to AIDS [28,29], was significantly protective against psoriasis in our dataset (p = 3.20610 26 , OR = 0.65 [0.54-0.78]). HLA-B*35 alleles can be segregated into B*35-Px and B*35-PY alleles, where Px alleles bind peptides with hydrophobic, nontyrosine residues at position 9 and PY alleles bind peptides with tyrosine at position 9. It has been shown that the influence of HLA-B*35 in accelerating progression to AIDS is mostly attributable to HLA-B*35-Px alleles [30]. In our psoriasis dataset, the B* 35

Author Summary
Individuals with autoimmune disease generally demonstrate excessive immune system activation, leading to inflammation and damage of specific target organs. However, in some cases the detrimental effects of an overactive immune system might be counterbalanced by a beneficial effect in protecting against certain infections. In this study, we investigated whether patients with psoriasis, a common autoimmune disease of the skin, harbor genetic variants that are associated with an enhanced ability to limit replication of the HIV-1 virus. We profiled the HLA (human leukocyte antigen) immune genes located on chromosome 6 in 1,727 Caucasian psoriasis cases and 3,581 healthy controls and found that psoriasis patients are significantly more likely than controls to have gene variants that are protective against HIV-1 disease. We found that this enrichment for HIV-1 protective variants was unique to psoriasis and largely absent in patients with other autoimmune or inflammatory diseases such as rheumatoid arthritis, Crohn's disease, type 1 diabetes, type 2 diabetes, and coronary artery disease. Our results suggest the possibility that the excessive skin inflammation in psoriasis may be associated with activation of antiviral immune pathways that were important to human ancestors who encountered viruses similar to HIV-1.
Amino Acids within HLA-B That Are Predictive of HIV-1 Viral Load Are Concordant between HIV-1 Controllers and Psoriasis Patients Specific amino acid positions within the peptide binding groove of HLA class I molecules have been shown to serve as important mediators for the protective and risk effects of individual HLA alleles on HIV-1 control [15]. Namely, amino acid residues at Table 1. Top ten classical HLA alleles associated with psoriasis, and comparison to HIV-1 controllers as published in [15]. positions 97, 67, and 70 within HLA-B were found to be more highly associated with HIV-1 control than HLA-B*57:01 and each of these amino acid positions was found to serve as a strong predictor of HIV-1 viral load levels in an independent cohort [15].
To determine whether psoriasis is associated with the groups of alleles that are marked by specific amino acids within HLA proteins, we used the official protein sequences [31] assigned to each four-digit HLA allele to perform association testing at each amino acid position within HLA-A, -B, -C, -DQA1, -DQB1, and -DRB1. As before, the association testing was adjusted for gender, ancestry, and cohort. We found that the five most  (Table S4). Each of these 5 positions is located within the peptide-binding groove of the HLA molecule and directly contacts the bound peptide [32]. At each of these positions, we investigated whether the direction of the association signal was consistent between psoriasis patients and HIV-1 controllers. We confirmed that for each position examined, the amino acid residues associated with psoriasis susceptibility were associated with HIV-1 virologic control, and the amino acids associated with a protective effect on psoriasis risk were associated with HIV-1 progression ( Figure 1). For example, alleles marked by Asn 97 , Thr 97 , and Val 97 in HLA-B were associated with psoriasis susceptibility and HIV-1 control while those marked by Arg 97 and Ser 97 in HLA-B were associated with psoriasis protection and HIV-1 progression.
To rule out the possibility that the observed similarities between psoriasis and HIV-1 control were the result of systematic bias of the imputation process or general amino acid variability at those positions, we examined whether these 5 amino acids positions were associated with other autoimmune or inflammatory diseases. We examined GWAS data from 5 diseases studied by the Wellcome Trust Case Control Consortium [33]-rheumatoid arthritis, Crohn's disease, type 1 diabetes, type 2 diabetes, and coronary artery disease-and used the same imputation process as performed with psoriasis. We found that none of these diseases displayed the degree of similarity between psoriasis and HIV-1 control when considering the direction and magnitude of the association signal at these amino acid positions ( Figure 1). Crohn's disease, which shares some pathophysiological features with psoriasis [34] and is also slightly enriched for HLA-C*06:02 (p = 4.2610 25 , OR 1.32) and HLA-B*57:01 (p = 3.68610 24 , OR 1.40), showed some similarity to psoriasis and HIV-1 control at  these positions, but the magnitude of the association was smaller and several important residues such as Asn 70 and Asn 97 in HLA-B and Gln 156 in HLA-C were not concordant. Interestingly, type 1 diabetics showed the opposite effect at many of these residues (i.e. patients with type 1 diabetes lack HIV-1 control alleles and have an excess of HIV-1 progression alleles), which could support the theory that type 1 diabetes is triggered by a viral infection. Another important amino acid within HLA-B that may be relevant for HIV-1 progression is position 116, which not only interacts with the carboxy-terminal residues of peptides in the F pocket, but also strongly influences the interaction of HLA class I molecules with the peptide-loading complex [35,36]. Studies of HLA-B*44:05 and HLA-B*44:02, which only differ at position 116, have shown that B*44:05 (containing tyrosine at position 116, ''116Y'') utilizes a tapasin-independent pathway that leads to a less optimal peptide repertoire compared to the tapasin-dependent HLA-B*44:02 (116D) [36]. 116Y is strongly associated with lack of HIV-1 control, p = 1.6610 210 , OR = 0.57 [15]. Our data show that 116Y is strongly associated with decreased susceptibility to psoriasis, p = 7.96610 217 , OR = 0.66 [0.60-0.73] (Table S4) Together, our data indicate that the genetic similarity between psoriasis patients and HIV-1 controllers extends to specific amino acid residues within class I molecules that mediate peptide binding, influence peptide loading, and which mark viral control or progression. We additionally performed stepwise regression to identify amino acid residues that were independently associated with psoriasis and found Trp 156 and Ala 24 in HLA-C; Val 97 , Leu 145 , Cys 67 , and Tyr 99 in HLA-B; and Gly 107 in HLA-A to be markers independently associated with psoriasis (Table S5). Similar to HIV-1 control [15], we found that HLA-B positions 97 and 67 remained in the model, while position 70 dropped out due to linkage disequilibrium with positions 97 and 67. However, we caution against an interpretation that the amino acids identified here as independent are necessarily the functional ones. Due to the complex LD patterns between the amino acids, the final output of the stepwise regression model is affected by the starting variable, and potentially functionally significant amino acids can be lost because they are tagged by other residues. For example, HLA-B Gly 62 , part of the a1 helix located in the B-pocket of the peptide binding groove, shows strong independent association with HIV-1 control (p = 4.6610 227 , OR = 5.03) [15]. Gly 62 is also strongly associated with susceptibility to psoriasis (p = 2.03610 239 , OR = 3.20) but is in high LD with HLA-B Val 97 and thus drops out of the final psoriasis model.
The HIV-1 Protective HLA-C 39UTR Deletion Polymorphism rs67384697 Is Enriched in Psoriasis Expression levels of HLA-C are modulated by the G/polymorphism rs67384697 located within the 39 UTR of HLA-C, where the presence of the deletion inhibits the binding of the microRNA hsa-miR-148 to the 39UTR and results in higher HLA-C surface expression [37]. In a multivariate model of HIV-1 control, the deletion allele of rs67384697 has a strong effect on viral control independent of the classical alleles HLA-B*57:01 and HLA-B*27:05, although the high LD of rs67384697 with other HLA-B and HLA-C alleles makes it difficult to determine whether rs67384697 (high HLA-C expression) is directly mediating this effect, or whether the HLA alleles themselves are causal [37]. Nevertheless, rs67384697 has been proposed to be the functional variant that explains the previously identified protective effect of rs9264942 on HIV-1 control, where rs9264942 is located -35 kb upstream of HLA-C and is in moderately high LD with rs67384697 (r 2 = 0.74). We investigated whether rs67384697 was associated with psoriasis by imputing the deletion genotype of all psoriasis cases and controls. This was made possible by the near perfect LD between HLA-C four digit classical alleles and presence or absence of the deletion [Supplementary Table 2 in [37]]. To confirm the validity of our imputation, we sequenced the region of the HLA-C 39UTR containing rs67384697 in a subset of our samples (n = 70) and found a concordance rate of 138/140 alleles (98.6%), indicating the imputation was robust. Using logistic regression and adjusting for sex, ancestry, and cohort, we found that deletion allele of rs67384697 was significantly associated with psoriasis (p = 1.02610 229 , OR = 1.72) (Table 5), again confirming the similarity between psoriasis patients and HIV-1 controllers. We found that the association of rs67384697 with psoriasis was largely driven by HLA-C*06:02, since conditioning on HLA-C*06:02 resulted in only a marginally significant p-value for rs67384697 (p = 0.044, OR = 1.12). We note, however, that among all HLA-C allotypes, HLA-C*06:02 shows the highest level of cell surface expression, which could explain, to some extent, its strong association with psoriasis.   [15]. (C-G) Comparison of odds ratios to five other autoimmune or inflammatory diseases from the Wellcome Trust Case-Control Consortium. None of these demonstrate the same degree of similarity as between psoriasis and HIV-1 control. To demonstrate that this similarity cannot be entirely explained by the association of psoriasis with B*57:01 and C*06:02, we conditioned each of these residues on both B*57:01 and C*06:02 and found that Cys 67 , Ser 67 , Lys 70 , Asn 97 , Arg 97 in HLA-B remained independently associated with psoriasis (all p,5610 24 ) in the same direction as the association with HIV-1 control. This suggests that the repertoire of peptides bound at HLA-B may be similar between psoriasis patients and HIV-1 controllers. doi:10.1371/journal.pgen.1002514.g001 Our HLA data revealed that psoriasis is associated with HLA-B alleles carrying the Bw4 epitope (p = 1.28610 225 , OR = 1.66, Table S4). The association was stronger for Bw4-80I [isoleucine] (p = 8.28610 222 , OR = 1.80) than for Bw4-80T [threonine] (p = 3.41610 24 , OR = 1.22), which is interesting because Bw4-80I is thought to have a higher binding affinity for its KIR receptor than Bw4-80T [44]. We therefore hypothesized that psoriasis susceptibility might be mediated through activation of NK cells through KIR3DS1 and its putative partner HLA-B Bw4-80I. We genotyped KIR3DS1 in a subset of our psoriasis samples (n = 397) and compared the results to a healthy control cohort with available KIR3DS1 and HLA genotypes (n = 282). We found that the presence of the compound genotype KIR3DS1+Bw4-80I was a strong risk factor for psoriasis (frequency 22.7% in cases vs 6.9% in controls, p = 1.54610 27 , OR = 3.92, Table 6). Individuals positive for KIR3DS1 but lacking Bw4-80I had no increased risk for psoriasis (p = 0.63, OR = 0.91), and individuals positive for Bw4-80I but lacking KIR3DS1 had only a borderline increased risk of psoriasis (p = 0.058, OR = 1.53). To our knowledge, this is the first report that the compound genotype KIR3DS1+Bw4-80I is a strong risk factor for psoriasis susceptibility. This finding is again consistent with our observation that there is significant genetic similarity between psoriasis patients and HIV-1 controllers; however, replication of the KIR3DS1+Bw4-80I association in additional psoriasis cohorts is warranted.

Discussion
In this study, we followed up on the observation that several of the top SNPs from genome-wide association studies of psoriasis were identical to the top SNPs from genome wide association studies of HIV-1 control. Using imputation of HLA alleles, we found that psoriasis patients are enriched for several of the most significant known genetic variants associated with HIV-1 control: HLA-B*57 and HLA-B*27, which are associated with decreased viral load and increased time to AIDS [14,15,[24][25][26][27]; specific amino acid residues at HLA-B positions 97, 67, and 70 that are strong markers of HIV-1 controller status and viral load [15]; the deletion SNP rs67384697 which is associated with decreased viral load independent of HLA-B*57 and HLA-B*27 [37]; and the activating KIR3DS1 allele in combination with HLA-B Bw4-80I [42]. Psoriasis patients also demonstrate a significant paucity of HLA alleles and variants associated with HIV-1 disease progression [15,28], including HLA-B*35 (especially B*35-Px), B*07:02, B*40, C*04:01, C*07, and tyrosine 116 in HLA-B associated with sub-optimal peptide loading ( Table S2, Table S4). These effects were consistent between the 3 psoriasis cohorts examined in this study, demonstrating that the effects observed were real (Table  S6).
An important question to address, however, is whether the structural similarity between HLA alleles in psoriasis and HIV-1 control reflects the same underlying causal variants, or merely a coincidental association due to linkage disequilibrium. Our data suggest that some, but not all, of the observed similarity can be attributed to linkage disequilibrium. Our haplotype analysis shows that both psoriasis patients and HIV-1 controllers are enriched for the same extended haplotype, B*57:01-C*06:02-DQA1*02:01-DQB1*03:03-DRB1*07:01. In the HIV-1 controllers, this haplotype is likely primarily driven by selection for B*57:01, since previous studies have shown that the association of C*06:02 with HIV-1 control is dependent on B*57:01 in Europeans [15] and B*5801 in Africans [45], although one indirect benefit of C*06:02 for HIV-1 control is its high LD (D9 = 1) with the rs67384697 deletion polymorphism. In psoriasis, the haplotype association appears to be driven more by C*06:02 than B*57:01, since rs67384697, a G/del SNP located in the 39UTR of HLA-C, has been shown to modulate HLA-C expression levels by affecting the binding of microRNA hsa-miR-148 to its target site [37]. rs67384697 is associated with HIV-1 control independent of the classical HLA-alleles HLA-B*5701 and HLA-B*27:05. doi:10.1371/journal.pgen.1002514.t005  (Table 3). In addition, the association of the deletion allele of rs67384697 with psoriasis appears to be largely driven by LD with C*06:02. However, it should be noted that variants in high LD with C*0602 may be contributing to the observed association signal for C*0602. Interestingly, one can take the association signal for C*06:02 in psoriasis and perform stepwise conditioning on all coding amino acids within C*06:02 to demonstrate that the coding residues do not account for the entire association signal (Table S7). Therefore, the association of C*06:02 with psoriasis reflects, in part, other variants in high LD with C*06:02.
Despite the effects of linkage disequilibrium, our data suggest that some HIV-1 control variants indeed contribute independently to psoriasis susceptibility.  Table 2). A previous analysis of the HLA region in psoriasis by Feng et al. [46] also found that B*57:01 was independent of C*06:02 in Caucasians; moreover, in this study B*57:01 was also found to be independent of C*06:02 in a Chinese psoriasis cohort, which is notable because the LD between C*06:02 and B*57:01 is lower in Asians (D9 = 0.41) compared to Europeans (D9 = 0.90) [47]. Prior studies have also shown that B*27 is a strong risk factor for psoriasis in the subset of patients with psoriatic arthritis, especially those with axial disease [48][49][50][51]. Second, linkage disequilibrium with C*06:02 does not explain the lower frequency of the HIV-1 progression allele B*35 in psoriasis, nor can it account for the concordance of amino acid residues at HLA-B positions 67, 70, and 97 whose association with psoriasis was shown to be independent of HLA-B*57:01 and HLA-C*06:02 (Cys 67 , Ser 67 , Lys 70 , Asn 97 , Arg 97 , see Figure 1. legend). Finally, the provisional association of KIR3DS1+HLA-B Bw4-80I with psoriasis cannot be due to linkage disequilibrium, because KIR3DS1 is located on chromosome 19 which segregates independently of chromosome 6.
Although B*57:01, B*27:05, and possibly B*35 may have independent effects in psoriasis, additional studies are needed to clarify the precise mechanism(s) by which these and other psoriasis-associated HLA alleles contribute to psoriasis susceptibility or protection. The observation that psoriasis patients and HIV-1 controllers display concordant amino acids within the peptide binding groove of HLA-B suggests the possibility that an unknown psoriasis antigen shares homology with HIV-1 epitopes. An alternative possibility is that B*57:01, B*27:05, and B*35 do not restrict antigen presentation in psoriasis, but primarily function through their ability or inability to activate NK cells. We have provisionally shown a strong effect of KIR3DS1+Bw4-80I on psoriasis susceptibility, and B*57:01 contains the Bw4-80I epitope. The second strongest HLA allele in our stepwise regression model, B*38:01, also contains the Bw4-80I epitope. B*27:05 contains the Bw4-80T epitope, while protective alleles B*35 and B*40 contain the Bw6 epitope, which do not serve as ligands for KIR. Previous studies have shown that the activating KIR allele KIR2DS1 also contributes to psoriasis or psoriatic arthritis susceptibility [52][53][54][55], supporting the notion that NK cells may play a role in the pathogenesis of psoriasis. Finally, we have discussed the potential role of peptide processing on susceptibility to psoriasis, with the presence of tyrosine at HLA-B position 116 associated with protection against psoriasis, where position 116 is located near the C terminus of the bound peptide. A role for peptide processing influencing psoriasis risk has been previously identified for the gene ERAP1, an amino peptidase which regulates the quality of peptides bound to MHC class I molecules through trimming the peptide N terminus [10].
The genetic similarity between psoriasis patients and HIV-1 controllers has interesting implications. On a population level, the data would predict that Caucasian individuals with psoriasis are more likely than Caucasian individuals without psoriasis to be HIV-1 controllers, and HIV-1 controllers are more likely than non-controllers to develop psoriasis. This does not imply that every individual with psoriasis will be an HIV-1 controller, since only a fraction of psoriasis patients will harbor, for example, B*57, and even the presence of B*57 does not guarantee HIV-1 control, as this allele is present in some HIV-1 progressors. Nevertheless, one would expect an enrichment of HIV-1 controllers in the psoriasis population relative to a non-psoriatic population.
Our data also suggest a hypothesis that the existence of psoriasis may represent aberrant activation of evolutionarily-derived viral control alleles [56]. Barreiro et al. have shown that several of the top MHC SNPs associated with both psoriasis and HIV-1 control reside on haplotypes which show strong evidence of recent positive selection in the genome, as evidenced by long haplotypes indicative of rapid expansion of an advantageous allele in the population [57]. Psoriasis could subsequently result from the activation of viral control alleles due to the presence of a psoriasis antigen with sequence homology to HIV-1, or due to other environmental triggers. Although this study has focused on HLA and KIR alleles, other non-MHC psoriasis genes are plausibly associated with host response to viral infection. ERAP1 is involved in class I peptide processing and demonstrates epistasis with C*06:02 [10] and IFIH1 encodes a cytoplasmic helicase that mediates induction of interferon response to viral RNA [10,58]. TNIP1, TNFAIP3, TRAF3IP2, NFKBIA, and REL are associated with the TNF-a pathway and activation of NF-kB; while IL23R, IL12B, IL23A, and TYK2 are associated with activation of the Th17 pathway. Psoriasis is characterized by the upregulation of the cytokines IFN-a, IFN-c, TNF-a, IL-17, IL-22, and IL-23 [59], while TNF-a, IFN-c, and Th17+ T cells have been shown to be important in HIV-1 controllers [60][61][62].
The enrichment of viral control alleles in psoriasis patients may also help explain the psoriasis HIV-1 paradox (Figure 2). Psoriasis patients are more likely to harbor alleles such as HLA-B*57, HLA-B*27, HLA-C*06 (in high LD with the HLA-C 39UTR deletion polymorphism), and KIR3DS1+Bw4-80I, theoretically resulting in vigorous cytotoxic T cell and NK cell responses upon infection with HIV-1 virus. The pro-inflammatory environment created by these anti-viral responses, resulting in the production of cytokines such as TNF-a and IFN-c, would worsen the psoriasis. In addition, if the psoriasis antigen had sequence homology to HIV-1, then antigen specific immune responses directed against HIV-1 might cross-react with the psoriasis antigen and also flare the psoriasis. In either case, reduction of viral load and removal of the antigenic trigger through treatment with anti-retroviral therapies would improve the psoriasis, which is indeed seen clinically [17]. This explanatory model is consistent with several observational studies that patients with severe psoriasis and HIV-1 infection tend to carry the HLA-C*06 and HLA-B*27 alleles [63,64], because such alleles would trigger the vigorous immune response associated with exacerbation of psoriasis.
The data presented here with psoriasis and HIV-1 control illustrate the delicate balance of the human immune response, in which processes that safeguard the body against pathogens may also engage deleterious inflammatory responses. A similar example occurs with a genetic variant in the SH2B3 gene which may be protective against bacterial infection but which increases suscep-tibility to celiac disease, an autoimmune disease of the gut resulting from gluten intolerance [65]. Another example can be seen with the identification of genetic variants in immune function genes that increase the risk of sepsis, a systemic inflammatory response to infection which can lead to death [66,67].
In summary, using a large dataset of psoriasis cases and controls, we have shown that psoriasis patients and HIV-1 controllers share a high degree of similarity at their HLA loci. While some of this similarity is attributable to linkage disequilibrium, we present evidence that much of the similarity may be attributable to shared biological mechanisms including activation of natural killer cells, specificity of antigen presentation, and use of optimal MHC class I peptide processing. The genetic similarity between psoriasis and HIV-1 control suggests the possibility that psoriasis represents aberrant activation of pathways associated with anti-viral immunity. If this hypothesis is true, then the study of the biological pathways active in psoriasis may provide new therapeutic insights for the treatment of HIV-1.

Study Subjects
The study population and source are shown in Table S1. Two independent genome-wide association scan datasets were used as cohort 1 and cohort 2 in the present study. All cases and controls were of European descent. More details on subject characteristics and recruitment can be found in Liu et al. [8] and Nair et al. [9]. Only the subjects whose HLA alleles were successfully imputed (see below) were included in our analysis. In cohort 3, 169 psoriasis cases recruited from Washington University, St. Louis were directly typed for the class I HLA alleles by combining locus-specific amplification with hybridization of sequence-specific oligonucleotide probes as described in [68]; 1,711 control samples of European ancestry were obtained from studies 66 and 67 of illumina iControlDB. There was no overlap between the subjects among the three cohorts. Informed consent was obtained from each participant.

HIV Cohorts
Data generated by this study were primarily compared against a published genome-wide association study of HIV-1 control involving 516 HIV controllers of European ancestry (viral load ,2,000 RNA copies/ml by three measurements over at least 12 months without antiviral therapy) and 1,196 controls (treatmentnaïve chronically infected individuals with advanced disease, median viral load 61,698 copies/ml) [15]. HLA haplotype analysis was performed on 214 Caucasian HIV-1 infected individuals (52 HIV-1 controllers, 162 non-controllers) from the SCOPE cohort (Study of the Consequences of the Protease Inhibitor Era), whose HLA class I and II alleles had been previously directly genotyped. SCOPE HIV-1 controllers were antiretroviral therapy-naïve subjects who had at least one year duration of documented plasma HIV RNA below 2,000 copies/ml, while SCOPE noncontrollers were subjects who had at least one documented viral load above 10,000 copies/ml.

WTCCC Data
The Wellcome Trust Case-Control Consortium data were obtained from the WTCCC official website (http://www.wtccc. org.uk/). In this study, we used Affymetrix 500 K genotyping data from approximately 2,000 samples from each of five diseases (rheumatoid arthritis, Crohn's disease, type 1 diabetes, type 2 diabetes, and coronary artery disease) and 3,000 shared control samples from the 1958 birth cohort (58C) and the National Blood Service (NBS). More details about these samples are described elsewhere [33].

KIR Analysis Cohorts
KIR3DS1 typing was performed on 397 psoriasis subjects from cohorts 2 and 3 described above. Control HLA and KIR3DS1 data were obtained from 282 healthy Caucasian blood donors from the Carrington laboratory.

Imputation and Validation of Classical HLA Alleles
The program HLA*IMP [23] was used to impute HLA loci -A, -B, -C, -DQA1, -DQB1 and -DRB1 to 4-digit resolution in our genome-wide SNP datasets. Individuals or SNPs with a missing data frequency above 0.20 were excluded as recommended in the software manual. A call threshold of 0.7 on the modes of the posterior HLA type distributions was employed, which represents a good compromise between accuracy and call rate as suggested by the author. Imputation accuracy was assessed by comparing the imputed HLA alleles to directly typed HLA class I alleles in a subset of our samples (n = 98), comprising 42 samples imputed from Illumina 300 K SNP data and 56 samples imputed from Perlegen SNP data. We found that HLA*IMP produces highly accurate HLA type imputations at HLA class I loci at the 4-digit level. The concordance for the Illumina 300 K platform was 244/252 alleles (96.8%) and the concordance for the Perlegen platform was 322/ 329 alleles (97.9%), for an overall concordance rate of 566/581 alleles (97.4%). To examine the imputation accuracy of infrequent/ rare HLA alleles, we identified all HLA class I alleles with a population frequency of less than 5% in individuals of European descent, according to the online database: http://www.allelefrequencies.net. We then examined the accuracy of imputation at the 4 digit level for these infrequent/rare HLA alleles in our subjects for whom we had both directly genotyped HLA alleles and imputed HLA alleles. The concordance of HLA alleles with frequency ,5% was 95/102 alleles (93.1%) for the Illumina 300 K platform and 116/122 alleles (95.1%) for the Perlegen platform, for an overall concordance rate of 211/224 alleles (94.2%). However, our manuscript excludes HLA alleles with frequency less than 1%. For HLA alleles with frequency greater than 1% but less than 5%, the concordance was 78/79 alleles (98.7%) for Illumina 300 K and 99/100 alleles (99.0%) for Perlegen, for an overall concordance rate of 177/179 alleles (98.9%).

Association Testing and Adjustment for Covariates
Additive logistic regression models in PLINK [69] were used for most of the association tests, except for the HLA haplotype association tests. To account for potential population stratification or admixture in these samples, principal component analyses (PCA) was performed using the EIGENSTRAT [70]. Seven PCs in cohort 1 and ten PCs in cohort 2 were used for ancestry adjustment, based on leveling off of the PCA scree plot. The principal component score for each individual was included as a covariate in all models along with cohort and gender in logistic regression models. Multivariate logistic regression was performed in R software package (http://www.r-project.org/). To examine the consistency of association signals seen in the 3 cohorts used, a heterogeneity index was calculated using the meta-analysis module in PLINK.

Stepwise Regression Modeling of Independent HLA Alleles
Conditional and stepwise logistic regression was performed using the 'condition' function in PLINK to determine whether independent effects existed. The method begins with an empty model to which variables are added in an iterative process as described by Barcellos et al [71]. Briefly, starting with HLA-C*06:02 which exhibits the strongest association with psoriasis, we conditioned candidate HLA alleles on C*06:02 to determine the next most significant independent effect. For the model including both class I and class II alleles, the iterative process completes when no candidate allele demonstrates p,0.0006, which corresponds to the Bonferroni correction for the 88 HLA candidate alleles with MAF.1% tested. For the model including only class I alleles, the iterative process completes when no candidate allele demonstrates p,0.00096, which corresponds to the Bonferroni correction for the 52 class I candidate alleles with MAF.1% tested.

Association Testing of Amino Acid Positions (Including WTCCC Diseases)
The amino acid sequence of all HLA alleles is completely determined by the HLA type at four-digit resolution. We used the official amino acids sequences defined for known HLA alleles [31] and our imputed HLA allele data to determine the frequency of amino acid residues in cases and controls. HLA amino acids residues were tested for association using a logistic regression model that corrects for population substructure, gender and cohort using PLINK. For amino acid positions with .2 alleles, the omnibus test in the conditional haplotype analysis module in PLINK was used to determine a single p-value for all alleles at that position.

Stepwise Regression Modeling of Amino Acid Residues
We performed stepwise logistic regression in PLINK to determine the amino acid residues that were independently associated with psoriasis, using the same approach as done with the HLA alleles. For the amino acid analysis, the algorithm completes when no remaining candidate residue has p,0.0001, which corresponds to the Bonferroni correction for the 480 HLA amino acid residues with MAF.1% tested. For the stepwise regression modeling of HLA-C 06:02 association signal (Table S7), the algorithm completes when no candidate residue has p, 0.0006(0.05/87).

Imputation and Validation of rs67384697 Deletion Polymorphism
Since there is strong linkage disequilibrium between specific HLA-C alleles and rs67384697 [Supplementary Table 2 in [37]], we were able to determine the rs67384697 genotype of all subjects using their HLA-C four digit classical alleles. To ensure the accuracy of our imputation, we directly sequenced (ABI 3730 DNA analyzer, Quintara Biosciences, Berkeley, CA) the region of the HLA-C 39UTR containing rs67384697 in a subset of our samples (n = 70) and determined a high concordance rate of 138/140 alleles (98.6%). The following primers were used for sequencing of genomic DNA samples: forward 59-gtgagattctggggagctga and reverse 59-gaacagcaactaggcacagg as specified in [37].

Haplotype Analysis
Arlequin V3.5, based on the EM algorithm, was used to estimate the frequency of HLA allele haplotypes in our psoriasis cohorts and the SCOPE HIV cohort. To ensure the accuracy of haplotype construction, we compared the haplotype frequencies generated by Arlequin to the frequencies obtained by direct counting of the phased HLA alleles output from HLA*IMP, and found the two methods to yield nearly identical results. Haplotype frequencies were tested for statistically significant differences between case and control groups using the Chi Square test or Fisher's exact test in the R software package. cycles of 15 s at 94uC, 1 min at 55uC, 2 min at 72uC followed by a final 7 min extension step at 72uC. To confirm the accuracy of the results, samples were replicated using a second set of KIR3DS1 and KIR3DL1 primers from [73]. Phenotype frequencies for the presence of each gene were estimated by direct counting. Frequency differences between psoriasis and control groups were tested for significance by two-sided Fisher's exact test.

Supporting Information
Table S1 Study population and source. GAIN: The Genetic Association Information Network. WashU: Washington University in St. Louis. UCSF: University of California San Francisco. Illumina iControlDB: an online database of genotype and phenotype data from individuals that can be used as controls in association studies. (DOC)

Table S2
Association results for the imputed classical HLA alleles at 4-digit resolution. P values and ORs were adjusted for ancestry, gender, and cohort. The last two columns show the p values and ORs after conditioning on HLA-C*06:02. Only alleles with frequency greater than 1% in the control group were analyzed. (DOC)

Table S3
Association results for the imputed classical HLA alleles at 2-digit resolution. P values and ORs were adjusted for ancestry, gender, and cohort. The last two columns show the p values and ORs after conditioning on HLA-C*06. Only alleles with frequency greater than 1% in the control group were analyzed. (DOC)

Table S4
Association results for the imputed amino acid residues in each of the classical HLA loci. P values and ORs were adjusted for ancestry, gender, and cohort. Only residues with frequency greater than 1% in the control group were analyzed. (DOC)  Comparison of p-values and odds ratios (ORs) of the top HLA alleles, amino acid residues, and HLA-C 39 UTR deletion SNP in the 3 psoriasis cohorts used for this study. I 2 represents the proportion of variation that is due to heterogeneity. Roughly, I 2 values of 25, 50, and 75 indicate low, moderate, and high heterogeneity. Most variants demonstrate low heterogeneity with a consistent direction of association across the 3 cohorts. (DOC)

Table S7
Stepwise conditioning of HLA-C*06:02 association signal on the coding amino acids of HLA-C*06:02. The significant residual association indicates that the association of HLA-C*06:02 with psoriasis cannot be entirely explained by coding differences between HLA-C*06:02 and other HLA-C alleles. Therefore, the association of HLA-C*06:02 with psoriasis likely reflects the contribution of other variants in LD with HLA-C*06:02 such as HLA-B*57:01. (DOC)