Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations

Multiple measures for self-identification improve matching donors with patients in unrelated hematopoietic stem cell transplant

BACKGROUND

Questions persist around whether and how to use race or geographic ancestry in biomedical research and medicine, but these forms of self-identification serve as a critical tool to inform matching algorithms for human leukocyte antigen (HLA) of varying levels of resolution for unrelated hematopoietic stem cell transplant in large donor registries.

METHODS

Here, we examined multiple self-reported measures of race and ancestry from a survey of a cohort of over 100,000 U.S. volunteer bone marrow donors alongside their high-resolution HLA genotype data.

RESULTS

We find that these self-report measures are often non-overlapping, and that no single self-reported measure alone provides a better fit to HLA genetic ancestry than a combination including both race and geographic ancestry. We also found that patterns of reporting for race and ancestry appear to be influenced by participation in direct-to-consumer genetic ancestry testing.

CONCLUSIONS

While these data are not used directly in matching for transplant, our results demonstrate that there is a place for the language of both race and geographic ancestry in the critical process of facilitating accurate prediction of HLA in the donor registry context.

Variation within the non-coding genome influences genetic and epigenetic regulation of the human leukocyte antigen genes

Variation within the non-coding genome may influence the regulation and expression of important genes involved in immune control such as the human leukocyte antigen (HLA) system. Class I and Class II HLA molecules are essential for peptide presentation which is required for T lymphocyte activation. Single nucleotide polymorphisms within non-coding regions of HLA Class I and Class II genes may influence the expression of these genes by affecting the binding of transcription factors and chromatin modeling molecules. Furthermore, an interplay between genetic and epigenetic factors may also influence HLA expression. Epigenetic factors such as DNA methylation and non-coding RNA, regulate gene expression without changing the DNA sequence. However, genetic variation may promote or allow genes to escape regulation by epigenetic factors, resulting in altered expression. The HLA system is central to most diseases, therefore, understanding the role of genetics and epigenetics on HLA regulation will tremendously impact healthcare. The knowledge gained from these studies may lead to novel and cost-effective diagnostic approaches and therapeutic interventions. This review discusses the role of non-coding variants on HLA regulation. Furthermore, we discuss the interplay between genetic and epigenetic factors on the regulation of HLA by evaluating literature based on polymorphisms within DNA methylation and miRNA regulatory sites within class I and Class II HLA genes. We also provide insight into the importance of the HLA non-coding genome on disease, discuss ethnic-specific differences across the HLA region and provide guidelines for future HLA studies.

Impact of HLA-DQA1*05 Genotype in Immunogenicity and Failure to Treatment with Tumour Necrosis Factor-alpha Antagonists in Inflammatory Bowel Disease: A Systematic Review and Meta-analysis

BACKGROUND

HLA-DQA1*05 carriage has been associated with an increased risk of immunogenicity in patients with immune-mediated inflammatory diseases treated with tumour necrosis factor-alpha [TNF-a] antagonists. Results have shown an inconsistent association with a loss of response [LOR] in patients with inflammatory bowel disease [IBD], which could be modified when using proactive optimisation and association with immunomodulatory drugs.

AIMS

To define the association of HLA-DQA1*05 on anti-drug antibody development and loss of response [LOR] to anti-TNF-a in IBD.

METHODS

We searched MEDLINE, EMBASE, and SCOPUS, for the period up to August 2023, to identify studies reporting the risk of immunogenicity and/or LOR in IBD patients with HLA-DQA1*05 genotype.

RESULTS

A total of 24 studies comprising 12 papers, 11 abstracts and one research letter, with a total of 5727 IBD patients, were included. In a meta-analysis of 10 studies [2984 patients; 41.9% with HLA-DQA1*05 genotype], HLA-DQA1*05 carriers had higher risk of immunogenicity compared with non-carriers (risk ratio, 1.54; 95% confidence interval [CI], 1.23 - 1.94; I2 = 62%) [low certainty evidence]. Lack of therapeutic drug monitoring [TDM] increased immunogenicity in the presence of risk human leukocyte antigen [HLA] [risk ratio 1.97; 95% CI, 1.35 - 2.88; I2 = 66%], whereas proactive TDM revoked this association [very low certainty of evidence]. A meta-analysis of six studies [765 patients] found that risk for secondary LOR was higher among HLA-DQA1*05 carriers [hazard ratio 2.21; 95% CI, 1.69 - 2.88; I2 = 0%] [very low certainty evidence], although definition and time to assessment varied widely among studies.

CONCLUSION

HLA-DQA1*05 carriage may be associated with an increased risk of immunogenicity and secondary LOR in IBD patients treated with TNF-a antagonists.

Contribution of HLA class I (A, B, C) and HLA class II (DRB1, DQA1, DQB1) alleles and haplotypes in exploring ethnic origin of central Tunisians

BACKGROUND

Estimation of HLA (Human leukocyte Antigen) alleles' frequencies in populations is essential to explore their ethnic origin. Anthropologic studies of central Tunisian population were rarely reported. Then, in this work, we aimed to explore the origin of central Tunisian population using HLA alleles and haplotypes frequencies.

METHODS

HLA class I (A, B, C) and HLA class II (DRB1, DQA1, DQB1) loci genotyping of 272 healthy unrelated organ donors was performed by Polymerase Chain Reaction-Sequence Specific Oligonucleotide (PCR-SSO). We compared central Tunisians with other populations (Arabs, Berbers, Mediterraneans, Europeans, Africans, etc.) using alleles and haplotypes frequencies, genetic distances, Neighbour-Joining dendrogram and correspondence analysis.

RESULTS

Among the 19 HLA A alleles, the 26 HLA B alleles, the 13 HLA C alleles, the 15 HLA DRB1 alleles, the 6 HLA DQA1 alleles and the 5 HLA DQB1 alleles identified in the studied population, HLA A*02 (22.8%), HLA B*50 (13.1%), HLA C*06 (21.8%), HLA DRB1*07 (17.8%), HLA DQA1*01 (32.1%) and HLA DQB1*03 (31.6%) were the most frequent alleles. The extended haplotypes HLA A*02-B*50-C*06-DRB1*07-DQA1*02-DQB1*02 (1.97%) was the most frequent HLA six-loci haplotype.

CONCLUSION

Central Tunisians were very close to other Tunisian populations, to Iberians and North Africans. They were rather distant from sub-Saharan populations and eastern Mediterraneans especially Arabs although the strong cultural and religious impact of Arabs in this population.

Associative role of HLA-DRB1 as a protective factor for susceptibility and progression of Parkinson's disease: a Chinese cross-sectional and longitudinal study

Background

Previous genome-wide association studies investigating the relationship between the HLA-DRB1 and the risk of Parkinson's disease (PD) have shown limited racial diversity and have not explored clinical heterogeneity extensively.

Methods

The study consisted of three parts: a case-control study, a cross-sectional study, and a longitudinal cohort study. The case-control study included 477 PD patients and 477 healthy controls to explore the relationship between rs660895 and PD susceptibility. The cross-sectional study utilized baseline data from 429 PD patients to examine the correlation between rs660895 and PD features. The longitudinal study included 388 PD patients who completed a 3-year follow-up to investigate the effects of rs660895 on PD progression.

Results

In the case-control study, HLA-DRB1 rs660895-G allele was associated with a decreased risk of PD in allele model (adjusted OR=0.72, p = 0.003) and dominant model (AG + GG vs. AA: adjusted OR = 0.67, p = 0.003). In the cross-sectional analysis, there was no association between rs660895 and the onset age, motor phenotype, or initial motor symptoms. In the longitudinal analysis, PD patients with the G allele exhibited a slower progression of motor symptoms (MDS-UPDRS-III total score: β = -5.42, p < 0.001, interaction ptime × genotype < 0.001) and non-motor symptoms (NMSS score: β = -4.78, p = 0.030, interaction ptime × genotype < 0.001).

