PyPop update--a software pipeline for large-scale multilocus population genomics

PyPop: a mature open-source software pipeline for population genomics

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations. Code available at: https://zenodo.org/records/10080668 and https://github.com/alexlancaster/pypop.

HLA allelic diversity in the Waorani population of Ecuador: Its significance to their ancestry and migration

The Waorani, an isolated indigenous tribe in Ecuador, have long been characterized by limited genetic diversity, with few studies delving into their genetic background. Human Leukocyte Antigen (HLA) genes which are located in the human major histocompatibility complex (MHC) provides valuable insights into population evolution due to its highly polymorphic nature. However, little is known about the HLA diversity and ancestry of the Waorani population. In this study, we sequenced eight HLA genes using Next Generation Sequencing (NGS) from 134 Waorani individuals and obtained up to four-field HLA allele resolution. Cluster and phylogenetic analysis show that the Waorani are genetically distant from other Ecuador populations, but instead show genetic affinities with the Puyanawa and Terena tribes from Brazil, as well as the Mixe tribe from Mexico. The identification of alleles common within the Waorani population, previously linked to specific health conditions, notably paves the way for future association analyses. This extensive study, employing Next-Generation Sequencing (NGS) technology, significantly enriches the sparse and segmented understanding of HLA diversity in the South American region. Our findings enhance the global comprehension of human genetic diversity and underscore the value of studying indigenous populations. Such research is vital for deepening our insights into human migration patterns and evolutionary processes.

Association of HLA class I and II genes with cervical cancer susceptibility in a Han Chinese population

Cervical cancer (CC) is one of the leading causes of cancer-related death in females worldwide. Genome-wide association studies (GWASs) have identified CC-related susceptibility loci in HLA regions. To investigate the associations between HLA genes and cervical intraepithelial neoplasia (CIN) or cervical cancer (CC), six loci of HLA class I (HLA-A, -B, and -C) and II (HLA-DRB1, -DPB1, and -DQB1) were selected for genotyping, and the associations between these alleles or their haplotypes with CIN or CC risk or protection from disease were evaluated. In total, 2193 participants, including 909 healthy individuals in the control group, 769 patients with CC, and 515 patients with CIN2+ (CIN II and III), were enrolled in the current study. HLA genes were genotyped using the NGSgo Illumina MiSeq workflow, and the associations between these loci and CIN2+ or CC at the allele and haplotype levels were analyzed. The allele frequencies of HLA-A*33:03, B*58:01, C*03:02, DPB1*05:01, and DRB1*12:01 were lower in both the CC and CIN2+ groups than in the control group, whereas those of B*55:02, C*04:03, and DPB1*03:01 were higher in the CC group than in the control group. In the histologic CC type analysis, the differences in the frequencies of these alleles in squamous cell carcinoma (SCC) of the cervix and stage I CC showed a consistent trend. In the haplotype analysis, the frequency of A*33:03-C*03:02-B*58:01 was lower in the CC and CIN2+ groups than in the control group, and that of A*24:02-C*04:03-B*15:25 was higher in the CC group than in both the control and CIN2+ groups. These three different haplotype frequencies were also identified in the FIGO CC stage analysis. In addition, in human papilloma virus (HPV) genotype analyses, the frequencies of HLA-C*03:02 and DPB1*05:01 were significantly lower in the CC and CIN2+ groups than in the control group, and in SCC subgroup, the frequencies of HLA-DQB1*04:01 and DRB1*04:05 were higher in the HPV other genotype infection group than in the HPV16 infection group. In both HPV16 single infection and coinfection with other HPVs, the frequency of haplotype A*33:03-C*03:02-B*58:01 was lower in both CC and CIN2+ than in the control group, while the frequencies of A*11:01-C*14:02-B*51:01 and A*24:02-C*03:04-B*13:01 were higher in the CIN2+ than in CC and the control group. In the HPV16 and other HPV infection comparisons, the frequencies of DRB1*04:05-DQB1*04:01-DPB1*02:01 and DRB1*11:01-DQB1*03:01-DPB1*05:01 were lower in the HPV16 infection group than in the other HPV infection group. Our results suggest that the HLA class I and II genes may affect the risk of CIN and CC as well as the histologic CC types and FIGO stages of CC in the Han Chinese population. In addition, HLA genes were associated with HPV16 infection at both the allelic and haplotype levels.

HLA Alleles, Genotype and Haplotype Analyzes from Central Anatolia Region of Turkey

Background

Although human leukocyte antigen (HLA) data for the Turkish population has been reported, there are no statistics on the HLA-DPB1 locus, which has recently received significant attention, particularly in hematopoietic stem cell transplantation. In addition, there is no study that has reported the 2-6 loci HLA haplotype distribution, 8-digit HLA allele frequency, and genotype frequency in the Turkish population.

Aims

To evaluate the low and high resolution (2-4-8 digits) HLA-A, -B, -C, -DRB1, -DQB1, -DPB1 allele data using the data of 6100 healthy individuals from the Central Anatolian region of Turkey.

Study Design

Retrospective cross-sectional study.

Methods

All tests were performed using molecular HLA techniques: low-resolution DNA-based sequence-specific oligonucleotides, low/high-resolution DNA-based sequence-specific primer, and high-resolution next generation sequencing. A total of 6100 healthy donors with a minimum of 3 loci (HLA-A, -B, -DRB1) were analyzed for their HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 data. Pypop and HLA-net GENE[RATE] were used to analyze the data.

Results

Among the HLA class I alleles, the following were the most frequently observed alleles: for HLA-A, A*02, A*24, A*03, and A*01; for HLA-B, B*35, B*51, and B*44; and for HLA-C, C*07, C*04, and C*12. Among the HLA class II alleles, the following alelles were the most frequently observed: for HLA-DRB1, DRB1*11, DRB1*04, and DRB1*13; for HLA-DQB1, DQB1*03, DQB1*05, and DQB1*06; and for HLA- DPB1, DPB1*04, DPB1*02, and DPB1*03. The most common alleles among HLA-DPB1 in the 4-digit evaluation were DPB1*04:01, DPB1*02:01, and DPB1*04:02. Among the HLA classes I and II, the following were the most frequently observed 8-digit alleles in HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 were A*02:01:01:01, B*49:01:01:01, C*04:01:01:06, DRB1*07:01:01:01, DQB1*03:01:01:02, and DPB1*02:01:02:05, respectively. The most common 6 loci haplotype was A*02~B*35~C*04~DRB1*11~DQB1*03~DPB1*04 (2.71%).

Conclusion

In this study, low and high resolution HLA-DPB1 allele frequency, 6 locus haplotype frequency and genotype frequency were reported for the first time in Turkish population. These new data can be used to map HLA in our country and may provide ideas for potential future studies.

Time trees and clock genes: a systematic review and comparative analysis of contemporary avian migration genetics

Timing is a crucial aspect for survival and reproduction in seasonal environments leading to carefully scheduled annual programs of migration in many species. But what are the exact mechanisms through which birds (class: Aves) can keep track of time, anticipate seasonal changes, and adapt their behaviour? One proposed mechanism regulating annual behaviour is the circadian clock, controlled by a highly conserved set of genes, collectively called 'clock genes' which are well established in controlling the daily rhythmicity of physiology and behaviour. Due to diverse migration patterns observed within and among species, in a seemingly endogenously programmed manner, the field of migration genetics has sought and tested several candidate genes within the clock circuitry that may underlie the observed differences in breeding and migration behaviour. Among others, length polymorphisms within genes such as Clock and Adcyap1 have been hypothesised to play a putative role, although association and fitness studies in various species have yielded mixed results. To contextualise the existing body of data, here we conducted a systematic review of all published studies relating polymorphisms in clock genes to seasonality in a phylogenetically and taxonomically informed manner. This was complemented by a standardised comparative re-analysis of candidate gene polymorphisms of 76 bird species, of which 58 are migrants and 18 are residents, along with population genetics analyses for 40 species with available allele data. We tested genetic diversity estimates, used Mantel tests for spatial genetic analyses, and evaluated relationships between candidate gene allele length and population averages for geographic range (breeding- and non-breeding latitude), migration distance, timing of migration, taxonomic relationships, and divergence times. Our combined analysis provided evidence (i) of a putative association between Clock gene variation and autumn migration as well as a putative association between Adcyap1 gene variation and spring migration in migratory species; (ii) that these candidate genes are not diagnostic markers to distinguish migratory from sedentary birds; and (iii) of correlated variability in both genes with divergence time, potentially reflecting ancestrally inherited genotypes rather than contemporary changes driven by selection. These findings highlight a tentative association between these candidate genes and migration attributes as well as genetic constraints on evolutionary adaptation.

