key: cord-0874741-t49zg8vw
authors: Sanchez, Anthony; Wagoner, Kent E.; Rollin, Pierre E.
title: Sequence-Based Human Leukocyte Antigen—B Typing of Patients Infected with Ebola Virus in Uganda in 2000: Identification of Alleles Associated with Fatal and Nonfatal Disease Outcomes
date: 2007-11-15
journal: J Infect Dis
DOI: 10.1086/520588
sha: f1ed6cf07443899783a1d364fbd7b44c467bccf9
doc_id: 874741
cord_uid: t49zg8vw

The Sudan species of Ebola virus (SEBOV) causes severe, often fatal infection in ∼50% of infected humans. We sought to determine whether the human leukocyte antigen—B (HLA-B) locus has a role in the outcome of SEBOV disease by typing 77 cases from an outbreak in northern Uganda in 2000–2001. Sequence-based HLA-B typing was performed using leukocytes isolated from 77 patients. Statistical analysis and a predictive discriminant analysis (PDA) were applied to typing data. Epitope prediction software was also applied to SEBOV sequences. Statistically significant associations were found between certain sets of alleles and either fatal or nonfatal disease outcomes. Alleles B*67 and B*15 were associated with fatal outcomes, whereas B*07 and B*14 were associated with nonfatal outcomes. The PDA-derived functions that were produced were 81.8% accurate in classifying patients into their correct outcome group. Several epitopes predicted to bind strongly to HLA-B*07 molecules were identified in the viral polymerase, nucleoprotein, and VP35 protein. HLA-B alleles associated with either fatal or nonfatal outcomes of SEBOV disease were identified and can be used in a predictive model. Studies of HLA-B—restricted epitopes could contribute to characterization of protective host responses and to vaccine development.

. Surviving these acute infections depends on the ability to efficiently clear virus and infected cells from the body before virus replication becomes unrestricted. Studies of EBOV infections in humans and animal models indicate that an important component of the immune response is cell-mediated immunity [2] [3] [4] . Humoral immunity appears to play only a minor role, if any, in surviving EBOV disease, because (1) antibody responses are characteristically delayed (i.e., they are absent early and usually are undetected in fatal cases) [5] , (2) passive transfer of antiserum has not been effective in protecting primates [6, 7] , and (3) antibody may actually enhance entry of EBOV into host cells [8] . NK cells may contribute to clearance, but there is no evidence suggesting that they play a significant role in controlling EBOV infections [9] . A study of infected macaques indicated a progressive depletion of these cells during the course of the disease [9] .

In early September 2000, a large outbreak of SEBOV disease was identified in the northern region of Uganda and centered around the city of Gulu [10] . A team from the Centers for Disease Control and Prevention (CDC; Atlanta), participating in an international effort to control this outbreak of SEBOV disease, provided on-site diagnostic testing of blood samples and observed the progression of disease in many patients. All persons infected with SEBOV showed signs and symptoms of disease (no asymptomatic cases were identified), yet some patients were observed to have disease of milder severity and shorter duration. Analysis of the acute-phase peripheral blood mononuclear cells (PBMCs) from SEBOV-infected patients showed no significant difference in expression of cytokines, Fas, or Fas ligand when patients with fatal and nonfatal cases were compared with each other or were compared with uninfected, hospitalized patients [3] . However, persons with nonfatal outcomes differed from persons with fatal outcomes in that they had lower viral loads, as well as greater numbers of CD8 + T cells and activated (HLA-DR + ) CD8 + T cells [3, 11] . These findings suggested that the production of cytotoxic T lymphocytes (CTLs) is important in limiting SE-BOV replication and reducing disease severity. Because cellmediated immunity is likely to be a critical component in surviving SEBOV infections, it was postulated that the HLA class I composition of infected patients could influence the course of disease. To investigate this possibility, blood specimens collected from a variety of SEBOV-infected patients were used for sequence-based typing of the HLA-B locus (exons 2 and 3), the most divergent of the major histocompatibility complex (MHC) class I genes. Fortunately, the large number of cases sampled and the ∼50% fatality rate associated with the outbreak allowed statistical analysis to be performed on the data generated. The results of this investigation are the subject of the present article and implicate certain alleles in either nonfatal or fatal outcomes of disease due to SEBOV.