Conclusion

Our findings support HLA-DRB1 rs660895-G allele is a protective genetic factor for PD risk in Chinese population. Furthermore, we also provide new evidence for the protective effect of rs660895-G allele in PD progression.

HLA frequency distribution of the Portuguese bone marrow donor registry

Introduction

The Portuguese donor Registry of CEDACE was the fifth largest per capita bone marrow donor Registry of the WMDA as of 2019 and has yet to be thoroughly analyzed. We aimed to characterize its various aspects, including demographics and HLA allele and haplotype frequencies, to evaluate the genetic matching propensity score and ultimately further develop it.

Methods

We described and compared characteristics of the donor population with census data and used an Expectation-Maximization algorithm and analyses of molecular variance to assess haplotype frequencies and establish phylogenetic distances between regions and districts within the country.

Results

We identified 396545 donors, corresponding to 3.85% of the Portuguese population; the median donor age was 39 years, with 60.4% of female donors. Most donors were Portuguese nationals, although 40 other nationalities were present, with a significant proportion of donors from Brazil and Portuguese-speaking African Countries; almost all donors self-reported as Western, with the second largest group reporting African ancestry. There was an asymmetric contribution of donors from different districts and regions, with most coming from coastal districts and few from the southern districts and autonomous regions; foreign and self-declared non-Western donors were mainly located in the Metropolitan Area of Lisbon and the South. Although most donors were typed in three loci (HLA-A, HLA-B and HLA-DRB1), only 44% were also typed in HLA-C, 1.28% in HLA-DQB1 and only 0.77% in all five loci and in high-resolution. There were varying allele and haplotype frequencies across districts and regions, with the most common three loci, low-resolution haplotypes, being HLA-A*01~B*08~DRB1*03, A*29~B*44~DRB1*07 and HLA-A*02~B*44~DRB1*04; some haplotypes were more prevalent in the South, others in the North and a few in the autonomous regions; African and foreign donors presented relevant differences in haplotype frequency distributions, including rare haplotypes of potential interest. We also report on four loci, low-resolution frequency distributions. Using AMOVA, we compared genetic distances between districts and regions, which recapitulated the country's geography.

Discussion

Our analysis showed potential paths to optimization of the Registry, including increasing the male donor pool and focusing on underrepresented districts and particular populations of interest, such as donors from Portuguese-speaking African countries.

HLA-DQB1*06 and Select Neighboring HLA Variants Predict Chlamydia Reinfection Risk

Associations of HLA class II alleles with genital chlamydial infection outcomes have been reported, especially HLA DQB1*06. However, the potential role of DQB1*06 in influencing reinfection risk has still not been established. The purpose of this study was to determine whether the association of DQB1*06 with chlamydia reinfection was impacted by any other nearby HLA class II variants that were also associated with reinfection. We used next-generation sequencing to map HLA class II variants spanning the HLA-DQ and -DR loci. DQB1*06 as well as DQB1*04 were confirmed as significant predictors of chlamydia reinfection, when controlling for age and percent African ancestry. SKAT analysis revealed one region each in DRB1, DRB5, DQA2, and three intergenic regions that had variants associated with reinfection. Further analyses of these variants revealed that rs112651494 within DRB5 and an intergenic SNP rs617058 in DRB1:DQA1 were significantly associated with reinfection, but this did not impact the significance of the association of DQB1*06 or DQB1*04 with reinfection.

Chemotherapy-induced toxic epidermal necrolysis in a patient with multiple myeloma, a case report and literature review

Rationale and patient concerns

Toxic epidermal necrolysis (TEN) and Stevens-Johnson syndrome (SJS) are severe drug-induced skin reactions associated with a high mortality rate. The patient in this case report developed TEN after receiving the Velcade-lenalidomide-dexamethasone (VRD) regimen for the treatment of multiple myeloma (MM). The patient's concerns included the progression of the rash, pain, itching, and potential long-term complications. TEN is a life-threatening condition that requires prompt medical intervention and hospitalization.

Interventions

The treatment approach for the patient included discontinuation of the causative medication (lenalidomide) and comprehensive supportive therapy. Supportive measures included the administration of systemic corticosteroids (methylprednisolone), intravenous immunoglobulin infusion, pain relief medication (ebastine), antibiotic prophylaxis, laminar bed use, and regular dressing changes. The goal was to alleviate symptoms, promote skin and mucous membrane healing, and prevent complications such as infection.

Diagnosis

The patient was diagnosed with stage III A DS and stage III ISS MM, specifically of the immunoglobulin G (λ) type. Diagnostic procedures included CT and MRI scans, bone marrow testing through flow cytometry and morphology analysis, and laboratory tests to assess blood markers. The diagnosis of TEN was made based on the clinical presentation, skin biopsy, and exclusion of other potential causes.

Outcomes

With the implemented interventions, the patient's condition gradually improved, and the rash resolved without any residual scarring. The patient's skin and mucosa healed, blood markers improved, and bone pain was relieved. The patient was discharged within a month of receiving the final treatment with bortezomib and dexamethasone. The patient got partial response(PR) of multiple myeloma.

Lessons

Drug-induced SJS/TEN is more prevalent in Asian populations, potentially due to differences in human leukocyte antigen (HLA) alleles. The use of systemic corticosteroid therapy in SJS/TEN cases is controversial due to the potential risks of immune suppression and complications. Balancing the immune response to prevent SJS/TEN while maintaining an effective cytotoxic immune response for tumor control remains a challenge. Lenalidomide, an immunomodulatory agent, can enhance antitumor immune responses but also contribute to the pathogenesis of SJS/TEN. Increased awareness of HLA variations and frequently mutated genes in different malignancies can help prevent SJS/TEN and improve patient outcomes.

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.

Congenital Adrenal Hyperplasia and Human Leukocyte Antigen B: A Meta-Analysis

The link between specific human leukocyte antigen (HLA)-B genes and congenital adrenal hyperplasia (CAH) has been a subject of interest. This study investigates the association between specific HLA-B haplotypes and CAH through a meta-analysis. Google Scholar was used as a database. Articles were included if the research was conducted between 1970 and 2022, was not a meta-analysis, and had odds ratios or enough data points to calculate an odds ratio. The National Institutes of Health (NIH) quality assessment tool of case-control studies was used to evaluate the risk of bias in individual studies, and MetaXL was used to generate data and create a forest plot for analysis. Twelve studies met the selection criteria and were included in the study (641 patients and 3,614 controls). Two HLA-B haplotypes showed increased odds of CAH compared to controls: B14 (OR=3.81; 95%CI=2.88, 5.05; I2=3%) and B35 (OR=1.88; 95%CI=1.22, 2.90; I2=25%). All other HLAs either showed no significant effect or had high heterogeneity. The results suggest that specific HLA-B haplotypes have increased odds of developing CAH, specifically B14 and B35. These findings may prove helpful in the pre- and post-natal diagnosis of CAH as well as the identification of carriers and prediction of patient prognosis.

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.

The HLA profile and genetic affinities of three primitive Tamil-speaking endogamous groups: Kallars of Thanjavur, Piramalai Kallar and Vanniyar

Background

The present study was aimed to study the frequencies of HLA-DRB1/-DQB1 alleles and haplotypes of three endogamous groups of Tamil Nadu state, South India. PCR-SSP typing of HLA-DRB1 and -DQB1 alleles were performed on 111 Kallars of Thanjavur, 80 Piramalai Kallar of Madurai and 119 Vanniyar. Genetic distances, neighbor-joining phylogenetic dendrograms and correspondence analysis have been performed.

