Computer simulation of human leukocyte antigen genes supports two main routes of colonization by human populations in East Asia

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

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

Comprehensive research on the distribution of HLA-DRB1 in Chinese populations

By presenting antigen peptides, HLA-DRB1 plays an important role in the immune system. However, the allele frequency of HLA-DRB1 exon 2 across China has not been comprehensively studied, especially in minority populations. We sampled 3757 individuals from 59 population. The HLA-DRB1 region from 212 to 463 bp (NM_002124.4 exon 2) in each population was sequenced by Sanger sequencing and genotyped via SBTengine® software, and the allele frequency was calculated by GenAlEx 6.5. Eighty-two DRB1 alleles were identified. The expected heterozygosity of DRB1 was lower in the south than in the north, which was inconsistent with the Y chromosome and mitochondrial DNA results. The Mantel test and nonparametric correlation analysis showed that the correlations of the genetic distance with geographical distance and of DRB1 allele frequencies with latitude weakened after the southern and northern groups were considered separately. Principal coordinate analysis showed that populations speaking the same languages were not codistributed. Compared with other genetic markers, the distribution of DRB1 seems less affected by geographic distance and ethnic origin. Local factors such as gene flow with neighbouring populations, geographic isolation or natural selection are important forces shaping the DRB1 gene pool of local populations.

Challenging ancient DNA results about putative HLA protection or susceptibility to Yersina pestis

In a recent article, Immel et al. (Immel A, Key FM, Szolek A, Barquera R, Robinson MK, Harrison GF, Palmer WH, Spyrou MA, Susat J, Krause-Kyora B, et al. 2021. Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes. Mol Biol Evol. 38:4059–4076) extracted DNA from 36 individuals dead from plague in Ellwangen, Southern Germany, during the 16th century. By comparing their human leukocyte antigen (HLA) genotypes with those of 50 present-day Ellwangen inhabitants, the authors reported a significant decrease of HLA-B*51:01 and HLA-C*06:02 and a significant increase of HLA-DRB1*13:01/13:02 frequencies from ancient to modern populations. After comparing these frequencies with a larger sample of 8,862 modern Germans and performing simulations of natural selection, they concluded that these changes had been driven by natural selection. In an attempt to provide more evidence on such stimulating results, we explored the HLA frequency patterns over all of Europe, we predicted binding affinities of HLA-B/C/DRB1 alleles to 106,515 Yersinia pestis-derived peptides, and we performed forward simulations of HLA genetic profiles under neutrality. Our analyses do not sustain the conclusions of HLA protection or susceptibility to plague based on ancient DNA.

Spatially explicit paleogenomic simulations support cohabitation with limited admixture between Bronze Age Central European populations

The Bronze Age is a complex period of social, cultural and economic changes. Recent paleogenomic studies have documented a large and rapid genetic change in early Bronze Age populations from Central Europe. However, the detailed demographic and genetic processes involved in this change are still debated. Here we have used spatially explicit simulations of genomic components to better characterize the demographic and migratory conditions that may have led to this change. We investigated various scenarios representing the expansion of pastoralists from the Pontic steppe, potentially linked to the Yamnaya cultural complex, and their interactions with local populations in Central Europe, considering various eco-evolutionary factors, such as population admixture, competition and long-distance dispersal. Our results do not support direct competition but rather the cohabitation of pastoralists and farmers in Central Europe, with limited gene flow between populations. They also suggest occasional long-distance migrations accompanying the expansion of pastoralists and a demographic decline in both populations following their initial contact. These results link recent archaeological and paleogenomic observations and move further the debate of genomic changes during the early Bronze Age.

Molecular epidemiologic study of citrin deficiency by screening for four reported pathogenic SLC25A13 variants in the Shaanxi and Guangdong provinces, China

BACKGROUND

Citrin deficiency (CD) is an autosomal recessive disease resulting from biallelic mutations of the SLC25A13 gene. This study aimed to investigate the molecular epidemiological features of CD in the Guangdong and Shaanxi provinces of China.

METHODS

A total of 3,409 peripheral blood samples from Guangdong and 2,746 such samples from Shaanxi province were collected. Four prevalent SLC25A13 mutations NG_012247.2 (NM_014251.3): c.852_855del, c.1638_1660dup, c.615+5G>A, and c.1751-5_1751-4ins(2684) were screened by using the conventional polymerase chain reaction (PCR)/PCR-restriction fragment length polymorphism and newly-developed multiplex PCR methods, respectively. The mutated SLC25A13 allele frequencies, carrier frequencies, and CD morbidity rates were calculated and then compared with the Chi-square and Fisher's exact tests.

RESULTS