In accordance with CDC institutional review board policies, all samples were anonymized and identified only as fatal or nonfatal cases. Sequence information for patients infected with EBOV was generated from reverse-transcribed mRNA (cDNA) produced in another study [3] . In brief, total RNA was isolated from PBMCs and cDNA produced by priming with oligo-d(T) [12] [13] [14] [15] [16] [17] [18] and random hexamers that were extended with SuperScript II reverse transcriptase and RNA removed with RNase H (Invitrogen; catalog no. 11904-018).

HLA-B locus sequence-based typing. Human HLA-B sequences were polymerase chain reaction (PCR) amplified from cDNA by use of Platinum PCR SuperMix (Invitrogen) or the Expand High Fidelity System (Roche Applied Science). To decrease background amplifications, 5% (v/v) dimethyl sulfoxide (D-2650; Sigma-Aldrich) was included in all reactions. Amplification of exons 2 through 3 was accomplished using the primers F1 (5 -CCGAACCGTCCTCCTGCTGCTC-3 ; end of exon 1) and R1 (5 -GTGGTGGGTCACGTGTGTCTTTG-3 ; beginning of exon 4). Separate amplification of exon 2 was accomplished using primers F1 and R2 (5 -CCCGCGGAGGA-GGCGCCCGT-3 ; beginning of exon 3), and amplification of exon 3 with primers F2 (5 -GCTACTACAACCAGAGCGAG-GCC-3 ; end of exon 2) and R1 (amplicon lengths of 640, 329, and 386 bp, respectively) (figure 1). Thermocycling conditions for amplification were as follows: 30-35 rounds of denaturing at 94ЊC for 15 s, annealing at 60ЊC for 15 s, and extension at 72ЊC for 40 s. PCR amplification generally required only the use of F1 and R1 primers to produce enough template DNA for direct sequencing; however, at times, the use of F1 and R2 and/or F2 and R1 primer combinations was needed to increase amplicon synthesis. Amplicons were separated from reaction components by electrophoresis in 1.3% agarose and ethidium bromide-stained bands isolated using Qiagen Qiex II (with DNA eluted from binding matrix using nuclease-free water). Cycle sequencing of amplicons was performed using Prism BigDye chemistry (version 1; Applied Biosystems) in 20-mL reactions containing 100 ng of primer (4 separate reactions were performed using primers F1, R1, F2, and R2). Reactions were prepared in 0.2-mL thin-walled PCR tubes and were heated to 96ЊC for 1 min, followed by 25 cycles of denaturing at 96ЊC for 10 s, annealing at 50ЊC for 5 s, and extension at 60ЊC for 4 min. Extension products were isolated by passing reactions through Centri-Sep columns (Princeton Separations) and thoroughly drying them in a speed vacuum. Products were analyzed on an automated sequencer (ABI 377 or ABI Prism 3100; Applied Biosystems) according to the manufacturer's protocols. Sequencher software (version 3.1.1; Gene Codes) was used to analyze chromatograph files in the generation of consensus HLA-B sequences and their exportation as ASCII files. HLA-B allele assignments were generated using the Assign software program (version 1.5.1) for sequence-based typing (developed by D. Sayer and J. Meier, Department of Clinical Immunology, Royal Perth Hospital, Perth, Australia). A description of Assign, version 2.0, has been published elsewhere [12] , and the most recent version is available from Conexio Genomics. HLA-B assignment was limited to the allele type (i.e., B*07 and not B*0713), which provided sample sizes large enough to perform statistical analyses.

Epitope prediction. The encoded amino acid sequences for all genes of Sudan EBOV Gulu [13] were analyzed using Internet-based software/tools to identify epitopes potentially involved in cell-mediated immune responses. Analyses used CTLPred [14] , a direct method for predicting CTL epitopes, as well as the following indirect methods, which were based on MHC antigen-binding properties: BIMAS [15] , Rankpep [16] , EpiJen [17] , NetCTL [18] , and ARB Matrix [19] . It is of note that some of these methods also factor in proteasomal C-terminal cleavage and TAP binding.