Results

The HLA class II alleles, DRB1*07 (25.2%), DRB1*15 (15.7%), DRB1*14 (11.7%) and DRB1*12 (9.90%) among Kallars of Thanjavur; DRB1*15 (28.7%), DRB1*04 (15.6%), DRB1*10 (14.3%), DRB1*13 (11.2%) and DRB1*03 (9.37%) among Piamalai Kallar and DRB1*15 (24.7%), DRB1*04 (15.9%), DRB1*07 (11.7%), DRB1*12 (11.3%) and DRB1*10 (10.0%) among Vanniyar were more frequent. Similarly, alleles DQB1*06 (31.0%), DQB1*02 (26.5%) and DQB1*05 (24.7%) among Kallars of Thanjavur; DQB1*05 (32.5%), DQB1*06 (31.8%), DQB1*02 (16.2%) and DQB1*03:02 (12.5%) among Piramalai Kallar and DQB1*05 (52.9%), DQB1*06 (22.6%) and DQB1*02 (11.3%) among Vanniyar were more frequent. We genotyped the two most frequent two-locus haplotypes, such as DRB1*15-DQB1*06 and DRB1*07-DQB1*02 for HLA-A/-B/–C alleles to identify the 5-locus extended haplotypes to extrapolate global affinities. We identified a number of five locus extended haplotypes among south Indian population with stronger global affinities. Further, we identified the presence of a highly unique extended haplotypes such as A*11-B*35-C*12-DRB1*07-DQB1*02 (HF:0.1458) in Kallars of Thanjavur, A*03-B*35-C*04-DRB1*15-DQB1*06 (HF:0.1833) in Piramalai Kallar and A*03-B*07-C*07-DRB1*15-DQB1*06 (HF: 0.1800) in Kallars of Thanjavur and (HF: 0.1081) in Vanniyar population.

Conclusions

Allele distribution and haplotype analysis have demonstrated that the Kallars of Thanjavur, Piramalai Kallar and Vanniyar populations shared HLA alleles with other ethnic and other Indian populations, while showing population specific haplotypes. Analysis of population-specific distribution of HLA alleles is proved to be important in finding out the relatedness of the ethnic groups across continents. The extensive polymorphism of the HLA system also has useful application in the study of the origin, evolution and migration patterns of human populations.

A comprehensive human leukocyte antigen analysis of 36 782 cord blood units stored in the Australian Public Cord Blood Banking Network

BACKGROUND AND AIMS

The network of public cord blood banks (CBBs) in Australia, known as AusCord, comprises CBBs located in Brisbane, Sydney and Melbourne. A novel comprehensive analysis has been performed to determine whether the cryopreserved, searchable cord blood unit (CBU) inventory of approximately 36 000 units share similar tissue types or haplotypes.

METHODS

Human leukocyte antigen (HLA) data was analysed using Microsoft Excel following standardisation of typing data.

RESULTS

The analysis has found that the majority of stored, searched and released CBU exhibit a tissue type that is unique within and between the CBBs. Therefore, each collection performed by the CBBs is likely to comprise a tissue type that is not already stored among the total AusCord inventory. HLA alleles (HLA-A*34, HLA-B*56, HLA-DRB1*08:03), which are uncommon in European populations, were associated with Pacific Islander and/or Indigenous Australian populations and confirmed to be more frequent among donors who, when screened, self-identified as these ethnicities.

CONCLUSIONS

These data indicate that (i) continued addition of CBU to existing inventories is likely to further increase the HLA diversity and (ii) screening donors for ethnicity or strategically locating collection sites where ethnic minorities reside can successfully result in collection of rare HLA associated with ethnic minority groups for whom finding donors might otherwise be more difficult.

Distribution of HLA-A, -C, -B, -DRB1, and -DQB1 polymorphisms in the Korean minority in Kazakhstan, and relatedness to neighboring and distant populations

BACKGROUND

Among the 125 ethnicities and linguistic groups in Kazakhstan, Koreans represent the eighth minority numerically. HLA class I and class II allele frequencies and extended haplotypes were studied for the first time and were compared to related and distant populations worldwide.

METHODS

HLA class I (A, B, C) and class II (DRB1, DQB1) profile was investigated in Kazakhstani Koreans, and were compared to other populations using standard genetic distances (SGD), neighbor-joining dendrograms, correspondence, and haplotype analysis.

RESULTS

One hundred and thirty-one HLA alleles were identified in Koreans living in Kazakhstan, with A*02:01 (23.08%), B*35:01 (8.24%), C*01:02 (15.38%), DRB1*08:03 (9.89%), and DQB1*03:01 (21.98%) being the most frequent alleles. A*03:01 ∼ B*07:02 (3.85%), B*08:01 ∼ DRB1*03:01 (3.85%), B*07:02 ∼ C*07:02 (7.14%), and DRB1*08:03 ∼ DQB1*06:01 (9.34%) were the most frequent two-locus haplotypes, while A*02:01 ∼ B*18:01 ∼ C*07:01 ∼ DRB1*11:04 ∼ DQB1*03:01 and A*33:03 ∼ B*44:03 ∼ C*14:03 ∼ DRB1*13:02 ∼ DQB1*06:04 (2.2% each) were the most frequent five-locus haplotypes.

CONCLUSION

Korean minority in Kazakhstan was closely related to East Asians, including Mongolians (SGD, 0.044), Tuvans (East Siberia; SGD, 0.081), Burayts (Siberia; SGD, 0.094), but distant from East Mediterranean such as Lebanese (SGD, 0.367), Greek (SGD, 0.377), and Saudi (SGD, 0.414), and most Siberians (SGD, 0.473-0.699). This relatedness could be mainly attributed to massive migration of Koryo Saram to Kazakhstan in the 20th century.

How HLA diversity is apportioned: influence of selection and relevance to transplantation

In his 1972 paper ‘The apportionment of human diversity’, Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low F_ST values. However, when we analyse haplotypes defined by these SNPs (which define ‘HLA alleles’), we find marked differences in frequencies between geographic regions. These differences are not reflected in the F_ST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting F_ST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals.

Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA).

Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies.

Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

Global Meta-Analysis on the Association between Behcet Syndrome and Polymorphisms from the HLA Class I (A, B, and C) and Class II (DRB1, DQB1, and DPB1) Genes

Behcet syndrome (BS) is a multisystemic perivasculitis whose genetic susceptibility is linked to HLA region. We first meta-analysed all HLA class I and II genes involved in BS susceptibility in all ethnic groups worldwide. We identified 1141 articles and finally included 31 case-control studies after multiple rounds of selection. We analysed frequencies for 24 HLA-A alleles (3 alleles for HLA-A∗26 at four digits), 50 HLA-B alleles (11 alleles for HLA-B∗51 at four digits), 15 HLA-C alleles, 16 HLA-DRB1 alleles, 6 HLA-DQB1 alleles, and 15 HLA-DPB1 alleles. We meta-analysed only HLA allelic frequencies from at least three studies; therefore, we investigated 21 alleles out of 140. Going from 7.00 to 1.6 OR, we found 11 class I alleles conferring risk for BS: B∗51 : 08, B∗51, B∗51 : 01, B∗51 : 02, DQB1∗03, A∗26 : 01, Cw∗14, Cw∗15, Cw∗16, B∗15, and A∗26. Overall, the studies included populations from Europe (Greece, Spain, Italy, Germany, and Ireland), Asia (Korea, China, China Han, and Thailand), Middle East (Israel, Saudi Arabia, and Iran), and Morocco (as no other North-African population was included). We collected a number of ethnical groups sufficient to conduct an ethnic-specific meta-analysis where Europeans showed 11.25 OR for B∗51:08 and Japan 3.50 OR for A∗26 : 01. A remarkable result was that the most frequent HLA - B∗51 two-digit alleles associated with BS were different among populations: HLA - B∗51 : 08 in Europe, HLA - B∗51 : 01 in Turkey, and HLA - B∗51 : 02 in Japan. Overall, we discussed our real-world results with other imputation studies.

