Heterogeneity of dN/dS Ratios at the Classical HLA Class I Genes over Divergence Time and Across the Allelic Phylogeny

Evidence for Pathogen-Driven Selection Acting on HLA-DPB1 in Response to Plasmodium falciparum Malaria in West Africa

African populations remain underrepresented in studies of human genetic diversity, despite a growing interest in understanding how they have adapted to the diverse environments they live in. In particular, understanding the genetic basis of immune adaptation to pathogens is of paramount importance in a continent such as Africa, where the burden of infectious diseases is a major public health challenge. In this study, we investigated the molecular variation of four Human Leukocyte Antigens (HLA) class II genes (DRB1, DQA1, DQB1 and DPB1), directly involved in the immune response to parasitic infections, in more than 1000 individuals from 23 populations across North, East, Central and West Africa. By analyzing the HLA molecular diversity of these populations in relation to various geographical, cultural and environmental factors, we identified divergent genetic profiles for several (semi-)nomadic populations of the Sahel belt as a signature of their unique demography. In addition, we observed significant genetic structuring supporting both substantial geographic and linguistic differentiations within West Africa. Furthermore, neutrality tests suggest balancing selection has been shaping the diversity of these four HLA class II genes, which is consistent with molecular comparisons between HLA genes and their orthologs in chimpanzees (Patr). However, the most striking observation comes from linear modeling, demonstrating that the prevalence of Plasmodium falciparum, the primary pathogen of malaria in Africa, significantly explains a large proportion of the nucleotide diversity observed at the DPB1 gene. DPB1*01:01, a highly frequent allele in Burkinabé populations, is identified as a potential protective allele against malaria, suggesting that strong pathogen-driven positive selection at this gene has shaped HLA variation in Africa. Additionally, two low-frequency DRB1 alleles, DRB1*08:06 and DRB1*11:02, also show significant associations with P. falciparum prevalence, supporting resistance to malaria is determined by multigenic and/or multiallelic combinations rather than single allele effects.

HLA Genetic Diversity and Chronic Hepatitis B Virus Infection: Effect of Heterozygosity Advantage

This research aims to determine whether HLA heterozygosity confers a protective effect against hepatitis B virus infection by analyzing the relationship between HLA diversity and the risk of hepatitis B virus (HBV) infection. A total of 327 hepatitis B patients were selected and categorized based on their clinical status: 284 patients with chronic HBV infection and 43 patients with HBV-related liver cirrhosis (LC). The control group included 304 healthy individuals. HLA genotyping for 11 loci, including HLA class I and class II, was conducted using next-generation sequencing. The results of this study indicate a statistically significant negative correlation between HLA class II heterozygosity and the risk of HBV infection. Specifically, heterozygosity in HLA-DQB1 (OR = 0.49, 95% CI = 0.31-0.76, p = 0.01277) and HLA-DRB1 (OR = 0.42, 95% CI = 0.24-0.77, p = 0.01855) were significantly associated with protection. Subgroup analysis was conducted to explore the effect of HLA diversity among pathological subtypes (chronic hepatitis B and control group, liver cirrhosis and control group). For liver cirrhosis, compared with the control group, a decreased risk of LC was possibly associated with the heterozygosity of HLA class I locus B (OR = 0.24, 95% CI = 0.09-0.65, p = 0.0591), but this hypothesis was not confirmed by other studies. The diversity of HLA, measured by HLA heterozygosity, was associated with a protective effect against HBV infection.

The Gallus gallus RJF reference genome reveals an MHCY haplotype organized in gene blocks that contain 107 loci including 45 specialized, polymorphic MHC class I loci, 41 C-type lectin-like loci, and other loci amid hundreds of transposable elements