Analysis of the Origin of Emiratis as Inferred from a Family Study Based on HLA-A, -C, -B, -DRB1, and -DQB1 Genes

In this study, we investigated HLA class I and class II allele and haplotype frequencies in Emiratis and compared them to those of Asian, Mediterranean, and Sub-Saharan African populations.

METHODS

Two-hundred unrelated Emirati parents of patients selected for bone marrow transplantation were genotyped for HLA class I (A, B, C) and class II (DRB1, DQB1) genes using reverse sequence specific oligonucleotide bead-based multiplexing. HLA haplotypes were assigned with certainty by segregation (pedigree) analysis, and haplotype frequencies were obtained by direct counting. HLA class I and class II frequencies in Emiratis were compared to data from other populations using standard genetic distances (SGD), Neighbor-Joining (NJ) phylogenetic dendrograms, and correspondence analysis.

RESULTS

The studied HLA loci were in Hardy-Weinberg Equilibrium. We identified 17 HLA-A, 28 HLA-B, 14 HLA-C, 13 HLA-DRB1, and 5 HLA-DQB1 alleles, of which HLA-A*02 (22.2%), -B*51 (19.5%), -C*07 (20.0%), -DRB1*03 (22.2%), and -DQB1*02 (32.8%) were the most frequent allele lineages. DRB1*03~DQB1*02 (21.2%), DRB1*16~DQB1*05 (17.3%), B*35~C*04 (11.7%), B*08~DRB1*03 (9.7%), A*02~B*51 (7.5%), and A*26~C*07~B*08~DRB1*03~DQB1*02 (4.2%) were the most frequent two- and five-locus HLA haplotypes. Correspondence analysis and dendrograms showed that Emiratis were clustered with the Arabian Peninsula populations (Saudis, Omanis and Kuwaitis), West Mediterranean populations (North Africans, Iberians) and Pakistanis, but were distant from East Mediterranean (Turks, Albanians, Greek), Levantine (Syrians, Palestinians, Lebanese), Iranian, Iraqi Kurdish, and Sub-Saharan populations.

CONCLUSIONS

Emiratis were closely related to Arabian Peninsula populations, West Mediterranean populations and Pakistanis. However, the contribution of East Mediterranean, Levantine Arab, Iranian, and Sub-Saharan populations to the Emiratis' gene pool appears to be minor.

Association Between HLA Polymorphisms and Sympathetic Ophthalmia in Han Chinese

PURPOSE

Sympathetic ophthalmia (SO) is considered as an autoimmune disease with unclear mechanisms. This study investigated the relationship between HLA polymorphisms and SO.

METHODS

HLA typing was performed using the LABType reverse SSO DNA typing method. The allele and haplotype frequencies were assessed using the PyPop software. Statistical significance of genotype distributions between 116 patients and 84 healthy individuals (control) was determined using Fisher's exact test or Pearson's chi-squared test.

RESULTS

The SO group had a higher frequency of HLA-DRB1 * 04:05, HLA-DQB1 * 04:01, DRB1 * 04:05-DQB1 * 04:01 haplotype   as compared to the control group (Pc < 0.001 for all).

CONCLUSION

This study revealed that DRB1 * 04:05 and DQB1 * 04:01 alleles, as well as DRB1 * 04:05-DQB1 * 04:01 haplotye could be potential risk factors for SO.

HLA allele and haplotype frequencies in Kazakhstani Russians and their relationship with other populations

HLA class I and class II genotypes from 947 Kazakhstani individuals of Russian origin were analyzed for investigating their most likely origin. The results were compared with similar data from other Russians (East and West), and also Worldwide populations, using standard genetic distances, neighbor-joining dendrograms, correspondence and haplotype analysis. Of the five HLA loci analyzed (HLA-A, HLA-C, HLA-B, HLA-DRB1, and HLA-DQB1) genotyped, 216 HLA alleles were identified. The most frequent alleles were A*02:01 (26.5%), B*07:02 (11.1%), C*04:01 (13.5%) and C*06:02 (12.1%), DRB1*07:01 (13.8%) and DRB1*15:01 (12.2%), and DQB1*03:01 (19.7%). Significant linkage disequilibrium was noted between all HLA pairs. DRB1*15:01 ~ DQB1*06:02 (10.5%), B*07:02 ~ C*07:02 (10.0%), B*07:02 ~ DRB1*15:01 (6.3%), and A*01:01 ~ B*08:01 (4.5%) were the most frequent two-locus haplotypes identified. Subsequent analyses showed that Kazakhstani Russians were closely related to West Russia-residing populations (Northwest Slavic, Vologda, Chelyabinsk, Moscow), East Europeans (Belarus Brest, Ukraine, Poland) and Scandinavians (Swedish, Finns), but distinct from East Russia-residing populations (Tuvians, Siberians from Chukotka, Kamchatka, and Ulchi) and East Mediterraneans (Levantines, Turks, North Macedonians, Albanians), and East Asians (Koreans, Japanese, Taiwanese, Mongolians). These results are in accordance with historical data indicating that the Russians of central Asia originate mainly from European Russia during the migratory flow of 18th and 19th centuries.

Human leukocyte antigens class I and class II alleles associated with vertical human immunodeficiency virus transmission - an exploratory study from Mumbai, India

Background Human leukocyte antigens (HLA) an important host genetic factor is responsible for influencing human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS) transmission and disease progression. Contributions of HLA I and II alleles have not been reported in the Indian population with respect to vertical HIV transmission. Aim In the current study we determined the frequencies of HLA class I and class II alleles in a cohort of children exposed to HIV through their mothers. Method In this exploratory study children perinatally exposed to HIV-1 who fit the study criteria and had completed 18 month follow-up were typed for HLA class I and class II alleles using polymerase chain reaction combined with sequence-specific oligonucleotides probes (PCR-SSOP) and sequence-specific primer (SSP) method. HLA typing was done in 30 positive and 60 HIV negative children along with confounding factors such as treatment regimens, viral load and CD4 count of the mother, feeding option, etc. SPSS software was used for statistical analysis and online docking tools for in-silico analysis. Results HLA-B*40 (p = 0.018) was significantly higher in negative children and was associated with protection, whereas HLA-A*01 (p = 0.05), HLA-B*37 (p = 0.032) and HLA-DRB1*09 (p = 0.017) were associated with transmission. Known protective allele HLA-B*27 was only present in negative children. Many specific haplotypes were exclusively present in the negative children or the positive ones. In-silico analysis was performed to predict the ability of HLA-B*40 to bind to antigenic peptides obtained from HIV-1 sequences in our study group. Limitations Small sample size is a concerning limitation of the study. Nonetheless this is a comprehensive study on HLA alleles in HIV exposed Indian children Conclusion Our study highlights the contribution of HLA class I and II alleles in the Indian children and further adds to understanding the immunogenetic mechanisms. These can be developed as markers for prediction of infection transmission. The observations also contribute to the database of genetic makeup of our population and can help in designing vaccine strategies.

Transporter Associated with Antigen Processing 1 Gene Polymorphisms Increase the Susceptibility to Tuberculosis

Purpose

Tuberculosis (TB) is known to result from a complex interaction between the host immune response and Mycobacterium infection. The transporter associated with antigen processing (TAP) plays an important role in the processing and presentation pathways for the Mycobacterium tuberculosis (M. tb) antigen. To investigate the possible association of the TAP1 and TAP2 genes with TB.

Patients and Methods

A total of 449 TB patients and 435 control subjects were included in this study, and single nucleotide polymorphisms (SNPs) in the TAP gene, as well as TAP1 and TAP2 alleles, were genotyped.

Results

TAP gene association analysis of TB diseases showed that rs41551515-T in the TAP1 gene was significantly associated with susceptibility to TB (P=7.96E-04, OR=4.124, 95% CI: 1.683-10.102), especially pulmonary TB (PTB, P=6.84E-04, OR=4.350, 95% CI: 1.727-10.945), and the combination of rs1057141-T-rs1135216-C in the TAP1 gene significantly increased the risk of TB susceptibility (P=5.51E-05, OR=10.899, 95% CI: 2.555-46.493). Five novel TAP1 alleles were detected in Yunnan Han people, and the allele frequency of TAP1*unknown_3 (rs41555220-rs41549617-rs1057141-rs1135216-rs1057149-rs41551515: C-A-T-C-C-T) was notably increased in all TB patients, including in the PTB and EPTB subgroups, and was significantly associated with the risk of susceptibility to TB. However, no association between the TAP2 gene and TB was found in this study.

