HLA class I haplotype diversity is consistent with selection for frequent existing haplotypes

Frequency distribution of HLA class I and II alleles in Greek population and their significance in orchestrating the National Donor Registry Program

Human leukocyte antigens (HLA) represent one of the most polymorphic systems in humans, responsible for the identification of foreign antigens and the presentation of immune responses. Therefore, HLA is considered to play a major role in human disorders, donor-recipient matching and transplantation outcomes. This study aimed to determine the HLA class I and II alleles and haplotypes in the Greek population. Moreover, a comparative analysis of HLA alleles and haplotype frequencies found in Greek and pooled European populations was also performed to acquire a better knowledge about the HLA alleles distribution. A total number of 1896 healthy individuals were typed for their HLA alleles in the National Tissue Typing Center of Greece. High-resolution HLA typing for the HLA-A, -B, -C and -DR, -DQ, -DP with the use of the next-generation sequencing analysis was performed, followed by data analysis for establishing the HLA allele and haplotype differences. The results of this study showed that the most frequent alleles for the HLA-A were the A*02:01:01 (27.1%), *24:02:01 (14.4%), *01:01:01 (9.3%), for the HLA-B were the B*51:01:01 (15.3%), *18:01:01 (9.7%), *35:01:01 (6.8%) and for the HLA-C were the C*04:01:01 (15.4%), *07:01:01 (13.1%), *12:03:01 (9.6%). For the HLA class II, the most frequent alleles for the HLA-DRB1 were the DRB1*11:04:01 (16.4%), *16:01:01 (11.3%), *11:01:01 (9.5%), for the HLA-DQB1 were the DQB1*03:01:01 (30.5%), *05:02:01 (15.1%), *05:01:01 (10.6%) and for the HLA-DPB1 were the DPB1*04:01:01 (34.8%), *02:01:01 (11.6%), *04:02:01 (7.3%). Additionally, the most frequent haplotypes were the A*02:01:01∼C*07:01:01-B*18:01:01∼DRB1*11:04:01 (2.3%), followed by the A*01:01:01∼C*07:01:01∼B*08:01:01∼DRB1*03:01:01 (2.2%), A*24:02:01∼C*04:01:01∼B*35:02:01∼DRB1*11:04:01 (1.4%) and A*02:01:01∼C*04:01:01∼B*35:01:01-DRB1*14:01:01 (1.2%). The results herein were comparable to those obtained from the pooled European populations. Moreover, these results can be used for the improvement of the donor-recipient matching procedure and to understand better the disease association in Greece.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Graph-Based Imputation Methods and Their Applications to Single Donors and Families

The outcome of Hematopoietic Stem Cell (HSCT) and organ transplant is strongly affected by the matching of the HLA alleles of the donor and the recipient. However, donors and sometimes recipients are often typed at low resolution, with some alleles either missing or ambiguous. Thus, imputation methods are required to detect the most probably high-resolution HLA haplotypes consistent with a typing. Such imputation algorithms require predefined haplotype frequencies. As such, the phasing of the typing is required for both imputation and frequency generation.We have developed a new approach to HLA haplotype and genotype imputation, where first all candidate phases of a typing are explicated, and then the ambiguity within each phase is solved. This ambiguity is solved through a graph structure of all partial haplotypes and the haplotypes consistent with them.This phasing approach was used to produce an imputation algorithm (GRIMM-Graph Imputation and Matching). GRIMM was then combined with the possibility of combining information from multiple races to produce MR-GRIMM (Multi-Race GRIMM). When family information is available, the phasing of each family member can be restricted by the others. We propose GRAMM (GRaph-bAsed faMily iMputation) to phase alleles in family pedigree HLA typing data and in mother-cord blood unit pairs. Finally, we combined MR-GRIMM with an expectation-maximization (EM) algorithm to estimate haplotype frequencies sharing information between races to produce MR-GRIMME (MR-GRIMM EM).We have shown that these algorithms naturally combine information between races and family members. The accuracy of each of these algorithms is significantly better than its current parallel methods. MR-GRIMM leads to high accuracy in matching predictions. GRAMM better imputes family members than either MR-GRIMM or any existing algorithm and has practically no phasing errors. MR-GRIMME obtains a higher likelihood than existing algorithms.MR-GRIMM, MR-GRIMME, and GRAMM are available as servers or through stand-alone versions in GITHUB and PyPi, as detailed in the appropriate sections.