The mutations were detected in 68 out of 6,818 SLC25A13 alleles in Guangdong and 29 out of 5,492 alleles in the Shaanxi population. The carrier frequencies were subsequently calculated to be 1/51 and 1/95, while the CD morbidity rates were 1/10,053 and 1/35,865, in the 2 populations, respectively. When compared with the Shaanxi population, Guangdong exhibited a higher frequency of mutated SLC25A13 allele (68/6,818 vs. 29/5,492, χ²=8.570, P=0.003) in general, with higher c.852_855del (54/6,818 vs. 13/5,492, χ²=17.328, P=0.000) but lower c.1751-5_1751 -4ins(2684) (2/6,818 vs. 9/5,492, P=0.015) allele frequencies. The distribution of c.615+5G>A and c.1638_1660dup between the 2 provinces, as well as all 4 prevalent mutations among different geographic regions within the 2 provinces, did not differed significantly.

CONCLUSIONS

Our findings depicted the CD molecular epidemiological features in Guangdong and Shaanxi populations, providing preliminary but significant laboratory evidences for the subsequent CD diagnosis and management in the 2 provinces of mainland China.

Influence of Paleolithic range contraction, admixture and long-distance dispersal on genetic gradients of modern humans in Asia

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long-distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

Deciphering the fine nucleotide diversity of full HLA class I and class II genes in a well-documented population from sub-Saharan Africa

With the aim to understand how next‐generation sequencing (NGS) improves both our assessment of genetic variation within populations and our knowledge on HLA molecular evolution, we sequenced and analysed 8 HLA loci in a well‐documented population from sub‐Saharan Africa (Mandenka). The results of full‐gene NGS‐MiSeq sequencing compared with those obtained by traditional typing techniques or limited sequencing strategies showed that segregating sites located outside exon 2 are crucial to describe not only class I but also class II population diversity. A comprehensive analysis of exons 2, 3, 4 and 5 nucleotide diversity at the 8 HLA loci revealed remarkable differences among these gene regions, notably a greater variation concentrated in the antigen recognition sites of class I exons 3 and some class II exons 2, likely associated with their peptide‐presentation function, a lower diversity of HLA‐C exon 3, possibly related to its role as a KIR ligand, and a peculiar molecular diversity of HLA‐A exon 2, revealing demographic signals. Based on full‐length HLA sequences, we also propose that the most frequent DRB1 allele in the studied population, DRB1*13:04, emerged from an allelic conversion involving 3 potential alleles as donors and DRB1*11:02:01 as recipient. Finally, our analysis revealed a high occurrence of the DRB1*13:04‐DQA1*05:05:01‐DQB1*03:19 haplotype, possibly resulting from a selective sweep due to protection to Onchorcerca volvulus, a prevalent pathogen in West Africa. This study unveils highly relevant information on the molecular evolution of HLA genes in relation to their immune function, calling for similar analyses in other populations living in contrasting environments.

Positive selection on schizophrenia-associated ST8SIA2 gene in post-glacial Asia

A number of loci are associated with highly heritable schizophrenia and the prevalence of this mental illness has had considerable negative fitness effects on human populations. Here we focused on one particular schizophrenia-associated gene that encodes a sialyltransferase (ST8SIA2) and is expressed preferentially in the brain with the level being largely determined by three SNPs in the promoter region. It is suggested that the expression level of the ST8SIA2 gene is a genetic determinant of schizophrenia risk, and we found that a geographically differentiated non-risk SNP type (CGC-type) has significantly reduced promoter activity. A newly developed method for detecting ongoing positive selection was applied to the ST8SIA2 genomic region with the identification of an unambiguous sweep signal in a rather restricted region of 18 kb length surrounding the promoter. We also found that while the CGC-type emerged in anatomically modern humans in Africa over 100 thousand years ago, it has increased its frequency in Asia only during the past 20-30 thousand years. These findings support that the positive selection is driven by psychosocial stress due to changing social environments since around the last glacial maximum, and raise a possibility that schizophrenia extensively emerged during the Upper Paleolithic and Neolithic era.

Discerning the Origins of the Negritos, First Sundaland People: Deep Divergence and Archaic Admixture

