A novel framework for human leukocyte antigen (HLA) genotyping using probe capture-based targeted next-generation sequencing and computational analysis.

Human leukocyte antigen (HLA) genes play pivotal roles in numerous immunological applications. Given the immense number of polymorphisms, achieving accurate high-throughput HLA typing remains challenging. This study aimed to harness the human pan-genome reference consortium (HPRC) resources as a potential benchmark for HLA reference materials. We meticulously annotated specific four field-resolution alleles for 11 HLA genes (HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3, -DRB4 and -DRB5) from 44 high-quality HPRC personal genome assemblies. For sequencing, we crafted HLA-specific probes and conducted capture-based targeted sequencing of the genomic DNA of the HPRC cohort, ensuring focused and comprehensive coverage of the HLA region of interest. We used publicly available short-read whole-genome sequencing (WGS) data from identical samples to offer a comparative perspective. To decipher the vast amount of sequencing data, we employed seven distinct software tools: OptiType, HLA-VBseq, HISAT genotype, SpecHLA, T1K, QzType, and DRAGEN. Each tool offers unique capabilities and algorithms for HLA genotyping, allowing comprehensive analysis and validation of the results. We then compared these results with benchmarks derived from personal genome assemblies. Our findings present a comprehensive four-field-resolution HLA allele annotation for 44 HPRC samples. Significantly, our innovative targeted next-generation sequencing (NGS) approach for HLA genes showed superior accuracy compared with conventional short-read WGS. An integrated analysis involving QzType, T1K, and DRAGEN was developed, achieving 100% accuracy for all 11 HLA genes. In conclusion, our study highlighted the combination of targeted short-read sequencing and astute computational analysis as a robust approach for HLA genotyping. Furthermore, the HPRC cohort has emerged as a valuable assembly-based reference in this realm.

Profiling genes encoding the adaptive immune receptor repertoire with gAIRR Suite.

Adaptive immune receptor repertoire (AIRR) is encoded by T cell receptor (TR) and immunoglobulin (IG) genes. Profiling these germline genes encoding AIRR (abbreviated as gAIRR) is important in understanding adaptive immune responses but is challenging due to the high genetic complexity. Our gAIRR Suite comprises three modules. gAIRR-seq, a probe capture-based targeted sequencing pipeline, profiles gAIRR from individual DNA samples. gAIRR-call and gAIRR-annotate call alleles from gAIRR-seq reads and annotate whole-genome assemblies, respectively. We gAIRR-seqed TRV and TRJ of seven Genome in a Bottle (GIAB) DNA samples with 100% accuracy and discovered novel alleles. We also gAIRR-seqed and gAIRR-called the TR and IG genes of a subject from both the peripheral blood mononuclear cells (PBMC) and oral mucosal cells. The calling results from these two cell types have a high concordance (99% for all known gAIRR alleles). We gAIRR-annotated 36 genomes to unearth 325 novel TRV alleles and 29 novel TRJ alleles. We could further profile the flanking sequences, including the recombination signal sequence (RSS). We validated two structural variants for HG002 and uncovered substantial differences of gAIRR genes in references GRCh37 and GRCh38. gAIRR Suite serves as a resource to sequence, analyze, and validate germline TR and IG genes to study various immune-related phenotypes.

Response to hepatitis B vaccination is co-determined by HLA-DPA1 and -DPB1.

BACKGROUND AND AIMS: No report explored the combined effects of HLA-DPA1 and -DPB1 with long-term response to hepatitis B (HB) vaccination (HBVac). The specific aims of the study were to assess the combined effects and relative contributions of DPA1 and DPB1 genes. METHODS: The cases were 152 adolescents who had undetectable (<1.0 mIU/mL) post-booster anti-HBs titers and the controls were adolescents who had residual anti-HBs ≥ 10 mIU/mL at aged 16 years (n = 207) or had detectable (≥1.0 mIU/mL) anti-HBs titers after booster HBVac (n = 481). HLA-DPA1 and -DPB1 genotypes were determined by sequence-based typing. RESULTS: HLA-DPA1*01:03:01 was correlated with lower ORs of undetectable anti-HBs titers, while -DPA1*02:02:02 and -DPB1*05:01:01 were correlated with higher ORs. The ORs for HLA-DPA1*01:03:01-DPB1*05:01:01 and DPA1*02:02:02-DPB1*protective combinatory types were significantly less than 1.0. As compared with subjects who had no protective allele, the adjusted ORs (95% CI) were 0.545 (0.328-0.906), 0.350 (0.174-0.702), and 0.122 (0.058-0.257), for subjects who had protective alleles on DPA1only, DPB1 only, and both genes, respectively. Analyses of amino acid polymorphisms showed that subjects who carried Arg81-Pro158-Val191-Pro259α + Met234ß and Gln62-Arg82α + Met234ß combinations had 4.3-to-4.6 folds of risks. CONCLUSION: Both DPA1 and DPB1 genes contribute to the persistence of immunological response to primary infantile HBVac. The effects of HLA-DP risk alleles were dominated by the protective alleles and there were significant gene-gene interactions. Our findings provide evidences for the design of more potent HB vaccine.

SNP rs7770370 in HLA-DPB1 loci as a major genetic determinant of response to booster hepatitis B vaccination: results of a genome-wide association study.

BACKGROUND AND AIM: Hepatitis B (HB) vaccination is highly effective in reducing the risk of hepatitis B virus infection. However, breakthrough and chronic hepatitis B virus infections in vaccinated subjects raised concern about its long-term efficacy. The specific aim of the study was to explore the host genetic determinants of long-term immunological memory against HB vaccination. METHODS: We conducted a case-control study nested in a cohort of HB booster recipients who had received primary HB vaccination during infancy but failed to reside an anti-HBs titers ≥ 10 mIU/mL at the age of 15-18 years. We used a genome-wide single nucleotide polymorphism (SNP) array plate to scan autosomal chromosomes and assayed the human leukocyte antigen (HLA)-DPB1 genotype by sequence-based techniques. RESULTS: We found that 10 of the 112 candidate SNPs (P-value < 5.0 × 10(-5) ) clustered within a 47-Kb region of the HLA-DP loci. All the minor alleles of these HLA-DP candidate SNPs were correlated with lower likelihoods of nonresponse to HB vaccine. There was a significant linkage disequilibrium between these HLA-DP candidate SNPs and HLA-DPB1 protective alleles. Multivariate analyses showed that rs7770370 was the most significant genetic factor. As compared with rs7770370 GG homozygotes, adjusted odds ratios were 0.524 (95% confidence interval, 0.276-0.993) and 0.095 (95% confidence interval, 0.030-0.307) for AG heterozygotes and AA homozygotes, respectively. CONCLUSION: Our results showed that rs7770370 was the most significant genetic factor of response to HB booster. The rs7770370 and nearby SNPs may also contribute to the long-term immunological memory against HB vaccination.