Combined imputation of HLA genotype and self-identified race leads to better donor-recipient matching

Allogeneic Hematopoietic Cell Transplantation (HCT) is a curative therapy for hematologic disorders and often requires human leukocyte antigen (HLA)-matched donors. Donor registries have recruited donors utilizing evolving technologies of HLA genotyping methods. This necessitates in-silico ambiguity resolution and statistical imputation based on haplotype frequencies estimated from donor data stratified by self-identified race and ethnicity (SIRE). However, SIRE has limited genetic validity and presents a challenge for individuals with unknown or mixed SIRE. We present MR-GRIMM "Multi-Race Graph IMputation and Matching" that simultaneously imputes the race/ethnic category and HLA genotype using a SIRE based prior. Additionally, we propose a novel method to impute HLA typing inconsistent with current haplotype frequencies. The performance of MR-GRIMM was validated using a dataset of 170,000 donor-recipient pairs. MR-GRIMM has an average 20 % lower matching error (1-AUC) than single-race imputation. The recall metric (sensitivity) of the race/ethnic category imputation from HLA was measured by comparing the imputed donor race with the donor-provided SIRE. Accuracies of 0.74 and 0.55 were obtained for the prediction of 5 broad and 21 detailed US population groups respectively. The operational implementation of this algorithm in a registry search could help improve match predictions and access to HLA-matched donors.

Bw4 ligand and direct T-cell receptor binding induced selection on HLA A and B alleles

Introduction

The HLA region is the hallmark of balancing selection, argued to be driven by the pressure to present a wide variety of viral epitopes. As such selection on the peptide-binding positions has been proposed to drive HLA population genetics. MHC molecules also directly binds to the T-Cell Receptor and killer cell immunoglobulin-like receptors (KIR).

Methods

We here combine the HLA allele frequencies in over six-million Hematopoietic Stem Cells (HSC) donors with a novel machine-learning-based method to predict allele frequency.

Results

We show for the first time that allele frequency can be predicted from their sequences. This prediction yields a natural measure for selection. The strongest selection is affecting KIR binding regions, followed by the peptide-binding cleft. The selection from the direct interaction with the KIR and TCR is centered on positively charged residues (mainly Arginine), and some positions in the peptide-binding cleft are not associated with the allele frequency, especially Tyrosine residues.

Discussion

These results suggest that the balancing selection for peptide presentation is combined with a positive selection for KIR and TCR binding.

Genetic swamping of the critically endangered Scottish wildcat was recent and accelerated by disease

The European wildcat population in Scotland is considered critically endangered as a result of hybridization with introduced domestic cats,1,2 though the time frame over which this gene flow has taken place is unknown. Here, using genome data from modern, museum, and ancient samples, we reconstructed the trajectory and dated the decline of the local wildcat population from viable to severely hybridized. We demonstrate that although domestic cats have been present in Britain for over 2,000 years,3 the onset of hybridization was only within the last 70 years. Our analyses reveal that the domestic ancestry present in modern wildcats is markedly over-represented in many parts of the genome, including the major histocompatibility complex (MHC). We hypothesize that introgression provides wildcats with protection against diseases harbored and introduced by domestic cats, and that this selection contributes to maladaptive genetic swamping through linkage drag. Using the case of the Scottish wildcat, we demonstrate the importance of local ancestry estimates to both understand the impacts of hybridization in wild populations and support conservation efforts to mitigate the consequences of anthropogenic and environmental change.

GRAMM: A new method for analysis of HLA in families

Recently, haplo-identical transplantation with multiple HLA mismatches has become a viable option for stem cell transplants. Haplotype sharing detection requires the imputation of donor and recipient. We show that even in high-resolution typing when all alleles are known, there is a 15% error rate in haplotype phasing, and even more in low-resolution typings. Similarly, in related donors, the parents' haplotypes should be imputed to determine what haplotype each child inherited. We propose graph-based family imputation (GRAMM) to phase alleles in family pedigree HLA typing data, and in mother-cord blood unit pairs. We show that GRAMM has practically no phasing errors when pedigree data are available. We apply GRAMM to simulations with different typing resolutions as well as paired cord-mother typings, and show very high phasing accuracy, and improved allele imputation accuracy. We use GRAMM to detect recombination events and show that the rate of falsely detected recombination events (false-positive rate) in simulations is very low. We then apply recombination detection to typed families to estimate the recombination rate in Israeli and Australian population datasets. The estimated recombination rate has an upper bound of 10%-20% per family (1%-4% per individual).

