Cost-efficient high-throughput HLA typing by MiSeq amplicon sequencing
- Vinzenz Lange†1,
- Irina Böhme†1Email author,
- Jan Hofmann2,
- Kathrin Lang1,
- Jürgen Sauter2,
- Bianca Schöne1,
- Patrick Paul1,
- Viviane Albrecht1,
- Johanna M Andreas1,
- Daniel M Baier2,
- Jochen Nething3,
- Ulf Ehninger2, 3,
- Carmen Schwarzelt1,
- Julia Pingel2,
- Gerhard Ehninger4 and
- Alexander H Schmidt1, 2
© Lange et al.; licensee BioMed Central Ltd. 2014
Received: 23 July 2013
Accepted: 7 January 2014
Published: 24 January 2014
A close match of the HLA alleles between donor and recipient is an important prerequisite for successful unrelated hematopoietic stem cell transplantation. To increase the chances of finding an unrelated donor, registries recruit many hundred thousands of volunteers each year. Many registries with limited resources have had to find a trade-off between cost and resolution and extent of typing for newly recruited donors in the past. Therefore, we have taken advantage of recent improvements in NGS to develop a workflow for low-cost, high-resolution HLA typing.
We have established a straightforward three-step workflow for high-throughput HLA typing: Exons 2 and 3 of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 are amplified by PCR on Fluidigm Access Array microfluidic chips. Illumina sequencing adapters and sample specific tags are directly incorporated during PCR. Upon pooling and cleanup, 384 samples are sequenced in a single Illumina MiSeq run. We developed “neXtype” for streamlined data analysis and HLA allele assignment. The workflow was validated with 1140 samples typed at 6 loci. All neXtype results were concordant with the Sanger sequences, demonstrating error-free typing of more than 6000 HLA loci. Current capacity in routine operation is 12,000 samples per week.
The workflow presented proved to be a cost-efficient alternative to Sanger sequencing for high-throughput HLA typing. Despite the focus on cost efficiency, resolution exceeds the current standards of Sanger typing for donor registration.
KeywordsHuman leukocyte antigen HLA typing NGS Dual indexing 4-primer approach Amplicon-based sequencing Fluidigm Access Array Illumina MiSeq
The success of unrelated hematopoietic stem cell transplantation is closely linked to a good match of the HLA genes between donor and recipient . To overcome the odds of finding a matching donor – given the high diversity of the HLA system (>8000 described alleles) – many organizations worldwide endeavor to characterize the HLA types of volunteers. In 2012 alone, 1.6 million new potential donors were registered. Typing potential donors upfront for 5 loci at high resolution would be highly beneficial in speeding up the search process, thereby saving precious time for patients in need of a transplantation . However, due to the cost of high resolution HLA typing based on Sanger sequencing, many organizations restrict the number of loci typed to HLA-A,-B and-DRB1 or even restrict the number of potential donors to be recruited. After revolutionizing genome sequencing, NGS is expected to also markedly change the diagnostic market. HLA typing by NGS has been demonstrated by performing long range PCR in combination with shotgun sequencing [3–5]. However, the required sophisticated sample preparation steps and associated costs make that approach currently more attractive for patient samples than for cost-sensitive high-throughput stem cell donor registry typing. Several groups have demonstrated the feasibility of HLA typing by short amplicon sequencing on the 454 system [6, 7]. Moonsamy et al. streamlined the workflow considerably by taking advantage of the Fluidigm Access Array technology for performing up to 48 × 48 (2048) PCR reactions on a single chip . Given the improvements in read length to 2 × 250 base pairs, the Illumina MiSeq system seemed compatible for an amplicon sequencing approach as well. We have therefore set out to transfer the strategy of Moonsamy et al. to the MiSeq. In contrast to the 454 system, the MiSeq features on-board clonal amplification, thereby greatly simplifying sample preparation and rendering the workflow more amendable for high-throughput routine operation. In addition, the operational cost is greatly reduced, with a more than 50 fold reduced cost per read compared with the 454 GS FLX system. Taken together, those features allowed us to develop a high-throughput low-cost HLA sequencing workflow that would render typing of the transplantation-relevant genes and exons affordable for donor centers.
Results and discussion
Optimization of primers and PCR conditions
In contrast to Sanger based sequencing, using NGS unspecific amplification does not result in noisy sequencing data. Due to the clonal sequencing, the unspecific amplification products are easily identified by the analysis software. Nevertheless, specificity is preferred to maximize the relevant information of the available reads.
Likewise, with regard to unbalanced amplification, the NGS approach is considerably more sensitive in detecting suppressed alleles than Sanger sequencing which requires at least 5% to 20% of the weaker allele for detection . In contrast, given sufficient read depth, our analysis software “neXtype” identifies suppressed alleles down to 2% read count relative to the dominant allele. In addition, to avoid false homozygous results, neXtype blocks homozygous calls with less than 100 reads (20 fold above the detection limit). Despite the increased sensitivity for suppressed alleles, we went through several rounds of primer optimization to come up with a set of primers with well-balanced amplification independent of the given allele combination (Figure 5 b,c). As part of our quality control, every new batch of primers is tested against a set of 95 samples to confirm well-balanced amplification, i.e. no second allele detected with less than 20% of the reads of the dominant allele. This set of samples was carefully selected to represent at least 2 samples for each allele group targeted by a particular primer in the primer group.
We also monitor the occurrence of artifacts previously reported in the literature: The generation of artificial “recombinant” or “chimeric” PCR products from the two alleles present in a sample. Such “crossover” events have been described as potential error source of HLA-DRB1 typing, as particular crossover events of HLA-DRB1 and HLA-DRB3/4/5 are identical to named HLA-DRB1 alleles [14–16]. We quantified the rate of crossover products (Figure 5a). Depending on the targeted exon, an average of up to 25% of the matched reads were identified as crossover products (Figure 5d). While we are not aware of previous quantitative analyses in the context of HLA, our findings are nonetheless in accordance with reports in the context of 16S sequencing . Since most of the crossover products are generated during the late PCR cycles, the artifact may be reduced by decreasing the number of PCR cycles  (Figure 5d). However, in the interest of maximum robustness with regard to samples with low DNA concentration, we have chosen not to reduce PCR cycles below 32 cycles but have rather extended neXtype to handle the artifacts properly.
HLA allele calling
We developed the NGS HLA typing software “neXtype” to match the throughput of the workflow. NeXtype takes full advantage of Q values and high coverage to automate data analysis as far as possible without compromising accuracy. Statistical algorithms are implemented to a) identify new alleles, and b) distinguish between sequencing artifacts and closely related alleles. Currently, more than 95% of the typing assignments do not require any manual corrections. Taking advantage of the high number of reads per amplicon and the high quality of base-calls, we have chosen a set of highly conservative parameters for “autotyping”: Those typing assignments (82%) do not even need to be reviewed by a human analyst. To qualify for this category (among other restrictions), homozygous calls are accepted only if the read count is high enough to safely exclude a potentially suppressed second allele, taking imbalanced amplification of up to 20:1 into account. The 18% of typings not meeting the autotyping criteria remain for analyst review and are divided into three categories: 13% require user confirmation, 2.7% requiring user editing and 2.3% failed typings. Due to the high level of automation, analysis of 384 samples (2304 loci) can be performed in less than 3 hours by one trained analyst.
We validated the workflow including neXtype analysis software by analyzing 1140 samples of known HLA type. All common alleles with frequency >0.001 in the German population were represented by at least 2 samples . Using 24 Fluidigm Access Array chips we amplified 12 HLA targeted exons per sample and sequenced them in 3 MiSeq runs. In the case of 160 sample-loci (2.3%), data quality was insufficient and these were therefore excluded from further analysis. All alleles of the 160 sample-loci that failed were successfully typed by this approach in other samples. Therefore, these failures were classified as technical failures independent of HLA type. In routine operation those technical failures need to be repeated. Those failures do however not compromise HLA typing accuracy. Subsequent to analysis using neXtype, the results were compared with the typing results obtained by Sanger sequencing: All NGS neXtype results were concordant with the Sanger sequences, demonstrating error-free typing for more than 6000 sample-loci.
Resolution in comparison with Sanger based sequencing
Comparing the proposed workflow with conventional Sanger sequencing, it is apparent that the NGS workflow is more cost-efficient and easier to set up. In contrast to Sanger sequencing, the individual sample needs to be handled only when it is combined with the barcode sequence, loaded onto the Access Array chip and pooled with other samples after PCR. To exclude sample mix up by design, we have set up 2 liquid handlers with barcode reading capabilities. However, the process as such could easily be handled by manual pipetting. Therefore, it appears to be an attractive alternative to Sanger sequencing even for labs with far lower throughput. Given the combined chip and sequencing run costs of about $1000 (list price Access Array chip and MiSeq Reagent Nano Kit 500 cycles), HLA typing by NGS as proposed here seems to be more cost-efficient, starting with as few as 24 samples, even when compared with Sanger sequencing costs alone. At the same time, it was possible to further increase resolution by sequencing more exons and introns of up to 48 targeted amplicons at no significant additional cost. For the purpose of typing for HLA registries, we have restricted the number of amplicons to the exons most relevant for the donor search to routinely multiplex 384 samples. In high throughput operation, we achieved a cost reduction of more than 50% compared with cost-optimized Sanger sequencing. These realized savings lower the cost for high resolution typing of volunteer donors, thus increasing the number of registered donors with the same budget and simultaneously improving typing resolution. Both factors should facilitate and accelerate the search process and thus improve the outcome of unrelated hematopoietic stem cell transplantation.
All samples used were blood samples, collected in venosafe 4 ml EDTA tubes (Terumo, Tokyo, Japan). Validation samples were chosen randomly from newly recruited DKMS donors in Germany. Ethnic background is self-assigned at recruitment by country code of origin. 85% of the donors described themselves as German, 7% as Turkish, 1% as polish donors, 5% came from 33 different countries and from 2% no ethnic information is reported. Upon registration, donors sign an informed consent approving HLA typing and storage of samples for future typings related to stem cell transplantation. No ethics committee approval was obtained as sequencing-based HLA typing of potential stem cell donors is standard practice for stem cell donor centers and is covered by the donor consent form signed at recruitment. DKMS Life Science Lab is an affiliated company with DKMS German Bone Marrow Donor Center and performs HLA typing on request of DKMS.
All samples have been typed for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 by Sanger sequencing using inhouse primers on a 3730xl capillary sequencer (Life Technologies, Carlsbad, USA) at the time of donor registration. These typings have been used for validation. The Sanger sequencing strategy involved amplification spanning exons 2 and 3 for HLA-A, -B, -C; two separate amplicons for exons 2 and 3 for HLA-DQB1 and -DPB1 and one amplicon for Exon 2 for HLA-DRB1. All primers were located in the introns. Sequencing was carried out separately for all exons, and group specific sequencing primers (GSSP) were used to resolve ambiguities.
DNA was isolated from 150 μl whole blood (EDTA) based on magnetic bead technology with “chemagic DNA Blood kit 150 special” (PerkinElmer, Baesweiler, Germany) using “chemagic Magnetic Separation Module I” (PerkinElmer, Baesweiler, Germany). DNA was eluted in 100 μl elution buffer (10 mM Tris–HCl pH8.0) yielding concentrations of between 35 and 100 ng/μl as determined by UV on NanoQuant plate using TECAN infinite 200Pro (Tecan, Männedorf, Switzerland) plate reader.
PCR amplification with the access array system
PCR amplification was performed using a 48.48 Fluidigm Access Array (Fluidigm Corporation, South San Francisco, USA) in combination with the Roche High Fidelity Fast Start Kit (Roche, Basel, Switzerland).
A 10 μl PCR master mix was prepared for each sample and contained 1 μl 10x buffer mix without MgCl2 (Roche High Fidelity Fast Start Kit), 1.8 μl 25 mM MgCl2, 0.5 μl DMSO, 0.5 μl 20xAA Loading Reagent (Fluidigm Corporation, South San Francisco, USA), 0.2 μl 10 mM dNTPs each (Roche High Fidelity Fast Start Kit), 0.1 μl Fast Start High Fidelity Taq Polymerase (5 U/μl) (Roche High Fidelity Fast Start Kit) and 1.9 μl PCR grade water. 6 μl PCR master mix, 2 μl DNA and 2 μl barcode primers (2 μM equimolar mix of index 1 and index 2, Additional file 1) were mixed in a 384 well plate using a Tecan EVO robotic station (Tecan, Männedorf, Switzerland). 4.5 μl of the mix were transferred to the sample inlets of a primed 48.48 Fluidigm Access Array by the Tecan EVO. 4 μl target specific primers (1 μM in 1×AA loading reagent) were transferred to the primer inlets. Target specific and index primer were obtained from metabion (metabion international AG, Martinsread, Germany).
After loading, the chip was placed in an IFC Controller AX (Fluidigm Corporation, South San Fransisco, USA) for loading the PCR reagents before the PCR was performed on a FC1 cycler. Thermal profile was 50°C for 2 min, 70°C for 20 min, 95°C for 10 min, followed by 32 cycles at 95°C for 25 s, 60°C for 30 s and 72°C for 90 s, and a finishing step at 72°C for 5 min. The PCR products were harvested from the PCR chambers by an IFC Controller AX (Fluidigm Corporation, South San Francisco, USA).
Dual indexing primers
Please refer to the supplementary information for the concept, indices and sample sheet information.
Amplicon pooling, purification and quantification
The 48 samples of an Access Array chips were pooled. 8 pools were purified in parallel using AMPure XP (Beckman Coulter, Brea, USA) with a ratio of 0.7:1 beads to DNA. Purified amplicon pools were diluted 1:4000 for quantification by qPCR using an ECO Real-Time PCR cycler (Illumina, San Diego, USA) and the Library Quant Illumina Kit (KAPA Biosystems, Boston, USA) with standards in a range from 0.2 fM to 20 pM. Pooling, purification and dilution for quantification were performed on a Biomek 3000 workstation (Beckman Coulter, Brea, USA).
Library preparation and MiSeq run
8 purified and quantified amplicon pools were mixed in equimolar amounts. The library was prepared as recommended by Illumina (MiSeq Reagent Kit v2 – Reagent Preparation Guide) and was loaded at 7.5 pM on one MiSeq flowcell with 10% phiX spiked in. Paired end sequencing was performed with 251 cycles.
Prior to the update of Real-Time Analysis Software (Illumina, San Diego, USA) to version 18.104.22.168, hardcoded matrix and phasing values were used for sequencing runs.
Classification of reads for primer optimization
Primer sets to be tested were combined with samples of known HLA types on the Fluidigm chip and processed as described above. To assess the quality of the primer sets, the resulting MiSeq reads were classified in two ways: 1) Finding the most similar HLA allele(s) amongst all known HLA alleles by simple 10mer comparison. Up to 10 mismatches were allowed when considering the read as “matched” to a reference allele. Reads with more than 10 mismatches were classified as “not matched”. 2) Base by base comparison of reads with all known HLA alleles. The classification was based on the Q values of the mismatches. If the sum of the Q values of non-matching bases (Qsum) was >80, the read was regarded as “not matching”. If the read matched considerably worse (Qsum difference >50) to the reference allele than to an allele as determined by approach 1, the read was regarded as “other HLA allele”. If the read matched to both alleles on the basis of the above criteria, the read was assigned to the allele with the lower Qsum. Reads fulfilling the matching criteria were tested to determine whether existing mismatches could be explained by PCR crossover, i.e. artificial recombinations of the two alleles during PCR. Those reads were classified as “Allele n with crossover” if at least two bases with good Q value (>25) matched to the alternate allele. Depending on PCR conditions and primer set, a considerable portion of the reads could be identified as “crossover” products. Reads fulfilling the matching criteria were tested to determine whether existing mismatches could be explained by PCR crossover. Those reads were classified as “Allele n with crossover” if at least two bases with good Q value (>25) matched to the alternate allele. “Not matched” reads were tested base by base against allele1 or allele2. If matching criteria were met, those reads were classified as “crossover”. In order to assess allele balancing, we calculated the ratio of the two alleles based on the classification for each primer set.
neXtype is designed for HLA typing from raw data provided by the Illumina MiSeq system. The methods implemented within neXtype take advantage of the large number of reads per sequence and the associated quality values (Q values). The software is specifically adapted to a paired-end amplicon sequencing approach. NeXtype is implemented as an Oracle PL/SQL application and scalable to any number of CPU cores and Oracle instances.
In a pre-typing run, neXtype assembles the input data that is independent of each specific run, such as HLA reference sequences from the IMGT/HLA database and primer information. The typing run itself is organized into several consecutive steps: primer recognition, allele matching, result classification, exon combination, result rating, and user interaction.
Primer recognition is carried out for each MiSeq read to assign HLA locus and exon. Primer recognition involves three levels:
Firstly, a string comparison between primer library and read is carried out.
If no match can be found then the read is compared to known artifacts of the primer set collected in a look-up table.
In typical cases, for each locus and exon more than two EAG are identified. NeXtype ranks these EAGs by score, i.e. number of “best match” reads in descending order. The possible classification of an EAG is as follows: result, potential new allele, noise, cross-over, and co-amplification. A base-by-base comparison of all “best match” reads to the EAG allows identification of potentially new alleles by identifying systematically occurring SNPs and evaluating the number of reads with this alteration. Likewise, the consensus sequence is determined based on the majority of base calls at each position. Using a binomial distribution, EAGs are checked pairwise to determine whether their occurrence can be explained as noise of higher scoring EAGs, given the number of reads and their Q values. EAGs are classified as “cross-over” if their consensus sequence can be obtained by combination of any two higher scoring EAGs. The consensus sequence is matched against the HLA system to identify off target PCR products (e.g. HLA-H, -DRB3, …). Those are marked as “co-amplificate”. EAGs not scoring at least 1/50 of the best scoring EAG are discarded. All EAGs not classified in the above categories are considered “results” by neXtype. If more than two EAGs are classified as “result”, neXtype indicates a warning or error depending on the relative number of “best matches” of the EAG.
The final HLA assignments are obtained by combining the “results” from each exon. Only alleles that can be found in the two EAGs classified as “result” with the highest number of “best matches” in each exon become part of the final typing result. In addition to EAGs classified as “result”, EAGs classified as “co-amplificate” are considered if they are needed to generate a valid result. Exons with less than 20 reads in sum over all “result” EAGs are ignored. The collected alleles are coded into NMDP multi-allele or G codes [20, 21]. In cases where results of exon 2 and exon 3 lead to more than one possible combination of exons (exon-shuffling), not all ambiguities can be resolved, thus resulting in intermediate resolution typing results.
To rate the quality of HLA assignments, neXtype uses several parameters, including the percentage of successful primer matches, the share of reads with a successful association to an EAG, and the percentage of reads in an exon which contribute to the final result in this specific exon. Each parameter has two predefined thresholds allowing classification of the HLA assignment’s quality. If a parameter exceeds the upper threshold, it scores 1, whereas if it falls below the lower threshold, it scores 0 and between the two thresholds the score is linearly interpolated between 0 and 1. The final rating is the product over all parameter scores. Furthermore, exons are checked for possible suppressed alleles using different thresholds depending on whether a homozygous or heterozygous result is retained.
All classified EAGs are displayed for each exon. The user is permitted to include or exclude EAGs from each exon based on his assessment and to recombine the exons for a final typing result. Any user interaction triggers recomputing of result ratings. The user can further mark individual EAGs as known artifacts that should be ignored in future runs. When the user agrees with the HLA assignment by neXtype, a status “finished” can be set and the result is submitted electronically. HLA assignments achieving an overall rating of 1 are automatically submitted.
Human leukocyte antigen
Next generation sequencing
Solid phase reversible immobilization
Group specific sequencing primer.
- Petersdorf EW, Anasetti C, Martin PJ, Gooley T, Radich J, Malkki M, Woolfrey A, Smith A, Mickelson E, Hansen JA: Limits of HLA mismatching in unrelated hematopoietic cell transplantation. Blood. 2004, 104: 2976-2980. 10.1182/blood-2004-04-1674.PubMedView ArticleGoogle Scholar
- Tu B, Cha N, Yang R, Ng J, Hurley CK: A one-step DNA sequencing strategy to HLA type hematopoietic stem cell donors at recruitment - rethinking typing strategies. Tissue Antigens. 2013, 81: 150-160. 10.1111/tan.12072.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang C, Krishnakumar S, Wilhelmy J, Babrzadeh F, Stepanyan L, Su LF, Levinson D, Fernandez-Viña MA, Davis RW, Davis MM, Mindrinos M: High-throughput, high-fidelity HLA genotyping with deep sequencing. Proc Natl Acad Sci USA. 2012, 109: 8676-8681. 10.1073/pnas.1206614109.PubMed CentralPubMedView ArticleGoogle Scholar
- Shiina T, Suzuki S, Ozaki Y, Taira H, Kikkawa E, Shigenari A, Oka A, Umemura T, Joshita S, Takahashi O, Hayashi Y, Paumen M, Katsuyama Y, Mitsunaga S, Ota M, Kulski JK, Inoko H: Super high resolution for single molecule-sequence-based typing of classical HLA loci at the 8-digit level using next generation sequencers. Tissue Antigens. 2012, 80: 305-316. 10.1111/j.1399-0039.2012.01941.x.PubMedView ArticleGoogle Scholar
- Lind C, Ferriola D, Mackiewicz K, Heron S, Rogers M, Slavich L, Walker R, Hsiao T, McLaughlin L, D’Arcy M, Gai X, Goodridge D, Sayer D, Monos D: Next-generation sequencing: the solution for high-resolution, unambiguous human leukocyte antigen typing. Hum Immunol. 2010, 71: 1033-1042. 10.1016/j.humimm.2010.06.016.PubMedView ArticleGoogle 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 CentralPubMedView ArticleGoogle Scholar
- Gabriel C, Danzer M, Hackl C, Kopal G, Hufnagl P, Hofer K, Polin H, Stabentheiner S, Pröll J: Rapid high-throughput human leukocyte antigen typing by massively parallel pyrosequencing for high-resolution allele identification. Hum Immunol. 2009, 70: 960-964. 10.1016/j.humimm.2009.08.009.PubMedView ArticleGoogle Scholar
- Moonsamy PV, Williams T, Bonella P, Holcomb CL, Höglund BN, Hillman G, Goodridge D, Turenchalk GS, Blake LA, Daigle DA, Simen BB, Hamilton A, May AP, Erlich HA: High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ system for simplified amplicon library preparation. Tissue Antigens. 2013, 81: 141-149. 10.1111/tan.12071.PubMedView ArticleGoogle Scholar
- Kircher M, Sawyer S, Meyer M: Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res. 2012, 40: e3-10.1093/nar/gkr771.PubMed CentralPubMedView ArticleGoogle Scholar
- Bochtler W, Maiers M, Bakker JNA, Oudshoorn M, Marsh SGE, Baier D, Hurley CK, Müller CR: World Marrow Donor Association framework for the implementation of HLA matching programs in hematopoietic stem cell donor registries and cord blood banks. Bone Marrow Transplant. 2011, 46: 338-343. 10.1038/bmt.2010.132.PubMedView ArticleGoogle Scholar
- Hurley CK, Oudshoorn M, Setterholm M: Donor registries and search strategies. Methods Mol Biol. 2012, 882: 531-547. 10.1007/978-1-61779-842-9_30.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiao Y, Lazaro AM, Masaberg C, Haagenson M, Vierra-Green C, Spellman S, Dakshanamurthy S, Ng J, Hurley CK: Evaluating the potential impact of mismatches outside the antigen recognition site in unrelated hematopoietic stem cell transplantation: HLA-DRB1*1454 and DRB1*140101. Tissue Antigens. 2009, 73: 595-598. 10.1111/j.1399-0039.2009.01245.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Davidson CJ, Zeringer E, Champion KJ, Gauthier M-P, Wang F, Boonyaratanakornkit J, Jones JR, Schreiber E: Improving the limit of detection for Sanger sequencing: a comparison of methodologies for KRAS variant detection. Biotechniques. 2012, 53: 182-188.PubMedGoogle Scholar
- Danzer M, Niklas N, Stabentheiner S, Hofer K, Pröll J, Stückler C, Raml E, Polin H, Gabriel C: Rapid, scalable and highly automated HLA genotyping using next-generation sequencing: a transition from research to diagnostics. BMC Genomics. 2013, 14: 221-10.1186/1471-2164-14-221.PubMed CentralPubMedView ArticleGoogle Scholar
- Erlich H: HLA DNA typing: past, present, and future. Tissue Antigens. 2012, 80: 1-11. 10.1111/j.1399-0039.2012.01881.x.PubMedView ArticleGoogle Scholar
- Rastrou MV, Holcomb CL, Williams TC, Goodridge D, Lazaro AM, Tilanus M, Erlich HA: 12-OR: In vitro generated PCR crossover products identified by next generation sequencing correspond to some alleles in the IMGT database. Hum Immunol. 2012, 73: 10-View ArticleGoogle Scholar
- Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, Methé B, DeSantis TZ, Petrosino JF, Knight R, Birren BW: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21: 494-504. 10.1101/gr.112730.110.PubMed CentralPubMedView ArticleGoogle Scholar
- Judo MSB, Wedel AB, Wilson C: Stimulation and suppression of PCR-mediated recombination. Nucleic Acids Res. 1998, 26: 1819-1825. 10.1093/nar/26.7.1819.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmidt AH, Baier D, Solloch UV, Stahr A, Cereb N, Wassmuth R, Ehninger G, Rutt C: Estimation of high-resolution HLA-A, -B, -C, -DRB1 allele and haplotype frequencies based on 8862 German stem cell donors and implications for strategic donor registry planning. Hum Immunol. 2009, 70: 895-902. 10.1016/j.humimm.2009.08.006.PubMedView ArticleGoogle Scholar
- Robinson J, Halliwell JA, McWilliam H, Lopez R, Parham P, Marsh SGE: The IMGT/HLA database. Nucleic Acids Res. 2013, 41 (Database issue): D1222-D1227.PubMed CentralPubMedView ArticleGoogle Scholar
- Marsh SGE, Albert ED, Bodmer WF, Bontrop RE, Dupont B, Erlich HA, Fernández-Viña M, Geraghty DE, Holdsworth R, Hurley CK, Lau M, Lee KW, Mach B, Maiers M, Mayr WR, Müller 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 CentralPubMedView ArticleGoogle 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.