Analysis of HLA gene polymorphisms in East Africans reveals evidence of gene flow in two Semitic populations from Sudan

Sudan, a northeastern African country, is characterized by high levels of cultural, linguistic, and genetic diversity, which is believed to be affected by continuous migration from neighboring countries. Consistent with such demographic effect, genome-wide SNP data revealed a shared ancestral component among Sudanese Afro-Asiatic speaking groups and non-African populations, mainly from West Asia. Although this component is shared among all Afro-Asiatic speaking groups, the extent of this sharing in Semitic groups, such as Sudanese Arab, is still unknown. Using genotypes of six polymorphic human leukocyte antigen (HLA) genes (i.e., HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1), we examined the genetic structure of eight East African ethnic groups with origins in Sudan, South Sudan, and Ethiopia. We identified informative HLA alleles using principal component analysis, which revealed that the two Semitic groups (Gaalien and Shokrya) constituted a distinct cluster from the other Afro-Asiatic speaking groups in this study. The HLA alleles that distinguished Semitic Arabs co-exist in the same extended HLA haplotype, and those alleles are in strong linkage disequilibrium. Interestingly, we find the four-locus haplotype "C*12:02-B*52:01-DRB1*15:02-DQB1*06:01" exclusively in non-African populations and it is widely spread across Asia. The identification of this haplotype suggests a gene flow from Asia, and likely these haplotypes were brought to Africa through back migration from the Near East. These findings will be of interest to biomedical and anthropological studies that examine the demographic history of northeast Africa.

An Immunoinformatics Approach for SARS-CoV-2 in Latam Populations and Multi-Epitope Vaccine Candidate Directed towards the World's Population

The coronavirus pandemic is a major public health crisis affecting global health systems with dire socioeconomic consequences, especially in vulnerable regions such as Latin America (LATAM). There is an urgent need for a vaccine to help control contagion, reduce mortality and alleviate social costs. In this study, we propose a rational multi-epitope candidate vaccine against SARS-CoV-2. Using bioinformatics, we constructed a library of potential vaccine peptides, based on the affinity of the most common major human histocompatibility complex (HLA) I and II molecules in the LATAM population to predict immunological complexes among antigenic, non-toxic and non-allergenic peptides extracted from the conserved regions of 92 proteomes. Although HLA-C, had the greatest antigenic peptide capacity from SARS-CoV-2, HLA-B and HLA-A, could be more relevant based on COVID-19 risk of infection in LATAM countries. We also used three-dimensional structures of SARS-CoV-2 proteins to identify potential regions for antibody production. The best HLA-I and II predictions (with increased coverage in common alleles and regions evoking B lymphocyte responses) were grouped into an optimized final multi-epitope construct containing the adjuvants Beta defensin-3, TpD, and PADRE, which are recognized for invoking a safe and specific immune response. Finally, we used Molecular Dynamics to identify the multi-epitope construct which may be a stable target for TLR-4/MD-2. This would prove to be safe and provide the physicochemical requirements for conducting experimental tests around the world.

The HLA diversity of the Anthony Nolan register

While the success of allogeneic stem cell transplantation depends on a high degree of HLA compatibility between donor and patient, finding a suitable donor remains challenging due to the hyperpolymorphic nature of HLA genes. We calculated high-resolution allele, haplotype and phenotype frequencies for HLA-A, -C, -B, -DRB1 and -DQB1 for 10 subpopulations of the Anthony Nolan (AN) register using an in-house expectation-maximisation (EM) algorithm run on mixed resolution HLA data, covering 676 155 individuals. Sample sizes range from 599 410 for British/Irish North West European (BINWE) individuals, the largest subpopulation in the United Kingdom to 1105 for the British Bangladeshi population. Calculation of genetic distance between the subpopulations based on haplotype frequencies shows three broad clusters, each following a major continental group: European, African and Asian. We further analysed the HLA haplotype and phenotype diversity of each subpopulation, and found that 35.52% of BINWE individuals ranging to 98.34% of Middle Eastern individuals on the register had a unique phenotype within their subpopulation. These analyses and the allele, haplotype and phenotype frequency data of the subpopulation on the AN register are a valuable resource in understanding the HLA diversity in the United Kingdom and can be used to improve the accuracy of match likelihoods and to inform future donor recruitment strategies.

How Ancestry Influences the Chances of Finding Unrelated Donors: An Investigation in Admixed Brazilians

A match of HLA loci between patients and donors is critical for successful hematopoietic stem cell transplantation. However, the extreme polymorphism of HLA loci - an outcome of millions of years of natural selection - reduces the chances that two individuals will carry identical combinations of multilocus HLA genotypes. Further, HLA variability is not homogeneously distributed throughout the world: African populations on average have greater variability than non-Africans, reducing the chances that two unrelated African individuals are HLA identical. Here, we explore how self-identification (often equated with "ethnicity" or "race") and genetic ancestry are related to the chances of finding HLA compatible donors in a large sample from Brazil, a highly admixed country. We query REDOME, Brazil's Bone Marrow Registry, and investigate how different criteria for identifying ancestry influence the chances of finding a match. We find that individuals who self-identify as "Black" and "Mixed" on average have lower chances of finding matches than those who self-identify as "White" (up to 57% reduction). We next show that an individual's African genetic ancestry, estimated using molecular markers and quantified as the proportion of an individual's genome that traces its ancestry to Africa, is strongly associated with reduced chances of finding a match (up to 60% reduction). Finally, we document that the strongest reduction in chances of finding a match is associated with having an MHC region of exclusively African ancestry (up to 75% reduction). We apply our findings to a specific condition, for which there is a clinical indication for transplantation: sickle-cell disease. We show that the increased African ancestry in patients with this disease leads to reduced chances of finding a match, when compared to the remainder of the sample, without the condition. Our results underscore the influence of ancestry on chances of finding compatible HLA matches, and indicate that efforts guided to increasing the African component of registries are necessary.

Demographic history and selection at HLA loci in Native Americans

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

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.

High-resolution characterization of 12 classical and non-classical HLA loci in Southern Brazilians

The human leukocyte antigen (HLA) are the most polymorphic genes in the human genome. Because of their importance for antigen recognition, HLA molecules play a central role in host defense and graft rejection upon transplantation. The aim of this study was to characterize allelic diversity of the classical HLA genes HLA-A, -B, -C, -DRA, -DRB1, -DQA1, -DQB1, -DPA1, -DPB1, and the non-classical class I genes HLA-E, -F and -G at high-resolution for a population of predominantly European ancestry from Curitiba, Brazil. Genotyping of 108 individuals was performed by next-generation sequencing on the MiSeq platform and also by Sanger sequencing. The genotype distributions of all loci were in accordance with Hardy-Weinberg equilibrium (P > 0.05) and a total of 202 HLA variants at second field resolution were observed for the 12 loci. The strongest linkage disequilibrium (r² = 1.0, P < 10^-5 ) was observed for the following pairs of alleles: HLA-B*42:01:01 ~ HLA-DRB1*03:02:01; HLA-B*14:02:01 ~ HLA-C*08:02:01; B*42:01:01 ~ HLA-C*17:01:01; HLA-DRB1*03:01:01 ~ HLA-DQB1*02:01:01 ~ DRB1*03:01:01 ~ HLA-DQB1*02:01:01; DRB1*13:01:01~ HLA-DQB1*06:03:01 and HLA-DRB1*09:01:02 ~ HLA-DQA1*03:02. This is the first study to characterize all 12 HLA genes at high resolution in a single population. On the basis of the allelic frequencies of worldwide populations and principal component analysis, we confirmed the similarity of the study population to European and other Euro-descendant populations.