Expanding the MAPPs Assay to Accommodate MHC-II Pan Receptors for Improved Predictability of Potential T Cell Epitopes

A critical step in the immunogenicity cascade is attributed to human leukocyte antigen (HLA) II presentation triggering T cell immune responses. The liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based major histocompatibility complex (MHC) II-associated peptide proteomics (MAPPs) assay is implemented during preclinical risk assessments to identify biotherapeutic-derived T cell epitopes. Although studies indicate that HLA-DP and HLA-DQ alleles are linked to immunogenicity, most MAPPs studies are restricted to using HLA-DR as the dominant HLA II genotype due to the lack of well-characterized immunoprecipitating antibodies. Here, we address this issue by testing various commercially available clones of MHC-II pan (CR3/43, WR18, and Tü39), HLA-DP (B7/21), and HLA-DQ (SPV-L3 and 1a3) antibodies in the MAPPs assay, and characterizing identified peptides according to binding specificity. Our results reveal that HLA II receptor-precipitating reagents with similar reported specificities differ based on clonality and that MHC-II pan antibodies do not entirely exhibit pan-specific tendencies. Since no individual antibody clone is able to recover the complete HLA II peptide repertoire, we recommend a mixed strategy of clones L243, WR18, and SPV-L3 in a single immunoprecipitation step for more robust compound-specific peptide detection. Ultimately, our optimized MAPPs strategy improves the predictability and additional identification of T cell epitopes in immunogenicity risk assessments.

COVID-19 annual update: a narrative review

Three and a half years after the pandemic outbreak, now that WHO has formally declared that the emergency is over, COVID-19 is still a significant global issue. Here, we focus on recent developments in genetic and genomic research on COVID-19, and we give an outlook on state-of-the-art therapeutical approaches, as the pandemic is gradually transitioning to an endemic situation. The sequencing and characterization of rare alleles in different populations has made it possible to identify numerous genes that affect either susceptibility to COVID-19 or the severity of the disease. These findings provide a beginning to new avenues and pan-ethnic therapeutic approaches, as well as to potential genetic screening protocols. The causative virus, SARS-CoV-2, is still in the spotlight, but novel threatening virus could appear anywhere at any time. Therefore, continued vigilance and further research is warranted. We also note emphatically that to prevent future pandemics and other world-wide health crises, it is imperative to capitalize on what we have learnt from COVID-19: specifically, regarding its origins, the world's response, and insufficient preparedness. This requires unprecedented international collaboration and timely data sharing for the coordination of effective response and the rapid implementation of containment measures.

Allele and haplotype frequencies of 11 HLA loci in Koreans by next-generation sequencing

Data on HLA genotype distribution, including DQA1 and DPA1, in the Korean population are limited. We aimed to investigate the allele and haplotype frequencies of 11 HLA loci in 339 Korean subjects using next-generation sequencing (NGS)-based HLA typing. A total of 339 samples from unrelated healthy subjects were genotyped for HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQB1, -DQA1, -DPB1, and -DPA1 using two different NGS-based HLA typing kits (166 tested using the NGSgo-MX11-3 kit [GenDx, Netherlands] and 173 by the AllType NGS 11 Loci Amplification kit [One Lambda, USA]). PyPop software was used to estimate allele and haplotype frequencies and linkage disequilibrium between the loci. Additionally, a principal component analysis was performed to compare the allele distribution of Koreans with that of other populations. A total of 214 HLA alleles (97 class I and 117 class II alleles) were assigned. The most frequent alleles for each locus were A*24:02:01 (24.78%), B*15:01:01 (10.18%), C*01:02:01 (18.44%), DRB1*04:05:01 (9.59%), DRB3*02:02:01 (13.72%), DRB4*01:03:01 (25.81%), DRB5*01:01:01 (9.0%), DQA1*01:02:01 (16.96%), DQB1*03:01:01 (14.31%), DPA1*01:03:01 (44.4%), and DPB1*05:01:01 (35.1%), respectively. The most frequent haplotypes were A*33:03:01-C*03:02:02-B*58:01:01 for HLA class I (5.01%) and DRB1*04:05:01-DQA1*03:03:01-DQB1*04:01:01-DPA1*02:02:02-DPB1*05:01:01 for HLA class II (6.23%). The total allelic ambiguities by NGS were estimated to be minimal and considerably decreased compared with those by Sanger sequencing. The Japanese population had the most similar allele distribution to Koreans, followed by the Chinese population. Frequency data of 11 HLA loci in Koreans can provide essential data for population genetics and disease association studies.