Statistical data analysis. Predictive discriminant analysis (PDA) techniques [20] were used to predict survival among Fisher's exact test was used as the initial screening device to determine which of the 26 HLA-B alleles were associated with patient survival. An allele was considered for inclusion in the PDA model if its univariate association with patient survival resulted in a P value of р.20. Once an initial set of predictor variables was selected from the univariate analyses, an all-possible subsets approach was used to select a set of HLA-B alleles that maximized classification accuracy.

Accuracy of the PDA model was assessed by examining the overall and individual group hit rates. The hit rates were estimated using the leave-one-out (L-O-O) method, a 1-sample cross-validation procedure that produces classification estimates that are "nearly" unbiased estimates of the true population hit rates. The overall hit-rate estimate was statistically compared with chance expectation by use of proportional chance criterion and maximum chance criterion methods. The group hit-rate estimates were statistically evaluated using the proportional chance criterion. If the overall and group hit rates were deemed to be statistically better than chance expectation, a descriptive measure (I) was calculated to denote percent improvement over chance expectation.

In this data set, a higher percentage of patients died (54.5%) than lived (45.5%). Because HLA-B allele information was not available for all patients, the data set does not include all patients who were infected with EHF in the outbreak in Uganda. Because the exact probability of surviving EHF is not exactly known, the prior probability of group membership was set at , which is a value close to the sample estimates in the P p .5 current data set.

All PDA analyses were performed using the Discriminant procedure in SSPS software (version 12.0; SPSS), and all PDA analyses/interpretations were done according to established guidelines [20] . Because the HLA-B predictor variables were dichotomous variables, and because the discriminant analysis procedure assumes multivariate normality, the classification analyses were also done using Proc Logistic in SAS software (version 9.1; SAS Institute), as well as various kernel discriminant analysis options, by use of Proc Discrim in SAS software (version 9.1). Identical or near identical results were obtained using all statistical procedures.

Results of HLA-B allele assignment for the 77 patients typed are presented in table 1 and show that all 11 patients with the B*07 allele survived SEBOV disease. Fortunately, this allele has been included in most MHC binding/CTL prediction programs, which allowed the analysis of the predicted amino acid sequences among all the genes of Sudan EBOV Gulu for potential epitopes. Table 2 lists potential HLA-B*07 epitopes identified by prediction software, and, of the top 3 ranked epitopes, 2 are present in the L protein and 1 is present in the nucleoprotein. Lower-ranking epitopes are located in the L protein and VP35 protein. The highest-ranked epitope (1207-KPRCPSAAL) is highly conserved in all EBOVs and is only partially conserved in Marburg virus (MARV), a related filovirus. The nucleoprotein epitope 313-YPQLSAIAL, which is the second-ranked epitope, is highly conserved for both EBOV and MARV. The third-ranked epitope (1879-MPKTHRGEL) is not conserved. Of the other sequences, only the fourth-and fifthranked epitopes (1871-TPRMLLPIM and 292-IPKACQKSL, respectively) show any conservation (highly conserved for EBOV only).

Univariate associations between patient survival and the By chance alone, it was expected that 38.5 (50%) of the patients would be classified into their correct outcome group. On the basis of results from the L-O-O method, the overall classification accuracy rate was 81.8% (63 of 77 patients), which is a marked improvement over chance expectation (I p ). By use of the proportional chance criterion (table 5) , 63.6% the overall classification accuracy of 81.8% was statistically better than chance, with and . On the basis of z p 20.11 P ! .0001 the maximum chance criterion (i.e., classification of all subjects into the largest group), 42 (54.5%) of 77 patients would be correctly classified by chance alone. On the basis of the maximum chance criterion, the overall hit rate (81.8%) was better than chance, with , , and . By use z p 16.63 P ! .0001 I p 60.0% of the proportional chance criterion, the group hit rates were also statistically better than chance, with , , z p 7.76 P ! .0001 and for the group of patients who died of EHF and I p 52. . Insertion of a "0" value for an allele to denote genetic absence and a "1" value to denote genetic presence of an allele allows the calculation of LCFDead and LCFAlive values. The patient is assigned to the group with the higher LCF value.

