A registry-based population study of the HLA in Québec, Canada

As part of the worldwide effort to better characterize HLA diversity in populations, we have studied the population of Québec in Canada. This province has been defined by a complex history with multiple founder effects and migration patterns. We analyzed the typing data of 3806 individuals registered in Héma-Québec's Registry, which covered most administrative regions in Québec. Typing information was resolved at the second field level of resolution by next-generation sequencing (NGS) or by Sanger sequencing. We used the HLA-net.eu GENE[RATE] tools to estimate allele and two-locus haplotype frequencies for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1, as well as Hardy-Weinberg equilibrium (HWE), selective neutrality, and linkage disequilibrium. The chord genetic distance was also calculated between administrative regions and was visualized using non-metric multidimensional scaling (NMDS) analysis. While most individual regions were in HWE, HWE was rejected for the province considered as a whole. Some regions exhibited signatures of selection, mostly toward an excess of heterozygotes. Allele and haplotype frequencies revealed outlier regions that strongly differed from the other regions. NMDS plots also showed differences between regions. The administrative regions of the province of Québec displayed heterogeneity in their HLA profiles. This heterogeneity was attributable to differing allele and haplotype specificities by region. In particular, regions 02-Saguenay-Lac-Saint-Jean and 01-Bas-St-Laurent diverged from the rest of the regions. The urban regions 06-Montréal and 13-Laval were very diversified in their HLA profiles. Together, these results will help optimize donor recruitment strategies in Québec.

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.

Binding affinities of 438 HLA proteins to complete proteomes of seven pandemic viruses and distributions of strongest and weakest HLA peptide binders in populations worldwide

We report detailed peptide‐binding affinities between 438 HLA Class I and Class II proteins and complete proteomes of seven pandemic human viruses, including coronaviruses, influenza viruses and HIV‐1. We contrast these affinities with HLA allele frequencies across hundreds of human populations worldwide. Statistical modelling shows that peptide‐binding affinities classified into four distinct categories depend on the HLA locus but that the type of virus is only a weak predictor, except in the case of HIV‐1. Among the strong HLA binders (IC₅₀ ≤ 50), we uncovered 16 alleles (the top ones being A*02:02, B*15:03 and DRB1*01:02) binding more than 1% of peptides derived from all viruses, 9 (top ones including HLA‐A*68:01, B*15:25, C*03:02 and DRB1*07:01) binding all viruses except HIV‐1, and 15 (top ones A*02:01 and C*14:02) only binding coronaviruses. The frequencies of strongest and weakest HLA peptide binders differ significantly among populations from different geographic regions. In particular, Indigenous peoples of America show both higher frequencies of strongest and lower frequencies of weakest HLA binders. As many HLA proteins are found to be strong binders of peptides derived from distinct viral families, and are hence promiscuous (or generalist), we discuss this result in relation to possible signatures of natural selection on HLA promiscuous alleles due to past pathogenic infections. Our findings are highly relevant for both evolutionary genetics and the development of vaccine therapies. However they should not lead to forget that individual resistance and vulnerability to diseases go beyond the sole HLA allelic affinity and depend on multiple, complex and often unknown biological, environmental and other variables.