Type 1 diabetes in diverse ancestries and the use of genetic risk scores

Over 75 genetic loci within and outside of the HLA region influence type 1 diabetes risk. Genetic risk scores (GRS), which facilitate the integration of complex genetic information, have been developed in type 1 diabetes and incorporated into models and algorithms for classification, prognosis, and prediction of disease and response to preventive and therapeutic interventions. However, the development and validation of GRS across different ancestries is still emerging, as is knowledge on type 1 diabetes genetics in populations of diverse genetic ancestries. In this Review, we provide a summary of the current evidence on the evolutionary genetic variation in type 1 diabetes and the racial and ethnic differences in type 1 diabetes epidemiology, clinical characteristics, and preclinical course. We also discuss the influence of genetics on type 1 diabetes with differences across ancestries and the development and validation of GRS in various populations.

High fraction of silent recombination in a finite-population two-locus neutral birth-death-mutation model

A precise estimate of allele and haplotype polymorphism is of great interest in theoretical population genetics, but also has practical applications, such as bone marrow registries management. Allele polymorphism is driven mainly by point mutations, while haplotype polymorphism is also affected by recombination. Current estimates treat recombination as mutations in an infinite site model. We here show that even in the simple case of two loci in a haploid individual, for a finite population, most recombination events produce existing haplotypes, and as such are silent. Silent recombination considerably reduces the total number of haplotypes expected from the infinite site model for populations that are not much larger than one over the mutation rate. Moreover, in contrast with mutations, the number of haplotypes does not grow linearly with the population size. We hence propose a more accurate estimate of the total number of haplotypes that takes into account silent recombination. We study large-scale human leukocyte antigen (HLA) haplotype frequencies from human populations to show that the current estimated recombination rate in the HLA region is underestimated.

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.

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.

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.

The changing landscape of HLA typing: Understanding how and when HLA typing data can be used with confidence from bench to bedside

Human leukocyte antigen (HLA) genes are extraordinary for their extreme diversity and widespread impact on human health and disease. More than 30,000 HLA alleles have been officially named and more alleles continue to be discovered at a rapid pace. HLA typing systems which have been developed to detect HLA diversity have advanced rapidly and are revolutionizing our understanding of HLA's clinical importance. However, continuous improvements in knowledge and technology have created challenges for clinicians and scientists. This review explains how differences in HLA typing systems can impact the HLA types that are assigned. The consequences of differences in laboratory testing methods and reference databases are described. The challenges of using HLA types that are not equivalent are illustrated. A fundamental understanding of the continual expansion of our understanding of HLA diversity and limitations in some of the typing data is essential for using typing data appropriately in clinical and research settings.

Comprehensive Analysis of CD4+ T Cell Responses to CMV pp65 Antigen Restricted by Single HLA-DR, -DQ, and -DP Allotype Within an Individual