Conclusion

Host genetic variants of rs41551515-T and the combination rs1057141-T-rs1135216-C, as well as TAP1*unknown_3 may play a critical role in susceptibility to TB disease.

A Novel Finding of an HLA Allele's and a Haplotype's Relationship with SARS-CoV-2 Vaccine-Associated Subacute Thyroiditis

Subacute thyroiditis (SAT) is a thyroid disease associated with viral infections. Its relationship with major histocompatibility complex (MHC) antigens was shown before. SAT cases triggered by different types of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines have been reported. In this study, human leukocyte antigen (HLA) genotypes of 27 SAT patients (13 vaccine-associated (V-SAT) and 14 non-SARS-CoV-2-infection non-vaccine-associated (non-V-SAT)) were compared with those of 362 healthy donors. HLA analyses were performed with low-resolution DNA-based sequence-specific oligonucleotide or sequence-specific primer methods. Statistical analyses were performed using IBM SPSS Statistics 25 and Stata/MP 14.1 with the hapipf function. Allele and haplotype frequencies were estimated by PyPop and gene[RATE] tool programs. The allele frequencies of HLA-A*11, HLA-B*35, and HLA-C*04 were higher in the patient groups. Both the allele frequency of HLA-A*11 and the haplotype frequency of A*11-B*35-C*04 were higher in the V-SAT group. The A*11-B*35-C*04 haplotype, including all three loci of MHC class I genes, is shown to be associated with the disease for the first time, especially in the V-SAT group. This finding will contribute to a better understanding of the etiopathogenesis of vaccine-associated SAT and the role of HLA genotypes in the functioning mechanisms of the SARS-CoV-2 vaccines.

Distribution of 22 Single Nucleotide Polymorphisms in 13 Cytokine Genes in Malays, Chinese, and Indians in Peninsular Malaysia

Introduction: Cytokines are cell signaling glycoproteins that are particularly important in immunity and inflammatory responses. Therefore, variations, such as single nucleotide polymorphisms (SNPs), in genes encoding for cytokines may have important consequences for their roles in health. Materials and Methods: A total of 222 unrelated, healthy, and un-admixed Malays (n = 97), Chinese (n = 77), and Indians (n = 48) with a median age of 30 years old (range 21-50) were typed for 22 cytokine gene SNPs: IL-1α -889 T/C, IL-1β (-511 T/C, +3962 T/C), IL-1R pst1 1970 T/C, IL-1RA mspa1 11100 T/C, IL-4Rα +1902 G/A, IL-12 - 1188 C/A, IFN-γ +874 A/T, TGF-β (cdn 10 C/T, cdn 25 G/C), TNF-α (-308 A/G, -238 A/G) IL-2 (+166 G/T, -330 T/G), IL-4 (-1098 T/G, -590 T/C, -33 T/C), IL-6 (-174 C/G, nt565 G/A), and IL-10 (-1082 G/A, -819 C/T, -592 A/C). This involved using well-established polymerase chain reaction procedures with sequence-specific primers and restriction fragment length polymorphism methods. Results: The majority of the screened cytokine gene SNPs are polymorphic in all three ethnicities. Exceptions include TGF-β cdn 25 (G/C), IL-1β +3962 (T/C), and TNF-α -238 (A/G), which were all observed to be monomorphic in Malays, Chinese and Indians. Many of the analyzed cytokine gene SNP genotypes deviated from Hardy-Weinberg equilibrium and the three ethnic study groups were all well-separated from reference Asian, African and European populations in a principal component analysis plot. Conclusion: We successfully typed 22 SNPs in 13 cytokine genes from genetic material collected from unrelated and un-admixed Malay, Chinese and Indian individuals in Peninsular Malaysia. These new cytokine gene population datasets reveal interesting contrasts with other populations.

HLA-A, -B, -C, -DRB1 and -DQB1 allele and haplotype frequencies in Lebanese and their relatedness to neighboring and distant populations

BACKGROUND

This study examined the origin of present-day Lebanese using high-resolution HLA class I and class II allele and haplotype distributions. The study subjects comprised 152 unrelated individuals, and their HLA class I and class II alleles and two-locus and five-locus haplotypes were compared with those of neighboring and distant communities using genetic distances, neighbor-joining dendrograms, correspondence, and haplotype analyses. HLA class I (A, B, C) and class II (DRB1, DQB1) were genotyped at a high-resolution level by PCR-SSP.

RESULTS

In total, 76 alleles across the five HLA loci were detected: A*03:01 (17.1%), A*24:02 (16.5%), B*35:01 (25.7%), C*04:01 (25.3%), and C*07:01 (20.7%) were the most frequent class I alleles, while DRB1*11:01 (34.2%) and DQB1*03:01 (43.8%) were the most frequent class II alleles. All pairs of HLA loci were in significant linkage disequilibrium. The most frequent two-locus haplotypes recorded were DRB1*11:01 ~ DQB1*03:01 (30.9%), B*35:01-C*04:01 (20.7%), B*35:01 ~ DRB1*11:01 (13.8%), and A*24:02 ~ B*35:01 (10.3%). Lebanese appear to be closely related to East Mediterranean communities such as Levantines (Palestinians, Syrians, and Jordanians), Turks, Macedonians, and Albanians. However, Lebanese appear to be distinct from North African, Iberian, and Sub-Saharan communities.

CONCLUSIONS

Collectively, this indicates a limited genetic contribution of Arabic-speaking populations (from North Africa or the Arabian Peninsula) and Sub-Saharan communities to the present-day Lebanese gene pool. This confirms the notion that Lebanese population are of mixed East Mediterranean and Asian origin, with a marked European component.

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.

Characterizing the diversity of MHC conserved extended haplotypes using families from the United Arab Emirates

Aside from its anthropological relevance, the characterization of the allele frequencies of genes in the human Major Histocompatibility Complex (MHC) and the combination of these alleles that make up MHC conserved extended haplotypes (CEHs) is necessary for histocompatibility matching in transplantation as well as mapping disease association loci. The structure and content of the MHC region in Middle Eastern populations remain poorly characterized, posing challenges when establishing disease association studies in ethnic groups that inhabit the region and reducing the capacity to translate genetic research into clinical practice. This study was conceived to address a gap of knowledge, aiming to characterize CEHs in the United Arab Emirates (UAE) population through segregation analysis of high-resolution, pedigree-phased, MHC haplotypes derived from 41 families. Twenty per cent (20.5%) of the total haplotype pool derived from this study cohort were identified as putative CEHs in the UAE population. These consisted of CEHs that have been previously detected in other ethnic groups, including the South Asian CEH 8.2 [HLA- C*07:02-B*08:01-DRB1*03:01-DQA1*05:01-DQB1*02:01 (H.F. 0.094)] and the common East Asian CEH 58.1 [HLA- C*03:02-B*58:01-DRB1*03:01- DQA1*05:01-DQB1*02:01 (H.F. 0.024)]. Additionally, three novel CEHs were identified in the current cohort, including HLA- C*15:02-B*40:06-DRB1*16:02-DQB1*05:02 (H.F. 0.035), HLA- C*16:02-B*51:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.029), and HLA- C*03:02-B*58:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.024). Overall, the results indicate a substantial gene flow with neighbouring ethnic groups in the contemporary UAE population including South Asian, East Asian, African, and European populations. Importantly, alleles and haplotypes that have been previously associated with autoimmune diseases (e.g., Type 1 Diabetes) were also present. In this regard, this study emphasizes that an appreciation for ethnic differences can provide insights into subpopulation-specific disease-related polymorphisms, which has remained a difficult endeavour.

Association study of TAP and HLA-I gene combination with chronic hepatitis C virus infection in a Han population in China

