Ultra-high resolution HLA genotyping and allele discovery by highly multiplexed cDNA amplicon pyrosequencing
© Lank et al.; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 7 July 2012
Published: 6 August 2012
High-resolution HLA genotyping is a critical diagnostic and research assay. Current methods rarely achieve unambiguous high-resolution typing without making population-specific frequency inferences due to a lack of locus coverage and difficulty in exon-phase matching. Achieving high-resolution typing is also becoming more challenging with traditional methods as the database of known HLA alleles increases.
We designed a cDNA amplicon-based pyrosequencing method to capture 94% of the HLA class I open-reading-frame with only two amplicons per sample, and an analogous method for class II HLA genes, with a primary focus on sequencing the DRB loci. We present a novel Galaxy server-based analysis workflow for determining genotype. During assay validation, we performed two GS Junior sequencing runs to determine the accuracy of the HLA class I amplicons and DRB amplicon at different levels of multiplexing. When 116 amplicons were multiplexed, we unambiguously resolved 99%of class I alleles to four- or six-digit resolution, as well as 100% unambiguous DRB calls. The second experiment, with 271 multiplexed amplicons, missed some alleles, but generated high-resolution, concordant typing for 93% of class I alleles, and 96% for DRB1 alleles. In a third, preliminary experiment we attempted to sequence novel amplicons for other class II loci with mixed success.
The presented assay is higher-throughput and higher-resolution than existing HLA genotyping methods, and suitable for allele discovery or large cohort sampling. The validated class I and DRB primers successfully generated unambiguously high-resolution genotypes, while further work is needed to validate additional class II genotyping amplicons.
KeywordsHLA Genotyping Roche/454 Pyrosequencing Galaxy Tissue typing Cellular immunity Multiplexing
The human leukocyte antigen (HLA) complex on the short arm of chromosome 6 encodes the most polymorphic genes in the human genome, and is extensively studied due to the central role its associated proteins play in the etiology of infectious diseases, autoimmune disorders, and organ transplantation outcomes [1–5]. The “classical” HLA class I and II gene products are involved in endogenous and exogenous antigen presentation respectively, mediating adaptive immune responses to pathogens and tissue transplants. As sequencing methods of characterizing human genetic diversity have proliferated, the number of known classical HLA alleles has exploded, with over 6,800 distinct variants described to date, 1,633 of which were identified in 2010 alone . The majority of observed allele variants occur within the class I genes HLA-A, -B and -C, as well as the class II DRB genes. These highly polymorphic loci are the most commonly studied for their differential effects on directing adaptive immunity [7, 8]. The rapid proliferation of reference alleles in the curated IMGT database demonstrates that coding HLA polymorphism is vast, with certain serologically distinct families containing many hundreds of individual alleles. For instance, in May of 2011, four new HLA-A*02 alleles were described  with a total of 334 protein-distinct alleles of this lineage described to date.
Non-synonymous allelic variants of HLA gene products can bind distinct antigenic peptides [10, 11], be subject to differential regulation [12, 13], and have varied interactions with T-Cell Receptors , Killer Immunoglobulin-like Receptors and viral proteins . The direct link between HLA polymorphism and various diseases is a subject of intensive investigation, with hundreds of published studies every year. High-resolution HLA typing (to the specific protein, rather than the serological group level) can elucidate the differential consequences of small protein changes within an allele family, such as the alleles B*57:02 and B*57:03, which can differ by a single nucleotide yet lead to differential outcomes during HIV infection . The ability to type individuals at true high-resolution may uncover further links between individual alleles and disease.
Traditional HLA typing methods are losing their effectiveness at accurately predicting high-resolution (four-digit) genotype due to the large number of highly similar alleles. Molecular sequence specific PCR (SSP) and sequence specific oligonucleotide hybridization (SSO) methods use unique sequence signatures to detect the presence or absence of an allele; however, with more alleles being described every month, many with only one novel single nucleotide polymorphism (SNP), it is impossible for SSP panels to rule out the presence of uncommon but similar alleles. For higher accuracy, these methods rely on an ever-greater number of reactions to distinguish potentially related alleles, or use complementary methods to resolve ambiguities . More sophisticated (and expensive) sequence-based methods improve resolution; however, even “gold-standard” commercial kits fail to produce the entire coding sequence of individually typed alleles. Allele ambiguity is an increasing problem as many new alleles have been described with SNP variants outside of the traditional exons examined by sequence-based typing (SBT). This has led to the proliferation of methods for imputing four-digit genotype from limited sequencing data by utilizing known allele frequencies in various populations [19–21].
Four-digit ambiguous alleles with standard SBT methods and cSBT
Class I exon 2, 3, 4 ambiguous to 4-digit (high freq)
DRB Exon 2 Ambiguous to 4-digit (high freq)
All cSBT ambiguous to 4-digit
Recently, we described a novel HLA class I genotyping method (cSBT), utilizing a PCR amplicon that universally amplifies a 581 bp diagnostic template from cDNA, encompassing all of exon 3 plus portions of exons 2 and 4 . This method allows unprecedented throughput for class I HLA typing (through the use of a single universal PCR amplification) but lacks the ability to fully distinguish class I exons 2 and 3 (the minimal requirements for novel allele discovery), producing statically ambiguous genotyping results.
The primary aim of this study was to expanded the area of sequence coverage using novel primer binding sites in exons 1 and 7 of the HLA class I open reading frame, pan-locus DRB primers, and preliminary use of universal primers for the other class II loci (DRA, DPA, DPB, DQA, and DQB) to unambiguously identify HLA genotypes in large patient cohorts with unprecedented resolution and throughput without relying on the use of population-specific allele frequency inferences [20, 23]. We also designed an open-source analysis workflow for rapidly imputing high-resolution genotype. In the process, several novel HLA alleles were identified and validated with direct PCR Sanger sequencing. Thus, the improved cSBT method allows rapid screening of large cohorts for potential novel allele characterization. This general approach will also be adaptable to improving second- and third-generation sequencing technologies that may allow longer amplicons to be completely sequenced [25, 26].
HLA Class I primer design
HLA Class II primer design
Library preparation and pyrosequencing
We generated cDNA from total cellular RNA, and PCR-amplified the three primary diagnostic HLA amplicons described above (HLA-1, HLA-2 and DRB) from 98 separate samples, including 40 International Histocompatibility Working Group reference cell lines http://ihwg.org. For a limited set of four reference samples we also generated amplicons from the additional class II loci. The quality and purity of the DNA amplicon sequencing library is important to ensure successful sequencing , so we performed multiple solid-phase reversible immobilization (SPRI) clean-ups on our samples to remove excess primer-dimer and shorter PCR products. We used 96 distinct multiplex identifier (MID)-tagged primer pairs for each amplicon to generate sequencing libraries Additional file 1: Table S3.
The presence of differently sized amplicons in the sequencing library is an important consideration. The final amplicon sizes of the primary products (including MIDs and adapter sequences) were 746 bp (HLA-1), 995 bp (HLA-2), and 737 bp (DRB Amplicon). Although HLA-2 is longer than typically recommended for Roche/454 pyrosequencing, the use of an augmented pooling strategy accounting for the differential sizes of amplicons (longer products at a higher molar ratio than the shorter ones) produced normalized sequence data. This was necessary since the Roche/454 emPCR process preferentially amplifies shorter PCR products.
Read metrics for both sequencing runs described
Ave Length (st. dev)
Reads per Sample (st. dev)
Amp 1 Reads per Sample (st.dev)
Amp 2 Reads per Sample (st. dev)
DRB Reads (st. dev)
Analysis design and implementation overview
The different length of HLA reference alleles relative to the sequence positions necessitates an iterative approach to determining genotypes. To improve data regularity and quality, all reads were subjected to stringent length filtering and trimming (between 300 and 350 bp trimmed reads, depending on the amplicon), which resulted in the loss of 25% to 60% of total filterpass sequences. The remaining normalized sequence reads were mapped at 100% identity to reference alleles, and the pattern and number of sequence matches for each individual was evaluated to deduce genotype.
For samples that sequenced or amplified at lower numbers, and for our second highly multiplexed experiment, we employed a less stringent, modified pipeline that allowed allele-specific read groups of any number. These calls are noted in orange on the genotyping tables and are of intermediate confidence due to the limited number of sequencing reads.
Class I genotyping analysis
For class I genotyping, allele mapping began with the HLA-2 forward reads (2 F) to determine the constellation of potential alleles present in an individual Figure 3c. This was possible since trimmed 2 F reads encompass only part of exons 2 and 3, and all named alleles have complete exon 2 and 3 sequences. Perfectly matched 2 F reads thus dictated the reference set for subsequent matches, reducing the number of alleles evaluated for subsequent mapping from several thousand to several hundred and dramatically decreasing computation time for subsequent BLAT mapping.
Using this restricted subset of reference alleles, the other trimmed read groups were evaluated for identity. Based on the positioning of these reads relative to the reference, we established match conditions for each potential allele. These conditions included fully matched reads (condition #1), partially matching reads where the query read contained nucleotide positions not known for the reference, such as exon 1 sequence for exon 2/3-only alleles (condition #0), and unmatched reads (condition null). For alleles of different lengths, a distinct matrix of read match conditions is expected Figure 3d. The matrix is generated by performing inner joins of the 2 F dataset to the 1 F and 1R matches, and an outer join with the 2R matches to account for short reference alleles with null 2R match conditions.
Class II genotyping analysis
For class II data, the same length and quality filtering were applied, and similar issues arose when mapping the forward and reverse reads, as many alleles are known only from exon 2 sequence. In this case forward reads are used as the key to limit the potentially matching allele set, considering matches with extension sequences at both the 5’ and 3’ end valid due to the positioning of the reads relative to exon 2. Reverse reads are then evaluated against a reference subset of only forward matching alleles. A condition matrix, similar to the one used for class I reads, is established. However, due to the prevalence of exon-2-only reference alleles, an additional condition is added for forward reads when the query reads contain additional 5’ and 3’ sequence relative to the reference.
Class I genotyping results
Our attempt to push the assay to 92 samples (271 amplicons, both class I and DRB) in a single GS Junior run for the second experiment missed many alleles (not seen, rather than incorrectly called) without using the modified, lower-confidence pipeline. Still, at this multiplexing level our assay yielded high-resolution, accurate genotyping for 93% of expected class I alleles (Additional file 1: Table S6).
Our sequence-based results overturned several incorrect typings in the reference samples that were subsequently confirmed with traditional typing methods. With the additional exon coverage in this experiment, we were able to demonstrate deficiencies in the original SSOP reference cell line typing. This included correcting the original HLA-C*07:01 calls for samples HLA-Ref22 andHLA-Ref28 to HLA-C*07:18. The distinguishing SNP between these alleles occurs at nucleotide 1043, and was not assayed with either previous SSP or SBT methods. We also identified a putative novel HLA-C*07 allele in sample HLA-Ref31.
Class II genotyping results
In the third experiment of additional class II loci amplicons we identified alleles from each locus sequenced for all four samples Additional file 1: Table S8. Comparison of this data with the reference class II typing previous done on the samples by SSOP was 72% concordant (13 of 18 alleles matching). However, in several cases our results were discordant with original typing, either identifying different variants of the same serological group, or identifying alleles not seen with SSOP typing. It is unclear whether these differences result from our analysis methods or further deficiencies in the SSOP reference typing, as was the case for several class I alleles in the other experiments.
Novel allele characterization
As a demonstration of the ability of the assay to detect novel alleles, we chose one of the putative novel sequences (Additional file 3: File 2) identified in sample KTS02 and confirmed its sequence using allele-specific PCR and Sanger sequencing. All alleles present in KTS02 were interrogated, and primers specific for the putative novel allele were identified. Following primary PCR, six allele-specific primers were used to Sanger sequence this product (Additional file 1: Table S9). Analysis of the resultant data confirmed the novel SNP (T587C), resulting in a proline to leucine substitution. This sequence was submitted for approval and the name HLA-C*08:56 was officially assigned by the WHO Nomenclature Committee in October 2011. This follows the agreed policy that, subject to the conditions stated in the most recent Nomenclature Report that names will be assigned to new sequences as they are identified . Lists of such new names will be published in the following WHO Nomenclature Report.
We present an improved version of a previously described cSBT assay with dramatic improvements in resolution through the use of additional universal primer binding sites in the HLA class I open reading frame, the addition of a universal DRB amplicon, and preliminary data from additional class II loci amplicons. Our results demonstrate a scalable technique for use on either the Roche/454 GS Junior or the FLX platform with Titanium chemistry.
The use of a cDNA template, combined with clonal second-generation pyrosequencing, allows comprehensive genotyping from only two PCR amplifications for all HLA class I loci, instead of the nine or more exon- and locus-specific amplifications required for genomic DNA template methods (both Sanger and Roche/454-based). Other amplicons, such as those presented for DRB and other class II loci, can be used to achieve high-resolution genotyping as needed. This updated method allows for dozens or hundreds of patient samples to be HLA class I (and DRB or other loci) genotyped to true high-resolution (allele-specific) simultaneously using fewer amplicons per sample, at a low per-sample cost for both large scale screening projects (1,000 s of samples) using a GS-FLX instrument, or smaller (dozens to hundreds of samples) on a GS Junior instrument. Per-sample cost varies depending on the level of multiplexing and other factors, but can be estimated at ~ $70 - $100 for a 48 to 96 sample multiplexed experiment on a GS Junior instrument (Additional file 1: Table S10). Our results also demonstrate that cSBT has higher throughput and higher resolution than existing methods, and does not rely on population-specific allele frequency inferences.
Higher throughput than reported here (>92 samples typed at a time) is possible for the entire pre-emPCR process in a lab setting utilizing liquid handling robotics for RNA extraction, cDNA synthesis, PCR amplification, SPRI purification, quantification, and pooling, similar to the methods employed by Erlich et al. . Despite the need to generate 288 amplicons per 96 samples typed with the cSBT method for class I and DRB loci, this is a dramatic reduction over the >960 amplicons needed to perform comparably accurate typing from genomic DNA starting material using locus- and exon-specific amplification primers [22, 23, 28].
An important caveat to the presented data is the different levels of multiplexing used in the first two genotyping experiments. This limited data is insufficient to establish an “ideal” multiplexing level for cSBT. The first GS Junior sequencing run, which multiplexed 116 amplicons (equivalent to 58 class I samples, 39 class I + DRB, or 15 class I + all class II loci) was highly accurate, and is the current recommended multiplexing level. Investigators who consistently recover abundant data may wish to increase the multiplexing level nearer to the second sequencing run (271 amplicons multiplexed), which, in the presented data, performed nearly as well as the lower multiplexed run when lowering sequence number requirements.
Another consideration is the use of and availability of high-quality RNA, which is necessary for cSBT. Traditionally, HLA typing facilities have collected DNA for use with established assays. DNA has the advantage of being easier to collect and store, whereas RNA requires more robust preservation and extraction methods. The limited access to RNA may not make cSBT ideal for all situations, particularly clinical settings that only have access to poor-quality starting material. In such cases a genomic DNA-based method may be more appropriate . Current studies, underway with collaborators, are evaluating the use of saliva, buccal swabs, or tissue biopsies as RNA starting material. Initial results with saliva starting material have led to successful amplicon generation. However, class II genotyping likely requires the use of whole blood or PBMC as a starting material, since class II transcripts are unlikely to be expressed in other tissues.
This general approach will be adaptable to improving second- and third- generation sequencing techniques that may allow even longer amplicons to be completely sequenced [25, 26]. As longer read-length clonal sequencing technologies become available (such as GS-FLX + from Roche/454 or SMRT sequencing from Pacific Biosciences), it should be possible to use the exon 1 and exon 7 primer sites described here to amplify a single PCR product, instead of the two tiled amplicons that are required currently. We have already successfully amplified and sequenced this ~1 kb amplicon with Roche/454 Titanium sequencing, although the high-quality read lengths were insufficient to provide complete coverage (data not shown). However, the use of a single amplicon increases the likelihood that a subset of alleles with primer incompatibilities could be missed, since a two-tiled amplicon approach produces redundant data with two separate amplifications. Using two class I amplicons acts as an important control to ensure that any allele with primer incompatibilities for a given amplicon is detected with the other, such as the detection of HLA-B*18:01:01:01 in sample HLA-Ref20 despite a primer mismatch for one of the amplicons. In our highly-multiplexed experiment, the loss of one amplicon read direction was more common, even for non-mismatched alleles, and this impacted our ability to distinguish allele resolution, though we were still able to detect signatures of alleles missing from a single read group. The use of a second class I amplicon also facilitates analysis of apparently homozygous loci, which can be artifacts of poorly amplified or potentially novel alleles. It should be noted that although the sequencing platform and read length may change, the cSBT primers and much of the analysis methods presented here would be applicable to any cDNA-based assay.
Methods with more limited coverage than cSBT can lead to incorrect genotype associations. For instance, in our reference panel we were able to determine that two allele calls originally labeled as HLA-C*07:01 were in fact HLA-C*07:18, which encodes a distinct protein. This underscores the power of our approach, which does not assume ambiguously amplified sequences to be the most likely allele in a given population. Additionally, our genotyping results were usually accurate to six-digit resolution, an unprecedented achievement for a high-throughput HLA genotyping assay. While it is currently unknown whether synonymous coding six-digit variants have differential effects on allele expression or disease correlations, it should now be possible to examine this question.
Our preliminary investigations into typing other class II loci have show that the presented primers generate sequencable products, which map to known alleles. However, more work is needed to validate these amplicons properly as we have not determined the cause of observed differences between reference SSOP typing and the presented cSBT typing. We provide this data, and the primers used, so that investigators can carry out further work on these loci if desired.
Two other potential applications of expanded cSBT are novel allele discovery and extending the reference sequence length of known alleles. In the first case, we have shown that cSBT can effectively screen large cohorts for novel alleles, which can then be characterized with traditional and accepted sequencing methods. In the second application, the majority of named HLA class I alleles are known only from exon 2 and 3 sequence (3,439 alleles). Thus, if cSBT was applied to samples with these rare alleles, exons 4, 5, and 6 and parts of exons 1 and 7 could be quickly deduced to improve the length and quality of the IMGT reference database.
The rapidly expanding number of described classical HLA alleles is undermining the utility of previously unambiguously high-resolution typing assays to provide unambiguous results. The cSBT method presented here offers one solution to overcome the difficulties associated with the ever-expanding HLA allele library: sequence a larger fraction of the open reading frames, and use cDNA starting material to eliminate the problems of exon-phase matching. Combining this with clonal, long-read-length Roche/454 pyrosequencing has produced the ability to unambiguously resolve alleles for large cohorts at a fraction of the cost and time of previous methods. This method also does not rely on population-specific allele frequency inferences, and only uses the a priori knowledge of the conserved primer binding sites to amplify both known and novel alleles. The modified cSBT method for high-throughput, high-resolution HLA typing should accelerate studies involving large cohort genotyping, and allow a deeper survey of worldwide HLA polymorphism in historically understudied regions .
Human cell line or patient-derived cellular RNA was collected from a variety of sources and in accordance with the University of Wisconsin and Dana Farber Cancer Research Center’s policy on human research subjects. Samples were anonymized prior to typing, and no genotyping results were released to physicians or used to inform clinical treatments. Samples were processed for RNA either as cell pellets, cell-culture homogenates, peripheral blood mononuclear cells, or whole blood. Samples included 15 patient samples from an unrelated HIV study at the University of Wisconsin, 40 HLA reference cell lines provided by the International Histocompatibility Working Group, 5 human embryonic cell-lines provided by WiCell, 40 samples derived from cell-lines or patients collected by the Dana Farber Cancer Center, and 10 blood samples from female Korean patients with HPV-16-induced cervical cancers from an unrelated study .
RNA isolation and cDNA synthesis
Methods for isolation of RNA and cDNA synthesis were identical to those presented previously , involving the MagaNA Pure LC RNA isolation platform (Roche Applied Science, Indianapolis, IN), quantification using a NanoDrop 1000 spectrometer (ThermoFisher, Waltham, MA) and oligo(dT)20-primed SuperScript III cDNA synthesis (Invitrogen, Carlsbad, CA). For all experiments an input of 50 ng RNA was used to seed the cDNA synthesis.
Universal primer design
The design of the 581 bp internally diagnostic amplicon has been previously described . The method for finding additional conserved sites for HLA class I and HLA DRB universal primers was similar, and is briefly described here. All known HLA class I and DRB reference sequences from the IMGT/HLA database were aligned using MUSCLE  or CodonCode Aligner (CodonCode Corporation, Dedham, MA). All analysis was done with IMHT/HLA Release 3.4 (April 2011) . In addition to the previously described HLA class I 581 bp amplicon spanning consensus positions 145–725  we identified novel conserved primer binding sites in exons 1 and 7. Due to current read length limitations of GS-FLX Titanium reagents we could not fully sequence this amplicon. However, by combining this 1020 bp amplicon with the 581 bp amplicon acting as a bridge, we were able to tile nucleotides 50–1069 of the HLA class I open reading frame, including all of exons 2 through 6. To further normalize the distribution of reads, we altered these original 581 and 1 kb amplicons by mixing their forward and reverse primers. The class II DRB and other loci genotyping amplicon was identified with similar methods.
Initial sequencing with the HLA-1 forward primer revealed a sequence incompatibility with most HLA-A alleles, leading to the generation of a short, chimeric DNA bead which interfered with genotyping Additional file 4: Figure S1. To correct this issue, we switched the Roche/454 adapter positions, using adapter “B” with the forward gene-specific primer, and adapter “A” with the gene-specific reverse primer.
Amplicon and library preparation
Fusion primers with Roche/454 adapters, MIDs, and gene-specific sequences were generated [Additional file 1: Table S3] and used to amplify PCR products using Phusion High-Fidelity DNA polymerase (New England Biolabs, Ipswich, MA) as previously described . PCR products were purified with SPRI Ampure-XP paramagnetic beads (Beckman Coulter Genomics, Danvers, MA). DNA concentrations were normalized, and products were pooled into a single amplicon library. Due to the different sizes of our pooled amplicons, and the expected number of alleles sequenced, the three amplicons (HLA-1, HLA-2 and DRB amplicon) were pooled at a 3:6:1 ratio.
Emulsion PCR and pyrosequencing were performed with the amplicon (Lib-A) kit with GS Junior reagents from Roche/454 using standard, manufacturer provided protocols (454 Life Sciences, Branford, CT). An input DNA molecule-to-bead ratio of 1.5 was used. DNA bead enrichment levels were expected (between 5 and 20%) for all experiments, and sequencing yielded >70,000 filterpass reads per GS Junior run.
Novel allele Sanger sequencing
One novel allele identified in sample KTS02 was confirmed through the use of Sanger sequencing. Primers specific to the novel allele (and not to the other alleles present in this individual) were used to amplify a primary PCR product from previously generated cDNA using High Fidelity Platinum Taq (Invitrogen, Carlsbad, CA). Products were purified with SPRI Ampure-XP paramagnetic beads, and sequenced on an Applied Biosystems 3730xl capillary sequencer (Applied Biosystems, Carlsbad, CA) using six different sequencing primers. Data was analyzed, and the novel SNP confirmed using CodonCode Aligner (CodonCode Corporation, Dedham, MA).
Alignment algorithm, matrix creation and Galaxy implementation
After standard read filtering using the Roche/454 RunProcessor software, the primary sff file was parsed by MID and amplicon sequence using the sfffile utility (454 Life Sciences, Branford, CT). Individual sff files were converted into fastq files using BioPerl, and then submitted to a local Galaxy server for further analysis using the Galaxy API to write a data library and recursively execute mapping workflows using Perl and Python. Reads were quality- and length-filtered to a standardized length (between 300 and 350 bp depending on amplicon/read direction) and then mapped to the set of HLA reference alleles using BLAT . Perfectly matching sequences were grouped, and matched alleles from each read direction were joined as described in the Results section. Positive allele rows were exported using a custom Galaxy tool, and the entire specific Galaxy workflow used is provided Additional file 2: File 1. The resulting matrix of allele matches was compared to expected match patterns to establish genotype. In some cases of low read numbers, potential homozygous loci, or novel alleles, manual interrogation of individual reads was necessary. RNA samples with poor amplifications or aberrantly sequenced reads were excluded from analysis.
Human Leukocyte Antigen
BLAST-like Alignment Tool
cDNA Sequence-Based Typing
International ImMunoGeneTics information system
Polymerase Chain Reaction
Sequence-Specific Primer assay
Sequence-Specific Oligonucleotide assay
Single Nucleotide Polymorphism
Solid Phase Reversible Immobilization
Human Immunodeficiency Virus.
Work at UW-Madison was supported in part by University of Wisconsin Foundation-Wisconsin Partnership-Medical Education and Research Committee grant number 20070709. Work at Dana Farber was supported by grant NIH-U01 A190043. Additional funding was provided by an award from the NIH, distributed by the Human Immunology Project Consortium, number 4 U01 AI08959-02.
- Caillat-Zucman S: Molecular mechanisms of HLA association with autoimmune diseases. Tissue Antigens. 2009, 73: 1-8. 10.1111/j.1399-0039.2008.01167.x.View ArticlePubMedGoogle Scholar
- Capittini C, Martinetti M, Cuccia M: MHC variation, mate choice and natural selection: the scent of evolution. Riv Biol. 2008, 101: 463-480.PubMedGoogle Scholar
- Kulpa DA, Collins KL: The emerging role of HLA-C in HIV-1 infection. Immunology. 2011, 134: 116-122. 10.1111/j.1365-2567.2011.03474.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez-Novoa S, Cuenca L, Morello J, Cordoba M, Blanco F, Jimenez-Nacher I, Soriano V: Use of the HCP5 single nucleotide polymorphism to predict hypersensitivity reactions to abacavir: correlation with HLA-B*5701. J Antimicrob Chemother. 2010, 65: 1567-1569. 10.1093/jac/dkq204.View ArticlePubMedGoogle Scholar
- Seliger B, Stoehr R, Handke D, Mueller A, Ferrone S, Wullich B, Tannapfel A, Hofstaedter F, Hartmann A: Association of HLA class I antigen abnormalities with disease progression and early recurrence in prostate cancer. Cancer Immunol Immunother. 2010, 59: 529-540. 10.1007/s00262-009-0769-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson J, Mistry K, McWilliam H, Lopez R, Parham P, Marsh SG: The IMGT/HLA database. Nucleic Acids Res. 2011, 39: D1171-D1176. 10.1093/nar/gkq998.PubMed CentralView ArticlePubMedGoogle Scholar
- Balazs I, Beekman J, Neuweiler J, Liu H, Watson E, Ray B: Molecular typing of HLA-A, -B, and DRB using a high throughput micro array format. Hum Immunol. 2001, 62: 850-857. 10.1016/S0198-8859(01)00274-9.View ArticlePubMedGoogle Scholar
- Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PI, Walker BD, Ripke S, Brumme CJ, Pulit SL, Carrington M, Kadie CM, Carlson JM, Heckerman D, Graham RR, Plenge RM, Deeks SG, Gianniny L, Crawford G, Sullivan J, Gonzalez E, Davies L, Camargo A, Moore JM, Beattie N, Gupta S, Crenshaw A, Burtt NP, Guiducci C, Gupta N, Gao X, Qi Y, et al: The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science. 2010, 330: 1551-1557.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh SG: Nomenclature for factors of the HLA system, update June 2011. Tissue Antigens. 2011, 78: 314-319. 10.1111/j.1399-0039.2011.01760.x.View ArticlePubMedGoogle Scholar
- Liao WW, Arthur JW: Predicting Peptide binding affinities to MHC molecules using a modified semi-empirical scoring function. PLoS One. 2011, 6: e25055-10.1371/journal.pone.0025055.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedrich D, Jalbert E, Dinges WL, Sidney J, Sette A, Huang Y, McElrath MJ, Horton H: Vaccine-induced HIV-specific CD8+ T cells utilize preferential HLA alleles and target specific regions of HIV-1. J Acquir Immune Defic Syndr. 2011Google Scholar
- Greene JM, Wiseman RW, Lank SM, Bimber BN, Karl JA, Burwitz BJ, Lhost JJ, Hawkins OE, Kunstman KJ, Broman KW, Wolinsky SM, Hildebrand WH, O'Connor DH: Differential MHC class I expression in distinct leukocyte subsets. BMC Immunol. 2011, 12: 39-10.1186/1471-2172-12-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim OY, Monsel A, Bertrand M, Coriat P, Cavaillon JM, Adib-Conquy M: Differential down-regulation of HLA-DR on monocyte subpopulations during systemic inflammation. Crit Care. 2010, 14: R61-10.1186/cc8959.PubMed CentralView ArticlePubMedGoogle Scholar
- Khanna R, Silins SL, Weng Z, Gatchell D, Burrows SR, Cooper L: Cytotoxic T cell recognition of allelic variants of HLA B35 bound to an Epstein-Barr virus epitope: influence of peptide conformation and TCR-peptide interaction. Eur J Immunol. 1999, 29: 1587-1597. 10.1002/(SICI)1521-4141(199905)29:05<1587::AID-IMMU1587>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Parham P: MHC class I molecules and KIRs in human history, health and survival. Nat Rev Immunol. 2005, 5: 201-214. 10.1038/nri1570.View ArticlePubMedGoogle Scholar
- Sampah ME, Ceccato CM, Blankson JN: HIV type 1-mediated downregulation of HLA-B*57/B*5801 proteins on elite suppressor CD4+ T cells. AIDS Res Hum Retroviruses. 2011, 27: 183-186. 10.1089/aid.2010.0144.PubMed CentralView ArticlePubMedGoogle Scholar
- Kloverpris HN, Stryhn A, Harndahl M, van der Stok M, Payne RP, Matthews PC, Chen F, Riddell L, Walker BD, Ndung'u T, Buus S, Goulder P: HLA-B*57 Micropolymorphism shapes HLA allele-specific epitope immunogenicity, selection pressure, and HIV immune control. J Virol. 2012, 86: 919-929. 10.1128/JVI.06150-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Lebedeva TV, Mastromarino SA, Lee E, Ohashi M, Alosco SM, Yu N: Resolution of HLA class I sequence-based typing ambiguities by group-specific sequencing primers. Tissue Antigens. 2011, 77: 247-250. 10.1111/j.1399-0039.2010.01616.x.View ArticlePubMedGoogle Scholar
- Bettencourt BF, Santos MR, Fialho RN, Couto AR, Peixoto MJ, Pinheiro JP, Spinola H, Mora MG, Santos C, Brehm A, Bruges-Armas J: Evaluation of two methods for computational HLA haplotypes inference using a real dataset. BMC Bioinformatics. 2008, 9: 68-10.1186/1471-2105-9-68.PubMed CentralView ArticlePubMedGoogle Scholar
- Listgarten J, Brumme Z, Kadie C, Xiaojiang G, Walker B, Carrington M, Goulder P, Heckerman D: Statistical resolution of ambiguous HLA typing data. PLoS Comput Biol. 2008, 4: e1000016-10.1371/journal.pcbi.1000016.PubMed CentralView ArticlePubMedGoogle Scholar
- Nunes JM, Riccio ME, Tiercy JM, Sanchez-Mazas A: Allele frequency estimation from ambiguous data: using resampling schema in validating frequency estimates and in selective neutrality testing. Hum Biol. 2011, 83: 437-447. 10.3378/027.083.0307.View ArticlePubMedGoogle Scholar
- Bentley G, Higuchi R, Hoglund B, Goodridge D, Sayer D, Trachtenberg EA, Erlich HA: High-resolution, high-throughput HLA genotyping by next-generation sequencing. Tissue Antigens. 2009, 74: 393-403. 10.1111/j.1399-0039.2009.01345.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Erlich RL, Jia X, Anderson S, Banks E, Gao X, Carrington M, Gupta N, DePristo MA, Henn MR, Lennon NJ, de Bakker PI: Next-generation sequencing for HLA typing of class I loci. BMC Genomics. 2011, 12: 42-10.1186/1471-2164-12-42.PubMed CentralView ArticlePubMedGoogle Scholar
- Lank SM, Wiseman RW, Dudley DM, O'Connor DH: A novel single cDNA amplicon pyrosequencing method for high-throughput, cost-effective sequence-based HLA class I genotyping. Hum Immunol. 2010, 71: 1011-1017. 10.1016/j.humimm.2010.07.012.PubMed CentralView ArticlePubMedGoogle Scholar
- Korlach J, Bjornson KP, Chaudhuri BP, Cicero RL, Flusberg BA, Gray JJ, Holden D, Saxena R, Wegener J, Turner SW: Real-time DNA sequencing from single polymerase molecules. Methods Enzymol. 2010, 472: 431-455.View ArticlePubMedGoogle Scholar
- Min SK, Kim WY, Cho Y, Kim KS: Fast DNA sequencing with a graphene-based nanochannel device. Nat Nanotechnol. 2011, 6: 162-165. 10.1038/nnano.2010.283.View ArticlePubMedGoogle Scholar
- Gilles A, Meglecz E, Pech N, Ferreira S, Malausa T, Martin JF: Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing. BMC Genomics. 2011, 12: 245-10.1186/1471-2164-12-245.PubMed CentralView ArticlePubMedGoogle Scholar
- Holcomb CL, Hoglund B, Anderson MW, Blake LA, Bohme I, Egholm M, Ferriola D, Gabriel C, Gelber SE, Goodridge D, Hawbecker S, Klein R, Ladner M, Lind C, Monos D, Pando MJ, Proll J, Sayer DC, Schmitz-Agheguian G, Simen BB, Thiele B, Trachtenberg EA, Tyan DB, Wassmuth R, White S, Erlich HA: A multi-site study using high-resolution HLA genotyping by next generation sequencing. Tissue Antigens. 2011, 77: 206-217. 10.1111/j.1399-0039.2010.01606.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, Nekrutenko A, Taylor J: Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol. 2010, 19: 1-21.Google Scholar
- Kent WJ: BLAT–the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh SG, Albert ED, Bodmer WF, Bontrop RE, Dupont B, Erlich HA, Fernandez-Vina M, Geraghty DE, Holdsworth R, Hurley CK, Lau M, Lee KW, Mach B, Maiers M, Mayr WR, Muller CR, Parham P, Petersdorf EW, Sasazuki T, Strominger JL, Svejgaard A, Terasaki PI, Tiercy JM, Trowsdale J: Nomenclature for factors of the HLA system, 2010. Tissue Antigens. 2010, 75: 291-455. 10.1111/j.1399-0039.2010.01466.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Koehler RN, Walsh AM, Sanders-Buell EE, Eller LA, Eller M, Currier JR, Bautista CT, Wabwire-Mangen F, Hoelscher M, Maboko L, Kim J, Michael NL, Robb ML, McCutchan FE, Kijak GH: High-throughput high-resolution class I HLA genotyping in East Africa. PLoS One. 2010, 5: e10751-10.1371/journal.pone.0010751.PubMed CentralView ArticlePubMedGoogle Scholar
- Keskin DB, Reinhold B, Lee SY, Zhang G, Lank S, O'Connor D, Berkowitz RS, Brusic V, Kim SJ, Reinherz EL: Direct identification of an HPV-16 tumor antigen from cervical cancer biopsy specimens. Frontiers in Immunology. 2011, 2:Google Scholar
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.