Within an individual, six different HLA class II heterodimers are expressed co-dominantly by two alleles of HLA-DR, -DQ, and -DP loci. However, it remained unclear which HLA allotypes were used in T cell responses to a given antigen. For the measurement of the CD4⁺ T cell responses restricted by a single HLA allotype, we established a panel of artificial antigen-presenting cells (aAPCs) expressing each single HLA allele of 20 HLA-DRB1, 16 HLA-DQ, and 13 HLA-DP alleles. CD4⁺ T cell responses to cytomegalovirus (CMV) pp65 restricted by single HLA class II allotype defined in 45 healthy donors. The average magnitude of CD4⁺ T cell responses by HLA-DR allotypes was higher than HLA-DQ and HLA-DP allotypes. CD4⁺ T cell responses by DRA*01:01/DRB1*04:06, DQA1*01:02/DQB1*06:02, DPA1*02:02/DPB1*05:01 were higher among the other alleles in each HLA-DR, -DQ, and -DP locus. Interestingly, the frequencies of HLA-DR alleles and the positivity of specific allotypes showed an inverse correlation. One allotype within individuals is dominantly used in CD4⁺ T cell response in 49% of donors, and two allotypes showed that in 7% of donors, and any positive response was detected in 44% of donors. Even if one individual had several dominant alleles, CD4⁺ T cell responses tended to be restricted by only one of them. Furthermore, CD8⁺ and CD4⁺ T cell responses by HLA class I and class II were correlated. Our results demonstrate that the CD4⁺ T cell preferentially use a few dominant HLA class II allotypes within individuals, similar to CD8⁺ T cell response to CMV pp65.

Usefulness of the Hematopoietic Stem Cell Donor Pool as a Source of HLA-Homozygous Induced Pluripotent Stem Cells for Haplobanking: Combined Analysis of the Cord Blood Inventory and Bone Marrow Donor Registry

Induced pluripotent stem cells (iPSCs) have opened up unprecedented opportunities for novel therapeutic options for precision medicine. Hematopoietic stem cell (HSC) donor pools with previously determined HLA types may be ideal sources for iPSC production. Based on the HLA distribution of cryopreserved cord blood units (CBUs) and registered bone marrow (BM) donors, we estimated how much of the Korean population could be covered by HLA-homozygous iPSCs. We analyzed a total of 143,866 Korean HSC donors (27,904 CBUs and 115,962 BM donors). Each donor sample was typed for the HLA-A, -B, and -DRB1 alleles at low to intermediate resolution by DNA-based molecular techniques: PCR sequence-specific oligonucleotide (PCR-SSOP), PCR with sequence-specific primers (PCR-SSP) and PCR with sequence-based typing (PCR-SBT). We also identified individuals possessing homozygous HLA haplotypes by direct counting. The matching probabilities for zero-mismatch transplantation were calculated for 143,866 Koreans and 50 million potential Korean patients. Among the HSC donor pool, 17 HLA-A alleles, 41 HLA-B alleles, and 13 HLA-DRB1 alleles, as well as 128 homozygous HLA-A-B-DRB1 haplotypes, were identified at serologic equivalents, and those haplotypes cumulatively matched 93.20% of the 143,866 Korean donors as zero HLA-mismatch iPSC sources. Among the combinations of 2,056 haplotypes with frequencies ≥ 0.001% in a population of 50 million, those 128 homozygous haplotypes can provide 93.65% coverage for potential Korean recipients. Haplobanking of a reasonable number of HLA-A, -B, and -DRB1 homozygous iPSC lines derived from CBUs and cells of registered BM donors may be an efficient option for allogenic iPSC therapy.

Invasion Rate Versus Diversity in Population Dynamics with Catastrophes

A key question in the current diversity crisis is how diversity has been maintained throughout evolution and how to preserve it. Modern coexistence theories suggest that a high invasion rate of rare new types is directly related to diversity. We show that adding almost any mechanism of catastrophes to a stochastic birth, death, and mutation process with limited carrying capacity induces a novel phase transition characterized by a positive invasion rate but a low diversity. In this phase, new types emerge and grow rapidly, but the resulting growth of very large types decreases diversity. This model also resolves two major drawbacks of neutral evolution models: their failure to explain balancing selection without resorting to fitness differences and the unrealistic time required for the creation of the observed large types. We test this model on a classical case of genetic polymorphism: the HLA locus.

HLA class I alleles are associated with clinic-based migraine and increased risks of chronic migraine and medication overuse

OBJECTIVE

We aimed to evaluate associations of human leukocyte antigen variants with migraine or headache in hospital and population-based settings.

METHODS