MHCY is a second major histocompatibility complex-like gene region in chickens originally identified by the presence of major histocompatibility complex class I-like and class II-like gene sequences. Up to now, the MHCY gene region has been poorly represented in genomic sequence data. A high density of repetitive sequence and multiple members of several gene families prevented the accurate assembly of short-read sequence data for MHCY. Identified here by single-molecule real-time sequencing sequencing of BAC clones for the Gallus gallus Red Jungle Fowl reference genome are 107 MHCY region genes (45 major histocompatibility complex class I-like, 41 c-type-lectin-like, 8 major histocompatibility complex class IIβ, 8 LENG9-like, 4 zinc finger protein loci, and a single only zinc finger-like locus) located amid hundreds of retroelements within 4 contigs representing the region. Sequences obtained for nearby ribosomal RNA genes have allowed MHCY to be precisely mapped with respect to the nucleolar organizer region. Gene sequences provide insights into the unusual structure of the MHCY class I molecules. The MHCY class I loci are polymorphic and group into 22 types based on predicted amino acid sequences. Some MHCY class I loci are full-length major histocompatibility complex class I genes. Others with altered gene structure are considered gene candidates. The amino acid side chains at many of the polymorphic positions in MHCY class I are directed away rather than into the antigen-binding groove as is typical of peptide-binding major histocompatibility complex class I molecules. Identical and nearly identical blocks of genomic sequence contribute to the observed multiplicity of identical MHCY genes and the large size (>639 kb) of the Red Jungle Fowl MHCY haplotype. Multiple points of hybridization observed in fluorescence in situ hybridization suggest that the Red Jungle Fowl MHCY haplotype is made up of linked, but physically separated genomic segments. The unusual gene content, the evidence of highly similar duplicated segments, and additional evidence of variation in haplotype size distinguish polymorphic MHCY from classical polymorphic major histocompatibility complex regions.

The HLA diversity of the Anthony Nolan register

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

Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

Divergent Allele Advantage Provides a Quantitative Model for Maintaining Alleles with a Wide Range of Intrinsic Merits

The Major Histocompatibility Complex (MHC) is the most genetically diverse region of the genome in most vertebrates. Some form of balancing selection is necessary to account for the extreme diversity, but the precise mechanism of balancing selection is unknown. Due to the way MHC molecules determine immune recognition, overdominance (also referred to as heterozygote advantage) has been suggested as the main driving force behind this unrivalled diversity. However, both theoretical results and simulation models have shown that overdominance in its classical form cannot maintain large numbers of alleles unless all alleles confer unrealistically similar levels of fitness. There is increasing evidence that heterozygotes containing genetically divergent alleles allow for broader antigen presentation to immune cells, providing a selective mechanism for MHC polymorphism. By framing competing models of overdominance within a general framework, we show that a model based on Divergent Allele Advantage (DAA) provides a superior mechanism for maintaining alleles with a wide range of intrinsic merits, as intrinsically less-fit MHC alleles that are more divergent can survive under DAA. Specifically, our results demonstrate that a quantitative mechanism built from the DAA hypothesis is able to maintain polymorphism in the MHC. Applying such a model to both livestock breeding and conservation could provide a better way of identifying superior heterozygotes, and quantifying the advantages of genetic diversity at the MHC.

Semen quality is affected by HLA class I alleles together with sexually transmitted diseases

BACKGROUND

The human leukocyte antigen (HLA) locus includes several genes with key roles in antigen presentation and immune response, some of them inclusively found to be associated with non-obstructive azoospermia. Still, HLA connections to other infertility phenotypes such as semen hyperviscosity (SHV), asthenozoospermia (AST), and oligozoospermia (OLI) have been often neglected.

OBJECTIVES

In this work, we aimed to evaluate the association of HLA class I and II genes with SHV, AST, and OLI phenotypes while exploring a possible role in an adaptive immune response to sexually transmitted diseases (STD).

MATERIALS AND METHODS

Whole-exome sequencing was performed in a Portuguese cohort of 71 infertility cases and 68 controls, followed by HLA typing using a specific software-HLA*PRG:LA tool. Molecular screenings of seven STD were carried out in a subset of 72 samples (30 cases and 42 controls).

RESULTS

Statistical tests uncovered three protective alleles: HLA-A*11:01, associated with all forms of male infertility (p = 0.0006); HLA-DQB1*03:02 with SHV and OLI (P_SHV   = 0.0303, P_OLI   = 0.0153); and HLA-A*29:02 with OLI (p = 0.0355), which was found to interfere in sperm number together with HPV (p = 0.0313). Five risk alleles were also identified: two linked with SHV (HLA-B*50:01, p = 0.0278; and HLA-C*06:02, p = 0.0461), another one with both SHV and OLI (HLA-DQA1*05:01, P_SHV   = 0.0444 and P_OLI =0.0265), and two with OLI (HLA-C*03:03, p = 0.0480; and HLA-DQB1*03:01, p = 0.0499). Here, HLA-C*03:03 carriers tend to be HPV infected.

CONCLUSIONS

The application of HLA*PRG:LA tool to the study of male infertility provided novel insights for an HLA correlation with semen quality, namely among SHV and OLI phenotypes. The discovery of an HLA-A*29:02/HPV crosstalk, together with former reports of HLA alleles conferring resistance-susceptibility to diverse human pathogens, raises the hypothesis of a mechanistic link between male infertility, HLA polymorphism, and host response to STD.

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.

The variations of TRBV genes usages in the peripheral blood of a healthy population are associated with their evolution and single nucleotide polymorphisms