Human presence in Southeast Asia dates back to at least 40,000 years ago, when the current islands formed a continental shelf called Sundaland. In the Philippine Islands, Peninsular Malaysia, and Andaman Islands, there exist indigenous groups collectively called Negritos whose ancestry can be traced to the "First Sundaland People." To understand the relationship between these Negrito groups and their demographic histories, we generated genome-wide single nucleotide polymorphism data in the Philippine Negritos and compared them with existing data from other populations. Phylogenetic tree analyses show that Negritos are basal to other East and Southeast Asians, and that they diverged from West Eurasians at least 38,000 years ago. We also found relatively high traces of Denisovan admixture in the Philippine Negritos, but not in the Malaysian and Andamanese groups, suggesting independent introgression and/or parallel losses involving Denisovan introgressed regions. Shared genetic loci between all three Negrito groups could be related to skin pigmentation, height, facial morphology and malarial resistance. These results show the unique status of Negrito groups as descended from the First Sundaland People.

The HLA-B landscape of Africa: Signatures of pathogen-driven selection and molecular identification of candidate alleles to malaria protection

Human leukocyte antigen (HLA) genes play a key role in the immune response to infectious diseases, some of which are highly prevalent in specific environments, like malaria in sub-Saharan Africa. Former case-control studies showed that one particular HLA-B allele, B*53, was associated with malaria protection in Gambia, but this hypothesis was not tested so far within a population genetics framework. In this study, our objective was to assess whether pathogen-driven selection associated with malaria contributed to shape the HLA-B genetic landscape of Africa. To that aim, we first typed the HLA-A and -B loci in 484 individuals from 11 populations living in different environments across the Sahel, and we analysed these data together with those available for 29 other populations using several approaches including linear modelling on various genetic, geographic and environmental parameters. In addition to relevant signatures of populations' demography and migrations history in the genetic differentiation patterns of both HLA-A and -B loci, we found that the frequencies of three HLA alleles, B*53, B*78 and A*74, were significantly associated with Plasmodium falciparum malaria prevalence, suggesting their increase through pathogen-driven selection in malaria-endemic environments. The two HLA-B alleles were further identified, by high-throughput sequencing, as B*53:01:01 (in putative linkage disequilibrium with one HLA-C allele, C*04:01:01:01) and B*78:01 in all but one individuals tested, making them appropriate candidates to malaria protection. These results highlight the role of environmental factors in the evolution of the HLA polymorphism and open key perspectives for functional studies focusing on HLA peptide-binding properties.

High-resolution HLA-A, -B and -DRB1 allele and haplotype frequencies in 7823 Han marrow donors of Liaoning province, China

BACKGROUND

The human leukocyte antigen (HLA) system is the most polymorphic gene cluster in humans. High-resolution donor-recipient matching for HLA genes improves patient survival after unrelated hematopoietic stem cell transplantation.

MATERIALS AND METHODS

In this study, we analyzed the high-resolution allele and haplotype frequencies at the HLA-A, -B and -DRB1 loci in the Liaoning Han population and analyzed its relationships with other populations.

RESULTS

The 3 most frequent alleles at the HLA-A, -B and -DRB1 loci were A*24:02, A*02:01:01G, A*11:01; B*13:02, B*46:01, B*40:01:01G; DRB1*09:01, DRB1*15:01 and DRB1*07:01, respectively. The most frequent 2-locus haplotypes were A*30:01-B*13:02 and B*13:02-DRB1*07:01. A*30:01-B*13:02-DRB1*07:01 was determined to be the predominant 3-locus haplotype. Hot maps and multiple correspondence analyses based on the frequencies of HLA specificities, which allow statistical visualization of dependent and independent relationships among variables, indicate that the Liaoning Han population is closely related to Northern populations of China and shows relative close relationships with Asian populations.

CONCLUSION

These data will provide an outline of the HLA characteristics of healthy individuals in our region and help bone marrow transplantation patients find suitable HLA-matched donors.

Molecular diagnosis of pediatric patients with citrin deficiency in China: SLC25A13 mutation spectrum and the geographic distribution

Citrin deficiency (CD) is a Mendelian disease due to biallelic mutations of SLC25A13 gene. Neonatal intrahepatic cholestasis caused by citrin deficiency (NICCD) is the major pediatric CD phenotype, and its definite diagnosis relies on SLC25A13 genetic analysis. China is a vast country with a huge population, but the SLC25A13 genotypic features of CD patients in our country remains far from being well clarified. Via sophisticated molecular analysis, this study diagnosed 154 new CD patients in mainland China and identified 9 novel deleterious SLC25A13 mutations, i.e. c.103A > G, [c.329 - 154_c.468 + 2352del2646; c.468 + 2392_c.468 + 2393ins23], c.493C > T, c.755 - 1G > C, c.845_c.848 + 1delG, c.933_c.933 + 1insGCAG, c.1381G > T, c.1452 + 1G > A and c.1706_1707delTA. Among the 274 CD patients diagnosed by our group thus far, 41 SLC25A13 mutations/variations were detected. The 7 mutations c.775C > T, c.851_854del4, c.1078C > T, IVS11 + 1G > A, c.1364G > T, c.1399C > T and IVS16ins3kb demonstrated significantly different geographic distribution. Among the total 53 identified genotypes, only c.851_854del4/c.851_854del4 and c.851_854del4/c.1399C > T presented different geographic distribution. The northern population had a higher level of SLC25A13 allelic heterogeneity than those in the south. These findings enriched the SLC25A13 mutation spectrum and brought new insights into the geographic distribution of the variations and genotypes, providing reliable evidences for NICCD definite diagnosis and for the determination of relevant molecular targets in different Chinese areas.