The case-control study population, aged 30-70, included 605 clinic-based migraine patients in a medical center and 8449 population-based participants in Taiwan Biobank (TWB). Clinic-based cases were ascertained by neurologists. Participants in Taiwan Biobank were interviewed by a structured questionnaire including headache and migraine history; among them, 2394 had headache or migraine history while 6055 were free of headache and served as controls. All subjects were genotyped by Axiom Genome-Wide Single Nucleotide Polymorphism Arrays and imputed for eight classical human leukocyte antigen genes. Human leukocyte antigen frequencies were compared between clinic-based and self-reported patients and controls. We utilized likelihood ratio tests to examine human leukocyte antigen-disease associations and logistic regressions to estimate the effect of human leukocyte antigen alleles on migraine.

RESULTS

Human leukocyte antigen-B and C showed significant associations with clinic-based migraine (q-value < 0.05). Human leukocyte antigen-B*39:01, human leukocyte antigen-B*51:01, human leukocyte antigen-B*58:01 and human leukocyte antigen-C*03:02 were significantly associated with migraine, with age and sex-adjusted odds ratios (95% CIs) of 1.80 (1.28-2.53), 1.50 (1.15-1.97), 1.36 (1.14-1.62) and 1.36 (1.14-1.62), correspondingly. Clinic-based migraineurs carrying human leukocyte antigen-B*58:01 or human leukocyte antigen-C*03:02 had 1.63 (1.11-2.39) -fold likelihood to have chronic migraine with medication-overuse headache compared to episodic migraine. However, no human leukocyte antigen genes were associated with self-reported headache or migraine in the community.

CONCLUSIONS

Human leukocyte antigen class I genetic variants are positively associated with risk of clinic-based migraine but not self-reported migraine or headache and may contribute to migraine chronification and medication overuse.

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.

Next-generation sequencing reveals new information about HLA allele and haplotype diversity in a large European American population

The human leukocyte antigen (HLA) genes are extremely polymorphic and are useful molecular markers to make inferences about human population history. However, the accuracy of the estimation of genetic diversity at HLA loci very much depends on the technology used to characterize HLA alleles; high-resolution genotyping of long-range HLA gene products improves the assessment of HLA population diversity as well as other population parameters compared to lower resolution typing methods. In this study we examined allelic and haplotype HLA diversity in a large healthy European American population sourced from the UCSF-DNA bank. A high-resolution next-generation sequencing method was applied to define non-ambiguous 3- and 4-field alleles at the HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci in samples provided by 2248 unrelated individuals. A number of population parameters were examined including balancing selection and various measurements of linkage disequilibrium were calculated. There were no detectable deviations from Hardy-Weinberg proportions at HLA-A, HLA-DRB1, HLA-DQA1 and HLA-DQB1. For the remaining loci moderate and significant deviations were detected at HLA-C, HLA-B, HLA-DRB3/4/5, HLA-DPA1 and HLA-DPB1 loci mostly from population substructures. Unique 4-field associations were observed among alleles at 2 loci and haplotypes extending large intervals that were not apparent in results obtained using testing methodologies with limited sequence coverage and phasing. The high diversity at HLA-DPA1 results from detection of intron variants of otherwise well conserved protein sequences. It may be speculated that divergence in exon sequences may be negatively selected. Our data provides a valuable reference source for future population studies that may allow for precise fine mapping of coding and non-coding sequences determining disease susceptibility and allo-immunogenicity.

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

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

Multiplicative fitness, rapid haplotype discovery, and fitness decay explain evolution of human MHC

The major histocompatibility complex (MHC) is a central component of the vertebrate immune system and hence evolves in the regime of a host-pathogen evolutionary race. The MHC is associated with quantitative traits which directly affect fitness and are subject to selection pressure. The evolution of haplotypes at the MHC HLA (HLA) locus is generally thought to be governed by selection for increased diversity that is manifested in overdominance and/or negative frequency-dependent selection (FDS). However, recently, a model combining purifying selection on haplotypes and balancing selection on alleles has been proposed. We compare the predictions of several population dynamics models of haplotype frequency evolution to the distributions derived from 6.59-million-donor HLA typings from the National Marrow Donor Program registry. We show that models that combine a multiplicative fitness function, extremely high haplotype discovery rates, and exponential fitness decay over time produce the best fit to the data for most of the analyzed populations. In contrast, overdominance is not supported, and population substructure does not explain the observed haplotype frequencies. Furthermore, there is no evidence of negative FDS. Thus, multiplicative fitness, rapid haplotype discovery, and rapid fitness decay appear to be the major factors shaping the HLA haplotype frequency distribution in the human population.