Host immune system genes play key roles in the progression of chronic hepatitis C virus (HCV) infection. Transporters associated with antigen processing (TAP) play an important role in the loading of viral peptides onto MHC class I molecules. This study aimed to investigate the association between TAP gene polymorphisms and chronic HCV in a Chinese Han population. A total of 232 chronic hepatitis C (CHC) patients and 362 healthy individuals were recruited from the Han population in Yunnan province in southwest China, and a TaqMan SNP genotyping assay was used to detect six single nucleotide polymorphisms (SNPs) of TAP1 and three SNPs of TAP2 genes. The association of the TAP gene with CHC was analysed at the allele, genotype, and haplotype levels. There were no significant differences in the allele and genotype frequencies of these SNPs in the TAP gene between CHC patients and controls after Bonferroni correction. A novel TAP1 allele (TAP1*unknown_1: rs41555220-rs41549617-rs1057141-rs1135216-rs1057149-rs41551515: G-G-A-G-G-G) was only present in the CHC group, and this allele significantly increased susceptibility to CHC (p = .005, odds ratio [OR] = 11.105. 95% confidence interval [CI]: 1.362-90.558). Homozygous TAP1*03:01/TAP1*03:01 was observed only in the CHC group that exhibited an obvious risk for CHC (p = .002, OR = 9.637, 95% CI: 1.153-80.574). And the haplotype TAP1*unknown_1-TAP2*01:01 was only present in the CHC group and indicated a significant risk for CHC (p = .002, OR = 9.498, 95% CI: 1.140-79.149). We observed significant interactions among HLA-A, -B,C, TAP1, and TAP2 alleles, and combination analysis revealed that the combination of TAP1*01:01-TAP2*01:01-HLA-B*35:01 was only present in the control group (2.2%) and resulted in significantly increased resistance to CHC (p = .002, OR = 0.096, 95% CI: 0.012-0.759). Whereas, the combination of TAP1*01:01-TAP2*01:01-HLA-C*07:02 and TAP1*03:01-TAP2*01:01-HLA-C*01:02 increased the susceptibility to CHC significantly (p = .001, OR = 2.016, 95% CI: 1.309-3.106 and p = .002, OR = 8.070, 95% CI: 1.018-63.997, respectively). Our results indicated that TAP and HLA-I may exert a combined effect on CHC susceptibility in the Chinese Han population.

HLA-DRB1 and -DQB1 Alleles, Haplotypes and Genotypes in Emirati Patients with Type 1 Diabetes Underscores the Benefits of Evaluating Understudied Populations

Background: HLA class II (DR and DQ) alleles and antigens have historically shown strong genetic predisposition to type 1 diabetes (T1D). This study evaluated the association of DRB1 and DQB1 alleles, genotypes, and haplotypes with T1D in United Arab Emirates. Materials and Methods: Study subjects comprised 149 patients with T1D, and 147 normoglycemic control subjects. Cases and controls were Emiratis and were HLA-DRB1 and -DQB1 genotyped using sequence-based typing. Statistical analysis was performed using Bridging Immunogenomic Data-Analysis Workflow Gaps R package. Results: In total, 15 DRB1 and 9 DQB1 alleles were identified in the study subjects, of which the association of DRB1*03:01, DRB1*04:02, DRB1*11:01, DRB1*16:02, and DQB1*02:01, DQB1*03:02, DQB1*03:01, and DQB1*06:01 with altered risk of T1D persisted after correcting for multiple comparisons. Two-locus haplotype analysis identified DRB1*03:01∼DQB1*02:01 [0.44 vs. 0.18, OR (95% CI) = 3.44 (2.33-5.1), Pc = 3.48 × 10^-10]; DRB1*04:02∼DQB1*03:02 [0.077 vs. 0.014, OR = 6.06 (2.03-24.37), Pc = 2.3 × 10^-3] and DRB1*04:05∼DQB1*03:02 [0.060 vs. 0.010, OR = 6.24 (1.79-33.34), Pc = 0.011] as positively associated, and DRB1*16:02∼DQB1*05:02 [0.024 vs. 0.075, OR = 0.3 (0.11-0.74), Pc = 0.041] as negatively associated with T1D, after applying Bonferroni correction. Furthermore, the highest T1D risk was observed for DR3/DR4 [0.104 vs. 0.006, OR = 25.03 (8.23-97.2), Pc = 2.6 × 10^-10], followed by DR3/DR3 [0.094 vs. 0.010, OR = 8.72 (3.17-25.32), Pc = 3.18 × 10^-8] diplotypes. Conclusion: While DRB1 and DQB1 alleles and haplotypes associated with T1D in Emiratis showed similarities to Caucasian and non-Caucasian populations, several alleles and haplotypes associated with T1D in European, African, and Asian populations, were not observed. This underscores the contribution of ethnic diversity and possible diverse associations between DRB1 and DQB1 and T1D across different populations.

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.

Association of HLA-DM and HLA class II Genes with Antibody Response Induced by Inactivated Japanese Encephalitis Vaccine

HLA (HLA) class II molecules, HLA-DR, DP, and DQ, together with HLA II-like protein DM, play a dominant role in the processing and presentation of antigens, which may influence vaccine effectiveness. We previously demonstrated that variations in the HLA-DRB1, DPB1, and DQB1 genes may affect the neutralising antibody (NAb) response induced by the inactivated Japanese encephalitis vaccine (IJEV). In the present study, we genotyped HLA-DPA1, DQA1, DMA, and DMB genes and used previous HLA-DRB1, DPB1, and DQB1 data to evaluate the association of these genes with IJEV-induced NAbs, at both the seroconversion and geometric mean titres (GMTs). We confirmed the seropositive association of DQB1*02:01 and NAbs (0.156 vs. 0.075, Padj = 0.018; OR = 2.270; 95% CI = 1.285-3.999) and seronegative association of DQB1*02:02 (0.014 vs. 0.09, Padj = 0.0002; OR = 0.130; 95% CI = 0.047-0.400). Furthermore, the DMB*01:03-DMA*01:01-DPA1*01:03-DPB1*04:01 haplotype was associated with a negative response (0.020 vs. 0.074; Padj = 0.03; OR = 0.250; 95% CI = 0.097-0.649), whereas DRB1*15:02-DMB*01:01-DMA*01:01 was associated with a positive response (0.034 vs. 0; Padj = 0.044). In addition, DRB1*12:02, DRB1*13:02, DPB1*04:01, DPB1*05:01, DPB1*09:01, DQA1*06:01, and DQA1*01:02 were associated with a higher GMT of NAbs, whereas DRB1*11:01, DPB1*13:01, and DQA1*05:05 were associated with a lower GMT of NAbs. In conclusion, the present study suggests that variations in the HLA-DM and HLA class II genes, as well as their combined allotypes, may influence the IJEV NAbs at seroconversion and GMT levels. This article is protected by copyright. All rights reserved.

HLA-C*04:01 Affects HLA Class I Heterozygosity and Predicted Affinity to SARS-CoV-2 Peptides, and in Combination With Age and Sex of Armenian Patients Contributes to COVID-19 Severity

The novel SARS-CoV-2 coronavirus infection has become a global health concern, causing the COVID-19 pandemic. The disease symptoms and outcomes depend on the host immunity, in which the human leukocyte antigen (HLA) molecules play a distinct role. The HLA alleles have an inter-population variability, and understanding their link to the COVID-19 in an ethnically distinct population may contribute to personalized medicine. The present study aimed at detecting associations between common HLA alleles and COVID-19 susceptibility and severity in Armenians. In 299 COVID-19 patients (75 asymptomatic, 102 mild/moderate, 122 severe), the association between disease severity and classic HLA-I and II loci was examined. We found that the advanced age, male sex of patients, and sex and age interaction significantly contributed to the severity of the disease. We observed that an age-dependent effect of HLA-B*51:01 carriage [odds ratio (OR)=0.48 (0.28-0.80), Pbonf <0.036] is protective against severe COVID-19. Contrary, the HLA-C*04:01 allele, in a dose-dependent manner, was associated with a significant increase in the disease severity [OR (95% CI) =1.73 (1.20-2.49), Pbonf <0.021] and an advancing age (P<0.013). The link between HLA-C*04:01 and age was secondary to a stronger association between HLA-C*04:01 and disease severity. However, HLA-C*04:01 exerted a sex-dependent differential distribution between clinical subgroups [females: P<0.0012; males: P=0.48]. The comparison of HLA-C*04:01 frequency between subgroups and 2,781 Armenian controls revealed a significant incidence of HLA-C*04:01 deficiency in asymptomatic COVID-19. HLA-C*04:01 homozygous genotype in patients blueprinted a decrease in heterozygosity of HLA-B and HLA class-I loci. In HLA-C*04:01 carriers, these changes translated to the SARS-CoV-2 peptide presentation predicted inefficacy by HLA-C and HLA class-I molecules, simultaneously enhancing the appropriate HLA-B potency. In patients with clinical manifestation, due to the high prevalence of HLA-C*04:01, these effects provided a decrease of the HLA class-I heterozygosity and an ability to recognize SARS-CoV-2 peptides. Based on our observations, we developed a prediction model involving demographic variables and HLA-C*04:01 allele for the identification of potential cases with the risk of hospitalization (the area under the curve (AUC) = 86.2%) or severe COVID-19 (AUC =71%).

