Demographic history and selection at HLA loci in Native Americans

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

Honduras HIV cohort: HLA class I and CCR5-Δ32 profiles and their associations with HIV disease outcome

IMPORTANCE

We identify both canonical and novel human leukocyte antigen (HLA)-HIV associations, providing a first step toward improved understanding of HIV immune control among the understudied Honduras Mestizo population. Our results are relevant to understanding the protective or detrimental effects of HLA subtypes in Latin America because their unique HLA diversity poses challenges for designing vaccines against HIV and interpreting results from such vaccine trials. Likewise, the description of the HLA profile in an understudied population that shows a unique HLA immunogenetic background is not only relevant for HIV immunology but also relevant in population genetics, molecular anthropology, susceptibility to other infections, autoimmune diseases, and allograft transplantation.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

HLA allele and haplotype frequencies of registered stem cell donors in Chile

Patients in need of hematopoietic stem cell transplantation often rely on unrelated stem cell donors matched in certain human leukocyte antigen (HLA) genes. Donor search is complicated by the extensive allelic variability of the HLA system. Therefore, large registries of potential donors are maintained in many countries worldwide. Population-specific HLA characteristics determine the registry benefits for patients and also the need for further regional donor recruitment. In this work, we analyzed HLA allele and haplotype frequencies of donors of DKMS Chile, the first Chilean donor registry, with self-assessed "non-Indigenous" (n=92,788) and "Mapuche" (n=1,993) ancestry. We identified HLA alleles that were distinctly more abundant in the Chilean subpopulations than in worldwide reference populations, four of them particularly characteristic for the Mapuche subpopulation, namely B*39:09g, B*35:09, DRB1*04:07g, and DRB1*16:02g. Both population subsamples carried haplotypes of both Native American and European origin at high frequencies, reflecting Chile's complex history of admixture and immigration. Matching probability analysis revealed limited benefits for Chilean patients (both non-Indigenous and Mapuche) from donor registries of non-Chilean donors, thus indicating a need for ongoing significant donor recruitment efforts in Chile.

Putative role of HLA polymorphism among a Brazilian HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) population

Around ten million people are infected with HTLV-1 worldwide, and 1-4% develop HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP), characterized by an important degeneration of the spinal cord, which can lead to death. Distinct HLA alleles have been associated with either HAM/TSP susceptibility or protection. However, these HLA alleles set may change according to the population studied. Brazil is the second country in the number of HTLV-1-infected people and there are few reports addressing the HLA influence on HTLV-1 infection as well as on disease outcome. The objective of this study was to evaluate the influence of HLA alleles as a risk factor for HAM/TSP and the proviral load (PVL) levels, clinical progression, and death outcomes in an admixed Brazilian population. The HLA-A, -B, -C, and -DRB1 were genotyped in 375 unrelated HTLV-1-infected individuals divided into asymptomatic carriers (AC) (n = 165) and HAM/TSP (n = 210) in a longitudinal cohort from 8 to 22 years of follow-up. Because locus B deviated from Hardy-Weinberg Equilibrium for the study groups, the results represented for HLA-B alleles were inconclusive. The alleles HLA-A*68 and -C*07 were related to HAM/TSP risk in multivariate analysis. The alleles HLA-A*33, and -A*36 were associated with protection against disease progression in HAM/TSP patients, while -C*12, -C*14, and -DRB1*08 were associated with increased risk of death. In the AC group, the presence of, -C*06 and -DRB1*15 alleles influenced an increased PVL, in an adjusted linear regression model, while -A*30, -A*34, -C*06, -C*17 and -DRB1*09 alleles were associated with increased PVL in HAM/TSP group compared to HAM/TSP individuals not carrying these alleles. All these alleles were also related to increased PVL associated with clinical progression outcome. Increased PVL associated with the death outcome was linked to the presence of HLA-A*30. PVL has been associated with HLA, and several alleles were related in AC and HAM/TSP patients with or without interacting with clinical progression outcomes. Understanding the prognostic value of HLA in HAM/TSP pathogenesis can provide important biomarkers tools to improve clinical management and contribute to the discovery of new therapeutic interventions.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing D_ST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.

Protective HLA alleles against severe COVID-19: HLA-A*68 as an ancestral protection allele in Tapachula-Chiapas, Mexico

