High-resolution HLA A∼B∼DRB1 haplotype frequencies from the Ezer Mizion Bone Marrow Donor Registry in Israel

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.

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.

HLA molecular study of patients in a public kidney transplant program in Guatemala

Guatemala is a country located in Central America, and while it is one of the most populated countries in the region, the genetic diversity of the population has been poorly analyzed. Currently, there are no analyses of the distribution of human leukocyte antigen (HLA) system alleles in mixed ancestry (i.e., ladino) populations in Guatemala. The HLA system exhibits the most extensive polymorphism in the human genome and has been extensively analyzed in a large number of studies related to disease association, transplantation, and population genetics (with particular importance in the understanding of diversity in the human population). Here, we present HLA typing data from 127 samples of unrelated individuals from the kidney transplant program of the San Juan de Dios General Hospital (Guatemala City) using a PCR-SSOP-based (PCR-sequence specific oligonucleotide probes) typing method. We found 16 haplotypes that accounted for 39.76 % of the total haplotype diversity, of which thirteen have been reported previously in Native American populations and three have been reported in European populations. The analyses showed no deviations from Hardy-Weinberg equilibrium, and admixture estimates calculated with k = 3 ancestral components showed that Native American was the most represented component, followed by the European component. The African component was less prominent in the Guatemala mixed ancestry sample in comparison to samples from other countries in Central America. The HLA-based admixture results for Central America showed a continuum in the distribution of Native American, European and African ancestries throughout the region, which is consistent with the complex demographic history of the region.

MHC Haplotyping of SARS-CoV-2 Patients: HLA Subtypes Are Not Associated with the Presence and Severity of COVID-19 in the Israeli Population

HLA haplotypes were found to be associated with increased risk for viral infections or disease severity in various diseases, including SARS. Several genetic variants are associated with COVID-19 severity. Studies have proposed associations, based on a very small sample and a large number of tested HLA alleles, but no clear association between HLA and COVID-19 incidence or severity has been reported. We conducted a large-scale HLA analysis of Israeli individuals who tested positive for SARS-CoV-2 infection by PCR. Overall, 72,912 individuals with known HLA haplotypes were included in the study, of whom 6413 (8.8%) were found to have SARS-CoV-2 by PCR. A total of 20,937 subjects were of Ashkenazi origin (at least 2/4 grandparents). One hundred eighty-one patients (2.8% of the infected) were hospitalized due to the disease. None of the 66 most common HLA loci (within the five HLA subgroups: A, B, C, DQB1, DRB1) was found to be associated with SARS-CoV-2 infection or hospitalization in the general Israeli population. Similarly, no association was detected in the Ashkenazi Jewish subset. Moreover, no association was found between heterozygosity in any of the HLA loci and either infection or hospitalization. We conclude that HLA haplotypes are not a major risk/protecting factor among the Israeli population for SARS-CoV-2 infection or severity. Our results suggest that if any HLA association exists with the disease it is very weak, and of limited effect on the pandemic.

Single haplotype admixture models using large scale HLA genotype frequencies to reproduce human admixture

The human leukocyte antigen (HLA) is the most polymorphic region in humans. Anthropologists use HLA to trace populations' migration and evolution. However, recent admixture between populations can mask the ancestral haplotype frequency distribution. We present a statistical method based on high-resolution HLA haplotype frequencies to resolve population admixture using a non-negative matrix factorization formalism and validated using haplotype frequencies from 56 world populations. The result is a minimal set of source components (SCs) decoding roughly 90% of the total variance in the studied admixtures. These SCs agree with the geographical distribution, phylogenies, and recent admixture events of the studied groups. With the growing population of multi-ethnic individuals, or individuals that do not report race/ethnic information, the HLA matching process for stem-cell and solid organ transplants is becoming more challenging. The presented algorithm provides a framework that facilitates the breakdown of highly admixed populations into SCs, which can be used to better match the rapidly growing population of multi-ethnic individuals worldwide.

The association between HLA and non-Hodgkin lymphoma subtypes, among a transplant-indicated population

Several studies have implicated HLA in non-Hodgkin lymphoma (NHL) subtype etiology. However, NHL patients indicated for stem cell transplants are underrepresented in these reports. We therefore evaluated the association between HLA and NHL subtypes among a transplant-indicated population. One thousand three hundred and sixty-six NHL patients HLA-typed and indicated for transplant at the City of Hope National Medical Center (Duarte, CA) were compared to 10,271 prospective donors. Odds ratios and 95% confidence intervals were calculated for HLA haplotype and alleles, adjusted for sex and age. The HLA-A*0201∼C*0602∼B*1302∼DRB1*0701∼DQB1*0201 haplotype was significantly associated with follicular lymphoma (FL) risk among Caucasians. Several haplotypes were associated with diffuse large B-cell lymphoma (DLBCL) risk among Caucasians, including the previously implicated DLBCL risk loci, HLA-B*0801. The HLA-A*0101 allele was also observed to be associated with mantle cell lymphoma (MCL) risk. Our results support the association between previously reported susceptibility loci and FL and suggest potentially new DLBCL and MCL risk loci.

The role of HLA-A*33:01 in patients with cholestatic hepatitis attributed to terbinafine

BACKGROUND & AIMS

Terbinafine is an antifungal agent that has been associated with rare instances of hepatotoxicity. In this study we aimed to describe the presenting features and outcomes of patients with terbinafine hepatotoxicity and to investigate the role of human leukocyte antigen (HLA)-A*33:01.

METHODS

Consecutive high causality cases of terbinafine hepatotoxicity enrolled into the Drug Induced Liver Injury Network were reviewed. DNA samples underwent high-resolution confirmatory HLA sequencing using the Ilumina MiSeq platform.