Allelic and haplotypic data of MHC class II Alu insertions in Ngazidja (Comoros archipelago) and insight on its historical biology

This study aimed to analyse, for the first time, five MHC class II polymorphic Alu insertions in a population with a strong Sub-Saharan African genetic background: the Ngazidja islanders and compare its allelic and haplotypic data with Worldwide populations. The genotyping was performed in 80 individuals, using simple PCR and agarose gel electrophoresis methods. Allele and haplotype frequencies, genetic diversity, Hardy-Weinberg equilibrium deviations, normalized deviate of homozygotes and pairwise linkage disequilibrium were estimated. The phylogenetic analyses included the available population data. In Ngazidja, the MHC class II Alu insertion frequencies ranged from 0.119 for AluORF10 to 0.588 for AluDPB2. Concerning haplotypes, the most predominant were the ones with only the AluDPB2 insertion allele (AluDPB2*2-AluDQA2*1-AluDQA1*1-AluDRB1*1-AluORF10*1), followed by the theoretical ancestral haplotype with no Alu insertions (AluDPB2*1-AluDQA2*1-AluDQA1*1-AluDRB1*1-AluORF10*1) and finally the haplotype with the AluDPB2 and AluDQA1 Alu insertions (AluDPB2*2-AluDQA2*1-AluDQA1*2-AluDRB1*1-AluORF10*1) with frequencies of 19.2%, 15% and 12.9%, respectively. In the phylogenetic analyses, our results indicate that the Ngazidja people are genetically differentiated from the other populations of the analysis; we found also a new haplotype that can be probably characteristic of Sub-Saharans and finally confirm the usefulness of these markers as genetic and evolutionary tools for studying genetic variations among populations of different origins.

Genetic Analysis of Plant Pathogens Natural Populations

Population genetics allow to address fundamental questions about the biology of plant pathogens. By testing specific hypotheses, population genetics provide insights into the population genetic variability of pathogens across different geographical areas, time, and associated plant hosts, as well as on the structure and differentiation of populations, and on the possibility that a population is introduced and from where it has originated. In this chapter, basic concepts of population genetics are introduced, as well as the five evolutionary factors affecting populations, that is, mutations, recombination, variation in population size, gene flow, and natural selection. A step-by-step workflow, from sampling to data analysis, on how to perform a genetic analysis of natural populations of plant pathogens is discussed. Increased knowledge of the population biology of pathogens is pivotal to improve management strategies of diseases in agricultural and forest ecosystems.

HLA repertoire of 115 UAE nationals infected with SARS-CoV-2

The class I and class II Human Leucocyte Antigens (HLA) are an integral part of the host adaptive immune system against viral infections. The characterization of HLA allele frequency in the population can play an important role in determining whether HLA antigens contribute to viral susceptibility. In this regard, global efforts are currently underway to study possible correlations between HLA alleles with the occurrence and severity of SARS-CoV-2 infection. Specifically, this study examined the possible association between specific HLA alleles and susceptibility to SARS-CoV-2 in a population from the United Arab Emirates (UAE). The frequencies of HLA class I (HLA-A, -B, and -C) and HLA class II alleles (HLA-DRB1 and -DQB1); defined using Next Generation Sequencing (NGS); from 115 UAE nationals with mild, moderate, and severe SARS-CoV-2 infection are presented here. HLA alleles and supertypes were compared between hospitalized and non-hospitalized subjects. Statistical significance was observed between certain HLA alleles and supertypes and the severity of the infection. Specifically, alleles HLA-B*51:01 and HLA-A*26:01 showed a negative association (suggestive of protection), whilst genotypes HLA-A*03:01, HLA-DRB1*15:01, and supertype B44 showed a positive association (suggestive of predisposition) to COVID-19 severity. The results support the potential use of HLA testing to differentiate between patients who require specific clinical management strategies.

Approaching Genetics Through the MHC Lens: Tools and Methods for HLA Research

The current SARS-CoV-2 pandemic era launched an immediate and broad response of the research community with studies both about the virus and host genetics. Research in genetics investigated HLA association with COVID-19 based on in silico, population, and individual data. However, they were conducted with variable scale and success; convincing results were mostly obtained with broader whole-genome association studies. Here, we propose a technical review of HLA analysis, including basic HLA knowledge as well as available tools and advice. We notably describe recent algorithms to infer and call HLA genotypes from GWAS SNPs and NGS data, respectively, which opens the possibility to investigate HLA from large datasets without a specific initial focus on this region. We thus hope this overview will empower geneticists who were unfamiliar with HLA to run MHC-focused analyses following the footsteps of the Covid-19|HLA & Immunogenetics Consortium.

Graphical user interface for the haplotype frequency estimation software Hapl-o-Mat