Next generation sequencing of 11 HLA loci characterises a diverse UK cord blood bank

The introduction of next generation sequencing (NGS) for stem cell donor registry typing has contributed to faster identification of compatible stem cell donors. However, the successful search for a matched unrelated donor for some patient groups is still affected by their ethnicity. In this study, DNA samples from 714 National Health Service (NHS) Cord Blood Bank donors were typed for HLA-A, -B, -C, -DRB1, -DRB345, -DQA1, -DQB1, -DPA1 and -DPB1 by NGS. Analysis of the ethnic diversity showed a high level of diversity, with the cohort comprising of 62.3% European and 37.7% of either multi-ethnic or non-European donors, of which 12.3% were multi-ethnic. The HLA diversity was further confirmed using PyPop analysis, 405 distinct alleles were observed in the overall NHS-CBB cohort, of which 37 alleles are non-CWD, including A*31:14N, B*35:68:02, C*14:23 and DQA1*05:10. Furthermore, HLA-DQA1 and HLA-DPA1 analysis showed 12% and 10%, respectively, of the alleles currently submitted to IMGT, confirming further diversity of the NHS-CBB cohort. The application of 11 HLA loci resolution by NGS revealed a high level of diversity in the NHS-CBB cohort. The incorporation of this data coupled with ethnicity data could lead to improved donor selection, contributing to better clinical outcomes for patients.

Identification, molecular characterization and expression of aminopeptidase N-1 (APN-1) from Anopheles stephensi in SF9 cell line as a candidate molecule for developing a vaccine that interrupt malaria transmission

BACKGROUND

According to the World Health Organization reports, billions of people around the world are at risk for malaria disease and it is important to consider the preventive strategies for protecting the people that are living in high risk areas. One of the main reasons of disease survival is diversity of vectors and parasites in different malaria regions that have their specific features, behaviour and biology. Therefore, specific regional strategies are necessary for successful control of malaria. One of the tools that needs to be developed for elimination and prevention of reintroduction of malaria is a vaccine that interrupt malaria transmission (VIMTs). VIMT is a broad concept that should be adjusted to the biological characteristics of the disease in each region. One type of VIMT is a vector-based vaccine that affects the sexual stage of Plasmodium life cycle. According to recent studies, the aminopeptidase N-1 of Anopheles gambiae (AgAPN-1) is as a potent vector-based VIMT with considerable inhibition activity against the sexual stage of Plasmodium parasite.

METHODS