RESULTS

All 15 patients with terbinafine hepatotoxicity were more than 40 years old (median = 57 years), 53% were female and the median latency to onset was 38 days (range 24 to 114 days). At the onset of drug-induced liver injury, 80% were jaundiced, median serum alanine aminotransferase was 448 U/L and alkaline phosphatase was 333 U/L. One individual required liver transplantation for acute liver failure during follow-up, and 7 of the 13 (54%) remaining individuals had ongoing liver injury at 6 months, with 4 demonstrating persistently abnormal liver biochemistries at month 24. High-resolution HLA genotyping confirmed that 10 of the 11 (91%) European ancestry participants were carriers of the HLA-A*33:01, B*14:02, C*08:02 haplotype, which has a carrier frequency of 1.6% in European Ancestry population controls. One African American patient was also an HLA-A*33:01 carrier while 2 East Asian patients were carriers of a similar HLA type: A*33:03. Molecular docking studies indicated that terbinafine may interact with HLA-A*33:01 and A*33:03.

CONCLUSIONS

Patients with terbinafine hepatotoxicity most commonly present with a mixed or cholestatic liver injury profile and frequently have residual evidence of chronic cholestatic injury. A strong genetic association of HLA-A*33:01 with terbinafine drug-induced liver injury was confirmed amongst Caucasians.

LAY SUMMARY

A locus in the human leukocyte antigen gene (HLA-A*33:01, B*14:02, C*08:02) was significantly overrepresented in Caucasian and African American patients with liver injury attributed to the antifungal medication, terbinafine. These data along with the molecular docking studies demonstrate that this genetic polymorphism is a plausible risk factor for developing terbinafine hepatotoxicity and could be used in the future to help doctors make a diagnosis more rapidly and confidently.

Genetic Basis of Delayed Hypersensitivity Reactions to Drugs in Jewish and Arab Populations

Genetic variation can affect drug pharmacokinetics and pharmacodynamics and contribute to variability between individuals in response to medications. Specifically, differences in allele frequencies among individuals and ethnic groups have been associated with variation in their propensity to develop drug hypersensitivity reactions (HSRs). This article reviews the current knowledge on the genetic background of HSRs and its relevance to Jewish and Arab populations. The focus is on human leukocyte antigen (HLA) alleles and haplotypes as predictive markers of HSRs ("immunopharmacogenetics"), but other genes and alleles are described as well. Also discussed is the translation of the pharmacogenetic information to practice recommendations.

HLA typing using genome wide data reveals susceptibility types for infections in a psychiatric disease enriched sample

BACKGROUND

The infections Toxoplasma gondii (T. gondii), cytomegalovirus, and Herpes Simplex Virus Type 1 (HSV1) are common persistent infections that have been associated with schizophrenia and bipolar disorder. The major histocompatibility complex (MHC, termed HLA in humans) region has been implicated in these infections and these mental illnesses. The interplay of MHC genetics, mental illness, and infection has not been systematically examined in previous research.

METHODS

In a cohort of 1636 individuals, we used genome-wide association data to impute 7 HLA types (A, B, C, DRB1, DQA1, DQB1, DPB1), and combined this data with serology data for these infections. We used regression analysis to assess the association between HLA alleles, infections (individually and collectively), and mental disorder status (schizophrenia, bipolar disorder, controls).

RESULTS

After Bonferroni correction for multiple comparisons, HLA C∗07:01 was associated with increased HSV1 infection among mentally healthy controls (OR 3.4, p = 0.0007) but not in the schizophrenia or bipolar groups (P > 0.05). For the multiple infection outcome, HLA B∗ 38:01 and HLA C∗12:03 were protective in the healthy controls (OR ≈ 0.4) but did not have a statistically-significant effect in the schizophrenia or bipolar groups. T. gondii had several nominally-significant positive associations, including the haplotypes HLA DRB∗03:01 ∼ HLA DQA∗05:01 ∼ HLA DQB∗02:01 and HLA B∗08:01 ∼ HLA C∗07:01.

CONCLUSIONS

We identified HLA types that showed strong and significant associations with neurotropic infections. Since some of these associations depended on mental illness status, the engagement of HLA-related pathways may be altered in schizophrenia due to immunogenetic differences or exposure history.

East Meets West-Impact of Ethnicity on Donor Match Rates in the Ezer Mizion Bone Marrow Donor Registry

HLA haplotype frequencies in a volunteer bone marrow donor registry should reflect the frequencies of potential transplant recipients served by that registry, a challenge in a country with diverse subethnicities of immigrants from Eastern and Western cultures, such as Israel. We evaluated the likelihood of finding suitable donors for hypothetical patients drawn from defined subethnicities in the Ezer Mizion Bone Marrow Donor Registry (EM BMDR) from donors both within and outside the registry now and during the coming decade. On average, bioinformatics modeling predicts that, given current donor recruitment trends, 6/6 high-resolution HLA match rates for Israelis, which currently stand at 40% to 55% for most subethnicities, will rise by up to 1% per year over the next decade. Subethnicities with historically lower rates of interethnic admixture are less likely to find matches outside of their designated group but will benefit from expansion of the registry, whereas ethnically directed drives will enhance matching rates for currently underrepresented subethnicities. Donor searches for the same cohort using a large extramural registry was of only slight benefit for most of the 19 EM BMDR subethnicities evaluated, confirming that local donor registries that reflect the ethnic diversity of the community being served are best equipped to serve the needs of their respective communities. Contemporary trends of an increasingly multiethnic admixture in Israel may impact the effect of ethnic profiling in assessing future match rates for EM BMDR.