After the SEBOV outbreak in Uganda in 2000-2001, it was determined that persons with fatal outcomes of SEBOV disease had higher virus loads than did those who survived the disease [3, 11] . Survival is likely dependent on the rapid development of an effective CTL response that limits the level of virus replication before it can overwhelm the body, as is evidenced by increased numbers of activated (HLA-DR + ) CD8 + lymphocytes in patients with nonfatal cases [3] . As our study suggests, the MHC class I composition appears to be an important factor in disease outcome, and the presence of certain HLA-B alleles can confer either a protective advantage or a disadvantage in patients infected with SEBOV. Our findings indicate that ex- pression of the B*07 allele, which is highly associated with survival, confers a phenotype that efficiently drives T cell activation and CTL expansion early in the SEBOV infection, whereas alleles that are highly associated with fatalities (e.g., B*67 and B*15) contribute to an inefficient or unresponsive phenotype. One of the SEBOV-infected patients who survived had a B*07/B*15 profile (involving good and bad alleles), which suggests that, by itself, the B*07 allele may be sufficient to restrict SEBOV replication. With more data, one might be able to model the survival of patients with one good allele and one bad allele. However, in the present study, only 2 patients had such a profile, and, therefore, only univariate associations between allele type and patient survival were possible. Correlations between HLA composition and disease susceptibility or severity have been identified for other viral infections, including those due to HIV [21] [22] [23] [24] [25] , herpes viruses [26, 27] , hepatitis B and C viruses [28] [29] [30] , dengue virus [31] , hantaviruses [32, 33] , and severe acute respiratory syndrome coronavirus [34, 35] . It has also been reported that HLA-B-restricted epitopes exhibit a greater functional avidity (T cell receptor avidity for HLA-peptide complex) than do HLA-Aor HLA-C-restricted epitopes [36] , which suggests a significant role for HLA-B alleles in acquired antiviral immunity. We believe that this is also the case for SEBOV infections.

Although we have shown that the HLA-B profile of a patient with SEBOV infection is associated with disease outcome, additional investigation is needed to further characterize this find-ing. The application of MHC class I tetramer technology [37] could be used to analyze PBMCs isolated from persons who have survived SEBOV disease or future cases of SEBOV disease, to determine whether any of the epitopes predicted in table 2 are immunodominant in HLA-B*07-restricted responses. There is also a need to examine the HLA-B profiles of patients infected with other filoviruses, which would be extremely helpful in gauging the influence of the HLA-B locus on filovirus clearance and could be useful in the development and refinement of vaccines.

The ability to predict the capacity of an individual to control acute and life-threatening infections, such as infections caused by filoviruses, could be useful in outbreaks. Current and future capabilities for rapid HLA typing could easily allow this type of testing to be performed in the field alongside other clinical assays. The classification functions (LCFdead and LCFAlive) generated from the PDA model could be used to gauge the risks to those providing medical assistance to patients with SEBOV infection (should they themselves become infected), and high-risk individuals might be assigned duties involving less contact with infected patients and infectious materials. Similarly, patients with SEBOV infection might be managed more appropriately if a fatal outcome is likely. Future typing of patients with SEBOV infection would be useful in assessing the reliability of such tools, and, if validated, these tools could help refine the classification functions and expand their scope.

In conclusion, the present study provides, to our knowledge, the first evidence of genetic predispositions to fatal and nonfatal disease outcomes in patients with EBOV infection, and it indicates that these traits are linked to alleles of the MHC HLA-B locus. This association emphasizes the importance of the development of timely CTL responses to EBOV infections that lead to efficient virus clearance. It is hoped that this data set will be expanded to include predicted disease outcomes for ZEBOV and MARV, which would be an important resource for future investigations of human immune responses to filovirus infections.

Marburg and Ebola viruses

Defective humoral responses and extensive intravascular apoptosis are associated with fatal outcome in Ebola virus-infected patients

Analysis of human peripheral blood samples from fatal and nonfatal cases of Ebola (Sudan) hemorrhagic fever: cellular responses, virus load, and nitric oxide levels

Ebola virus pathogenesis: implications for vaccines and therapies

