Abstract
Drug-induced hypersensitivity reactions can significantly impact drug development and use. Studies to understand risk factors for drug-induced hypersensitivity reactions have identified genetic association with specific human leukocyte antigen (HLA) alleles. Interestingly, drug-induced hypersensitivity reactions can occur in nonhuman primates; however, association between drug-induced hypersensitivity reactions and major histocompatibility complex (MHC) alleles has not been described. In this study, tissue samples were collected from 62 cynomolgus monkeys from preclinical studies in which 9 animals had evidence of drug-induced hypersensitivity reactions. Microsatellite analysis was used to determine MHC haplotypes for each animal. A total of 7 haplotypes and recombinant MHC haplotypes were observed, with distribution frequency comparable to known MHC I allele frequency in cynomolgus monkeys. Genetic association analysis identified alleles from the M3 haplotype of the MHC I B region (B*011:01, B*075:01, B*079:01, B*070:02, B*098:05, and B*165:01) to be significantly associated (χ2 test for trend, p < 0.05) with occurrence of drug-induced hypersensitivity reactions. Sequence similarity from alignment of alleles in the M3 haplotype B region and HLA alleles associated with drug-induced hypersensitivity reactions in humans was 86% to 93%. These data demonstrate that MHC alleles in cynomolgus monkeys are associated with drug-induced hypersensitivity reactions, similar to HLA alleles in humans.
The major histocompatibility complex (MHC) gene loci encode proteins that play a critical role in immune response. Human MHC genes are divided into 3 classes: class I and II which encode complexes that present endogenous and exogenous peptide sequences, respectively (Vyas, Van der Veen, and Ploegh 2008), and class III which encodes specialized immune function genes (Milner and Campbell 2001). The genes encoded by MHC class I and II are involved in antigen presentation such as human leukocyte antigen (HLA)-A and HLA-B (class I) and HLA-DR, HLA-DP, and HLA-DQ (class II) genes (Traherne 2008). These genes are highly polymorphic and experience genomic recombination (HLA-DRβ duplication) to introduce additional sequence variation, which impacts peptide-binding specificity and receptor interaction (Parham 1988; Rollini, Mach, and Gorski 1985). In the cynomolgus monkey (Macaca fascicularis), MHC genes have been extensively studied based on the diverse geographic locations that harbor cynomolgus monkeys (Campbell et al. 2009; Pendley et al. 2008; Creager et al. 2011; Krebs et al. 2005). MHC class I genes are defined as either Mafa-A or Mafa-B (HLA homologues). Interestingly, the MHC genes in these animals have undergone extensive genetic rearrangement and duplication when compared to human MHC gene loci (O’Connor et al. 2007; Doxiadis et al. 2006; Blancher et al. 2006). The result is thousands of nonhuman primate MHC haplotypes, with cynomolgus monkey and the closely related rhesus monkey sharing MHC class I and II haplotypes across species (Doxiadis et al. 2006; de Groot et al. 2012; Robinson et al. 2013). Despite this high degree of variation, similarities to human MHC genes persist, with human and cynomolgus monkeys sharing MHC epitopes across species (Mothe et al. 2015).
Specific groups of linked genetic variants in HLA genes have been defined as HLA* allele haplotypes and are associated with variable immune responses (Price et al. 1999; Marsh et al. 2010). HLA* alleles have also been associated with drug-induced hypersensitivity reactions for patients receiving therapies such as abacavir, flucloxacillin, carbamazepine, and allopurinol (Daly et al. 2009; Mallal et al. 2002; Martin et al. 2004; Hung et al. 2005, 2006). Skin reaction associated with drug administration is one of the most common adverse drug reactions (ADRs) observed in humans (Bigby et al. 1986; Pichler 2003; Roujeau 2005). ADRs can be predictable in humans when due to the drug’s pharmacologic action (type A reactions) or unpredictable when not related to the drug’s pharmacologic action (type B reactions; Rawlins and Thompson 1977). It has been reported that the vast majority of drug-induced skin reactions in humans tend to be type B reactions and are idiosyncratic in nature (Uetrecht and Naisbitt 2013).
Standard preclinical toxicology testing has generally poorly predicted clinical drug-induced hypersensitivity reactions (Olson et al. 2000). However, it was recently reported that following treatment with metabotropic glutamate receptor 5 (mGluR5) negative allosteric modulators, clinical skin lesions were observed in Mauritius cynomolgus monkeys (Palanisamy et al. 2015). Further, similar skin lesions were observed during preclinical development of a second drug candidate (Pfizer internal data), and in total, 9 of the 62 cynomolgus monkeys exhibited clinical and histopathologic skin effects that were consistent with characteristics of drug-induced skin reactions observed in humans. Clinical effects in the monkeys appeared between 5 and 8 days of dosing and ranged from generalized to localized erythema often involving the genital region (scrotum) in males and/or the head or limbs in either sex. Localized skin reactions in the genital region were often edematous and occasionally associated with bulla formation. With one of the mGluR5 compounds, a dosing holiday followed by rechallenge resulted in more rapid clinical development of skin effects (day 3), which is consistent with a hypersensitivity reaction. Microscopically, the lesions were characterized by lymphocytic inflammation at the dermal–epidermal interface (interface dermatitis) with additional perivascular lymphocytic inflammation in the dermis that was sometimes associated with edema. Interface inflammation was accompanied by single cell necrosis of keratinocytes, particularly along the basal epidermal cell layer (Figure 1); in some instances, the degenerating keratinocytes were surrounded by lymphocytes. Bullae from the genital region were characterized by complete separation of the entire epidermis from the dermis (subepidermal bullae), sometimes with full thickness necrosis of the overlying epidermis. Bullae were filled with amphophilic fluid containing small numbers of inflammatory cells (predominantly lymphocytes). Immunohistochemical characterization of the skin lesions observed in the cynomolgus monkeys given mGluR5 inhibitors demonstrated a T lymphocyte–dominant immunophenotype. Collectively, the nature of the clinical effects (localized to generalized erythema with or without bullae formation), microscopic lesions (interface dermatitis with basal keratinocyte necrosis), the immunophenotype (T lymphocyte dominant), and the time of initial clinical onset with earlier onset upon rechallenge are consistent with an immune-mediated delayed type-IV hypersensitivity reaction (Palanisamy et al. 2015). These drug-induced clinical and microscopic skin effects in the monkeys have several overlapping features with human forms of drug-induced skin reactions, suggesting a similar underlying mechanism. While the mechanism may be similar, how drugs cause skin hypersensitivity reactions in humans and monkeys is not well understood and hypotheses are that the drug can form a hapten or superantigen-like structure that results in an immune response (Uetrecht 2007).