Systems for rapid amplification of cDNA ends (3'-RACE) and genome walking methods were used for sequence determination of apn-1 gene from Anopheles stephensi and distinct bioinformatics software were used for structural analysis. AsAPN-1 was expressed in Spodoptera frugiperda (Sf9) insect cell line using the baculovirus expression system. Recombinant AsAPN-1 was purified under the hybrid condition and its biological activity was assayed.

RESULTS

Asapn-1 gene and its coded protein from An. stephensi were characterized for the first time in this study. Subsequently, the structural features and immunological properties of its coded protein were evaluated by in silico approaches. Enzymatic activity of the recombinant AsAPN-1, which was expressed in Sf9 insect cell line, was equal to 6 unit/μl.

CONCLUSIONS

Results of this study revealed that AsAPN-1 is very similar to its counterpart in An. gambiae. In silico evaluation and fundamental data which are necessary for its evaluation as a VIMT-based vaccine in the next steps were acquired in this study and those could be useful for research groups that study on malaria vaccine for countries that An. stephensi is the main malaria vector there.

Review article: the genetics of the human leucocyte antigen region in inflammatory bowel disease

BACKGROUND

The human leucocyte antigen (HLA) complex, located at chromosome 6p21.3 is a highly polymorphic region containing the classical class I and II HLA genes. The region is highly associated with inflammatory bowel disease (IBD), largely through genome-wide association studies (GWAS).

AIMS

To review the role of HLA in immune function, summarise data on risk/protective HLA genotypes for IBD, discuss the role of HLA in IBD pathogenesis, treatment and examine limitations that might be addressed by future research.

METHODS

An organised search strategy was used to collate articles describing HLA genes in IBD, including Crohn's disease and ulcerative colitis.

RESULTS

All classical HLA genes with variation (including HLA-A, B, C, DRB1, DQA1, DQB1, DPA1 and DPB1) harbour IBD-associated genotypes. The most implicated gene is HLA-DRB1, with HLA-DRB1*03:01 the most associated risk allele in both Crohn's disease and ulcerative colitis. Elucidating precise disease associations is challenging due to high linkage disequilibrium between HLA genotypes. The mechanisms by which risk alleles cause disease are multifactorial, with the best evidence indicating structural and electrostatic alteration impacting antigen binding and downstream signalling. Adverse medication events have been associated with HLA genotypes including with thiopurines (pancreatitis) and anti-TNF agents (antibody formation).

CONCLUSIONS

The HLA complex is associated with multiple risk/protective alleles for IBD. Future research utilising long-read technology, ascertainment of zygosity and integration in disease modelling will improve the functional understanding and clinical translation of genetic findings.

Using next-generation sequencing for characterising HLA-DRB1 and DQB1 loci in a cohort of Colombian women

The Colombian population is characterised by a high genetic diversity, secondary to the ethnic mixture arising from colonisation. Unfortunately, few reports are available regarding HLA-DRB1 and DQB1 diversity in Colombia to date. HLA-DRB1 and DQB1 diversity was identified in this study using next-generating sequencing (NGS) on a cohort of Colombian women. Cervical samples taken from 276 women were used for typing DRB1 and DQB1 loci by Illumina MiSeq. Allele and haplotype frequencies were calculated using an expectation-maximisation algorithm. Hardy-Weinberg Equilibrium and linkage disequilibrium (LD) between loci were evaluated. Forty-seven DRB1 alleles and 14 DQB1 alleles were identified. DRB1*04:07:01G and DQB1*03:02:01G alleles occurred most frequently in the target population. Significant LD was found in 44 out of the 144 identified haplotypes, within which DRB1*04:07:01G-DQB1*03:02:01G occurred most frequently (6.56%). The alleles and haplotypes found with NGS agreed with that found in previous reports involving lower resolution for the Colombian population, and greater genetic variability was found, especially concerning DRB1. Comparing allele and haplotype frequency distribution in the target population to that of other populations denoted HLA system intra- and inter-population diversity.

High-resolution HLA allele and haplotype frequencies in majority and minority populations of Costa Rica and Nicaragua: Differential admixture proportions in neighboring countries

The HLA system shows the most extensive polymorphism in the human genome. Allelic and haplotypic frequencies of HLA genes vary dramatically across human populations. Due to a complex history of migration, populations in Latin America show a broad variety of admixture proportions, usually varying not only between countries, but also within countries. Knowledge of HLA allele and haplotype frequencies is essential for medical fields such as transplantation, but also serves as a means to assess genetic diversity and ancestry in human populations. Here, we have determined high-resolution HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a sample of 713 healthy subjects from three Mestizo populations, one population of African descent, and Amerindians of five different groups from Costa Rica and Nicaragua and compared their profiles to a large set of indigenous populations from Iberia, Sub-Saharan Africa, and the Americas. Our results show a great degree of allelic and haplotypic diversity within and across these populations, with most extended haplotypes being private. Mestizo populations show alleles and haplotypes of putative European, Amerindian, and Sub-Saharan African origin, albeit with differential proportions. Despite some degree of gene flow, Amerindians and Afro-descendants show great similarity to other Amerindian and West African populations, respectively. This is the first comprehensive study reporting high-resolution HLA diversity in Central America, and its results will shed light into the genetic history of this region while also supporting the development of medical programs for organ and stem cell transplantation.

HLA Class II Allele and Haplotype Diversity in Libyans and Their Genetic Relationships with Other Populations

Background: Libya witnessed the succession of many civilizations and ethnic groups throughout history, thereby questioning the origin of present-day Libyans. Indeed, they were considered Africans given the geographical position of the country, Arabs at the cultural level, and Berbers because of the notable presence of Berber tribes. Genetic anthropology studies investigating the origin of Libyans were rarely reported, and thus little was known about the population structure of current Libyans, particularly at autosomic markers level. Methods: We examined HLA class II (DRB1, DQB1) gene profiles of 101 unrelated Libyans, and compared them with Arab-speaking communities and with Sub-Saharan and Mediterranean populations using Neighbour-Joining dendrograms, genetic distances, correspondence, and haplotype analysis. Results: Of the 42 DRB1 alleles identified, DRB1*07:01 (14.36%), DRB1*03:01 (12.38%) were the most frequent, while DQB1*02:01 (24.17%), DQB1*02:02 (13.86%), and DQB1*03:01 (12.38%) were the most frequent of the 17 DQB1 alleles detected. DRB1*03:01-DQB1*02:01 (6.93%), DRB1*07:01-DQB1*02:02 (4.45%), and DRB1*04:03-DQB1*03:02 (3.46%) were the most frequent DRB1-DQB1 haplotypes. Conclusion: Libyans appear to be closely related to North Africans, Saudis, and Iberians, but distinct from Levantine Arabs, East Mediterraneans, and Sub-Saharan Africans. This indicates limited genetic contribution of Levantine Arabs and Sub-Saharans on the makeup of Libyan gene pool. Our study confirmed genetic heterogeneity among Arab populations, with three identified groups. The first comprises North Africans, Saudis, and Kuwaitis who were related to Iberians and West Mediterraneans, while the second consists of Levantine Arabs who were close to East Mediterraneans, and the third contained Sudanese and Comorians, with a close relatedness to Sub-Saharans.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

HLA class I allele lineages and haplotype frequencies in Arabs of the United Arab Emirates

The high degree of polymorphism of the HLA system provides suitable genetic markers to study the diversity and migration of different world populations and is beneficial for forensic identification, anthropology, transplantation and disease associations. Although the United Arab Emirates (UAE) population of about nine million people is heterogeneous, information is limited for the HLA class I allele and haplotype frequencies of the Bedouin ethnic group. We performed low-resolution PCR-SSP genotyping of three HLA class I loci at HLA-A, -B and -C for 95 unrelated healthy Bedouins from the cities of Al Ain and Abu Dhabi in the UAE. A total of 54 HLA allele lineages were detected; the most frequent low-resolution allele lineages at each HLA locus were A*02 (0.268), B*51 (0.163) and C*07 (0.216). The inferred estimates for the two most frequent HLA-A and HLA-B haplotypes were HLA-A*02 ~ HLA-B*50 (0.070) and HLA-A*02 ~ HLA-B*51 (0.051), and the most frequent 3-locus haplotype was HLA-A*02 ~ HLA-B*50 ~ HLA-C*06 (0.068). The HLA allele lineage frequencies of the UAE Arabs were compared to those previously reported for 70 other world populations, and a strong genetic similarity was detected between the UAE Arabs and the Saudi Arabians from the west with evidence of a limited gene flow between the UAE Arabs and Pakistani across the Gulf from the east, and the UAE Arabs and Omani from the south of the Gulf Peninsula.

High-resolution characterization of allelic and haplotypic HLA frequency distribution in a Spanish population using high-throughput next-generation sequencing

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS.

A meta-analysis of association of Human Leukocyte Antigens A, B, C, DR and DQ with Human Papillomavirus 16 infection

Human Papillomavirus (HPV) induced cervical cancer (CaCx) is a major health problem in women from both developing and developed regions of the world. This virus accounts for >95% of the CaCx cases with a preponderance of HPV type -16 (65%). Paradoxically HPV-16 is prevalent even in the cervix of healthier women and anti HPV-16 T-cell response is considered critical for the viral clearance. Studies on HLA association with HPV-16 infection and cervical cancer have yielded varied HLA associations in different epidemiological settings. To validate these associations, we performed a meta-analysis of HLA-A, B, C, DR and DQ association with HPV-16 infection. Of the 1409 studies retrieved, 26 qualified for meta-analysis based on stringent inclusion and exclusion criteria. HLA-B*47, B*57, DRB1*10, DRB1*15 and DQB1*0303 were significantly associated with HPV-16 infection (OR = 3.4, 1.8, 1.5, 1.1 and 1.5 respectively). HLA-B*49, B*39, A28 (serotype), C*04 and DRB1*13 were negatively associated with HPV-16 (OR = 0.5, 0.6, 0.7, 0.7, and 0.7 respectively). Certain HLA alleles such as B*07, DRB1*15, DRB1*11 and DRB1*07 showed weakly positive associations. A comprehensive analysis coupling HPV-16 antigenic diversity and the HLA variation in various global populations shall provide further insights into the immunogenetic predisposition to HPV-16 and shall help identify host-parasite co-evolution.

Genetics and genomic medicine in Tunisia

Genetics and genomic medicine in Tunisia.

Genetic polymorphism analysis of MICB gene in Jing ethnic minority of Southern China

In the present study, the polymorphism in the 5'-upstream regulation region (5'-URR), coding region (exons 2-4), and the 3'-untranslated region (3'-UTR) of MICB gene were investigated for 150 healthy unrelated Jing individuals in Guangxi Zhuang Autonomous Region, by using PCR-SBT method. A total of 14 variation sites in the 5'-URR, 9 in coding region, and 6 in the 3'-UTR were detected in the Jing population. The MICB gene seems to present two different lineages showing functional variations mainly in nucleotides of the promoter region. Nineteen different MICB extended haplotypes (EHs) encompassing the 5'-URR, exons 2-4, and 3'-UTR were found in this population, and the most frequent was EH2 (20.33%). The findings here are of importance for future studies on the potential role of regulation region of MICB gene in disease association, transplantation, viral infection, and tumor progression among Jing population.

Benchmarking the HLA typing performance of Polysolver and Optitype in 50 Danish parental trios

BACKGROUND

The adaptive immune response intrinsically depends on hypervariable human leukocyte antigen (HLA) genes. Concomitantly, correct HLA phenotyping is crucial for successful donor-patient matching in organ transplantation. The cost and technical limitations of current laboratory techniques, together with advances in next-generation sequencing (NGS) methodologies, have increased the need for precise computational typing methods.

RESULTS

We tested two widespread HLA typing methods using high quality full genome sequencing data from 150 individuals in 50 family trios from the Genome Denmark project. First, we computed descendant accuracies assessing the agreement in the inheritance of alleles from parents to offspring. Second, we compared the locus-specific homozygosity rates as well as the allele frequencies; and we compared those to the observed values in related populations. We provide guidelines for testing the accuracy of HLA typing methods by comparing family information, which is independent of the availability of curated alleles.

CONCLUSIONS

Although current computational methods for HLA typing generally provide satisfactory results, our benchmark - using data with ultra-high sequencing depth - demonstrates the incompleteness of current reference databases, and highlights the importance of providing genomic databases addressing current sequencing standards, a problem yet to be resolved before benefiting fully from personalised medicine approaches HLA phenotyping is essential.