Long Noncoding RNA HCP5, a Hybrid HLA Class I Endogenous Retroviral Gene: Structure, Expression, and Disease Associations

The HCP5 RNA gene (NCBI ID: 10866) is located centromeric of the HLA-B gene and between the MICA and MICB genes within the major histocompatibility complex (MHC) class I region. It is a human species-specific gene that codes for a long noncoding RNA (lncRNA), composed mostly of an ancient ancestral endogenous antisense 3′ long terminal repeat (LTR, and part of the internal pol antisense sequence of endogenous retrovirus (ERV) type 16 linked to a human leukocyte antigen (HLA) class I promoter and leader sequence at the 5′-end. Since its discovery in 1993, many disease association and gene expression studies have shown that HCP5 is a regulatory lncRNA involved in adaptive and innate immune responses and associated with the promotion of some autoimmune diseases and cancers. The gene sequence acts as a genomic anchor point for binding transcription factors, enhancers, and chromatin remodeling enzymes in the regulation of transcription and chromatin folding. The HCP5 antisense retroviral transcript also interacts with regulatory microRNA and immune and cellular checkpoints in cancers suggesting its potential as a drug target for novel antitumor therapeutics.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

Identification and Characterization of Positively Selected Mutations in Nef of Four HIV-1 Major Subtypes from Los Alamos National Laboratory

BACKGROUND

Human immunodeficiency virus-1 (HIV-1) mutates rapidly to escape host immune pressure. This results in the generation of positively selected mutations (PSM) throughout the viral genome. Escape mutations in Nef, one of the accessory proteins of HIV-1, which plays an important role in viral pathogenicity have previously been identified in several large cohort studies, but the evolution of PSMs overtime in various HIV-1 subtypes remains unknown.

METHODS

161 clade A1, 3093 clade B, 647 clade C and 115 clade D HIV-1 nef sequences were obtained from the HIV Database of Los Alamos National Laboratory and aligned using MEGA 6.0. The sequences from each clade were grouped based on the year of collection. Quasi analysis was used to identify PSMs and the number and locations of PSMs were compared among different subtypes.

RESULTS

PSMs for all four subtypes were distributed across the sequence of Nef, and conserved residues F90, W113, PxxPxR (a.a 72-77) remain unaltered overtime. The frequency of PSMs was stable among subtype B sequences but increased overtime for other subtypes. Phylogenetic analysis shows that sequences containing PSMs tend to cluster together at both inter and intra- subtype levels.

CONCLUSION

Identification of PSMs and their changes overtime within various subtypes of HIV-1 is important in defining global viral evolutionary patterns that can provide insights for designing therapeutic strategies.

Transplantation of photoreceptors into the degenerative retina: Current state and future perspectives

The mammalian retina displays no intrinsic regenerative capacities, therefore retinal degenerative diseases such as age-related macular degeneration (AMD) or retinitis pigmentosa (RP) result in a permanent loss of the light-sensing photoreceptor cells. The degeneration of photoreceptors leads to vision impairment and, in later stages, complete blindness. Several therapeutic strategies have been developed to slow down or prevent further retinal degeneration, however a definitive cure i.e. replacement of the lost photoreceptors, has not yet been established. Cell-based treatment approaches, by means of photoreceptor transplantation, have been studied in pre-clinical animal models over the last three decades. The introduction of pluripotent stem cell-derived retinal organoids represents, in principle, an unlimited source for the generation of transplantable human photoreceptors. However, safety, immunological and reproducibility-related issues regarding the use of such cells still need to be solved. Moreover, the recent finding of cytoplasmic material transfer between donor and host photoreceptors demands reinterpretation of several former transplantation studies. At the same time, material transfer between healthy donor and dysfunctional patient photoreceptors also offers a potential alternative strategy for therapeutic intervention. In this review we discuss the history and current state of photoreceptor transplantation, the techniques used to assess rescue of visual function, the prerequisites for effective transplantation as well as the main roadblocks, including safety and immune response to the graft, that need to be overcome for successful clinical translation of photoreceptor transplantation approaches.