Human Leukocyte Antigen (HLA) B27 Allotype-Specific Binding and Candidate Arthritogenic Peptides Revealed through Heuristic Clustering of Data-independent Acquisition Mass Spectrometry (DIA-MS) Data*
Expression of HLA-B27 is strongly associated with ankylosing spondylitis (AS) and other spondyloarthropathies. While this is true for the majority of HLA-B27 allotypes, HLA-B*27:06 and HLA-B*27:09 are not associated with AS. These two subtypes contain polymorphisms that are ideally positioned to influence the bound peptide repertoire. The existence of disease-inducing peptides (so-called arthritogenic peptides) has therefore been proposed that are exclusively presented by disease-associated HLA-B27 allotypes. However, we have recently demonstrated that this segregation of allotype-bound peptides is not the case and that many peptides that display sequence features predicted to favor binding to disease-associated subtypes are also capable of being presented naturally by protective alleles. To further probe more subtle quantitative changes in peptide presentation, we have used a combination of data-independent acquisition (DIA) and multiple reaction monitoring (MRM) mass spectrometry to quantify the abundance of 1646 HLA-B27 restricted peptides across the eight most frequent HLA-B27 allotypes (HLA-B*27:02-HLA-B*27:09). We utilized K means cluster analysis to group peptides with similar allelic binding preferences across the eight HLA-B27 allotypes, which enabled us to identify the most-stringent binding characteristics for each HLA-B27 allotype and further refined their existing consensus-binding motifs. Moreover, a thorough analysis of this quantitative dataset led to the identification of 26 peptides, which are presented in lower abundance by HLA-B*27:06 and HLA-B*27:09 compared with disease-associated HLA-B27 subtypes. Although these differences were observed to be very subtle, these 26 peptides might encompass the sought-after arthritogenic peptide(s).
Human leukocyte antigen (HLA) class I molecules present endogenous peptides derived from self-proteins or intracellular pathogens at the cell surface for scrutiny by CD8+ T lymphocytes. These peptide ligands are constrained within the peptide-binding cleft, which consists of six binding pockets (A-F) that accommodate the amino acid side chains of the peptides (1). Polymorphisms characterizing and differentiating HLA alleles and allotypes mostly cluster around these binding pockets, thereby shaping and defining their specific peptide cargo.
HLA-B27, which belongs to the family of HLA class I molecules, is an extremely polymorphic B-locus with more than 100 allelic variants known to date (2, 3). The most frequent allotypes with comprehensively described peptide repertoires are HLA-B*27:02–HLA-B*27:09 (4⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–18). Due to the architecture of their B pocket, the most dominant feature of all HLA-B27 bound peptides is the presence of an Arg residue at position P2 (19⇓–21). Additional primary and secondary anchor positions are P3, PΩ-2, and PΩ, which mostly carry hydrophobic and aromatic amino acid residues that bind into the D, E, and F pockets, respectively (18, 20, 21). All other peptide positions are considered as nonanchor positions and tolerate a large variety of amino acid residues with notable differences between the various allotypes (summarized in (18).
Ankylosing spondylitis (AS)1, a chronic inflammatory autoimmune disease, and other spondyloarthropathies are stronglylinked to the expression of HLA-B27 (reviewed in (22, 23). Even though this association has been known for more than 40 years (24, 25), the pathogenic property of HLA-B27 remains still unknown and many hypotheses are currently being investigated to elucidate its apparent role in AS (reviewed in (26, 27). Due to the peptide-presenting function of HLA-B27 and other HLA molecules, one of the most popular hypotheses is the arthritogenic peptide theory. This theory assumes that cytotoxic T lymphocytes (CTL) primed against pathogenic HLA-B27-restricted antigen(s) cross-react on molecularly similar self-peptide(s), which subsequently leads to tissue damage and onset of autoimmunity. In support of this theory, onset of AS is often preceded by infection with enteric bacteria and HLA-B27-restricted CTLs have been identified in the synovial fluid of AS patients, which recognize both bacterial epitopes and self-peptides (28, 29). Moreover, two allelic variants of HLA-B27 (HLA-B*27:06 and HLA-B*27:09), which contain polymorphisms ideally placed to alter the presented peptide repertoire, are not associated with AS (30⇓⇓⇓–34).
Based on this apparent protective role of HLA-B*27:06 and HLA-B*27:09 in disease pathogenesis, a number of analyses focused on comparing the peptide repertoires of disease-associated and nonassociated HLA-B27 subtypes with the goal to identify putative arthritogenic peptides, which were expected to be qualitatively absent from nondisease-associated subtypes (11⇓⇓–14, 17, 18). However, such a peptide has not been described yet, and a peptide feature that stratifies disease- and nondisease-associated subtypes does not appear to exist. In agreement of this notion, we could recently show that a peptide, which is capable of binding to each disease-associated HLA-B27 subtype, will also bind to the nondisease-associated allotypes HLA-B*27:06 and HLA-B*27:09 (18). This observation challenged the arthritogenic peptide theory in such that absolute binding preferences do not explain disease association. However, quantitative changes in peptide presentation between disease-associated and nonassociated subtypes may be relevant for disease pathogenesis by exceeding (or falling below) a threshold required for the activation of autoreactive T cells and, thus, for the onset of autoimmunity. Therefore, we sought to quantify the abundance of known HLA-B27 specific peptides across the eight most frequent HLA-B27 allotypes (HLA-B*27:02–HLA-B*27:09) with the goal to specifically identify candidate arthritogenic peptides that are expected to be presented in lower abundance by HLA-B*27:06 and HLA-B*27:09 compared with disease-associated HLA-B27 subtypes.
Data-independent acquisition (DIA) mass spectrometry—often referred to as SWATH-MS on SCIEX TripleTOF® instruments—is the state-of-the-art method to quantify thousands of peptides across multiple samples (35, 36). A mass spectrometer operated in DIA/SWATH mode uses wide ms1 isolation windows (typically 20–25 Da) to scan sequentially and repeatedly through a given mass range. All precursor ions isolated within each window are simultaneously fragmented in the collision cell resulting in multipeptide-derived ms2 spectra. Even though these spectra are too complex to be interrogated by conventional search engine algorithms, the continuous acquisition of all fragment ions over the chromatographic run time allows the extraction of quantifiable ms2 traces for virtually every peptide of interest. However, a faithful and accurate evaluation of DIA/SWATH-MS data requires prior knowledge of spectral and chromatographic information for all peptides of interest, although significant efforts are currently being made to develop software packages that extract sequence information directly from DIA/SWATH data (for example DIA Umpire: (37)).
Despite unique and compelling advantages of DIA/SWATH technology, multiple reaction monitoring (MRM)-MS still excels in unparalleled sensitivity (38, 39). This level of sensitivity is achieved by analyzing only one fragment ion (transition) at a time on a triple quadrupole mass spectrometer, which combines excellent ion transmission efficiencies with dual ion selectivity processes to significantly reduce background noise in favor of the analyte of interest. However, MRM-MS is severely restricted by a relatively low number of peptides that can be analyzed in one acquisition.
In this study, we used a combination of SWATH- and MRM-MS to quantify the abundance of 1646 epitopes across eight different HLA-B27 allotypes (HLA-B*27:02-HLA–B*27:09), including the non-AS associated subtypes HLA-B*27:06 and HLA-B*27:09. Computational algorithms were utilized to group all quantified peptides into distinct clusters with similar binding preferences across the HLA-B27 allotypes investigated. A detailed analysis of these cluster-specific consensus-binding motifs provided a means to attain the most stringent binding characteristics for each HLA-B27 allotype and further refined our understanding of the intricate mechanisms underlying HLA-peptide binding. Of note, a thorough comparison of this comprehensive quantitative dataset revealed 26 putative arthritogenic peptides, whose immunological potential will be analyzed in subsequent studies.
Previous SectionNext Section
Experimental Design and Statistical Rationale
We recently described the comprehensive analysis of the peptide repertoires of the eight most frequent HLA-B27 allotypes (HLA-B*27:02–HLA-B*27:09) by LC-MS/MS (18). In brief, HLA-B27 bound peptides were isolated from stable C1R transfectants by immunoaffinity capture, separated by reverse-phase high-performance liquid chromatography (RP-HPLC), and sequenced on a 5600+ TripleTOF® (SCIEX, Framingham, USA) mass spectrometer. At least four independent immunoprecipitations (biological replicates) were conducted to establish comprehensive immunopeptidomes for each HLA-B27 allomorph.
To obtain sufficient material for the eight SWATH-MS and 152 MRM-MS acquisitions presented in this study, the remaining aliquots of these subtype-specific peptide eluates were combined to form one biological pool for each of the eight distinct samples (HLA-B*27:02–HLA-B*27:09). By its nature, SWATH-MS eliminates technical variation allowing comprehensive interrogation of such pooled biological samples.
Data-Independent Acquisition (SWATH-MS)
Samples were loaded for 15 min with an isocratic flow (5 μl/min) of loading buffer (0.1% formic acid, 2% acetonitrile) onto a trap column (Eksigent, Dublin, USA; 200 μm × 0.5 mm ChromXP C18-CL 3 μm 120Å) using an Eksigent NanoUltra cHiPLC system. A microfluidic analytical column (Eksigent; 75 μm × 15 cm ChromXP C18-CL 3 μm 120Å) was switched in-line to separate the peptides by increasing concentrations of 80% acetonitrile /0.1% formic acid at a flow of 300 nl/min for 85 min.
The peptides were analyzed on a 5600+ TripleTOF mass spectrometer (SCIEX) operated in positive ion and SWATH-MS mode. In each cycle, a 200 ms precursor scan (range: 300–1800 Da) preceded the acquisition of 28 sequential SWATH-MS scans with an accumulation time of 100 ms each resulting in a total cycle time of 3 s. Each of the 28 SWATH-MS acquisitions scanned a 25 Da window between 300 to 1000 Da.
Raw SWATH data files have been deposited at PeptideAtlas (www.peptideatlas.org) with the dataset identifier PASS00763 (password: DT4743iu).
SWATH-MS data analysis
All SWATH-MS data were evaluated with the open source software Skyline (v126.96.36.19976; (40). The six most abundant b- or y- fragment ions of each peptide were selected for monitoring based on a spectral library, which contained information of all HLA-B27 peptides identified in our previous study (18). A custom-made retention time predictor based on retention time alignment of 15 intrinsic HLA-B27 peptides was used to correct retention time drift between analyses (Supplemental Fig. 1).
SWATH-MS data were imported into Skyline and the performance of the peak picking was visually assessed for each peptide and, if necessary, manually adjusted. A peptide identification was considered successful if all of the following criteria were met in at least one sample: (i) the deviation between the peptide's observed and predicted retention time was less than 5 min, (ii) the dot product, which delineates the similarity between the library spectrum and the observed transition chromatograms, was greater than 0.8, (iii) the mass error was smaller than 5 ppm, and (iv) the maximal peak height of a transition exceeded 500 arbitrary intensity units, which prevents quantification of low abundant background peaks.
A QTRAP® 5500 mass spectrometer (SCIEX) coupled to an Eksigent NanoUltra cHiPLC system was used for MRM detection as described previously (41, 42). In brief, the sample was loaded via a microfluidic trap column (Eksigent; 200 μm × 0.5 mm ChromXP C18-CL 3 μm 120Å) onto a microfluidic analytical column (Eksigent; 75 μm × 15 cm ChromXP C18-CL 3 μm 120Å), and the peptides were separated by increasing concentrations of 98% acetonitrile/0.1% formic acid at a flow rate of 300 nl/min for 90 min. The QTRAP® 5500 was operated in MRM mode (unit resolution for Q1 and Q3) analyzing each transition for 5 ms per cycle (dwell time). After each cycle, a single enhanced product ion scan (10,000 Da/s; rolling collision energy; unit resolution) was triggered to obtain sequence information of the most abundant precursor ion during this cycle.
Raw MRM data files have been deposited at PeptideAtlas (www.peptideatlas.org) with the dataset identifier PASS00762 (password: CS3537k).
MRM-MS Data Analysis
MRM data analysis and peptide quantification was done in Skyline (v188.8.131.5276; (40) using identical evaluation criteria as described for the SWATH-MS data analysis with the exception of the mass error criterion, which could not be used due to the modality of MRM-MS. Of note, instead of measuring the peptides of the iRT-B27 retention time predictor in each MRM-MS run, synthetic versions (GL Biochem, Shanghai, China) were separately analyzed in technical quadruplicates to determine the MRM-gradient-specific regression curve of the iRT-B27 peptides. The equation of this regression curve was applied to all MRM-MS runs for retention time alignment (Supplemental Fig. 2).
To correct for differences in sample density and loading, the peak areas for each peptide were normalized based on the assumption that a peptide with a universal HLA-B27 consensus-binding motif is presented in equal amounts across the eight HLA-B27 allotypes investigated. To minimize the intrinsic error of this assumption, individual normalization values were calculated for 428 universal HLA-B27 peptides and subsequently averaged (Supplemental Fig. 3; Supplemental Table 1). A peptide was considered a universal HLA-B27 peptide if it (i) was 9–12 amino acids in length, (ii) contained an Arg residue at P2, (iii) ended in a Phe, Leu, Ile, or Val residue, and (iv) did not carry a His, Arg, or Lys residue at P1 (18).
The k-means algorithm of the TM4:MeV software suite (http://www.tm4.org/mev.html; (43) was used to group the peptides according to their relative expression levels into 30 distinct clusters using the default parameters.
Previous SectionNext Section
Quantification of B27-Specific Peptides across the Eight Most Frequent HLA-B27 Allotypes
We utilized a combination of SWATH-MS, a DIA method accessible on TripleTOF® 5600+ instruments (SCIEX), and MRM-MS to comprehensively and accurately quantify HLA-B27-specific peptides across the eight most frequent HLA-B27 allotypes (HLA-B*27:02–HLA-B*27:09) including the nondisease-associated subtypes HLA-B*27:06 and HLA-B*27:09 (see Fig. 1 for experimental workflow). Equal amounts of subtype-specific HLA-B27 peptide eluates were pooled and analyzed on a TripleTOF® 5600+ mass spectrometer (SCIEX) operated in SWATH-MS mode. The data were evaluated in Skyline (40) based on a comprehensive library containing spectral information of peptides isolated from all HLA-B27 allotypes and a custom-made retention time predictor to align peak picking across each sample (iRT-B27 retention time predictor; Supplemental Fig. 1). In total, we analyzed 3306 peptides encompassing 3277 high confidence HLA-B27-bound peptides from our previous studies (18) and 29 selected HLA-Cw4-specific peptides, which are present in all the eluates due to low level expression of HLA-Cw4 on the parental C1R cell line (18, 44); Supplemental Table 2).
Workflow. We recently determined the peptide repertoires of the eight most frequent HLA-B27 allotypes (HLA-B*27:02 - HLA-B*27:09) by analyzing 44 independent DDA acquisitions on a TripleTOF 5600+ (SCIEX) mass spectrometer (18). Sequence information was obtained using the Paragon algorithm implemented in the ProteinPilot 4.5 software suite (SCIEX) and all ProteinPilot search result (.group) files were used to generate a spectral library in Skyline (MacLean et al. 2010). A custom-made retention time predictor (iRT-B27 retention time predictor) consisting of 15 intrinsic HLA-B27 and HLA-Cw4 peptides was created and used for retention time alignment (Supplemental Fig. 1). Residual, HLA-B27 allotype specific samples were pooled and analyzed by DIA/SWATH-MS on a TripleTOF 5600+ (SCIEX) mass spectrometer. Data were imported into Skyline and the performance of the peak picking was manually assessed for each peptide. Peptides, whose assignments failed our stringent evaluation criteria, were reanalyzed and reevaluated by MRM-MS. The quantitative data obtained by either SWATH-MS or MRM-MS were normalized for subsequent data analysis.
Because of the expected low abundance of certain peptides in some samples, the performance of the peak-picking algorithm was manually assessed for all peptides, and the peak boundaries were, if necessary, manually corrected to ensure faithful quantitative measurements. Only peptide assignments that fulfilled all of the following stringent evaluation criteria in at least one sample were considered for subsequent analyses: (i) the deviation between the observed and predicted retention time was less than 5 min, (ii) the dot product was greater than 0.8, (iii) the mass error was smaller than 5 ppm, and (iv) the maximal peak height exceeded 500 arbitrary intensity units. Using this strategy, quantitative measurements of 1095 HLA-B27-specific peptides were obtained across the eight different HLA-B27 subtypes and subsequently normalized to correct for differences in sample density and loading (Supplemental Table 1; Supplemental Table 3; Supplemental Fig. 3; Supplemental Fig. 4).
This comprehensive SWATH-MS dataset contained 92 peptides, which were recently quantified across all HLA-B27 allotypes by MRM-MS (18). To validate the accuracy of our SWATH-MS data, we compared the relative expression levels of these peptides to the published MRM-MS results (Supplemental Fig. 5A). Of note, 78 peptides (85%) showed an average variation of less than 20% across all HLA-B27 allotypes, confirming highly quantitative accuracy within the SWATH-MS dataset (Supplemental Fig. 5B). In addition, four peptides with both doubly and triply charged precursor ions were present in our SWATH-MS dataset. As expected, quantification of the differently charged precursor ions resulted in nearly identical relative expression levels across all HLA-B27 allotypes (Supplemental Fig. 5C).
In comparison to SWATH-MS, which results in MRM-like transition chromatograms, MRM-MS is a more focused and sensitive approach due to enhanced detection and use of target-optimized collision energies. Therefore, the 2211 peptides that did not pass our stringent SWATH-MS evaluation criteria were reanalyzed in 152 separate MRM-MS runs on a QTRAP® 5500 mass spectrometer (SCIEX) operated in MRM mode and evaluated in Skyline using the same criteria as described for the SWATH-MS acquisitions. Pursuing this strategy, quantitative data of an additional 551 HLA-B27-specific peptides was obtained (Supplemental Table 3: Quantification results).
Taken together, we were able to quantify the abundance of 1646 peptides across the eight most frequent HLA-B27 allotypes using a combination of SWATH- and MRM-MS.
Cluster Analyses Confirm and Extend the Allele Specific-Binding Motifs of the Different HLA-B27 Allotypes
We applied cluster analysis to group peptides with similar relative binding patterns across the different HLA-B27 allotypes. Based on adjusted “figure of merit” graphs (Supplemental Fig. 6), which propose appropriate input parameters for the k-Means algorithm (45), we opted to group the 1646 peptides into 30 distinct clusters, which appeared to be sufficient to adequately separate the peptides according to their allotype-specific binding patterns (Fig. 2; Supplemental Fig. 7). The resulting clusters contained between 17 and 128 peptides, and we determined the peptide length distributions and the consensus-binding motifs for each cluster (Supplemental Fig. 8, Supplemental Fig. 9).
Cluster analysis. The 1646 quantified peptides were grouped into 30 distinct clusters according to their relative expression levels across the eight HLA-B27 allotypes. The diagrams show the relative expression level of each individual peptide within a cluster (thin gray lines). The thick red lines indicate the average expression level of all peptides within a cluster and the number of peptides within a cluster is shown below the cluster number.
The first eight cluster groups contained peptides that are mainly presented by a single HLA-B27 allotype (Fig. 2;Supplemental Fig. 7). As expected, the binding motifs of these clusters were similar to previously determined canonical consensus-binding motifs (4⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–18). However, the consensus-binding motifs detailed in these studies are defined by analyzing all peptides eluted from a specific HLA-B27 allotype. Conversely, peptides within any of the first eight cluster groups are predominantly presented by only one HLA-B27 allotype and therefore, must contain feature(s) that are favored by this allotype but less tolerated by the remaining subtypes. Indeed, a thorough comparison of these cluster-specific binding motifs revealed several allotype-favored peptide features for each HLA-B27 subtype and provided a means to attain the most stringent binding characteristics for each allotype (Table I).
The remaining 22 clusters are composed of peptides that are abundantly presented by more than one HLA-B27 allotype (Fig. 2; Supplemental Fig. 7). Of note, the consensus-binding motifs of these clusters show a mixture of some but not all allotype-favored peptide features of the respective HLA-B27 subtypes (Table I). For example, cluster 29 represents peptides that are highly presented by HLA-B*27:03, HLA-B*27:04, and HLA-B*27:06. Their consensus-binding motif shows a preferred usage of His, Arg, or Lys at P1 (HLA-B*27:03-favored); Leu at PΩ-2 (HLA-B*27:04-favored); and Gly, Ser, or Ala at PΩ-1 (HLA-B*27:06-favored). These observations suggest that the combination of allotype-specific features present in a given peptide dictates (at least to some extent) its abundance levels across the different HLA-B27 allotypes. Moreover, these features can provide a means to estimate the relative abundance level of a given peptide bound by each HLA-B27 allotype.
Nested Sets of Peptides Suggest Central Bulging as the Preferred Mechanism to Accommodate Extended Peptides
Approximately 6% of all HLA-B27-specific peptides are part of nested sets of peptides, which differ only in their N- and/or C-terminal extensions but share an identical core sequence (15, 18). Our dataset contained quantitative information of 55 of such nested sets (Supplemental Table 4). Recent studies suggested that central bulging—rather than terminal protrusion—is the preferred mechanism to accommodate extended peptides into the binding cleft of HLA-B27 allotypes (15, 18). In support of this theory, we observed that most peptides within one nested set show significant differences in their abundance levels across the HLA-B27 allotypes. More importantly, their abundance level distribution often correlates with the presence of HLA-B27 allotype-specific feature(s) if these peptides would anchor with their termini into the binding cleft. For example, the peptide VRMNVLADAL is predominantly presented by HLA-B*27:08 (presumably because of the Ala residue at PΩ-1). The extended peptide VRMNVLADALK is virtually absent from HLA-B*27:08 but mainly presented by HLA-B*27:05 due to the HLA-B*27:05-favored C-terminal Lys. These observations confirm and extend the hypothesis that most peptides anchor with their termini into the binding cleft.
Identification of Candidate Arthritogenic Peptides
To identify possible arthritogenic peptides, we screened our quantitative dataset for peptides that are presented in less amounts by HLA-B*27:06 and HLA-B*27:09 compared with disease-associated HLA-B27 subtypes. A total of 26 peptides fulfilled this criterion (Table II; Supplemental Table 5). The majority of these peptides carry a Tyr or Phe residue at their C terminus but do not contain other obvious similarities in addition to the conserved Arg at P2. To assess the degree by which these peptides are less presented by HLA-B*27:06 and HLA-B*27:09, we calculated for each peptide the quantitative difference to the disease-associated subtype with the lowest relative expression level (Supplemental Fig. 10). On average, this difference was 3.9% (±4.4%) with individual values varying between less than 1% to more than 20% (Supplemental Fig. 10). This reflects the complexity of the dataset and is consistent with our view that qualitative differences are rarely seen between HLA-B27 allotypic peptide ligands but rather there are subtle to quite substantial quantitative changes in peptide abundance. To highlight this, we performed an identical calculation across all 1646 quantified peptides and obtained an average quantitative difference of 6.7% (±7.5%) between the three lowest relative expression levels on different HLA B27 allotypes. Only two peptides considered as candidate arthritogenic peptides (ARYVFQSENTF and ARVLLVPDNTF) exceeded this error region in such that the quantitative difference between HLA-B*27:06 or HLA-B*27:09 to a disease-associated subtype was significantly greater than the average quantitative difference between the three lowest relative expression levels of a random peptide (Supplemental Fig. 10). However, the noise from allotypic binding preferences almost certainly masks the candidate peptides and even subtle quantitative differences are sufficient to cause a differential T cell response leading to the onset and progression of AS. Therefore, we believe each of these 26 peptides might represent the potential arthritogenic self-peptide and should be analyzed in future studies.
Previous SectionNext Section
This study describes for the first time the quantification of >1500 epitopes across several closely related HLA class I allotypes (HLA-B*27:02–HLA-B*27:09) using a combination of SWATH- and MRM-MS. SWATH-MS is a novel DIA method enabled on TripleTOF® 5600+ mass spectrometers (SCIEX). Its potential lies in the ability to obtain quantitative MRM-like transitions for every peptide present in the sample and in the unique possibility to reinterrogate already acquired SWATH-MS data with a different set of peptides. However, due to the targeted nature of SWATH-MS, its data evaluation is strictly dependent on the availability of spectral and chromatographic information of all peptides of interest. This information is typically obtained by first analyzing the sample by standard DDA approaches in the presence of retention time alignment peptides (46). For this work, we generated such a spectral library in the open-source software Skyline (40) based on all HLA-B27 peptides identified in our previous study (18). Instead of spiking in peptides for retention time alignment, we exploited the presence of intrinsic peptides that were identified in each of the different HLA-B27 DDA acquisitions. Based on these peptides, we developed a custom-made “iRT-B27” retention time predictor (Supplemental Fig. 1) that allowed us to quantify 1095 peptides by SWATH-MS across all HLA-B27 allotypes.
MRM-MS - in comparison to SWATH-MS - is a more targeted approach and excels in unparalleled sensitivity but is limited by the number of peptides that can be analyzed in one run. We therefore opted to utilize MRM-MS to reanalyze only those peptides that were not identified by SWATH-MS. In addition, we did not measure the peptides of the iRT-B27 retention time predictor in each MRM-MS run but relied on a predetermined, MRM-gradient-specific regression curve of the iRT-B27 peptides for retention time alignment (correlation coefficient R2 = 0.9776; Supplemental Fig. 2). Pursuing this strategy, we were able to analyze 2211 peptides across eight samples in only 152 individual MRM-MS runs resulting in quantitative measurements of additional 551 peptides, which constitute ∼25% of the peptides that escaped quantification by SWATH-MS.
Quantifying peptide ligands across various HLA molecules is hampered by the low abundance of some shared HLA-associated peptides and that not all peptides are expected to be present in each sample. Despite significant advances in MS instrumentation during the last decade, faithful quantitative measurements are not possible for peptides whose abundance levels are very low and/or only marginally exceed the sensitivity threshold of detection. In those cases, the analysis software might pick an incorrect peak or an unspecific signal in the background region, which inevitably leads to erroneous quantification measurements and the inability to unambiguously determine the presence (or absence) of a peptide in a given sample. To attenuate this issue, we manually assessed the performance of the peak-picking algorithm of all peptides and corrected the peak boundaries as necessary. Moreover, we applied strict evaluation criteria to our comprehensive datasets and annotated all peptides, which did not fulfill these criteria in some samples, accordingly in Supplemental Table 3.
Through the utility of the TM4:MeV software suite (43), we grouped all quantified peptides into 30 different clusters according to their relative binding to the eight HLA-B27 allotypes. Instead of accepting missing values in our cluster analysis, we opted to use all quantified values irrespective of whether a peptide assignment in a given sample passed our stringent evaluation criteria or not. This decision is based on the assumption that a peptide, which failed our evaluation criteria in a sample, is most likely either absent or of extremely low abundance. Either way, the quantified (background) signal is sufficiently low to not perturb the results of the cluster analysis.
The presence of clusters that contained peptides, which were exclusively presented by only one HLA-B27 allotype, enabled us to identify the most stringent binding characteristics for each HLA-B27 allotype (Table I). As expected, most of these allotype-favored peptide features are present in previously determined canonical consensus-binding motifs (4⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–18). However, not all features identified by canonical consensus-binding motifs represent allotype-favored peptide features, which highlights the importance of allotype-specific cluster analyses. Of note, allotype-favored peptide features also suggest the preferred HLA-B27 subtype of generally tolerated features. For example, all HLA-B27 allotypes are known to permit binding of peptides carrying a C-terminal Leu residue. However, it appears that HLA-B*27:06 is the most-suited subtype to accommodate this peptide feature.
The remaining 22 clusters, which contained peptides that are abundantly presented by more than one HLA-B27 allotype, showed a mixture of allotype-favored peptide features of the respective HLA-B27 subtypes (Table I). This observation suggests that (allotype-favored) peptide features represent a key element in determining the abundance levels of peptides across different HLA allotypes. However, not all individual peptides within these 22 clusters show appropriate peptide features implying that other factors such as (i) secondary structure formations, (ii) conformational changes, (iii) availability and prevalence of source proteins, and (iv) differential interactions with components of the peptide-loading complex act in concert to regulate peptide binding and abundance. Thus, it is not surprising that several individual peptides within these 22 clusters show unexpected primary sequences, which are not in agreement with the allotype-favored peptide features of the respective cluster. Nevertheless, this cluster analysis describes and highlights an easy, empirical, and universal approach to attain the most stringent binding characteristics for a given HLA molecule and bypasses the requirement for complex models and prediction algorithms, which often fail to produce unambiguous results. Unraveling the determinant selection of closely related HLA molecules will lead to a better understanding of a panel of clinically relevant topics such as allograft rejections, hypersensitivity reactions and other autoimmune diseases.
A detailed analysis of nested sets of N- and C-terminally extended peptides suggested that most peptides are able to anchor with their termini into the binding cleft of HLA-B27 subtypes resulting in central bulging rather than terminal protrusion. This notion, which also has been proposed previously (15, 18), is based on the observation that the abundance levels of peptides across the HLA-B27 subtypes only correlate with the presence of allotype-specific feature(s) if these peptides anchor with their termini into the binding cleft. While this is true for nearly all peptides that are extended by only one or two amino acids, some unusually long peptides (extensions of more than three amino acids resulting in peptides longer than 15 amino acids) appear to have similar abundance distribution patterns as their base peptide. This observation may suggest the usage of identical anchor residues and implies that the termini of these peptides stick out from the binding cleft resulting in terminal protrusions. This is particularly evident in case of the nested set surrounding the base peptide ARLALSPVPSH, which is mainly presented by HLA-B*27:05 (Supplemental Table 4). As expected, ARLALSPVPSHW, which carries a one amino acid C-terminal Trp extension, is predominantly presented by HLA-B*27:02. An additional Met extension at the C terminus (ARLALSPVPSHWM), which is a HLA-B*27:08-favored feature, completely abolishes presentation by HLA-B*27:02 in favor of HLA-B*27:08. All further C-terminally extended peptides (ARLALSPVPSHWMVA and ARLALSPVPSHWMVAL) show a nearly identical distribution pattern across the HLA-B27 allotypes as ARLALSPVPSHWM. Although additional experiments such as binding and refolding assays have to be carried out to support this notion, we speculate that the vast majority of peptides bind with their termini into the peptide-binding groove resulting in central bulging.
The main goal of this study was the identification of peptides that are less presented by HLA-B*27:06 and HLA-B*27:09 compared with disease-associated HLA-B27 subtypes. We identified 26 peptides fulfilling this criterion and therefore might encompass the sought-after arthritogenic peptide(s). In analyzing their characteristics, 16 peptides were found to carry either a Tyr or Trp residue at PΩ (Supplemental Table 4, Supplemental Fig. 10), which explains at least to some extent their low abundance in HLA-B*27:06 or HLA-B*27:09 eluates (7, 9, 14, 18, 47). Conversely, the amino acid sequences of the remaining 10 peptides show no obvious features that would imply an impaired binding to HLA-B*27:06 or HLA-B*27:09. This observation suggests that candidate arthritogenic peptides cannot be predicted solely based on primary sequences and consensus-binding motifs. In fact, direct quantitative analyses—as presented in this study—are necessary and indispensable for their identification.
In summary, we utilized a combination of targeted mass spectrometry techniques to identify 26 candidate arthritogenic peptides whose immunogenic potential will be analyzed in further immunological and structural studies.
National Health and Medical Research Council http://dx.doi.org/10.13039/50110000092510225091085017
Previous SectionNext Section
Author contributions: R.B.S. and A.W.P. designed the research; R.B.S., S.S., T.C.L., and N.P.C. performed the research; R.B.S. and A.W.P. analyzed data; and R.B.S. and A.W.P. wrote the paper.
↵* This work was supported by project grants 1022509 and 1085017 from the National Health and Medical Research Council of Australia (NH&MRC). A.W.P. is an NH&MRC Senior Research Fellow. R.B.S. acknowledges support from the Swiss National Science Foundation.
↵ This article contains supplemental material.
↵1 The abbreviations used are:
- ankylosing spondylitis
- cytotoxic T lymphocytes
- data-dependent acquisition
- data-independent acquisition; t
- human leukocyte antigen
- multiple reaction monitoring
- reverse-phase high-performance liquid chromatography
- sequential window acquisition of all theoretical mass spectra.
- Received February 9, 2016.
- Revision received February 9, 2016.
- © 2016 by The American Society for Biochemistry and Molecular Biology, Inc.
HLA-B27 has a modifier effect on the phenotype of multiple diseases, both associated and non-associated with it. Among these effects, an increased frequency of clinical enthesitis in patients with Rheumatoid Arthritis (RA) has been reported but never explored again. We aimed to replicate this study with a sensitive and quantitative assessment of enthesitis by using standardized ultrasonography (US).
The Madrid Sonography Enthesitis Index (MASEI) was applied to the US assessment of 41 HLA-B27 positive and 41 matched HLA-B27 negative patients with longstanding RA. Clinical characteristics including explorations aimed to evaluate spondyloarthrtitis and laboratory tests were also done.
A significant degree of abnormalities in the entheses of the patients with RA were found, but the MASEI values, and each of its components including the Doppler signal, were similar in HLA-B27 positive and negative patients. An increase of the MASEI scores with age was identified. Differences in two clinical features were found: a lower prevalence of rheumatoid factor and a more common story of low back pain in the HLA-B27 positive patients than in the negative. The latter was accompanied by radiographic sacroiliitis in two HLA-B27 positive patients. No other differences were detected.
We have found that HLA-B27 positive patients with RA do not have more enthesitis as assessed with US than the patients lacking this HLA allele. However, HLA-B27 could be shaping the RA phenotype towards RF seronegativity and axial involvement.
Citation: Mera-Varela A, Ferreiro-Iglesias A, Perez-Pampin E, Porto-Silva M, Gómez-Reino JJ, Gonzalez A (2013) Ultrasonographic Assessment of Enthesitis in HLA-B27 Positive Patients with Rheumatoid Arthritis, a Matched Case-Only Study. PLoS ONE 8(3): e58616. https://doi.org/10.1371/journal.pone.0058616
Editor: Chi Zhang, University of Texas Southwestern Medical Center, United States of America
Received: November 2, 2012; Accepted: February 5, 2013; Published: March 7, 2013
Copyright: © 2013 Mera-Varela et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The study was supported by grants 10CSA918040PR from the Xunta de Galicia (http://www.sergas.e/MostrarContidos_N3_T01.aspx?IdPaxina=10142) and PI08/0744 of the Instituto de Salud Carlos III (http://www.isciii.es/) that are partially financed by the European Regional Development Fund of the European Union. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Research in the genetic component of RA etiology has experienced a marked progress in recent years with the identification of a large number of susceptibility loci , . However, there are still many unsolved questions. Among them the relative to genotype-phenotype relationships are very prominent. It is expected that genotypes will have an important role in shaping the RA phenotype and that identification of these relationships will help manage patients in a personalized manner. The most established relationship has been the observed between a large number of RA susceptibility factors and the production of anti-citrullinated protein antibodies (ACPA). The earlier discovered association of the shared epitope (SE) with severe RA, including progression of erosive arthritis, seems to be explained by SE association with ACPA production. Other lines of active research in this field include the search for genetic factors modifying progression of erosions , , or response to treatment ,  and of loci associated with ACPA negative patients , . Some of these studies are focused in RA susceptibility loci , , , which are good candidates because of their involvement in the disease mechanisms. Additional candidates are the genetic factors associated with related diseases. These loci are also of known functional relevance and they have the potential of being involved in subgroups of RA patients sharing clinical characteristics with the disease where the loci were identified.
One of these candidates is HLA-B27 , , , , . It is strongly associated with ankylosing spondylitis (AS) and with other spondyloarthritis (SpA), but not with RA , . HLA-B27 has also modifier effects on the phenotype that have been observed in its associated diseases and in some non-HLA-B27 associated diseases. Among the HLA-B27 associated diseases, there are multiple reports showing earlier disease onset, more involvement of sacroiliac and spine joints, a higher tendency to chronicity and more prevalence of back pain in HLA-B27 positive than in negative patients , , , , , . In addition, non-HLA-B27 associated diseases can be modified to show higher prevalence of archetypical AS clinical features. For example, HLA-B27 positive patients with inflammatory bowel disease show an increased prevalence of sacroiliitis, spondylitis and enthesitis , . HLA-B27 also increases the prevalence of sacroiliitis in patients with familial Mediterranean fever ; and of sacroiliitis, inflammatory back pain and enthesitis in children with Juvenile Idiopatic Arthritis , . Therefore, it is not surprising that some researchers have thought the presence of HLA-B27 could modify the phenotypes of patients with RA towards resembling SpA. There have been some reports in this direction , , , , , but all these studies were done more than a decade ago and none of their findings have become established. One of these studies reported an increased prevalence of clinical enthesitis in patients with early arthritis fulfilling RA criteria .
Enthesitis, the inflammation of the entheses, is very prevalent and relatively specific of all forms of SpA motivating its inclusion as part of the disease classification criteria , . Enthesitis is present in early SpA and, occasionally, it is the unique disease manifestation for some time. This has lead some authors to consider enthesitis the primary lesion in this group of diseases . Now we have the possibility to identify early enthesitis when it is not yet clinically apparent either to patients or to physical exploration thanks to magnetic resonance imaging and ultrasonography (US) , , , , . This property of US has been applied to a variety of clinical situations demanding sensitive evaluation. Examples include SpA with subclinical enthesopathy , recurrent acute anterior uveitis without SpA , long term dialysis , and patients with psoriasis without psoriatic arthritis , , showing in all these cases an increased frequency of abnormalities. US evaluation has also shown an increased frequency of abnormalities in patients with RA relative to healthy controls in spite of the complete absence of clinical enthesitis , . Only one of these studies assessed HLA-B27, concerning recurrent anterior uveitis, and it found that enthesitis was associated with this HLA allele .
Therefore, we have now the opportunity to analyze the prevalence of enthesitis in RA patients stratified for HLA-B27 status with the sensitivity and accuracy allowed by US. With this aim, we have selected all available RA patients positive for HLA-B27 in our hospital and matched RA patients negative for HLA-B27. The two groups were compared with the Madrid Sonography Enthesitis Index (MASEI) , focused anamnesis and physical exploration. No differences in the MASEI or any of its components were detected. This result excludes a role of HLA-B27 in the enthesal abnormalities of these patients. However, a modifier effect of HLA-B27 on the phenotype of the patients with RA was suggested by its association with RF negative status and a more frequent story of low back pain.
Materials and Methods
All patients gave their written informed consent to participate. Sample collection and the study protocol were approved by the Comite de Investigacion Clinica de Galicia (Spain).
Selection of patients
DNA and serum samples from patients with RA according to ACR classification criteria were used . All patients were of Spanish ancestry and were attending the Rheumatology Unit of the Hospital Clinico Universitario de Santiago. HLA-B27 positive patients were selected. Gender, age at disease onset and current age of these patients was used to select matched HLA-B27 negative patients at a 1∶1 ratio. Patients were invited to participate in the study and they gave their written informed consent. The rheumatologist and the nurse involved in recruitment and evaluation of the patients were blind to their HLA-B27 status.
Two complementary genotyping reactions were used to determine the HLA-B27 status (available from the authors upon request). The first PCR was based on the method of Bon et al. . It amplifies HLA-B*2701-*2706 alleles, which are the most common in the Caucasian population and the most associated with AS. The second PCR used a combination of 3 primers described by Faner et al. . It amplifies a larger number of alleles than the first, from HLA-B*2701 to *2724 except for HLA-B*2718 and HLA-B*2723. PCR amplification was done with a multiplex PCR system (KAPA2G fast HotStart, Kapa Biosystems, Woburn, MA). Detection of the PCR products was performed with single-base extension (SNaPshot Multiplex Kit from Applied Biosystems, Foster City, CA) using probes targeting a non polymorphic base. In this way, the most frequents alleles of HLA-B27 were covered by two PCR. Four samples were positive only in the second amplification, which is consistent with the low frequency of the rare alleles in Spain .They were excluded from analysis because of their more uncertain status.
HLA-DBR1 alleles were determined by a sequencing based typing method (SBT) using the AlleleSEQR HLA-DRB1 Typing kit (Abbott Diagnostics, Abbott Park, Germany), which includes bidirectional sequencing of the second exon of DRB1. Ambiguous samples were additionally sequenced with group-specific primers (AlleleSEQR HLA-DRB1 GSSP, Abbott). The anti-CCP status of the patients was determined using the EDIA ACPA Kit (Euro-Diagnostica, Arnhem, The Netherlands). Quantification and setting of the cut-off level at 5 units/ml were done according to the manufacturer's instructions. RF status was determined by rate nephelometry with the IMMAGE Immunochemistry Systems (Beckman Coulter, Ireland), which covers all Ig isotypes.
Clinical and Ultrasound evaluation
Patient history was reviewed and a specific anamnesis of symptoms and signs characteristic of HLA-B27 associated diseases was conducted. Physical exploration was done looking for signs of axial involvement with instruments developed for SpA as the Schober test, BASDAI and BASFI scores and assessment of the lateral trunk flexion. A modified Schober test according with Moll and Wright was used. Inflammatory back pain was defined as lumbar pain at night or with morning stiffness that does not improve with rest or that improves with exercise and that persists for more than a month. Plain antero-posterior pelvic radiographs were evaluated for the presence of sacroiliitis. This assessment was independently done by two rheumatologists in a blind form.
In addition, the presence of enthesitis was evaluated by a rheumatologist experienced in US focused on articular and periarticular locations and blind to the HLA-B27 status of patients. A General Electric LogiqQ7 US machine with a 10–14 MHz linear array transducer was used. The power Doppler setting was standardized with a pulse repetition frequency (PRF) of 500 Hz with a low wall filter and 35–40 dB of gain. The Madrid Sonography Enthesitis Index (MASEI) was used for quantitative and standardized assessment . This index evaluates 6 features in 6 entheses. A value ≥ 18 has been defined as characteristic of SpA. The 6 features are enthesis thickness, structure, calcifications, erosions, bursae and power Doppler signal. The 6 bilaterally assessed locations are proximal plantar fascia, distal Achilles tendon, distal and proximal patellar tendon insertion, distal quadriceps tendon and brachial triceps tendon.
HLA-B27 positive patients with RA were considered cases and the matched HLA-B27 negative patients were controls. Demographic, clinical, laboratory, radiographic and ultrasound characteristics were compared between the two groups using Student T test or the Fisher exact test for contingency tables depending on the quantitative or qualitative nature of the variables, respectively. Detailed US results were compared with the Man-Whitney U test and non-parametric statistics because they showed many zero values. All analyses were done with Statistica 7.0 (StatSoft, Tulsa, OK). Differences with P<0.05 were considered statistically significant.
Analysis of 672 patients with RA identified 65 that were HLA-B27+, the remaining 607 were HLA-B27- (11 additional patients showed an uncertain HLA-B27 status and were not included). The HLA-B27 positive subgroup of patients showed less RF positive subjects, 50.0 %, than the HLA-B27 negative subgroup 64.6 % (P = 0.025). Also the ACPA status showed a trend to decreased prevalence in this subgroup (Table 1). However, there were not more patients with low titers of ACPA (between 5 and 45 units) in the HLA-B27 positive (33.3 %) than in the HLA-B27 negative subgroup (32.8 %). The other characteristics, gender, age at disease onset, erosive arthritis and carrier status of the SE, were similar in the two subgroups of patients (Table 1).
A total of 41 HLA-B27 positive patients were available and willing to participate in the study. One HLA-B27- patient was selected for each HLA-B27+ patient trying to match them for gender, age at disease onset and current age. The resulting groups were very similar not only in the selected variables but also in other characteristics as height, weight, serology and SE status (Table 2). During the visit for evaluation one of the patients in the HLA-B27 positive subgroup was found to have SpA in spite of his previous classification as RA and was excluded from further analysis. Specific anamnesis and physical evaluation of the patients were done blind to their HLA-B27 status. They disclosed a higher prevalence of low back pain, mechanical or inflammatory, in the HLA-B27 positive subgroup, 27.5 %, than in the HLA-B27 negative subgroup, 7.5 %, P = 0.037. This difference was similarly distributed between inflammatory and mechanical pain (Table 3). No other symptom or evaluation showed differences between the two patient groups. Specifically there were not differences in modified Schober's test, BASDAI or BASFI scores or in the lateral flexion of the spine (Table 3). A few patients reported a story of past skin lesions that could correspond to psoriasis but without differences between the two subgroups.
Two rheumatologists that were blind to the HLA-B27 status of the patients evaluated plain pelvic radiographs for the presence of sacroiliitis. Their assessment was fully concordant identifying 2 patients with sacroiliitis. The two were positive for HLA-B27. This was not significantly different from the result in the HLA-B27 negative subgroup, but we investigated it further. A review of the two patients with sacroiliitis did not challenge their classification as RA: the two showed erosive RA and one of them was positive for ACPA and homozygous for the SE. The two patients had a story of inflammatory low back pain.
Systematic evaluation of the entheses following the MASEI procedure (representative images in Figure 1) yielded no differences between the two subgroups of patients (Table 4). Mean MASEI values were even nominally lower in the HLA-B27 positive subgroup. Also, the threshold score of 18 proposed as specific of AS , was not discriminating between the two subgroups: there were 13 patients above this value in the HLA-B27 positive subgroup and 10 in the HLA-B27 negative subgroup. In addition, no difference was detected with any of the components of the index (Table 4). This is specially relevant for the power Doppler signal (Figure 2) because it has been differentially associated with SpA relative to other diseases including RA . But no difference was detected comparing the whole distribution of values (Table 4), or the percentage of positive signals (9/40 in HLA-B-27 positive patients vs. 13/41 in HLA-B27 negative patients) or the mean values in the patients showing Doppler signals (4.7 vs. 4.6, in the HLA-B27 positive and negative patients, respectively). The only associations we identified were unrelated with the HLA-B27 status: a significant increase of the MASEI value in men relative to women (not shown), and an increase of the score with age due to calcifications (not shown). The two patients with sacroiliitis lacked signs of enthesitis (MASEI values of 4 and 8).
Figure 1. Representative images of features detected with US exploration.
A) Example of analysis of the Achilles tendon thickness measured between the two yellow crosses in one patient, and B) Erosion detected in the superior pole of the calcaneous in a different patient.
Figure 2. Positive power Doppler signal identifying tibial tuberosity enthesitis.
The signal (in red) is detected in the tibial insertion of the patellar ligament in one of the studied patients.
Our analyses have not shown any specific enthesal abnormality in the HLA-B27 positive patients with RA. The sensitive and quantitative evaluation allows us to exclude a significant HLA-B27 modifier effect towards enthesitis in the RA patients. However, an excess of patients in the HLA-B27 positive group referred a story of low back pain that in two of the patients was associated with radiographic sacroiliitis. These findings together with a lower prevalence of RF keep open the possibility of a modifier effect of the HLA-B27 allele on the phenotype of RA patients, but excluding the presence of enthesitis.
Our primary aim was to compare the frequency of enthesitis between HLA-B27 positive and negative patients with RA using systematic US examination. This imaging technology has shown sensitivity for detecting a high percentage of SpA patients with subclinical enthesitis , , , , , , , . It is applied with scoring protocols that help differentiate SpA patients from controls, from subjects that have suffered mechanical injury and from other forms of inflammatory arthritis , , . Specifically, the MASEI protocol used here has shown high sensitivity and specificity and it is very comprehensive because it includes assessment of 6 features in 6 enthesal sites and both grey scales and power Doppler , . The assessment of power Doppler is a distinct advantage of this method over the most commonly used GUESS procedure for our study because the Doppler signal is the most specific feature distinguishing SpA from mechanical and RA enthesal abnormalities .
We have found significant abnormalities in the entheses of patients with RA confirming findings of previous US studies , , . However, no differences were found between HLA-B27 positive and negative patients, even for components of the index that are more specific for SpA like the Doppler signal. Therefore, it is clear that the presence of HLA-B27 does not induce enthesal pathology in RA patients and that all the abnormalities found in these patients are produced with independence of this genetic factor.
Other interesting aspect of our MASEI results is that only patients older than 50 years of age showed values over the threshold identifying patients with AS or SpA , . Therefore, our results do not question the specificity of this threshold for early SpA, which starts most often below this age. An increase of enthesal abnormalities with age was already shown more than a decade ago . It was described as independent of the underlying disease and, as in our study, to be mostly due to calcification.
Two of the other clinical features we have analyzed were different between HLA-B27 positive and negative patients with RA. The first was the prevalence of RF, which was less common in HLA-B27 positive patients. It was accompanied by a trend to lower prevalence of ACPA. This result could imply either association of HLA-B27 with RF seronegative RA or misclassification of patients. The latter possibility can be reasonably excluded because we only found a patient with SpA among the 41 that were specifically revised. In addition, the HLA-B27+ RF- patients were not different in any respect from the HLA-B27- RF- subgroup. For example, the HLA-B27+RF- and B27-RF- patients were comparable at percentage of ACPA+, 26.7 % vs. 25.6 %, erosive arthritis, 51.7 % vs. 45.7 %, or percentage of carriers of SE, 44.4 % vs. 43.1 %, respectively. This leaves us with the possibility that HLA-B27 could be a susceptibility factor for RF seronegative RA. Only studies in additional sample collections will be needed to clarify this matter.
An increased prevalence of low back pain was also associated with HLA-B27 in the patients with RA. This outcome from the anamnesis was combined with the identification of two patients with radiographic sacroiliitis among the patients with inflammatory low back pain. The two results are in agreement with the idea that HLA-B27 could modify RA towards the axial joints. This idea is supported by previous studies showing association between sacroiliitis and HLA-B27, both in HLA-27-associated diseases , , , , and in some non-associated diseases , , ,  including RA , , . In fact, this is the unique association with HLA-B27 that has been replicated in different RA studies. However, two previous studies did not found differences in the sacroiliac joints between HLA-B27 positive and negative patients with RA , . These findings should motivate new studies directed to the analysis of low back pain and to the identification of sacroiliitis. MRI will be the technology of choice for these studies because it is very sensitive for incipient changes in the sacroiliac joints , and it could provide more definitive results than radiography given the small size of the HLA-B27 positive subgroup.
In relation with the low back pain association, it could be argued that its presence questions the classification of the patients as RA. However, this type of reasoning is incompatible with the aim of our study. When the aim is to identify patients with RA showing phenotypes resembling SpA, we need to consider classification of RA before doing any extra analysis and keep this classification constant along the study.
We have found that HLA-B27 positive patients with RA do not have more enthesitis as assessed with US than those lacking this allele. However, these patients referred a more prevalent story of low back pain and were more often seronegative for RF than the HLA-B27 negative patients indicating that HLA-B27 could be shaping the RA phenotype in other directions deserving further and more focused analysis.
Authors acknowledge the patients with RA participating in this study by their availability and generosity. AF-I has a pre-doctoral bursary from the Instituto de Salud Carlos III (Spain).
Conceived and designed the experiments: AM-V AG. Performed the experiments: AM-V AF-I EP-P MP-S JJG-R. Analyzed the data: AM-V AF-I AG. Wrote the paper: AM-V AF-I AG.
- 1. Gregersen PK (2010) Susceptibility genes for rheumatoid arthritis - a rapidly expanding harvest. Bull NYU Hosp Jt Dis 68: 179–182.
- 2. Plenge RM, Raychaudhuri S (2010) Leveraging human genetics to develop future therapeutic strategies in rheumatoid arthritis. Rheum Dis Clin North Am 36: 259–270.
- 3. Knevel R, Krabben A, Brouwer E, Posthumus MD, Wilson AG, et al. (2012) Genetic variants in IL15 associate with progression of joint destruction in rheumatoid arthritis: a multicohort study. Ann Rheum Dis
- 4. van der Linden MP, Feitsma AL, le Cessie S, Kern M, Olsson LM, et al. (2009) Association of a single-nucleotide polymorphism in CD40 with the rate of joint destruction in rheumatoid arthritis. Arthritis Rheum 60: 2242–2247.
- 5. Plant D, Bowes J, Potter C, Hyrich KL, Morgan AW, et al. (2011) Genome-wide association study of genetic predictors of anti-tumor necrosis factor treatment efficacy in rheumatoid arthritis identifies associations with polymorphisms at seven loci. Arthritis Rheum 63: 645–653.
- 6. Cui J, Saevarsdottir S, Thomson B, Padyukov L, van der Helm-van Mil AH, et al. (2010) Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy. Arthritis Rheum 62: 1849–1861.
- 7. Padyukov L, Seielstad M, Ong RT, Ding B, Ronnelid J, et al. (2011) A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis. Ann Rheum Dis 70: 259–265.
- 8. Seddighzadeh M, Gonzalez A, Ding B, Ferreiro-Iglesias A, Gomez-Reino JJ, et al. (2012) Variants Within STAT Genes Reveal Association with Anticitrullinated Protein Antibody-negative Rheumatoid Arthritis in 2 European Populations. J Rheumatol 39: 1509–1516.
- 9. El-Gabalawy HS, Goldbach-Mansky R, Smith D 2nd, Arayssi T, Bale S, et al. (1999) Association of HLA alleles and clinical features in patients with synovitis of recent onset. Arthritis Rheum 42: 1696–1705.
- 10. Jajic Z, Jajic I (1991) HLA-B27 antigen and rheumatoid arthritis. Acta Med Iugosl 45: 195–202.
- 11. Jaraquemada D, Ollier W, Awad J, Young A, Silman A, et al. (1986) HLA and rheumatoid arthritis: a combined analysis of 440 British patients. Ann Rheum Dis 45: 627–636.
- 12. Rantapaa Dahlqvist S, Nordmark LG, Bjelle A (1984) HLA-B27 and involvement of sacroiliac joints in rheumatoid arthritis. J Rheumatol 11: 27–32.
- 13. Rundback JH, Rosenberg ZS, Solomon G (1993) The radiographic features of rheumatoid arthritis in HLA-B27-positive patients. Skeletal Radiol 22: 263–267.
- 14. Reveille JD, Maganti RM (2009) Subtypes of HLA-B27: history and implications in the pathogenesis of ankylosing spondylitis. Adv Exp Med Biol 649: 159–176.
- 15. Dougados M, Baeten D (2011) Spondyloarthritis. Lancet 377: 2127–2137.
- 16. Feldtkeller E, Khan MA, van der Heijde D, van der Linden S, Braun J (2003) Age at disease onset and diagnosis delay in HLA-B27 negative vs. positive patients with ankylosing spondylitis. Rheumatol Int 23: 61–66.
- 17. Skare TL, Leite N, Bortoluzzo AB, Goncalves CR, Da Silva JA, et al. (2012) Effect of age at disease onset in the clinical profile of spondyloarthritis: a study of 1424 Brazilian patients. Clin Exp Rheumatol 30: 351–357.
- 18. Queiro R, Torre JC, Gonzalez S, Lopez-Larrea C, Tinture T, et al. (2003) HLA antigens may influence the age of onset of psoriasis and psoriatic arthritis. J Rheumatol 30: 505–507.
- 19. Chung HY, Machado P, van der Heijde D, D'Agostino MA, Dougados M (2011) HLA-B27 positive patients differ from HLA-B27 negative patients in clinical presentation and imaging: results from the DESIR cohort of patients with recent onset axial spondyloarthritis. Ann Rheum Dis 70: 1930–1936.
- 20. Kaarela K, Jantti JK, Kotaniemi KM (2009) Similarity between chronic reactive arthritis and ankylosing spondylitis.A 32-35-year follow-up study. Clin Exp Rheumatol 27: 325–328.
- 21. Queiro R, Sarasqueta C, Belzunegui J, Gonzalez C, Figueroa M, et al. (2002) Psoriatic spondyloarthropathy: a comparative study between HLA-B27 positive and HLA-B27 negative disease. Semin Arthritis Rheum 31: 413–418.
- 22. Rodriguez-Reyna TS, Martinez-Reyes C, Yamamoto-Furusho JK (2009) Rheumatic manifestations of inflammatory bowel disease. World J Gastroenterol 15: 5517–5524.
- 23. Orchard TR, Holt H, Bradbury L, Hammersma J, McNally E, et al. (2009) The prevalence, clinical features and association of HLA-B27 in sacroiliitis associated with established Crohn's disease. Aliment Pharmacol Ther 29: 193–197.
- 24. Kasifoglu T, Calisir C, Cansu DU, Korkmaz C (2009) The frequency of sacroiliitis in familial Mediterranean fever and the role of HLA-B27 and MEFV mutations in the development of sacroiliitis. Clin Rheumatol 28: 41–46.