Population-specific human leukocyte antigen (HLA) haplotype frequencies are an essential basis of advanced algorithms for donor selection in unrelated hematopoietic stem cell transplantation. In 2016, we introduced Hapl-o-Mat, a versatile tool for haplotype frequency estimation based on an expectation-maximization algorithm (https://github.com/DKMS/hapl-o-Mat). Hapl-o-Mat is specifically tailored to the analysis of HLA genes and able to cope with the heterogeneous genotyping data usually found in donor registries. To make Hapl-o-Mat accessible to a wider range of users, we designed a graphical user interface module that considerably facilitates the interaction with the application (https://github.com/DKMS/hapl-o-Mat_GUI). We further provide a precompiled version of Hapl-o-Mat that can be used on Windows personal computers without dependency on additional software libraries (https://github.com/DKMS/hapl-o-Mat_WinBin). We are confident that these new, user-oriented features will encourage more researchers to apply Hapl-o-Mat to their data, thereby increasing knowledge and public availability of population-specific HLA haplotype frequencies.

Allelic and haplotypic HLA diversity in indigenous Malaysian populations explored using Next Generation Sequencing

The heterogenous population of Malaysia includes more than 50 indigenous groups, and characterizing their HLA diversity would not only provide insights to their ancestry, but also on the effects of natural selection on their genome. We utilized hybridization-based sequence capture and short-read sequencing on the HLA region of 172 individuals representing seven indigenous groups in Malaysia (Jehai, Kintaq, Temiar, Mah Meri, Seletar, Temuan, Bidayuh). Allele and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DQA1, -DQB1, -DPA1, and -DPB1 revealed several ancestry-informative markers. Using SNP-based heterozygosity and pairwise Fst, we observed signals of natural selection, particularly in HLA-A, -C and -DPB1 genes. Consequently, we showed the impact of natural selection on phylogenetic inference using HLA and non-HLA SNPs. We demonstrate the utility of Next Generation Sequencing for generating unambiguous, high-throughput, high-resolution HLA data that adds to our knowledge of HLA diversity and natural selection in indigenous minority groups.

Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts.

Significance of regional population HLA immunogenetic datasets in the efficacy of umbilical cord blood banks and marrow donor registries: a study of Cretan HLA genetic diversity

BACKGROUND AIMS

The high genetic diversity of HLA across populations significantly confines the effectiveness of a donor or umbilical cord blood search for allogeneic hematopoietic stem cell transplantation (HSCT). This study aims to probe the HLA immunogenetic profile of the population of Crete, a Greek region with specific geographic and historical characteristics, and to investigate potential patterns in HLA distribution following comparison with the Deutsche Knochenmarkspenderdatei (DKMS) donor registry. It also aims to highlight the importance of regional public cord blood banks (PCBBs) in fulfilling HSCT needs, especially in countries with significant genetic diversity.

METHODS

A cohort of 1835 samples representative of the Cretan population was typed for HLA class I (HLA-A, HLA-B, HLA-C) and class II (HLA-DRB1, HLA-DQB1, HLA-DPB1) loci by high-resolution second field next-generation sequencing. Data were compared with the respective HLA profiles of 12 DKMS populations (n = 20 032). Advanced statistical and bioinformatics methods were employed to assess specific intra- and inter-population genetic indexes associated with the regional and geographic distribution of HLA alleles and haplotypes.

RESULTS

A considerable HLA allelic and haplotypic diversity was identified among the Cretan samples and between the latter and the pooled DKMS cohort. Even though the HLA allele and haplotype frequency distribution was similar to regions of close geographic proximity to Crete, a clinal distribution pattern from the northern to southern regions was identified. Significant differences were also observed between Crete and the Greek population of DKMS.

CONCLUSIONS

This study provides an in-depth characterization of the HLA immunogenetic profile in Crete and reveals the importance of demographic history in HLA heterogeneity and donor selection. The novel HLA allele and haplotype frequency comparative data between the Cretan and other European populations signify the importance of regional PCBBs in prioritizing HLA diversity to efficiently promote the HSCT program at the national level and beyond.

Exome Sequencing Reveals a Putative Role for HLA-C*03:02 in Control of HIV-1 in African Pediatric Populations

Human leucocyte antigen (HLA) class I molecules present endogenously processed antigens to T-cells and have been linked to differences in HIV-1 disease progression. HLA allelotypes show considerable geographical and inter-individual variation, as does the rate of progression of HIV-1 disease, with long-term non-progression (LTNP) of disease having most evidence of an underlying genetic contribution. However, most genetic analyses of LTNP have occurred in adults of European ancestry, limiting the potential transferability of observed associations to diverse populations who carry the burden of disease. This is particularly true of HIV-1 infected children. Here, using exome sequencing (ES) to infer HLA allelotypes, we determine associations with HIV-1 LTNP in two diverse African pediatric populations. We performed a case-control association study of 394 LTNPs and 420 rapid progressors retrospectively identified from electronic medical records of pediatric HIV-1 populations in Uganda and Botswana. We utilized high-depth ES to perform high-resolution HLA allelotyping and assessed evidence of association between HLA class I alleles and LTNP. Sixteen HLA alleles and haplotypes had significantly different frequencies between Uganda and Botswana, with allelic differences being more prominent in HLA-A compared to HLA-B and C allelotypes. Three HLA allelotypes showed association with LTNP, including a novel association in HLA-C (HLA-B∗57:03, aOR 3.21, Pc = 0.0259; B∗58:01, aOR 1.89, Pc = 0.033; C∗03:02, aOR 4.74, Pc = 0.033). Together, these alleles convey an estimated population attributable risk (PAR) of non-progression of 16.5%. We also observed novel haplotype associations with HLA-B∗57:03-C∗07:01 (aOR 5.40, Pc = 0.025) and HLA-B∗58:01-C∗03:02 (aOR 4.88, Pc = 0.011) with a PAR of 9.8%, as well as a previously unreported independent additive effect and heterozygote advantage of HLA-C∗03:02 with B∗58:01 (aOR 4.15, Pc = 0.005) that appears to limit disease progression, despite weak LD (r 2 = 0.18) between these alleles. These associations remained irrespective of gender or country. In one of the largest studies of HIV in Africa, we find evidence of a protective effect of canonical HLA-B alleles and a novel HLA-C association that appears to augment existing HIV-1 control alleles in pediatric populations. Our findings outline the value of using multi-ethnic populations in genetic studies and offer a novel HIV-1 association of relevance to ongoing vaccine studies.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

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.

Clinical features, biochemistry, and HLA-DRB1 status in youth-onset type 1 diabetes in Sudan

OBJECTIVE

To further understand clinical and biochemical features, and HLA-DRB1 genotypes, in new cases of diabetes in Sudanese children and adolescents.

RESEARCH DESIGN AND METHODS

Demographic characteristics, clinical information, and biochemical parameters (blood glucose, HbA1c, C-peptide, autoantibodies against glutamic acid decarboxylase 65 [GADA] and insulinoma-associated protein-2 [IA-2A], and HLA-DRB1) were assessed in 99 individuals <18 years, recently (<18 months) clinically diagnosed with T1D. HLA-DRB1 genotypes for 56 of these Arab individuals with T1D were compared to a mixed control group of 198 healthy Arab (75%) and African (25%) individuals without T1D.

RESULTS

Mean ± SD age at diagnosis was 10.1 ± 4.3 years (range 0.7-17.6 years) with mode at 9-12 years. A female preponderance was observed. Fifty-two individuals (55.3%) presented in diabetic ketoacidosis (DKA). Mean ± SD serum fasting C-peptide values were 0.22 ± 0.25 nmol/L (0.66±0.74 ng/ml). 31.3% were autoantibody negative, 53.4% were GADA positive, 27.2% were IA-2A positive, with 12.1% positive for both autoantibodies. Association analysis compared to 198 controls of similar ethnic origin revealed strong locus association with HLA-DRB1 (p < 2.4 × 10-14 ). Five HLA-DRB1 alleles exhibited significant T1D association: three alleles (DRB1*03:01, DRB1*04:02, and DRB1*04:05) were positively associated, while three (DRB1*10:01, DRB1*15:02, and DRB1*15:03) were protective. DRB1*03:01 had the strongest association (odds ratio = 5.04, p = 1.7 × 10-10 ).

CONCLUSIONS

Young Sudanese individuals with T1D generally have similar characteristics to reported European-origin T1D populations. However, they have higher rates of DKA and slightly lower autoantibody rates than reported European-origin populations, and a particularly strong association with HLA-DRB1*03:01.

HLA haplotype frequency estimation for heterogeneous populations using a graph-based imputation algorithm

HLA haplotype frequencies are estimated from ambiguous unphased HLA genotyping data using Expectation-Maximization (EM) algorithms. Current population genetics methods require independent EM frequency estimates for each population, and assume that each population is in Hardy-Weinberg Equilibrium (HWE). The HWE assumption of EM has thus far resulted in the exclusion of individuals from mixed or unknown ethnic backgrounds from reference datasets. Multi-region populations are currently poorly served by stem cell donor registry HLA imputation and matching implementations due to the inability of such algorithms to incorporate admixture into their population genetics models. To address this unmet need, we have expanded the imputation component of our GRaph IMputation and Matching (GRIMM) framework, where imputation becomes the expectation step in an iterative EM algorithm. Our novel multi-region EM implementation considers region as a Bayesian prior, enabling integration of HLA information from multiple single-region population groups, and for the first time including individuals with ambiguous or mixed ethnic backgrounds. We show that our multi-region EM produces much higher likelihood values and better haplotype recovery as measured by Kullback-Leibler divergence than all evaluated EM implementations when tested on real datasets of US donor registry HLA typings as well as simulated multi-region datasets of ambiguous HLA typings.

Origin of the Ukrainian minority of Kazakhstan as inferred from HLA-A, -B, -C, -DRB1, and -DQB1 alleles and haplotypes distribution

The present-day population of Kazakhstan comprises more than a hundred ethnic and linguistic groups, of which Ukrainians rates numerically as the fourth, following native Kazakhs, Russians, and Uzbeks. We investigated the HLA profile of Kazakhstani Ukrainians and compared them with those of other populations by applying genetic distances, correspondence analysis, admixture, and haplotype examination. Of the 128 HLA alleles identified, A*02:01 (28.64%), B*07:02 (9.71%), B*13:02 (9.71%), C*06:02 (15.05%), DRB1*07:01 (15.20%), and DQB1*03:01 (24.76%) were the utmost common alleles. The most common 2-locus haplotypes identified were A*02:01 ~ B*13:02 (6.66%), B*13:02 ~ DRB1*07:01 (6.57%), B*13:02 ~ C*06:02 (10.10%), and DRB1*01:01 ~ DQB1*05:01 (12.12%), while A*02:01 ~ B*13:02 ~ C*06:02 ~ DRB1*07:01 ~ DQB1*02:01 were the most frequent five-locus haplotype (4.04%) in Kazakhstani Ukrainians. Comparative analysis revealed that Ukrainians of Kazakhstan are closely related to Eastern Europeans (included Ukrainians), and European Russians, but distant from Asian populations. This supports the historical notion that Kazakhstani Ukrainians originated from Ukraine, following their massive migrations to central Asia in the 18th-20th centuries. This article is protected by copyright. All rights reserved.

Haplotypic Associations and Differentiation of MHC Class II Polymorphic Alu Insertions at Five Loci With HLA-DRB1 Alleles in 12 Minority Ethnic Populations in China

The analysis of polymorphic variations in the human major histocompatibility complex (MHC) class II genomic region on the short-arm of chromosome 6 is a scientific enquiry to better understand the diversity in population structure and the effects of evolutionary processes such as recombination, mutation, genetic drift, demographic history, and natural selection. In order to investigate associations between the polymorphisms of HLA-DRB1 gene and recent Alu insertions (POALINs) in the HLA class II region, we genotyped HLA-DRB1 and five Alu loci (AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10), and determined their allele frequencies and haplotypic associations in 12 minority ethnic populations in China. There were 42 different HLA-DRB1 alleles for ethnic Chinese ranging from 12 alleles in the Jinuo to 28 in the Yugur with only DRB1^∗08:03, DRB1^∗09:01, DRB1^∗12:02, DRB1^∗14:01, DRB1^∗15:01, and DRB1^∗15:02 present in all ethnic groups. The POALINs varied in frequency between 0.279 and 0.514 for AluDPB2, 0 and 0.127 for AluDQA2, 0.777 and 0.995 for AluDQA1, 0.1 and 0.455 for AluDRB1 and 0.084 and 0.368 for AluORF10. By comparing the data of the five-loci POALIN in 13 Chinese ethnic populations (including Han-Yunnan published data) against Japanese and Caucasian published data, marked differences were observed between the populations at the allelic or haplotypic levels. Five POALIN loci were in significant linkage disequilibrium with HLA-DRB1 in different populations and AluDQA1 had the highest percentage association with most of the HLA-DRB1 alleles, whereas the nearby AluDRB1 indel was strongly haplotypic for only DRB1^∗01, DRB1^∗10, DRB1^∗15 and DRB1^∗16. There were 30 five-locus POALIN haplotypes inferred in all populations with H5 (no Alu insertions except for AluDQA1) and H21 (only AluDPB2 and AluDQA1 insertions) as the two predominant haplotypes. Neighbor joining trees and principal component analyses of the Alu and HLA-DRB1 polymorphisms showed that genetic diversity of these genomic markers is associated strongly with the population characteristics of language family, migration and sociality. This comparative study of HLA-DRB1 alleles and multilocus, lineage POALIN frequencies of Chinese ethnic populations confirmed that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic and evolutionary markers for the identification of allele and haplotype lineages and genetic variations within the same and/or different populations.

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.

Allele and haplotype frequencies of human leukocyte antigen-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 by next generation sequencing-based typing in Koreans in South Korea

Allele frequencies and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 have been rarely reported in South Koreans using unambiguous, phase-resolved next generation DNA sequencing. In this study, HLA typing of 11 loci in 173 healthy South Koreans were performed using next generation DNA sequencing with long-range PCR, TruSight^® HLA v2 kit, Illumina MiSeqDx platform system, and Assign^™ for TruSight^™ HLA software. Haplotype frequencies were calculated using the PyPop software. Direct counting methods were used to investigate the association with DRB1 for samples with only one copy of a particular secondary DRB locus. We compared these allele types with the ambiguous allele combinations of the IPD-IMGT/HLA database. We identified 20, 40, 26, 31, 19, 16, 4, and 16 alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, respectively. The number of HLA-DRB3/4/5 alleles was 4, 5, and 3, respectively. The haplotype frequencies of most common haplotypes were as follows: A*33:03:01-B*44:03:01-C*14:03-DRB1*13:02:01-DQB1*06:04:01-DPB1*04:01:01 (2.89%), A*33:03:01-B*44:03:01-C*14:03 (4.91%), DRB1*08:03:02-DQA1*01:03:01-DQB1*06:01:01-DPA1*02:02:02-DPB1*05:01:01 (5.41%), DRB1*04:05:01-DRB4*01:03:01 (12.72%), DQA1*01:03:01-DQB1*06:01:01 (13.01%), and DPA1*02:02:02-DPB1*05:01:01 (30.83%). In samples with only one copy of a specific secondary DRB locus, we examined its association with DRB1. We, thus, resolved 10 allele ambiguities in HLA-B, -C (each exon 2+3), -DRB1, -DQB1, -DQA1, and -DPB1 (each exon 2) of the IPD-IMGT/HLA database. Korean population was geographically close to Japanese and Han Chinese populations in the genetic distances by multidimensional scaling (MDS) plots. The information obtained by HLA typing of the 11 extended loci by next generation sequencing may be useful for more exact diagnostic tests on various transplantations and the genetic population relationship studies in South Koreans.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7, and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from ‘hitchhiking’ alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

High-resolution HLA allele and haplotype frequencies in several unrelated populations determined by next generation sequencing: 17th International HLA and Immunogenetics Workshop joint report

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

HLA-A, -C, -B, -DRB1, -DQA1, and -DQB1 Allele and Haplotype Repertoires in the Albanian Population from Kosovo

To give new insight into the huge polymorphism of HLA system and supplement the existing data, an analysis of HLA alleles and HLA-A~C~B~DRB1~DQA1~DQB1 haplotype distribution in 124 Albanian individuals from Kosovo was performed. All samples were HLA-typed applying the polymerase chain reaction-sequence specific oligonucleotide probing (PCR-SSOP) method and all ambiguous HLA typing results were additionally confirmed by the standard PCR-Sequence Specific Primers (PCR-SSP) high-resolution protocol. Twenty-two HLA-A, 21 HLA-C, 37 HLA-B, 27 HLA-DRB1, 11 HLA-DQA1 and 14 HLA-DQB1 allele groups were detected. Sixteen out of 172 different six-locus estimated haplotypes were found at a frequency higher than 1.00% with a cumulative frequency of 28.82%. The most prevalent haplotype was found to be HLA-A*02:01~C*07:01~B*18:01~DRB1*11:04~DQA1*05:05~DQB1*03:0(5.2%).A total of 13 haplotypes were observed with higher frequency than in populations reported in HaploStats and The Allele Frequency Net Database. The proposed origin of the most frequent haplotypes reflects a basic Euro-Mediterranean background of Albanians in Kosovo. This is the first report of high-resolution HLA-A~C~B~DRB1~DQA1~DQB1 haplotype distribution among the Albanian population from Kosovo, which provides valuable anthropological data and confirms population-specific characteristics.

Genetic variability and biodiversity of Ukrainian Gray cattle by the BoLA-DRB3 gene

At the current stage of genetic studies of cattle, more and more attention is being drawn to autochthonous breeds. Native cattle have a number of prominent phenotypic traits and have preserved unique genes and their combinations lost by modern commercial breeds, which would be valuable to use in selective programs. We surveyed polymorphism of the Ukrainian autochthonous Gray breed according to alleles of exon 2 of the BoLA-DRB3 gene. The uniqueness of the gene lies in the broad variability of its allele variants. Significant informativeness at DNA level is quite important for genetic studies. We surveyed allele polymorphism using the PCR-RLFP method on DNA isolated from 88 samples of blood of cows and 5 samples of sperm. We identified 28 alleles, of which 23 variants were nomenclature ones and 5 (jba, *jab, *jbb, *nad and *nda) were “without established nomenclature”, their share accounting for 8.9%. Four alleles *06, *12, *16 and *jba had a frequency above 5% and occupied 69.9% of the breed’s allele fund overall. The commonest allele was BoLA-DRB3.2*16 (44.1%). In total, we found 40 genotypes. Considering the significant dominance of variant *16, as expected, 5 genotypes with its inclusion occurred: *16/*16, *12/*16, *06/*16, *16/*24 and *jba/*16. It was present in the genotype of two out three studied animals. Parameters of heterozygosity, effective number of alleles, Shannon and Pielou indices indicate that Ukrainian Gray cattle are characterized by lowest level of genetic variability and biodiversity according to the BoLA-DRB3 gene compared with other breeds. Due to significant dominance of allele *16, the breed has no inbred motifs. We noted deviation toward increase in homozygosity without deviations from the norm of the distribution according to Hardy-Weinberg equilibrium. The obtained results will be used for genetic-populational programs with the purpose of improving the genetic potential of cattle breeds in terms of economically beneficial traits and diseases of cattle.

HLA Profile of Kami Population Refutes the Earlier Proposition of Exclusive Closer Genetic Affinity of All the Gorkhas to Mongoloids

OBJECTIVE

Based on the HLA profile of Indian Gorkhas, Debnath and Chaudhuri (2006) proposed that Gorkhas are genetically closer to Mongoloids, and they may have originated from Mongolians or Tibetan stocks. However, the major limitation of the earlier study was that Gorkhas comprise 2 broad groups, i.e. Tibeto-Burmans and Indo-Aryans. Besides, Gorkhas have an assemblage of many sociocultural and linguistically distinct populations such as Rai, Magar, Limbu, Tamang, Newar, Bahun, Kami, and so on. Thus, the generalization of the findings on Gorkhas by considering them as a single homogenous population may not be free from biases. Therefore, the present study aims to understand the genetic affinity of a constituent population from the Gorkha community, i.e. Kami, based on HLA polymorphism.

METHODS

First field HLA typing was performed among 158 Kami individuals by PCR-SSP methods.

RESULTS

The most frequent genes observed were HLA-A*11, HLA-B*15, HLA-DRB1*15. The frequency of HLA-DRB1*15 reported here is the highest recorded among the North Indian population to date, which is a noteworthy finding of the study. The hierarchical cluster analysis and principal component analysis showed that the Kami population lies within the cluster of the Indian subcontinental population.

CONCLUSION

The study refutes the earlier proposition of exclusive belongingness of all the Gorkhas to Mongoloids.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Allogeneic hematopoietic stem cell transplant recipients in Spain: Human leukocyte antigen characteristics and diversity by high-resolution analysis

There are many studies on the polymorphism of the HLA system in healthy donor populations, such as registries of unrelated bone marrow donors. Investigations on the characterization of the HLA complex in hematopoietic stem cell transplant (HSCT) patients, however, are scarce, at least in the Spanish population. This study presents a large-scale analysis of allelic diversity and HLA distribution at a high-resolution level in 2886 patients undergoing HSCT in Spanish centres of the "Grupo Español de Trasplante Hematopoyético y Terapia Celular" during a period of 11 years. Allelic diversity analysis identified 67 HLA-A, 133 HLA-B, 60 HLA-C, 63 HLA-DRB1, 24 HLA-DQB1 and 27 HLA-DPB1 different alleles. Rare alleles were detected among which 33 alleles had not been reported in the European catalog of common and well-documented HLA alleles. Regarding the distribution of five genes-haplotypes, it was observed that the five most frequent extended haplotypes found in our population were between the most common in other Spanish populations, both in patients and in healthy subjects. However, some particular haplotypes were also detected. Bilocus associations HLA-C ~ B and -DRB1 ~ DQB1 were analyzed in order to predict the probability of finding 10/10 matched donors in registries. We found HLA-B alleles showing a great diversity of combinations with HLA-C alleles and unusual associations involving a negative predicting factor. In the field of adoptive therapies, our work supports the necessity to expand further research of TCR-engineered cells, adoptive transfer of virus-specific T-cells and vaccines to target HLA alleles other than A*02:01. HLA alleles such as A*01:01, A*03:01, A*24:02, B*44:03, B*07:02 or B*51:01, might be considered new targets due to its high frequency in our population.

Chronic idiopathic axonal polyneuropathy: Electrophysiological progression and human leukocyte antigen associations

BACKGROUND

We aimed to describe the electrophysiological progression rate of chronic idiopathic axonal polyneuropathy (CIAP) and look into the potential role of human leukocyte antigen (HLA) genetic susceptibility in its development.

METHODS

We recruited 57 patients with CIAP (mean age at diagnosis 67, mean follow-up 7 years). The assessments included clinical and electrophysiological data and HLA-DQ genotyping.

RESULTS

The DQA1*05 allele was found more frequently in patients than in healthy controls (odds ratio, 1.96, P = .011). In patients with length-dependent CIAP, a linear effect of time on the electrophysiological findings was found in the superficial radial (3.2% mean annual decrement, P < .001), sural (4.7% mean annual decrement, P = .002) and tibial nerve (6.1% mean annual decrement, P = .007) amplitudes, independently from age or gender.

CONCLUSIONS

Patients with length-dependent CIAP, show a linear progression over time. Interesting associations of HLA-DQA1*05 allele with length-dependent CIAP and non-DQ2/DQ8 with idiopathic sensory ganglionopathy were found. These merit further investigation in larger cohorts and may suggest a role of the immune system in the pathogenesis of CIAP.

Distributions of HLA-A, -B, and -DRB1 alleles typed by amplicon-based next generation sequencing in Korean volunteer donors for unrelated hematopoietic stem cell transplantation

HLA genes play a pivotal role for successful hematopoietic stem cell transplantation (HSCT). There is an increasing need for sophisticated screening of donor HLA genotypes for unrelated HSCT. Next generation sequencing (NGS) has emerged as an alternative for classical Sanger sequence for HLA typing. In this study, HLA-A, -B, and -DRB1 alleles were genotyped at the allelic (6-digit) level using MiSeqDx in 26,202 volunteers from the Korean Network for Organ Sharing. Exon 2 and 3 of HLA-A and -B and exon 2 of HLA-DRB1 were amplified by polymerase chain reaction (PCR) and each allele was determined by matching the targeted exons and the reference sequence consisting of the IPD-IMGT/HLA Database. Seventy alleles of HLA-A, 102 alleles of HLA-B, and 69 alleles of HLA-DRB1 were identified. According to common and well-documented catalogs, 34 alleles in HLA-A, 61 in HLA-B, and 45 in HLA-DRB1 locus were common alleles, and 12, 14, and 11 kinds, were well-documented alleles, respectively. Thirteen novel alleles including 3 alleles in HLA-A, 8 alleles in HLA-B, and 2 alleles in HLA-DRB1 loci were found. Ten haplotypes with a frequency of more than 1.0% accounted for 22.4% of the total haplotype frequencies. Cis/trans ambiguities of HLA-A and -B loci by combination of exons 2 and 3 were analyzed to be 0.17% of 3 and 3.95% of 22 genotypes, respectively. This information on rare and novel alleles found by accurate HLA typing with NGS may be helpful for unrelated HSCT among Koreans.

HLA allele and haplotype frequencies in the Panamanian population

In this study, we report for the first time HLA allele and haplotype frequencies in the modern Panamanian population at a two-field (four digits) resolution level. Reported frequencies were calculated from genotype data for the HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 loci of 462 healthy unrelated Panamanian adults of Hispanic ethnicity. In addition to providing new insights on the allelic structure of the Panamanian population and its origin, these data are critical for better planning of healthcare strategies in the country and for future research exploring the association with certain chronic and infectious diseases.

The role of HLA in Balkan endemic nephropathy

Balkan endemic nephropathy (BEN), a progressive chronic tubulointerstitial disease, occurs in the endemic focus of Croatia in a population of about 10,000 inhabitants. One of its most peculiar characteristics is a strong association with upper tract urothelial carcinoma (UTUC). Despite a high number of studies, currently there are insufficient data about the association of BEN and HLA genes. The aim of this study was to investigate the polymorphism of HLA-A, -B, and -DRB1 alleles and haplotypes among BEN patients and to determine whether an association between HLA and BEN exists. In this study, we investigated HLA-A, -B, and -DRB1 alleles and haplotypes in a population of patients with BEN (N = 111) and matched healthy controls (N = 190). All individuals were tested by PCR-SSO and PCR-SSP methods to assess the possible contribution of HLA alleles and haplotypes to the development of/protection from BEN. Our results showed a positive association between the presence of HLA-B*35:02 and DRB1*04:02 alleles and BEN (P = 0.0179 and P = 0.0151, respectively) in contrast to the protective effect of HLA-A*01:01, B*27:05 and B*57:01 alleles (P = 0.0111, P = 0.0330 and P = 0.0318, respectively). Moreover, when BEN patients' HLA haplotypes were compared to controls, two haplotypes were associated with BEN susceptibility among Croatians (HLA-A*02:01~B*08:01~DRB1*03:01 and HLA-A*02:01~B*27:02~DRB1*16:01, P = 0.0064 and P = 0.0023, respectively), while haplotypes HLA-A*02:01~B*27:05~DRB1*01:01 and HLA-A*02:01~B*38:01~DRB1*13:01 each showed a possible protective effect (P = 0.0495). Our results point toward genetic susceptibility to BEN and observed differences in both susceptible/protective HLA profiles indicate the necessity of further studies in order to elucidate the pathogenesis of this disease.