Full-length next-generation sequencing of HLA class I and II genes in a cohort from Thailand

Allele and haplotype frequencies of 11 HLA loci in Koreans by next-generation sequencing

Data on HLA genotype distribution, including DQA1 and DPA1, in the Korean population are limited. We aimed to investigate the allele and haplotype frequencies of 11 HLA loci in 339 Korean subjects using next-generation sequencing (NGS)-based HLA typing. A total of 339 samples from unrelated healthy subjects were genotyped for HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQB1, -DQA1, -DPB1, and -DPA1 using two different NGS-based HLA typing kits (166 tested using the NGSgo-MX11-3 kit [GenDx, Netherlands] and 173 by the AllType NGS 11 Loci Amplification kit [One Lambda, USA]). PyPop software was used to estimate allele and haplotype frequencies and linkage disequilibrium between the loci. Additionally, a principal component analysis was performed to compare the allele distribution of Koreans with that of other populations. A total of 214 HLA alleles (97 class I and 117 class II alleles) were assigned. The most frequent alleles for each locus were A*24:02:01 (24.78%), B*15:01:01 (10.18%), C*01:02:01 (18.44%), DRB1*04:05:01 (9.59%), DRB3*02:02:01 (13.72%), DRB4*01:03:01 (25.81%), DRB5*01:01:01 (9.0%), DQA1*01:02:01 (16.96%), DQB1*03:01:01 (14.31%), DPA1*01:03:01 (44.4%), and DPB1*05:01:01 (35.1%), respectively. The most frequent haplotypes were A*33:03:01-C*03:02:02-B*58:01:01 for HLA class I (5.01%) and DRB1*04:05:01-DQA1*03:03:01-DQB1*04:01:01-DPA1*02:02:02-DPB1*05:01:01 for HLA class II (6.23%). The total allelic ambiguities by NGS were estimated to be minimal and considerably decreased compared with those by Sanger sequencing. The Japanese population had the most similar allele distribution to Koreans, followed by the Chinese population. Frequency data of 11 HLA loci in Koreans can provide essential data for population genetics and disease association studies.

HLA-B∗46 associates with rapid HIV disease progression in Asian cohorts and prominent differences in NK cell phenotype

Human leukocyte antigen (HLA) alleles have been linked to HIV disease progression and attributed to differences in cytotoxic T lymphocyte (CTL) epitope representation. These findings are largely based on treatment-naive individuals of European and African ancestry. We assessed HLA associations with HIV-1 outcomes in 1,318 individuals from Thailand and found HLA-B∗46:01 (B∗46) associated with accelerated disease in three independent cohorts. B∗46 had no detectable effect on HIV-specific T cell responses, but this allele is unusual in containing an HLA-C epitope that binds inhibitory receptors on natural killer (NK) cells. Unbiased transcriptomic screens showed increased NK cell activation in people with HIV, without B∗46, and simultaneous single-cell profiling of surface proteins and transcriptomes revealed a NK cell subset primed for increased responses in the absence of B∗46. These findings support a role for NK cells in HIV pathogenesis, revealed by the unique properties of the B∗46 allele common only in Asia.

Associations of human leukocyte antigen with neutralizing antibody titers in a tetravalent dengue vaccine phase 2 efficacy trial in Thailand

The recombinant, live, attenuated, tetravalent dengue vaccine CYD-TDV has shown efficacy against all four dengue serotypes. In this exploratory study (CYD59, NCT02827162), we evaluated potential associations of host human leukocyte antigen (HLA) alleles with dengue antibody responses, CYD-TDV vaccine efficacy, and virologically-confirmed dengue (VCD) cases. Children 4-11 years old, who previously completed a phase 2b efficacy study of CYD-TDV in a single center in Thailand, were included in the study. Genotyping of HLA class I and II loci was performed by next-generation sequencing from DNA obtained from 335 saliva samples. Dengue neutralizing antibody titers (NAb) were assessed as a correlate of risk and protection. Regression analyses were used to assess associations between HLA alleles and NAb responses, vaccine efficacy, and dengue outcomes. Month 13 NAb log geometric mean titers (GMTs) were associated with decreased risk of VCD. In the vaccine group, HLA-DRB1*11 was significantly associated with higher NAb log GMT levels (beta: 0.76; p = 0.002, q = 0.13). Additionally, in the absence of vaccination, HLA associations were observed between the presence of DPB1*03:01 and increased NAb log GMT levels (beta: 1.24; p = 0.005, q = 0.17), and between DPB1*05:01 and reduced NAb log GMT levels (beta: -1.1; p = 0.001, q = 0.07). This study suggests associations of HLA alleles with NAb titers in the context of dengue outcomes. This study was registered with clinicaltrials.gov: NCT02827162.

Immunoinformatics Analysis of SARS-CoV-2 ORF1ab Polyproteins to Identify Promiscuous and Highly Conserved T-Cell Epitopes to Formulate Vaccine for Indonesia and the World Population