T cell receptors (TCRs) are a class of T cell surface molecules that recognize the antigen-derived peptides presented by the major histocompatibility complex (MHC) and are able to trigger a series of immune responses. TCRs are important members of the adaptive immune system that arose in the jawed fish 500 million years ago. T cell receptor beta variable (TRBV) genes have been widely used to characterize TCR repertoires. Studying the evolution of TRBV may help us to better understand the adaptive immune system. To investigate TRBV evolution and its impacts on the usages of TRBV genes in human populations, we compared the TRBV genes and their homologous sequences among humans, mouse, rhesus and chimpanzee, analyzed the single-nucleotide polymorphisms (SNPs) located at TRBV loci, and sequenced TCR repertoires in the peripheral blood of 97 healthy donors. We found that functional TRBVs are more evolutionarily conserved but possess more SNPs in human populations than do nonfunctional (pseudo) TRBVs. Based on the conservation levels in the four species, we classified the functional TRBVs into 2 groups: old (conserved between mouse and humans) and new (conserved only in primates). The new TRBVs evolve faster and possess more SNPs than the old TRBVs. The variations in TRBV genes frequencies in the peripheral blood of healthy donors are negatively correlated with SNP density. These observations suggest that TRBV usages may be influenced by TCR-MHC co-evolution.

A genomic perspective on HLA evolution

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

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

The above article originally published with an incomplete bibliographic information for Bitarello et al. (2016) and presented correctly in this article.

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

The major histocompatibility complex (MHC) contains the most polymorphic genetic system in humans, the human leukocyte antigen (HLA) genes of the adaptive immune system. High allelic diversity in HLA is argued to be maintained by balancing selection, such as negative frequency-dependent selection or heterozygote advantage. Selective pressure against immune escape by pathogens can maintain appreciable frequencies of many different HLA alleles. The selection pressures operating on combinations of HLA alleles across loci, or haplotypes, have not been extensively evaluated since the high HLA polymorphism necessitates very large sample sizes, which have not been available until recently. We aimed to evaluate the effect of selection operating at the HLA haplotype level by analyzing HLA A~C~B~DRB1~DQB1 haplotype frequencies derived from over six million individuals genotyped by the National Marrow Donor Program registry. In contrast with alleles, HLA haplotype diversity patterns suggest purifying selection, as certain HLA allele combinations co-occur in high linkage disequilibrium. Linkage disequilibrium is positive (D_ij'>0) among frequent haplotypes and negative (D_ij'<0) among rare haplotypes. Fitting the haplotype frequency distribution to several population dynamics models, we found that the best fit was obtained when significant positive frequency-dependent selection (FDS) was incorporated. Finally, the Ewens-Watterson test of homozygosity showed excess homozygosity for 5-locus haplotypes within 23 US populations studied, with an average Fnd of 28.43. Haplotype diversity is most consistent with purifying selection for HLA Class I haplotypes (HLA-A, -B, -C), and was not inferred for HLA Class II haplotypes (-DRB1 and—DQB1). We discuss our empirical results in the context of evolutionary theory, exploring potential mechanisms of selection that maintain high linkage disequilibrium in MHC haplotype blocks.

The adaptive immune system presents antigens derived from pathogenic and normal self proteins on the cell surface using human leukocyte antigen (HLA) molecules. The HLA loci coding for these molecules are found in major histocompatibility complex (MHC) region, the most polymorphic region in the human genome, with over 15,000 HLA alleles observed so far in the world population. A high frequency of many different HLA alleles is thought be sustained by balancing selection. New HLA alleles may have an advantage over existing frequent alleles since immune escape mutations in pathogens within a population are maintained primarily in epitopes presented on frequent HLA alleles. Host immune function is not determined by single HLA alleles, but by both copies of autosomal HLA genes together (genotypes). Complementarity in function across the two potentially-variant copies of HLA at each locus can result in overdominance and heterozygote advantage at the genotype level. Less explored are selection mechanisms that may be operating across combinations of HLA alleles across loci (haplotypes). Indeed, in addition to high allelic diversity, HLA also has distinctive patterns of haplotype diversity, as certain HLA alleles co-occur in high linkage disequilibrium across five classical HLA loci (HLA-A, -B, -C, -DRB1, -DQB1). We applied multiple population genetic models to a dataset of HLA haplotype frequencies derived from over six million individuals with the goal of determining what type of selection may impact HLA haplotype diversity. We found frequent haplotypes were preferentially maintained in the population across 23 US populations studied. Thus, balancing selection at the allele level and purifying selection at the haplotype level may together affect HLA diversity in human populations.

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.