HLA class I (-A, -B, -C) and class II (-DR, -DQ) polymorphism in the Mauritanian population

BACKGROUND

HLA antigens have been widely studied for their role in transplantation biology, human diseases and population diversity. The aim of this study was to provide the first profile of HLA class I and class II alleles in the Mauritanian population.

METHODS

HLA typing was carried in 93 healthy Mauritanian blood donors, using single specific primer amplification (PCR-SSP).

RESULTS

Occurrences of the main HLA class I (-A, -B, -C) and class II (-DR, -DQ) antigens in the general population showed that out of the 17 HLA-A allele groups detected, five main HLA-A allele groups: A*02 (18.42%), A*01 (14.04%), A*23 (14.04%), A*30 (13.16%) and A*29 (12.28%) were the most common identified along other 12 relatively minor allele groups. Twenty three allele groups were observed in the locus B of which B*07 (13.46%) was the most prevalent followed by B*15, B*35, B*08 and B*27 all, with a frequency between 7 to 8%. Three prevalent HLA-C allele groups (C*02: 35.09%, C*07: 20.19% and C*06: 13.6%) were detected. The main HLA class II observed allele groups were: DRB1*13 (27.42%), DRB1*03 (24.73%), DRB1*11 (13.98%), DQB1*03 (36.03%), DQB1*02 (22.06%) and DQB1*05 (18.8%). Except for few haplotype in class I (A*02-B*07: 4.45%, A*02-C02: 10%, A*23-C*02: 8.8%, B*07-C*02: 8.8%, B*15-C*02: 8.8%) and in class II (DRB1*13-DQB1*06: 11.94%, DRB1*03-DQB1*02:11.19% and DRB1*03-DQB1*03: 10.45%), the majority of locus combination were in the range of 2-3%. A single predominant haplotype C*02-DRB1*03 (16.67%) was found.

CONCLUSIONS

These results, in agreement with previous data using different tissues markers, underlined the ethnic heterogeneity of the Mauritanian population.

Human perforin gene variation is geographically distributed

Background

Deleterious mutations in PRF1 result in lethal, childhood disease, familial hemophagocytic lymphohistiocytosis type 2 (FHL 2). However, not all mutations in PRF1 are deleterious and result in FHL 2. Currently, these nondeleterious mutations are being investigated in the onset of numerous disorders, such as lymphomas and diabetes. Yet, there is still an overwhelmingly large amount of PRF1 mutations that are not associated with disease.

Methods

We conducted a post hoc analysis of the PRF1 mutations in the coding region using the recently published Exome Aggregation Consortium genomes, Leiden Open Variation Database, NCBI SNP database, and primary literature to better understand PRF1 variation in the human population.

Results

This study catalogs 460 PRF1 mutations in the coding region, and demonstrates PRF1 is more variant then previously predicted. We identify key PRF1 mutations with high allelic frequency and are only found in certain populations. Additionally, we define PRF1SNVs are geographically distributed.

Conclusions

This study concludes with a novel hypothesis that nondeleterious mutation in PRF1, which decreases perforin expression and/or activity, may be an example of selective advantage in the context of environmental stressors prevalent near the equator. Our studies illustrate how perforin deficiency can be protective from injuries resulting in blood–brain barrier (BBB) disruption.

Assessment of Spanish Panel Reactive Antibody Calculator and Potential Usefulness

Background and objectives

The calculated panel reactive of antibodies (cPRAs) necessary for kidney donor-pair exchange and highly sensitized programs are estimated using different panel reactive antibody (PRA) calculators based on big enough samples in Eurotransplant (EUTR), United Network for Organ Sharing (UNOS), and Canadian Transplant Registry (CTR) websites. However, those calculators can vary depending on the ethnic they are applied. Here, we develop a PRA calculator used in the Spanish Program of Transplant Access for Highly Sensitized patients (PATHI) and validate it with EUTR, UNOS, and CTR calculators.

Methods

The anti-human leukocyte antigen (HLA) antibody profile of 42 sensitized patients on waiting list was defined, and cPRA was calculated with different PRA calculators.

Results

Despite different allelic frequencies derived from population differences in donor panel from each calculator, no differences in cPRA between the four calculators were observed. The PATHI calculator includes anti-DQA1 antibody profiles in cPRA calculation; however, no improvement in total cPRA calculation of highly sensitized patients was demonstrated.

Interpretation and conclusion

The PATHI calculator provides cPRA results comparable with those from EUTR, UNOS, and CTR calculators and serves as a tool to develop valid calculators in geographical and ethnic areas different from Europe, USA, and Canada.

HLA-A, B, DRB1, DQA1, DQB1 alleles and haplotype frequencies in Dene and Cree cohorts in Manitoba, Canada

BACKGROUND

First Nations in the Canadian province of Manitoba have disproportionately high rates of epidemic and endemic TB. Gene polymorphisms that modulate HLA Class I and II antigens are among the risk markers for TB, along with other biologic, and social determinants of health. HLA-A, B, DRB1, DQA1, DQB1 were typed in two Manitoba First Nation indigenous groups to identify and compare the frequency of gene polymorphisms that may influence susceptibility or resistance to TB.

METHODS

Participants who self-identified as either Dene or Cree enrolled into the study from two First Nation communities in Manitoba, Canada. Genomic DNA was extracted from blood samples collected with informed consent from Dene (N=63) and Cree (N=42) First Nation study participants. Participants self-reported having treated active TB, treated latent TB or no TB. HLA Class I and II molecules were typed using sequence-specific oligonucleotide (SSO) probes from commercially available kits.

RESULTS

The rates of treated active and latent TB were marginally higher among the Dene than the Cree participants (p=0.112). Class I and II HLA loci were in Hardy-Weinberg equilibrium in both the Dene and Cree groups. In this exploratory analysis of TB and HLA allele frequencies in Dene and Cree cohorts HLA-A*03 and HLA-DQB1*05:03 were significantly associated with TB.

CONCLUSIONS

The high incidence of TB in both Dene and Cree populations in Canada requires both biomedical and socioeconomic prevention and control measures. Among the former, an understanding of HLA diversity among First Nations groups may aid the development of new effective vaccine and therapeutic modalities that depend on the interaction between small molecules and specific HLA epitopes.

Can MHC-assortative partner choice promote offspring diversity? A new combination of MHC-dependent behaviours among sexes in a highly successful invasive mammal

Sexual selection involving genetically disassortative mate choice is one of several evolutionary processes that can maintain or enhance population genetic variability. Examples of reproductive systems in which choosers (generally females) select mates depending on their major histocompatibility complex (MHC) genes have been reported for several vertebrate species. Notably, the role of MHC-dependent choice not in mating contexts, but in other kinds of social interactions such as in the establishment of complex social systems, has not yet drawn significant scientific interest and is virtually absent from the literature. We have investigated male and female MHC-dependent choice in an invasive population of North American raccoons (Procyon lotor) in Germany. Both male and female raccoons rely on olfaction for individual recognition. Males have an unusually complex social system in which older individuals choose unrelated younger ones to form stable male coalitions that defend territories and a monopoly over females. We have confirmed that females perform MHC-disassortative mate choice and that this behaviour fosters genetic diversity of offspring. We have also observed that males build coalitions by choosing male partners depending on their MHC, but in an assortative manner. This is the first observation of antagonistic MHC-dependent behaviours among sexes. We show that this is the only combination of MHC-dependent partner choice that leads to outbreeding. In the case of introduced raccoons, such behaviours can act together to promote the invasive potential of the species by increasing its adaptive genetic divergence.

Anthropological analysis of Tunisian populations as inferred from HLA class I and class II genetic diversity: A meta-analysis

