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.