HLA is a polymorphic antigen presenter which has provided valuable information on the susceptibility of populations to viruses. Therefore, the study of HLA can reveal specific susceptibility or resistance alleles to severe COVID-19 in an ethnically dependent manner. This pilot study investigated HLA alleles associated with COVID-19 severity in Tapachula, Chiapas, Mexico. A total of 146 Mexican Mestizos were typed for HLA class I and II using PCR-SSP. The patients were classified according to the outcome (death or improvement) and the infection's severity (mild or severe). In addition, a group of exposed uninfected individuals was included. HLA-A*68 was found to be a protective allele against the severe infection and fatal outcome; pC = 0.03, OR = 0.4, 95% CI =0.20-0.86, and pC =0.009, OR = 0.3, 95% CI =0.13-0.71 respectively. HLA-DRB1*03 also appears to be a protective factor against fatal outcome pC = 0.009, OR = 0.1, 95%IC = 0.01-0.66; however, the low frequency of this allele in the studied population limits the statistical power. The severity and fatal outcome of COVID-19 patients in Tapachula, Chiapas depend more on the lack of resistance than susceptibility HLA alleles.

Heterogeneity of Genetic Admixture Determines SLE Susceptibility in Mexican

Systemic Lupus Erythematosus (SLE) is an autoimmune inflammatory disorder for which Major Histocompatibility Complex (MHC) genes are well identified as risk factors. SLE patients present different clinical phenotypes, which are partly explained by admixture patterns variation among Mexicans. Population genetic has insight into the high genetic variability of Mexicans, mainly described through HLA gene studies with anthropological and biomedical importance. A prospective, case-control study was performed. In this study, we recruited 146 SLE patients, and 234 healthy individuals were included as a control group; both groups were admixed Mexicans from Mexico City. The HLA typing methods were based on Next Generation Sequencing and Sequence-Based Typing (SBT). The data analysis was performed with population genetic programs and statistical packages. The admixture estimations based on HLA-B and -DRB1 revealed that SLE patients have a higher Southwestern European ancestry proportion (48 ± 8%) than healthy individuals (30 ± 7%). In contrast, Mexican Native American components are diminished in SLE patients (44 ± 1%) and augmented in Healthy individuals (63 ± 4%). HLA alleles and haplotypes' frequency analysis found variants previously described in SLE patients from Mexico City. Moreover, a conserved extended haplotype that confers risk to develop SLE was found, the HLA-A^∗29:02∼C^∗16:01∼B^∗44:03∼DRB1^∗07:01∼DQB1^∗02:02, pC = 0.02, OR = 1.41. Consistent with the admixture estimations, the origin of all risk alleles and haplotypes found in this study are European, while the protection alleles are Mexican Native American. The analysis of genetic distances supported that the SLE patient group is closer to the Southwestern European parental populace and farthest from Mexican Native Americans than healthy individuals. Heterogeneity of genetic admixture determines SLE susceptibility and protection in Mexicans. HLA sequencing is helpful to determine susceptibility alleles and haplotypes restricted to some populations.

Bioinformatics Prediction of SARS-CoV-2 Epitopes as Vaccine Candidates for the Colombian Population

Coronavirus disease (COVID-19) pandemic caused by the coronavirus SARS-CoV-2 represents an enormous challenge to global public health, with thousands of infections and deaths in over 200 countries worldwide. The purpose of this study was to identify SARS-CoV-2 epitopes with potential to interact in silico with the alleles of the human leukocyte antigen class I (HLA I) and class II (HLA II) commonly found in the Colombian population to promote both CD4 and CD8 immune responses against this virus. The generation and evaluation of the peptides in terms of HLA I and HLA II binding, immune response, toxicity and allergenicity were performed by using computer-aided tools, such as NetMHCpan 4.1, NetMHCIIpan 4.0, VaxiJem, ToxinPred and AllerTop. Furthermore, the interaction between the predicted epitopes with HLA I and HLA II proteins frequently found in the Colombian population was studied through molecular docking simulations in AutoDock Vina and interaction analysis in LigPlot+. One of the promising peptides proposed in this study is the HLA I epitope YQPYRVVVL, which displayed an estimated coverage of over 82% and 96% for the Colombian and worldwide population, respectively. These findings could be useful for the design of new epitope-vaccines that include Colombia among their population target.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.