Representative microscopic image (original objective 20×) of a drug-induced skin lesion in a cynomolgus monkey treated with a small molecule metabotropic glutamate receptor (mGluR) 5–negative allosteric modulator. The lesion is localized to the dermal–epidermal junction and characterized by lymphocytic interface inflammation (asterisks) accompanied by single cell necrosis of basal keratinocytes (arrows).
In this study, the association between MHC haplotypes and the skin hypersensitivity reactions observed in Mauritius cynomolgus monkeys was evaluated. The Mauritius cynomolgus monkey subpopulation provides a unique opportunity to study MHC haplotypes and their relationship to immune response in monkeys, as geographic isolation of a small founder population on the island of Mauritius has reduced MHC class I and class II sequence diversity to only 7 haplotypes (O’Connor et al. 2007; Wiseman et al. 2007; Krebs et al. 2005; Bonhomme et al. 2008). This study characterized these haplotypes in a study population of cynomolgus monkeys and identified that alleles in the M3 haplotype of the MHC I B region are associated with drug-induced hypersensitivity reactions in cynomolgus monkeys.
Materials and Methods
Tissue Samples
Peripheral blood or formalin-fixed, paraffin-embedded (FFPE) spleen samples obtained from 62 cynomolgus monkeys of Mauritian origin (Charles River Laboratories, Houston, TX) were selected from preclinical toxicology studies conducted for small molecule drug candidates from two independent and unrelated programs at Pfizer, Inc. (5 studies over approximately 2 years). All procedures performed on animals were in accordance with regulations and guidelines reviewed and approved by the Pfizer Institutional Animal Care and Use Committee. In total, 9 of the 62 samples collected were from monkeys that had clinical evidence of drug-induced skin reactions (localized or generalized erythema, edema and/or bullae formation), and 53 were from monkeys that had no observed clinical signs following administration of the same test article.
Genomic DNA Isolation and Microsatellite Analysis
Genomic DNA was isolated from FFPE tissues using QIAamp DNA FFPE Tissue Kit (Qiagen, Valencia, CA), from peripheral blood using DNeasy Blood and Tissue Kit (Qiagen), according to manufacturer’s instructions, and was used as a template for microsatellite analysis. Microsatellite analysis was performed by the University of Wisconsin using published methods (Wiseman et al. 2007). Microsatellite profiles for each animal were used to infer the MHC class I and class II transcripts previously defined for each of the 7 ancestral MHC haplotypes in the Mauritius cynomolgus monkey population (Wiseman et al. 2013).
Sequence Alignment
Sequences for MHC M3 B region alleles and human HLA alleles were aligned using Clustal W (http://bio.lundberg.gu.se/edu/msf1.html) and sequence similarity scores were reported.
Statistical Analysis
Due to the clinical and histological similarity of the drug-induced skin reactions (reviewed in Introduction) that occurred with small molecule drugs from two independent drug programs, the samples were combined for genetic association analysis in order to increase statistical power. Genetic association analysis was performed using the χ2 trend test. Significance was considered at p < 0.05.
Results
MHC Haplotype Identification and Distribution
Using microsatellite analysis, there were 7 haplotypes and recombinant MHC class I haplotypes observed in the 62 monkeys. The distribution frequencies of these haplotypes in this population are shown in Table 1. Overall, the MHC class I haplotype frequencies were very similar to what has been published for the general Mauritius cynomolgus monkey population (Budde et al. 2010). Haplotypes in the MHC class II region were also identified in the study population (Table 1), and the distribution was similarly comparable to the general Mauritius monkey population (O’Connor et al. 2007), although to a lesser extent. Notably, the M4 haplotype was relatively enriched in the study population. Also notable was that no recombination was detected within the MHC A, B, or class II regions in the 62 cynomolgus monkeys analyzed in this study (data not shown).
Haplotype Distribution of MHC I and II Alleles in the Study Population and the General Mauritius Cynomolgus Monkey Population.
Note. MHC = major histocompatibility complex.
aData from Budde et al. (2010).
bData from O’Connor et al. (2007).
cHaplotypes not determined by microsatellite analysis.
Genetic Association Analysis
There were 4 MHC I haplotypes (M2, M3, M4, and M6) identified in the 9 animals that experienced drug-induced hypersensitivity reactions (Figure 2), suggesting that the alleles from these 4 haplotypes are more likely involved in drug-induced hypersensitivity reactions. Genetic analysis was conducted using the alleles from the 4 MHC haplotypes to identify possible association between these haplotypes and drug-induced hypersensitivity reactions in cynomolgus monkeys. Interestingly, the alleles from the M3 haplotype of the MHC class I B region (B*011:01, B*075:01, B*079:01, B*070:02, B*098:05, and B*165:01) were found to be significantly associated with drug-induced hypersensitivity reactions (Table 2, χ2 test for trend, p <0.05). In contrast, the remaining 11 MHC haplotype regions (A, B, and class II regions in the M2, M4, and M6 haplotypes, and A and class II regions in the M3 haplotype) had no statistical association with drug-induced hypersensitivity reactions in the 62 cynomolgus monkeys analyzed.

Alleles of the 4 haplotypes identified in cynomolgus monkeys experiencing drug-induced hypersensitivity reactions. Microsatellite analysis was used to determine the major histocompatibility complex haplotypes in 62 cynomolgus monkeys, and 4 haplotypes were identified in the 9 monkeys with drug-induced hypersensitivity reactions: M2 (red), M3 (blue), M4 (green), and M6 (gray).
Results of Genetic Association Analysis with Alleles from MHC Haplotypes and Drug-induced Hypersensitivity Reactions in Cynomolgus Monkeys.
Note. MHC = major histocompatibility complex.
Homology between HLA Alleles Associated with Drug-induced Hypersensitivity Reactions and MHC I M3 B Alleles
The DNA sequences for HLA alleles (HLA-B*57:01, B*15:02, B*58:01 and B*35:05) that have been reported to be associated with drug-induced hypersensitivity reactions (Pompeu et al. 2012) were compared with the DNA sequences of the cynomolgus monkey MHC I M3 B region alleles that were genetically associated with the skin reactions observed in the monkeys. Sequence similarity ranged from 86% to 93% (Table 3).
Sequence Similarity Scores from Alignments between Alleles in MHC I M3 Haplotype B Region (Mafa-B) and HLA Alleles.
Note. HLA = human leukocyte antigen; MHC = major histocompatibility complex.
aHomology between cynomolgus monkey MHC class I M3 B region alleles (Mafa-B* alleles) and human HLA alleles were aligned using Clustal W (http://bio.lundberg.gu.se/edu/msf1.html).
bHLA allele known to be associated with drug-induced hypersensitivity reaction (Pompeu et al. 2012); only subtype :01 was selected for each HLA allele.
Discussion
Given the extensively studied and well-known associations between certain HLA alleles and drug-induced hypersensitivity reactions in humans (Sukasem et al. 2014; Yip, Alfirevic, and Pirmohamed 2015), it was hypothesized that MHC haplotypes in cynomolgus monkeys may also be associated with drug-induced hypersensitivity reactions in this species. Drug-induced hypersensitivity reactions were observed in cynomolgus monkeys with small molecule drug candidates from two independent and unrelated programs, and had similar clinical and histological characteristics suggesting that a common immune-mediated pathway (such as MHC signaling) downstream of or unrelated to either pharmacological mechanism of action was involved in the pathogenesis. Indeed, similar to the observations in humans, 6 alleles from the M3 haplotype of the MHC B region were identified to be significantly associated with the drug-induced hypersensitivity reactions observed in these monkeys.
Interestingly, unlike drug-induced hypersensitivity reactions in humans, which typically manifest as idiosyncratic (low incidence), the skin reactions in the cynomolgus monkeys displayed a clear dose–response relationship and high incidence (Palanisamy et al. 2015). Given that certain populations of cynomolgus monkeys have relatively limited MHC diversity (Aarnink et al. 2010), it was initially suspected that the M3 haplotype was overrepresented in the population of Mauritian monkeys utilized in these preclinical studies, making what would typically appear as an idiosyncratic pattern in a genetically diverse human population appear more dose related in the less diverse study population, which would in turn skew the genetic association analysis. However, very comparable distribution of MHC I haplotypes was observed between the study population and the general cynomolgus monkey population of Mauritian origin, including the M3 haplotype which had a distribution frequency of 17% in the study population and 16% in the general population (Budde et al. 2010), excluding the possibility that the observed significant association between the M3 haplotype and drug-induced hypersensitivity reactions resulted from a statistical sampling artifact. Furthermore, no recombination was detected within MHC A or B, or MHC class II regions in the study population, confirming the relatively limited genetic diversity of MHC molecules in cynomolgus monkeys of Mauritius origin. Of note, alleles from the MHC II region were not found to be associated with the drug-induced hypersensitivity reactions in cynomolgus monkeys, which is surprising given that cluster of differentiation 4 (CD4)-positive T lymphocytes were predominant in the interface inflammation observed microscopically in the skin from these monkeys (Palanisamy et al. 2015). This might be explained by the less comparable distribution of MHC II haplotypes that were observed between the study population and the general Mauritian monkey population, which may have skewed the statistical analysis.
The mechanism of how small molecule drugs can cause skin hypersensitivity reactions is not well understood, but it is hypothesized that the parent drug or a reactive metabolite binds to a host protein forming a hapten and resulting in an autoimmune response against the drug-modified protein, or alternatively forming a superantigen-like cross-link between MHC and the T-cell receptor resulting in an autoimmune-like response (Uetrecht 2007). Abacavir-induced skin hypersensitivity in HLA B*57:01 carriers is thought to involve the drug binding within the antigen-binding cleft of HLA B*57:01 in an allele-specific manner and altering the HLA-bound peptide repertoire leading to subsequent T-cell responses (Ostrov et al. 2012). It has been reported that there are specific structural features of the antigen-binding cleft that are common or shared with different HLA subtypes in humans (Pompeu et al. 2012). In the present study, it is undetermined which specific allele(s) or residues of the M3 haplotype in the MHC B region are causative for the observed skin reactions in the cynomolgus monkeys, but it is likely that the molecular mechanisms are similar to what is reported in humans, given the high sequence conservation between alleles in the M3 haplotype and human HLA alleles.
In humans, immunological variations in T-cell receptor repertoire and development of tolerance have been reported to contribute to the idiosyncratic nature of drug-induced hypersensitivity reactions (Uetrecht and Naisbitt 2013). Although interanimal variability in T-cell receptor repertoire was not evaluated in the current study, CD4-positive T lymphocytes were predominant at the interface between the dermis and epidermis or within the epidermis itself (Palanisamy et al. 2015). T cells have previously been shown to contribute to skewed peripheral peripheral blood mononuclear cell (PBMC) repertoire in rhesus monkeys (Currier et al. 1999), suggesting it is possible that variability in T-cell receptor repertoire may be one of the factors that contributed to the observed drug-induced hypersensitivity reactions in the cynomolgus monkeys.
In conclusion, this is the first report showing that alleles from the M3 haplotype of the MHC I B region in cynomolgus monkeys are significantly associated with drug-induced hypersensitivity reactions in this species. These findings suggest that cynomolgus monkeys may be a good preclinical model for investigating the mechanisms underlying drug-induced hypersensitivity reactions. Moreover, because of the high degree of sequence similarity between these alleles in the M3 haplotype and HLA alleles associated with drug-induced hypersensitivity reactions in humans, identifying the specific alleles in the M3 haplotype associated with drug-induced hypersensitivity reactions in the monkeys may enable better human translation and will help in the identification of HLA alleles that may impart risk for patients administered these drugs.
Footnotes
Authors’ Note
This contents of this article are solely the responsibility of the authors and do not necessarily represent the official views of NCRR or NIH.
Acknowledgments
The authors would like to acknowledge Karrie Tartaro and Linda Nelms for assistance with DNA isolation.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was made possible in part by Grant Number P51 RR000167 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH), to the Wisconsin National Primate Research Center, University of Wisconsin–Madison. This research was conducted in part at a facility constructed with support from Research Facilities Improvement Program Grant Numbers RR15459-01 and RR020141-01.
