For conditional association analysis, we considered the regression model including the additional HLA variants as covariates. When conditioning on specific HLA amino acid position(s), we included multiallelic variants of the amino acid residues as covariates. When conditioning on specific HLA gene(s), we included all two- and four-digit classical alleles of the HLA gene(s) (but not alleles with strong correlations [R2 > 0.97]). We consecutively selected the HLA variants to be included as covariates for each HLA gene separately in a forward-type stepwise fashion until no variant satisfied the genome-wide significant threshold (p < 5.0 × 10−8). We tested a multivariate full regression model by including the HLA-C, HLA-B, HLA-A, and HLA-DQA1 risk variants identified by the stepwise regression analysis as covariates and excluding the most frequent allele (or residue) from each locus (or amino acid position) as a reference allele (Table 1). Assuming a PsV prevalence of 2.0%, we estimated phenotypic variance explained by the risk HLA alleles and amino acid polymorphisms on the basis of the effect sizes obtained from the multivariate regression analysis and a liability threshold model.2