SARS-CoV-2 and its variants caused the COVID-19 pandemic. Vaccines that target conserved regions of SARS-CoV-2 and stimulate protective T-cell responses are important for reducing symptoms and limiting the infection. Seven cytotoxic (CTL) and five helper T-cells (HTL) epitopes from ORF1ab were identified using NetCTLpan and NetMHCIIpan algorithms, respectively. These epitopes were generated from ORF1ab regions that are evolutionary stable as reflected by zero Shannon's entropy and are presented by 56 human leukocyte antigen (HLA) Class I and 22 HLA Class II, ensuring good coverage for the Indonesian and world population. Having fulfilled other criteria such as immunogenicity, IFNγ inducing ability, and non-homology to human and microbiome peptides, the epitopes were assembled into a vaccine construct (VC) together with β-defensin as adjuvant and appropriate linkers. The VC was shown to have good physicochemical characteristics and capability of inducing CTL as well as HTL responses, which stem from the engagement of the vaccine with toll-like receptor 4 (TLR4) as revealed by docking simulations. The most promiscuous peptide 899WSMATYYLF907 was shown via docking simulation to interact well with HLA-A*24:07, the most predominant allele in Indonesia. The data presented here will contribute to the in vitro study of T-cell epitope mapping and vaccine design in Indonesia.

HLA-A, -B, -C, -DRB1 allele and haplotype frequencies of the Korean population and performance characteristics of HLA typing by next-generation sequencing

INTRODUCTION

Human leukocyte antigen (HLA) identification at the allelic level is important for haematopoietic stem cell transplantation (HSCT). Next-generation sequencing (NGS) resolves ambiguous alleles by determining the phase of the polymorphisms. The aim of this study was to validate the software for HLA-SBT (sequence-based typing), assess Korean allele frequency, and characterise the performance of NGS-HLA typing.

METHODS

From the 2009 to 2016 registry, 1293 unrelated healthy donors with a complete dataset of previously characterised HLA-A, -B, -C, and -DRB1 loci were selected and assessed for frequency, haplotype inference, and relative linkage disequilibrium. For performance characteristics of NGS-HLA, alleles included in 1293 cases and ambiguous or alleles assigned as new by SBT-HLA software, or unassigned alleles were included. A total of 91 and 41 quality control samples resulted in 1056 alleles (132 samples × 4 loci × 2 diploid) for analysis. The GenDx NGSgo kit was used for NGS-HLA typing using the Illumina MiSeq platform.

RESULTS

A panel of 132 samples covered 231 alleles, including 53 HLA-A, 80 HLA-B, 43 HLA-C, and 55 HLA-DRB1 by HLA-SBT typing. Comparison of SBT-HLA and NGS-HLA typing showed 99.7% (1053/1056) concordance and discrepant cases were resolved by manual evaluation. Typing by NGS resulted in 67 HLA-A, 112 HLA-B, 71 HLA-C, and 72 HLA-DRB1 alleles. A total of 132 ambiguous, 4 new, and 1 unassigned alleles by HLA-SBT were resolved by NGS-HLA typing.

CONCLUSIONS

NGS-HLA typing provided robust and conclusive results without ambiguities, and its implementation could support HSCT in clinical settings.

HLA-B*27:04 associated with enthesitis and younger age of onset, and HLA-B allele profile in patients with ankylosing spondylitis in Thailand: A cross-sectional study

AIM

The aims of this study were to estimate human leukocyte antigen (HLA)-B allele frequency, to identify alleles associated with ankylosing spondylitis (AS), and to explore manifestations in various HLA-B*27 in Thai AS patients.

METHODS

This was a cross-sectional study. Thai patients older than 18 years with diagnosed AS according to modified New York criteria who visited Siriraj Hospital (Bangkok, Thailand) were consecutively enrolled. HLA-B alleles were determined by reverse sequence-specific oligonucleotide assays, and were assigned at a 4-digit resolution. HLA-B alleles of 334 unrelated healthy Thai donors who participated in a previous phase 2b dengue vaccine clinical trial were included as controls. Odds ratio (OR) and Fisher's exact test were used to estimate association between allele and AS. The P value significance threshold was calculated according to Bonferroni.

RESULTS

Among the 88 patients who were recruited, 34 HLA-B alleles were identified, and all patients were heterozygous. The prevalence of HLA-B*27 was 89.8%, and 4 alleles of HLA-B*27 were identified. HLA-B*27:04 (OR: 39.4, P < .0001) and HLA-B*27:05 (OR: 13.8, P = .0011) were associated with AS. In contrast, HLA-B*27:06 was not found to be associated with AS (OR: 0.4, P = .241). AS patients carrying HLA-B*27:04 were more likely to have enthesitis and younger age at onset than those carrying HLA-B*27:05.

CONCLUSIONS

HLA-B*27:04 and HLA-B*27:05 were both found to be strongly associated with Thai AS. HLA-B*27:06 showed a neutral allele for Thai AS. AS patients with HLA-B*27:04 had more enthesitis and younger age at onset than those with HLA-B*27:05.

Next-generation sequencing of 11 HLA loci in a large dengue vaccine cohort from the Philippines

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4-6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.