Clinical virology of Ebola hemorrhagic fever (EHF): virus, virus antigen, and IgG and IgM antibody findings among EHF patients in Kikwit, Democratic Republic of the Congo

Passive immunization of Ebola virus-infected cynomolgus monkeys with immunoglobulin from hyperimmune horses

Evaluation of immune globulin and recombinant interferon-a2b for treatment of experimental Ebola virus infections

Antibody-dependent enhancement of Ebola virus infection

Depletion of peripheral blood T lymphocytes and NK cells during the course of Ebola hemorrhagic fever in cynomolgus macaques

Outbreak of Ebola viral hemorrhagic fever Uganda

Rapid diagnosis of Ebola hemorrhagic fever by reverse transcription-PCR in an outbreak setting and assessment of patient viral load as a predictor of outcome

Assign 2.0: software for the analysis of Phred quality values for quality control of HLA sequencing-based typing

Complete genome sequence of an Ebola virus (Sudan species) responsible for a 2000 outbreak of human disease in Uganda

Prediction of CTL epitopes using QM, SVM and ANN techniques

Scheme for ranking potential HLA-A2 binding peptides based on independent binding of individual peptide side-chains

Enhancement to the RANKPEP resource for the prediction of peptide binding to MHC molecules using profiles

EpiJen: a server for multistep T cell epitope prediction

An integrative approach to CTL epitope prediction: a combined algorithm integrating MHC I binding, TAP transport efficiency, and proteasomal cleavage predictions

Automated generation and evaluation of specific MHC binding predictive tools: ARB matrix applications

Applied discriminant analysis

Associations between MHC class I and susceptibility to HIV-2 disease progression

Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA

Characterization of HLA-B57-restricted human immunodeficiency virus type 1 Gagand RT-specific cytotoxic T lymphocyte responses

HLA B*5701 is highly associated with restriction of virus replication in a subgroup of HIV-infected long term nonprogressors

Association of HLA types A1-B8-DR3 and B27 with rapid and slow progression of HIV disease

Heirarchy of Epstein-Barr virus-specific cytotoxic T-cell responses in individuals carrying different subtypes of an HLA allele: implications for epitope-based antiviral vaccines

Relative dominance of HLA-B*07 restricted CD8 + T-lymphocyte immune responses to human cytomegalovirus pp65 in persons sharing HLA-A*02 and HLA-B*07 alleles

Cytotoxic T lymphocyte response and viral load in hepatitis C virus infection

Human leukocyte antigen phenotypes and hepatitis C viral load

HLA phenotypes and outcomes of hepatitis B virus infection in Taiwan

HLA-A and -B allele associations with secondary dengue virus infections correlate with disease severity and the infecting viral serotype in ethnic Thais

Role of specific CD8 + T cells in the severity of a fulminant zoonotic viral hemorrhagic fever, hantavirus pulmonary syndrome

Association of HLA B27 with benign clinical course of nephropathia epidemica caused by Puumala hantavirus

Association of HLA class I with severe acute respiratory syndrome coronavirus infection

Association of human-leukocyte-antigen class I (B*0703) and class II (DRB1*0301) genotypes with susceptibility and resistance to the development of severe acute respiratory syndrome

Impact of HLA-B alleles, epitope binding affinity, functional avidity, and viral coinfection on the immunodominance of virus-specific CTL responses

Major histocompatibility complex (MHC) tetramer technology: an evaluation

We thank the Ugandan Ministry of Health and the World Health Organization Response Team for their support in facilitating the collection and processing of patient samples. We acknowledge the clinical leadership and counsel of the late Matthew Lukwiya; the encouragement of the late Piero Corti; the isolation ward support of Daniel Bausch, Maria Di Santo, and Yoti Zebulon; the laboratory and logistical support of Franca Cian and Elio Croce; and the many efforts and sacrifices of the staffs of the St. Mary's Lacor and the Gulu Hospitals. Supplement sponsorship. This article was published as part of a supplement entitled "Filoviruses: Recent Advances and Future Challenges," sponsored by the Public Health Agency of Canada, the National Institutes of Health, the Canadian Institutes of Health Research, Cangene, CUH2A, Smith Carter, Hemisphere Engineering, Crucell, and the International Centre for Infectious Diseases.