Despite their importance, anthropological meta-analyses which allow for comprehensive evaluation of the relationships of a given population were rare. This meta-analysis evaluates the origin of Tunisians using polymorphic profile of HLA class I (A, B), and class II (DRB1, DQB1) genes, in historical, social and cultural context, and is the only analysis in the Middle East-North Africa (MENA) region. A total of 20 eligible populations were selected from several databases, and included representing 2553 Tunisian individuals, who were compared with Mediterranean and sub-Saharan populations. In total, 204 HLA alleles were detected in Tunisians, which comprised 54 HLA-A, 76 HLA-B, 50 DRB1, and 24 DQB1 alleles. The most frequent alleles were A*02:01(24.72%) in Berbers of Zrawa, B*50:01 (13.90.11%) in Tunisian-So, DRB1*07:01 (28.66%) in Ghannouchians, and DQB1*02:01 (42.79%) in Tunisians-H. The A, B, DRB, and DQB1 genotypes of 420 individuals were further subjected to a selection study. Despite the relatively large sample size, the loci depicted non-significant negative Fnd values, an indication of overall trend to balancing selection or gene flow. Except for Berbers of Djerba, dendrograms, correspondence analyses, genetic distances and haplotype analysis demonstrated the close relatedness of Berbers, Southern and Northern Tunisians, and strong relatedness was evident to Western Mediterranean, North African and Iberian populations, but not Sub-Saharans and Eastern Mediterranean populations, including Arabs. Collectively, this suggests that the contribution of Arabs and sub-Saharans to the present Tunisian gene pool is low. In addition, all Mediterranean populations depict a typical Mediterranean substratum, except for Greeks.

Human leukocyte antigen allele linkage disequilibrium and haplotype structure in volunteer bone marrow donors of Paraná State

Background

Bone marrow transplantation has been used in the treatment of various diseases, especially hematologic diseases. The success of this treatment, among other factors, requires human leukocyte antigens (HLA) compatibility between patient and donor. Knowing the human leukocyte antigens allele group and haplotype frequencies as well as the linkage disequilibrium between alleles of different human leukocyte antigens loci can shorten the search time for a compatible bone marrow donor.

Objective

To assemble and analyze data on human leukocyte antigens frequencies available in the Laboratory of Immunogenetics and Histocompatibility (LIGH) database of the Universidade Federal do Paraná adding an estimation of the Hardy–Weinberg equilibrium and linkage disequilibrium.

Methods

The sample was composed of seven populations grouped by self-declared ancestry or inferred from the surname as follows: Laboratory of Immunogenetics and Histocompatibility database (all groups), descendants of Italians, Poles, and Asians, Afro-Brazilians, Mulattos (mixed ancestry) and Amerindians. Human leukocyte antigens genotyping was carried out using the polymerase chain reaction-sequence specific primers (PCR-SSP) and -sequence specific oligonucleotide (PCR-SSO) technologies.

Results

There were high frequencies of the HLA-A*02, HLA-B*35 and HLA-DRB1*13 allelic groups in all groups. The same was observed for the HLA-A*01-B*08-DRB1*03 haplotype except for Asian descendants. It was observed that the human leukocyte antigens Laboratory of Immunogenetics and Histocompatibility database and the Asian group are not in Hardy–Weinberg equilibrium. The Italian, Polish, Asian, Mulatto and Amerindian descendants showed haplotypes in complete linkage disequilibrium. Our results were compared with data on the human leukocyte antigens in the Paraná population available from the Brazilian Voluntary Bone Marrow Donor Registry (REDOME) and data published on the population of Curitiba and the northern region of Paraná.

Conclusions

Haplotypes frequent in the Asian group were not the most frequently observed in the Laboratory of Immunogenetics and Histocompatibility database and the National Bone Marrow Donor Registry for the state of Paraná. Linkage disequilibrium information may prove useful in the search for bone marrow donors for patients awaiting a suitable donor.

Determination of HLA-A, -B, and -DRB1 Allele and Haplotype Frequencies in the Croatian Population Based on a Family Study

In the present study, HLA allele and haplotype frequencies among Croatian families were investigated to evaluate valuable information about HLA genotypes and to compare them with data from the Croatian Bone Marrow Donor Registry (CBMDR). A total of 350 families which have been typed for the purpose of HSCT were included. All individuals were tested using PCR-SSO or PCR-SSP methods for HLA-A, -B, and -DRB1 alleles. The HLA-A-B-DRB1 haplotypes were determined by segregation and directly counted. A total of 30 HLA-A, 54 HLA-B, and 38 HLA-DRB1 alleles and 716 different HLA-A-B-DRB1 haplotypes were identified. Of these, the three most frequent alleles at HLA-A, -B, and -DRB1 loci, respectively, were A*02:01 (30.39%), A*11:01 (13.37%), A*24:02 (10.91%); B*51:01 (12.48%), B*18:01 (8.35%), B*08:01 (8.06%); DRB1*03:01 (11.20%), DRB1*01:01 (9.84%), DRB1*16:01 (9.63%). The following HLA alleles were detected only once: A*02:09, A*24:03, A*24:04, A*24:07; B*07:04, B*15:07, B*15:08, B*39:04, B*39:10, B*39:24, B*40:04; DRB1*08:03, DRB1*11:06, DRB1*13:32, DRB1*14:05. Five most frequent haplotypes were: A*01:01-B*08:01-DRB1*03:01 (5.34%), A*02:01-B*18:01-DRB1*11:04 (1.57%), A*02:01-B*27:02-DRB1*16:01 (1.50%), A*02:01-B*27:05-DRB1*01:01 (1.42%) and A*02:01-B*44:02:01G-DRB1*16:01 (1.28%). The haplotype frequencies based on the family study were compared with the frequencies from CBMDR, and similar results were obtained for all except for the HLA-A*26:01-B*38:01-DRB1*04:02 haplotype. A significantly higher frequency (P = 0.0017) of this haplotype was observed among family individuals. Nine haplotypes were unique and data about their frequencies do not exist in current databases. The data obtained in this study could be useful for anthropology, transplantation and disease association studies.

The investigation of the origin of Southern Tunisians using HLA genes

The south of Tunisia is characterized by marked ethnic diversity, highlighted by the coexistence of native Berbers with Blacks, Jews and Arab-speaking populations. Despite this heterogeneity, genetic anthropology studies investigating the origin of current Southern Tunisians were rarely reported. We examined human leukocyte antigen (HLA) class I (A, B) and class II (DRB1, DQB1) gene profiles of 250 unrelated Southern Tunisians, and compared them with those of Arab-speaking communities, along with Mediterranean and sub-Sahara African populations using genetic distances, neighbor-joining dendrograms, correspondence and haplotype analysis. In total, 137 HLA alleles were detected, which comprised 32 HLA-A, 52 HLA-B, 32 DRB1 and 21 DQB1 alleles. The most frequent alleles were HLA-A*02:01(18.02%), HLA-B*50:01 (9.11%), HLA-DRB1*07:01 (22.06%) and HLA-DQB1*02:01 (17.21%). All pairs of HLA loci show significant linkage disequilibrium. The four loci depict negative F_nd (the normalized deviate of the homozygosity) values indicating an overall trend to balancing selection. Southern Tunisians appear to be closely related to others Tunisian populations including Berbers, North Africans and Iberians. On the contrary, Southern Tunisians were distinct from Palestinian, Lebanese and Jordanian Middle Eastern Arab-speaking population, despite the deep Arab incursions and Arabization that affected Southern Tunisia. In addition, Southern Tunisians were distant from many sub-Saharan communities, evidenced by genetic distance analysis. Collectively, this indicates a limited genetic contribution of Arab invasion and Black caravans on the makeup of Southern Tunisian gene pool.