Cost-efficient multiplex PCR for routine genotyping of up to nine classical HLA loci in a single analytical run of multiple samples by next generation sequencing

Background

HLA genotyping by next generation sequencing (NGS) requires three basic steps, PCR, NGS, and allele assignment. Compared to the conventional methods, such as PCR-sequence specific oligonucleotide primers (SSOP) and -sequence based typing (SBT), PCR-NGS is extremely labor intensive and time consuming. In order to simplify and accelerate the NGS-based HLA genotyping method for multiple DNA samples, we developed and evaluated four multiplex PCR methods for genotyping up to nine classical HLA loci including HLA-A, HLA-B, HLA-C, HLA-DRB1/3/4/5, HLA-DQB1, and HLA-DPB1.

Results

We developed multiplex PCR methods using newly and previously designed middle ranged PCR primer sets for genotyping different combinations of HLA loci, (1) HLA-DRB1/3/4/5, (2) HLA-DQB1 (3.8 kb to 5.3 kb), (3) HLA-A, HLA-B, HLA-C, and (4) HLA-DPB1 (4.6 kb to 7.2 kb). The primer sets were designed to genotype polymorphic exons to the field 3 level or 6-digit typing. When we evaluated the PCR method for genotyping all nine HLA loci (9LOCI) using 46 Japanese reference subjects who represented a distribution of more than 99.5% of the HLA alleles at each of the nine HLA loci, all of the 276 alleles genotyped, except for HLA-DRB3/4/5 alleles, were consistent with known alleles assigned by the conventional methods together with relevant locus balance and no excessive allelic imbalance. One multiplex PCR method (9LOCI) was able to provide precise genotyping data even when only 1 ng of genomic DNA was used for the PCR as a sample template.

Conclusions

In this study, we have demonstrated that the multiplex PCR approach for NGS-based HLA genotyping could serve as an alternative routine HLA genotyping method, possibly replacing the conventional methods by providing an accelerated yet robust amplification step. The method also could provide significant merits for clinical applications with its ability to amplify lower quantity of samples and the cost-saving factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1514-4) contains supplementary material, which is available to authorized users.     

Mapping of susceptibility and protective loci for acute GVHD in unrelated HLA-matched bone marrow transplantation donors and recipients using 155 microsatellite markers on chromosome 22

Despite matching donors and recipients for the human leukocyte antigens (HLAs) expressed by the major histocompatibility genomic region of the short arm of chromosome 6, several recipients still develop acute graft-versus-host disease (aGVHD) after bone marrow transplantation (BMT). This is possibly due to non-HLA gene polymorphisms, such as minor histocompatibility antigens (mHas) and genes coding for cytokines. However, a detailed genetic background for aGVHD has not yet been established. To find novel susceptibility and/or protective loci for aGVHD, a whole genome-wide association study of donors and recipients needs to be performed. As the first step to such a study, we retrospectively analyzed polymorphisms of 155 microsatellite markers spread across the long arm of chromosome 22 in 70 pairs of HLA-matched unrelated BMT donors and recipients. We performed individual typing and then compared the markers’ allele frequencies (1) between all the aGVHD (grades III and IV GVHD) and GVHD-free (grade 0 GVHD) groups in donors and recipients and (2) between the aGVHD and aGVHD-free groups in donor/recipient pairs that were matched and mismatched for the microsatellite marker’s allele. Screening of the microsatellite markers revealed five loci with a significant difference between the aGVHD and GVHD-free groups and revealed eight loci on chromosome 22, where the microsatellite allele mismatched markers were associated with aGVHD. This screening analysis suggests that several aGVHD-associated susceptible and protective loci exist on chromosome 22, which may encompass novel gene regions that need to be elucidated for their role in aGVHD.     

Comparative genomic analysis of two avian (quail and chicken) MHC regions

We mapped two different quail Mhc haplotypes and sequenced one of them (haplotype A) for comparative genomic analysis with a previously sequenced haplotype of the chicken Mhc. The quail haplotype A spans 180 kb of genomic sequence, encoding a total of 41 genes compared with only 19 genes within the 92-kb chicken Mhc. Except for two gene families (B30 and tRNA), both species have the same basic set of gene family members that were previously described in the chicken "minimal essential" Mhc. The two Mhc regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated genes with 7 class I, 10 class IIB, 4 NK, 6 lectin, and 8 B-G genes. Comparisons between the quail and chicken Mhc class I and class II gene sequences by phylogenetic analysis showed that they were more closely related within species than between species, suggesting that the quail Mhc genes were duplicated after the separation of these two species from their common ancestor. The proteins encoded by the NK and class I genes are known to interact as ligands and receptors, but unlike in the quail and the chicken, the genes encoding these proteins in mammals are found on different chromosomes. The finding of NK-like genes in the quail Mhc strongly suggests an evolutionary connection between the NK C-type lectin-like superfamily and the Mhc, providing support for future studies on the NK, lectin, class I, and class II interaction in birds.     

Rhesus macaque class I duplicon structures, organization, and evolution within the alpha block of the major histocompatibility complex

The alpha block of the human and chimpanzee major histocompatibility complex (MHC) class I genomic region contains 10 to 11 duplicated MHC class I genes, including the HLA/Patr-A, -G, and -F genes. In comparison, the alpha block of the rhesus macaque (Macaca mulatta, Mamu) has an additional 20 MHC class I genes within this orthologous region. The present study describes the identification and analysis of the duplicated segmental genomic structures (duplicons) and genomic markers within the alpha block of the rhesus macaque and their use to reconstruct the duplication history of the genes within this region. A variety of MHC class I genes, pseudogenes, transposons, and retrotransposons, such as Alu and ERV16, were used to categorize the 28 duplicons into four distinct structural categories. The phylogenetic relationship of MHC class I genes, Alu, and LTR16B sequences within the duplicons was examined by use of the Neighbor-Joining (NJ) method. Two single-duplicon tandem duplications, two polyduplicon tandem duplications with an accompanying inversion product per duplication, eight polyduplicon tandem duplications steps, 12 deletions, and at least two recombinations were reconstructed to explain the highly complex organization and evolution of the 28 duplicons (nine inversions) within the Mamu alpha block. On the basis of the phylogenetic evidence and the reconstructed tandem duplication history of the 28 duplicons, the Mamu/Patr/HLA-F ortholog was the first MHC class I gene to have been fixed without further duplication within the alpha block of primates. Assuming that the rhesus macaque and the chimpanzee/human lineages had started with the same number of MHC class I duplicons at the time of their divergence approximately 24 to 31 MYA, then the number of genes within the alpha block have been duplicated at an approximately three times greater rate in the rhesus macaque than in either the human or chimpanzee.     

Alu polymorphism within the MICB gene and association with HLA-B alleles

We describe the finding of an Alu repeat dimorphism within the first intron of the MICB gene. The frequencies of the two AluyMICB alleles, AluyMICB*0(absence of insertion) and AluyMICB*1(presence of insertion), and their associations with the highly polymorphic HLA-B locus were determined for 51 human cell lines and for 109 and 200 Caucasians and northeastern Thais, respectively. Analysis of the AluyMICB and HLA-B allelic relationships revealed that AluyMICB*1 occurred at relatively low gene frequency (0.118-0.157) [corrected] but was strongly associated with HLA-B17 (HLA-B57,HLA-B58) and HLA-B13. The AluyMICB locus provides a useful dimorphic marker for investigations on the level of linkage disequilibrium between MICB, MICA, and HLA-B loci.