Haplotype analysis of the HFE gene among populations of Northern Eurasia, in patients with metabolic disorders or stomach cancer, and in long-lived people

BACKGROUND

Previously, it was shown that the HFE gene (associated with human hereditary hemochromatosis) has several haplotypes of intronic polymorphisms. Some haplotype frequencies are race specific and hence can be used in phylogenetic analysis. We assumed that analysis of Caucasoid patients-living now in Western Siberia and having diseases associated with dietary habits and metabolic rate-will allow us to understand the processes of possible selection during settling of the northern part of Asia.

RESULTS

Haplotype analysis of Northern Eurasian native and recently settled ethnic groups was performed on polymorphisms rs1799945, rs1800730, rs1800562, rs2071303, rs1800708, rs1572982, rs2794719, rs807209, and rs2032451 of this gene. The CCA haplotype of the rs2071303, rs1800708, and rs1572982 was found to be associated with HLA-A2 (39 %) in Asian populations. Haplotype analysis for the rs1799945, rs1800730, rs1800562, rs2071303, rs1800708, and rs1572982 was performed on Russian patients with some metabolic disorders or stomach cancer and among long-lived people. Decreased frequencies of the TTA haplotype (T in rs2071303, T in rs1800708, and A in rs1572982) were observed in the groups of patients with diseases associated with overweight (fatty liver disease, type 2 diabetes mellitus, or metabolic syndrome + arterial hypertension) as compared with the control sample. We detected significant differences in this haplotype's frequency between the patients with type 2 diabetes mellitus and Russian adolescents, elderly citizens, and long-lived people (χ(2) P value = 0.003, 0.010, and 0.015, respectively).

CONCLUSIONS

No significant differences in frequencies of the alleles with mutations in coding regions of the HFE gene (C282Y, H63D, and S65C) were detected between the analyzed patients (with stomach cancer, metabolic syndrome, fatty liver disease, or type 2 diabetes mellitus) and the control Caucasoid sample. Monophyletic origin of H63D (rs1799945) was confirmed in Caucasoids and Northern Asians. The reasons for a sharp increase in the frequency of CCA haplotype of HFE in the Asian race remain unclear.

HLA class I molecular variation and peptide-binding properties suggest a model of joint divergent asymmetric selection

The main function of HLA class I molecules is to present pathogen-derived peptides to cytotoxic T lymphocytes. This function is assumed to drive the maintenance of an extraordinary amount of polymorphism at each HLA locus, providing an immune advantage to heterozygote individuals capable to present larger repertories of peptides than homozygotes. This seems contradictory, however, with a reduced diversity at individual HLA loci exhibited by some isolated populations. This study shows that the level of functional diversity predicted for the two HLA-A and HLA-B genes considered simultaneously is similar (almost invariant) between 46 human populations, even when a reduced diversity exists at each locus. We thus propose that HLA-A and HLA-B evolved through a model of joint divergent asymmetric selection conferring all populations an equivalent immune potential. The distinct pattern observed for HLA-C is explained by its functional evolution towards killer cell immunoglobulin-like receptor (KIR) activity regulation rather than peptide presentation.

Forward-in-Time, Spatially Explicit Modeling Software to Simulate Genetic Lineages Under Selection

SELECTOR is a software package for studying the evolution of multiallelic genes under balancing or positive selection while simulating complex evolutionary scenarios that integrate demographic growth and migration in a spatially explicit population framework. Parameters can be varied both in space and time to account for geographical, environmental, and cultural heterogeneity. SELECTOR can be used within an approximate Bayesian computation estimation framework. We first describe the principles of SELECTOR and validate the algorithms by comparing its outputs for simple models with theoretical expectations. Then, we show how it can be used to investigate genetic differentiation of loci under balancing selection in interconnected demes with spatially heterogeneous gene flow. We identify situations in which balancing selection reduces genetic differentiation between population groups compared with neutrality and explain conflicting outcomes observed for human leukocyte antigen loci. These results and three previously published applications demonstrate that SELECTOR is efficient and robust for building insight into human settlement history and evolution.