# The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control

Myriam Rahmouni, Lorenzo de Marco, Jean-Louis Spadoni, Maxime Tison, Raissa Medina-Santos, Taoufik Labib, Josselin Noirel, Ryad Tamouza, Sophie Limou, Olivier Delaneau, et al.

► To cite this version:

Myriam Rahmouni, Lorenzo de Marco, Jean-Louis Spadoni, Maxime Tison, Raissa Medina-Santos, et al.. The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control. Frontiers in Immunology, 2023, 14, 10.3389/fimmu.2023.1305856. hal-04351478

HAL Id: hal-04351478

https://cnam.hal.science/hal-04351478v1

Submitted on 18 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY

Seik-Soon Khor,
Nanyang Technological University,
Singapore

REVIEWED BY

Vinod Kumar,
Radboud University Medical Centre,
Netherlands
Jürgen Sauter,
Deutsche Knochenmarkspenderdatei,
Germany

*CORRESPONDENCE

Jean-François Zagury

<figure><img src="image_3.png" /></figure>

zagury@cnam.fr

†These authors have contributed equally to this work

RECEIVED 02 October 2023

ACCEPTED 10 November 2023

PUBLISHED 11 December 2023

CITATION

Rahmouni M, De Marco L, Spadoni J-L, Tison M, Medina-Santos R, Labib T, Noirel J, Tamouza R, Limou S, Delaneau O, Fellay J, Bensussan A, Le Clerc S, McLaren PJ and Zagury J-F (2023) The *HLA-B*\*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control.
*Front. Immunol.* 14:1305856.
doi: 10.3389/fimmu.2023.1305856

COPYRIGHT

© 2023 Rahmouni, De Marco, Spadoni, Tison, Medina-Santos, Labib, Noirel, Tamouza, Limou, Delaneau, Fellay, Bensussan, Le Clerc, McLaren and Zagury. This is an open-access article distributed under the terms of the *Creative Commons Attribution License (CC BY)*. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control

Myriam Rahmouni¹, Lorenzo De Marco¹, Jean-Louis Spadoni¹,
Maxime Tison¹, Raissa Medina-Santos¹, Taoufik Labib¹,
Josselin Noirel¹, Ryad Tamouza², Sophie Limou³,
Olivier Delaneau⁴, Jacques Fellay⁵, Armand Bensussan⁶,
Sigrid Le Clerc¹, Paul J. McLaren⁷,⁸† and
Jean-François Zagury¹*†

¹Laboratoire Génomique, Bioinformatique, et Chimie Moléculaire, EA7528, Conservatoire National des Arts et Métiers, HESAM Université, Paris, France, ²Université Paris Est Créteil, INSERM U955, IMRB, Laboratoire Neuro-Psychiatrie translationnelle, Créteil, France, ³Nantes Université, Ecole Centrale Nantes, INSERM, Center for Research in Transplantation and Translational Immunology (CR2TI), Nantes, France, ⁴Department of Computational Biology, University of Lausanne, Lausanne, Switzerland, ⁵Global Health Institute, School of Life Sciences, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, ⁶Unité INSERM U976, Hôpital Saint-Louis, Paris, France, ⁷Sexually Transmitted and Blood-Borne Infections Division at JC Wilt Infectious Diseases Research Centre, National Microbiology Laboratory Branch, Public Health Agency of Canada, Winnipeg, MB, Canada, ⁸Department of Medical Microbiology and Infectious Diseases, University of Manitoba, Winnipeg, MB, Canada

**Introduction:** We have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers.

**Methods:** We performed a genome-wide Association Study comparing 543 HIV Elite Controllers with 3,272 uninfected controls of European descent. Using the latest database for imputation, we analyzed 35,552 Single Nucleotide Polymorphisms (SNPs) within the Major Histocompatibility Complex (MHC) region.

**Results:** Our analysis identified 2,626 SNPs significantly associated (p<5.10⁻⁸) with elite control of HIV-1 infection, including well-established MHC signals such as the rs2395029-G allele which tags *HLA-B*\*57:01. A thorough investigation of SNPs in linkage disequilibrium with rs2395029 revealed an extensive haploblock spanning 1.9 megabases in the MHC region tagging *HLA-B*\*57:01, comprising 379 SNP alleles impacting 72 genes. This haploblock contains damaging variations in proteins like NOTCH4 and DXO and is also associated with a strong differential pattern of expression of multiple MHC genes such as *HLA-B*, *MICB*, and *ZBTB12*. The study was expanded to include two cohorts of seropositive African-American individuals, where a haploblock tagging the *HLA-B*\*57:03 allele was similarly associated with control of viral load. The mRNA expression profile of this haploblock in African Americans closely mirrored that in the European cohort.
**Discussion:** These findings suggest that additional molecular mechanisms beyond the conventional antigen-presenting role of class I HLA molecules may contribute to the observed influence of HLA-B*57:01/B*57:03 alleles on HIV-1 elite control. Overall, this study has uncovered a large haploblock associated with HLA-B*57 alleles, providing novel insights into their massive effect on HIV-1 elite control.

KEYWORDS

elite controllers, HIV-1, AIDS, genetics, GWAS, viral load, MHC, haplotype

Introduction

Despite the availability of effective antiretroviral drugs, HIV
continues to be an important public health concern, with a large
number of new infections and deaths occurring each year,
particularly in low-income countries (1) In the early 2000s, a
subset of HIV-positive individuals who maintained consistently
low viral loads for several years without receiving any treatment
were identified and referred to as Elite Controllers (ECs). These
individuals represent approximately 0.2 to 0.5% of the HIV-positive
population (2), including in African cohorts (3), a finding that has
been confirmed in the GRIV cohort comprising Long-Term Non-
Progressors and ECs (see methods). Genetic association studies
initially revealed a strong association between EC status and the
presence of the HLA-B*57:01 allele (4, 5). Subsequently, large-scale
genome-wide association studies (GWAS) confirmed that genetic
variants within the Major Histocompatibility Complex (MHC)
region exerted the most significant genetic influence on viral load
control and disease progression (6, 7). In 2010, a GWAS focusing on
ECs identified four independent top signals, all localized within the
HLA region (8). In 2013, the ICGH consortium, a global
collaborative effort comprising multiple AIDS GWAS studies
published a first large genome-wide association study comparing
6,300 HIV-1 infected individuals and 7,200 uninfected controls and
found no difference except in the well-known CCR5 region (9). In
2015, the same ICGH consortium studied genetic associations with
the viral load at setpoint among the 6,300 HIV-1 infected
individuals and found signals mostly in the MHC and to a lesser
extent in the CCR5 region, reaffirming the prominent impact of the
HLA region on viral load and disease progression (10). However,
since then, no new findings regarding this region have emerged.
Given that EC subjects possess natural control over viral infection,
they represent a valuable population for investigating the molecular
mechanisms of protection. So far, the explanations for the biological
impact of the MHC region in HIV-1 control have been limited to
the role of classical HLA class I alleles, HLA-A, HLA-B, and HLA-C
(10–13). To delve deeper into the understanding of the MHC
region's biological impact on HIV-1 control, we reanalyzed the
genetic data of 543 ECs (with viral loads<1,000 copies/ml) and 3,272
uninfected controls from the ICGH consortium (references IHAC1/
2), with a specific focus on the MHC region. Considering that ECs
constitute 0.2-0.5% of the seropositive population (2), this group

encompasses an extreme phenotype representing over 100,000
seropositive individuals. Leveraging the latest bioinformatics
databases, we imputed 35,552 Single Nucleotide Polymorphisms
(SNPs) within the MHC region, and identified 2,626 significant
signals surpassing the genome-wide significance threshold (5.10-8).
The MHC region spans approximately 5 million base pairs (genomic
coordinates 28,477,797 to 33,448,354) on chromosome 6, according
to the latest GENCODE gene annotation for the GRCh38 reference
genome (14). Within this region, a total of 373 protein-coding genes,
18 pseudogenes, and 12 non-coding RNA genes were identified by
the GENCODE annotation.

To translate these 2,626 SNP associations into function, we then
conducted a comprehensive analysis of these positive signals with
the goal of addressing the following key biological questions,
namely, **a.** the reason for the massive genetic impact of the MHC
region in AIDS, **b.** the factors explaining why only a fraction of
HLA-B*57:01 individuals are Elite Controllers, and **c.** the
understanding of the molecular mechanisms of protection in
Elite Controllers.

Results

Huge number of MHC SNPS associating to
elite control in AIDS

The 543 ECs (VL< 1,000 copies/ml) and 3,272 controls (CTR) from ICGH originated from 2 series of case-control individuals genotyped on two Illumina chips (Table 1A). We could impute 35,552 SNPs with an info score test > 0.75 for all the case and control individuals in the MHC region (see Methods). We compared a total of 543 ECs with 3,272 control individuals using a meta-analysis on the 2 case-control comparisons (see Methods) and identified 2626 significant SNPs (p<5.10-8) out of the 35,552 SNPs tested. Among these SNPs, the top four signals previously reported by Pereyra et al. (8) were confirmed (Table 1B) and also replicated in the independent GRIV study involving 50 ECs (15). Interestingly, by stratifying our EC sample into those with viral loads below 100 copies/ml (referred to as EC2, N=217) and those with viral loads between 100 and 1,000 copies/ml (referred to as EC3, N=326), we noticed a strong enrichment (>30%) of the protective minor allele for rs2395029 corresponding to HLA-
<table><caption>TABLE 1 Description of the cohorts and main signals of Pereyra et al. (8).</caption><thead><tr><th colspan="9">A.</th></tr><tr><th>Groups</th><th>origins</th><th>Platform</th><th>N (EC)</th><th>HLA-B*57:01</th><th>EC2</th><th>EC2_HLA-B*57:01</th><th>EC3</th><th>EC3_HLA-B*57:01</th></tr></thead><tbody><tr><td>ill1</td><td>North America and non-Dutch</td><td>Illumina 550</td><td>418</td><td>119 (28%)</td><td>184</td><td>66 (36%)</td><td>234</td><td>53 (23%)</td></tr><tr><td>ill2</td><td>Non-French European</td><td>Illumina 650 Illumina 1M</td><td>125</td><td>27 (22%)</td><td>33</td><td>9 (27%)</td><td>92</td><td>18 (20%)</td></tr><tr><td colspan="2">TOTAL</td><td></td><td>543</td><td>146 (27%)</td><td>217</td><td>75 (35%)</td><td>326</td><td>71 (22%)</td></tr><tr><th>Groups</th><th>Origins</th><th>Platform</th><th>N (CTR)</th><th colspan="6">HLA-B*57:01</th></tr><tr><td>CTR1</td><td>North America and non-Dutch</td><td>Illumina 550</td><td>2759</td><td colspan="6">185 (6,7%)</td></tr><tr><td>CTR2</td><td>Non-French European</td><td>Illumina 650 Illumina 1M</td><td>513</td><td colspan="6">25 (4,8%)</td></tr><tr><td colspan="2">TOTAL</td><td></td><td>3272</td><td colspan="6">210 (6,4%)</td></tr><tr><th>Groups</th><th>Origins</th><th>Platform</th><th>N (HIV+)</th><th colspan="4">HLA-B*57:03</th><th>EC2</th><th>EC3</th></tr><tr><td>cs1</td><td>African American</td><td>Illumina 650</td><td>380</td><td colspan="4">65 (17%)</td><td>43</td><td>67</td></tr><tr><td>cs2</td><td>African American</td><td>Illumina 1M</td><td>379</td><td colspan="4">66 (17%)</td><td>45</td><td>43</td></tr><tr><td colspan="2">TOTAL</td><td></td><td>759</td><td colspan="4">131 (17%)</td><td>88</td><td>110</td></tr></tbody></table>

<table><thead><tr><th colspan="9">B.</th></tr><tr><th>SNP</th><th>Gene</th><th>MAF CTR</th><th>MAF EC</th><th>MAF EC3</th><th>MAF EC2</th><th>Enrichment</th><th>P; OR EC vs CTR</th><th>P; OR Pereyra et al.</th></tr></thead><tbody><tr><td>rs9264942</td><td>LINC02571</td><td>38.5%</td><td>58.1%</td><td>55.7%</td><td>61.7%</td><td>10.9%</td><td>2.91E-32; 2.2</td><td>2.81E-35; 2.9</td></tr><tr><td>rs4418214</td><td>MICA</td><td>7.2%</td><td>24.4%</td><td>21.2%</td><td>29%</td><td>37.1%</td><td>3.89E-63; 4.2</td><td>1.40E-34; 4.4</td></tr><tr><td>rs2395029</td><td>HCP5</td><td>3.3%</td><td>13.9%</td><td>11.5%</td><td>17.5%</td><td>52.3%</td><td>1.24E-44; 4.8</td><td>9.70E-26; 5.3</td></tr><tr><td>rs3131018</td><td>PSOR1C3</td><td>36.3%</td><td>22.7%</td><td>22.4%</td><td>22.8%</td><td>0%</td><td>9.84E-18; 1.9</td><td>4.210-16; 2.1</td></tr></tbody></table>

A. Description of the populations used in the study and the genotyping platforms, including the number of super-elite controllers with a viral load below 100 copies/ml (EC2) or regular elite controllers with a viral load comprised between 100 and 1000 copies/ml (EC3). This table also confirms the enrichment of subjects carrying the HLA-B*57:01 allele among the various EC populations. There is an enrichment of HLA-B*57:03 in the African-American Cs1 and Cs2 seropositive cohorts caused by the adjunction of additional Ecs. B. The 4 main independent SNP identified by Pereyra et al. 8 in MHC are also found in our EC-CTR case control study with similar p values. We have also specified the enrichment observed for these SNPs in the EC2 group compared to the EC3 group, one can see it is higher than 30% for SNPs rs4418214 and rs2395029.

B*57:01 and rs4418214 corresponding to MICA (Table 1B). To prioritize the most relevant SNPs among the 2,626 SNPs found significant in the MHC region, we thus decided to focus on those exhibiting a similar enrichment, i.e. those with a frequency difference >30% between EC2 and EC3. Out of the 2,626 SNPs, 379 showed an enrichment >30% between EC2 (VL<100 copies/ml) and EC3 (VL between 100 and 1,000 copies/ml) (Figures 1A, B). To enhance the precision of this analysis, we applied clumping to this group of 379 SNPs using an r2 threshold of 0.8, resulting in 54 representative SNPs ranging from position 30,374,976 to position 32,103,233 (Figure 1C). Among these representative SNPs, 44 showed a protective effect (higher minor allele frequency in EC2 compared to EC3 and controls) and 10 SNPs showed a negative effect (lower minor allele frequency in EC2 compared to EC3 and controls).

The distribution of these 54 SNPs was consistent across the two ICGH case groups, as well as in the two control groups.

## A large haploblock is in complete LD with the HLA-B*57:01 allele

For each individual, we then computed 2 scores as the number of carried “protective” alleles (from the 44 SNPs), or alternatively, the number of carried “risk” alleles (from the 10 negatively enriched SNPs). Figures 2A, B depict the box-plot representation of the scores obtained in the ECs (regrouped together as a single group) and in the CTR (regrouped together as a single group). As expected, the CTR individuals globally carry fewer protective alleles than the
<figure><img src="image_4.png" /><figcaption>FIGURE 1
Description of the SNPs associated with EC in the MHC region. Representation of the SNPs with their localization in chromosome 6 GRCh38 (x-axis) and their MAF (y-axis). (A) Representation of the 2626 SNPs significantly associated (p<5.10⁻⁸) with elite control in the EC vs CTR GWAS. In yellow, the SNPs with minor alleles highly enriched (increase>30%) in EC2 subjects (VL<100 copies/ml) compared to EC3 subjects (VL comprised between 100 and 1000 copies/ml) or highly decreased (decrease of the MAF higher than 30% in EC2 compared to EC3). In blue, SNPs with a moderate enrichment or decrease (between 15 to 30%). (B) Representation of the 379 SNPs significant in the GWAS and highly enriched (E>30%) (in blue) or highly decreased (in red). (C) Representation of the 54 SNPs remaining after clumping at r2>=0.8, representative of the previous 379 significant SNPs composing the haploblock.</figcaption></figure>

EC individuals, and also fewer than EC3 and much fewer than EC2 subjects.

To address the question of why only 1-2% of *HLA-B*\*57:01 individuals are Elite Controllers (see Methods), we compared the 210 CTR-*HLA-B*\*57:01 subjects with the 75 EC2-*HLA-B*\*57:01 and 71 EC3-*HLA-B*\*57:01 subjects with respect to each of the 54 variants (Figure 2). Figure 2A shows that the CTR-*HLA-B*\*57:01 subjects carry most of the protective variants (median is 43 out of 44 possible variants), as many as the EC2-*HLA-B*\*57:01 (or EC3-*HLA-B*\*57:01 subjects, not shown). We thus decided to check whether the distribution of the 54 variants was similar between CTR-*HLA-B*\*57:01 and EC2-*HLA-B*\*57:01. Figure 2C illustrates the proportion of individuals carrying the protective (or negative) alleles for each SNP. Remarkably, there was no significant difference observed for any SNP between CTR-*HLA-B*\*57:01 and EC2-*HLA-B*\*57:01. Overall, Figure 2 demonstrates the existence of a huge haploblock that encompasses almost all of the 44 positive SNPs and the major alleles of the 10 negative SNPs, which are carried concurrently by the majority of *HLA-B*\*57:01 individuals, both in the CTR and EC2 groups. As shown in Figure 2C, the carriers of these variants were higher than 59% on the left side of the figure and higher than 78% on the right side, indicating a slight loss of linkage disequilibrium (LD) compared with the middle region (peak at 100%) corresponding to SNP rs112515516 (position 31,361,323), a representative of *HLA-B*\*57:01.

Overall, the 379 enriched SNPs before the clumping form a haploblock ranging from position 30,273,759 to position 32,168,252

(from rs28894080 to rs1802036 SNPs). Figure 2D illustrates that the minor alleles of these SNPs within the haploblock are carried by at least 50% and often by 75% to 100% of *HLA-B*\*57:01 carrying individuals. Overall, *HLA-B*\*57:01 individuals carry a substantial proportion of the minor alleles of the 379 SNPs simultaneously, and there was no significant difference observed between CTR-*HLA-B*\*57:01 and EC-*HLA-B*\*57:01 individuals in terms of distribution of these alleles. In summary, we observe a remarkably large haploblock corresponding to the *HLA-B*\*57:01 allele in Europeans: it spans 1,894,493 base pairs and involves 379 SNPs in LD (Figure 2D). Interestingly, there was no significant signal when comparing by GWAS the ECs carrying the *HLA-B*\*57:01 allele (EC-*HLA-B*\*57:01) with the uninfected controls carrying the *HLA-B*\*57:01 allele (CTR-*HLA-B*\*57:01).

Detailed information about these variants and those in high LD (r²>0.8) can be found in Supplementary Table 1. These variants correspond to SNPs within or close to genes therefore having the potential to impact 72 genes (Supplementary Table 1A), SNPs in 22 long non-coding RNAs (Supplementary Table 1B), and SNPs in 19 intergenic regions (Supplementary Table 1C). Table 2 specifically highlights the exonic non-synonymous variations found in the various SNPs of the haploblock. Two variations found in the proteins NOTCH4 and DXO, respectively G534S and H261Q, are of particular interest due to their potential damaging effects, and they exhibit a strong LD with *HLA-B*\*57:01 (Table 2). NOTCH4 is a transmembrane receptor of 2003 amino-acids, expressed ubiquitously, which has been associated with increased apoptosis in breast cancer (16) and increased inflammation in HIV-1 associated nephropathy (17). DXO is a decapping exoribonuclease, of size 396 AA, which leads to an increase of HIV-1 infectivity (18), its inactivation may thus favor protection against HIV-1 replication.

We then decided to investigate the gene expression pattern associated with these SNP alleles. ECs, with their very low circulating virus, have a quasi-normal physiology. It is thus reasonable to use the GTex data (i.e. uninfected subjects) to deal with the level of mRNA expression in the PBMCs of ECs. Supplementary Table 2A presents the strongest differential expression patterns in PBMCs found for the 54 representative SNPs of the haploblock according to the GTEx database (19). As shown in Supplementary Figure 1, the ratio of the mean level of mRNA expression for the minor allele (which corresponds to the haploblock) over that of the major allele may reach extreme values under 0.5 or higher than 3, demonstrating the high impact of the haploblock on mRNA -and thus protein- expression. Twenty-one genes for which the differential mRNA expression linked to the SNP alleles of the haploblock yield strongly significant p values according to GTEx, are presented in Figure 3A. Among the overexpressed mRNAs, we observe *Mir6891* (p=10⁻⁸¹), *ZBTB12* (p=10⁻⁶⁴), *NOTCH4* (p=10⁻²²), *HLA-S* (p=10⁻³⁴), *HLA-B* (p=10⁻³¹), *HLA-L* (p=10⁻³¹), and *C4A* (p=10⁻²²). Conversely, *MICB* (p=10⁻⁴³), *C4B* (p=10⁻³⁰), and *ZFP57* (p=10⁻²⁶) exhibit underexpression. The overexpression of HLA-B and underexpression of MICB, two well-known immune proteins associated with the *HLA-B*\*57:01 haploblock, align well with expectations. The simultaneous
<figure><img src="image_5.png" /><figcaption>FIGURE 2</figcaption></figure>

Description of the representative SNPs of the European haploblock. (A) Box-plot presenting the numbers of minor alleles of the 44 enriched representative SNPs carried for each individual in a given group: the CTR population (N=2443), the EC group (N=543), the EC2 group (N=217), the EC3 group (N=326), the CTR-HLA-B*57:01 group (N=110), the EC2-HLA-B*57:01 group (N=75). (B) Box-plot presenting the numbers of minor alleles of the 10 decreased representative SNPs carried for each individual in a given group: the CTR population (N=2443), the EC group (N=543), the EC2 group (N=217), the EC3 group (N=326), the CTR-HLA-B*57:01 (CTR who are HLA-B*57:01) group (N=110), the EC2-HLA-B*57:01 (EC2 who are HLA-B*57:01) group (N=75). (C) Histogram presenting the percentage of individuals carrying each minor allele of the 54 representative SNPs: for the CTR-HLA-B*57:01 group (N=110) in blue, and for the EC2-HLA-B*57:01 group (N=75) in green. The SNPs are represented by their localization (GRCh38) on the x-axis. For all SNPs, there is no significant difference observed between the CTR-HLA-B*57:01 and EC-HLA-B*57:01 groups. (D) Histogram presenting the percentage of individuals carrying each minor allele of the 379 SNPs of the haploblock: for the CTR-HLA-B*57:01 group (N=110) in blue, and for the EC2-HLA-B*57:01 group (N=75) in green. The SNPs are represented by their localization in chromosome 6 (GRCh38) on the x-axis.

increase of C4A and decrease of C4B is supposed to promote inflammatory cells (20). NOTCH4 is overexpressed in its mutated form associated with the haploblock minor alleles (Table 2), potentially reinforcing its inhibitory effect. ZBTB12 is a transcriptional factor of 459 residues, specifically expressed in PBMCs, with a potential role in hematopoietic development (21), known to regulate the human endogenous retroviruses (22). ZFP57 is a 452 AA protein involved in DNA methylation and imprinting during development, and a study has also shown its activity in adult PBMCs (23). Overall, we see that the haploblock alleles may induce a differential impact on numerous proteins with potential anti-viral or immune regulatory.

activities. This is a likely source of potent anti-HIV-1 response observed in carriers of HLA-B*57:01, in addition to its known role for antigen recognition and presentation.

## A similar haploblock associated with HLA-B*57:03 in African Americans

In African populations, the Minor Allele Frequency (MAF) of rs1131446 that tags HLA-B*57:03 is 0.037 in 1000Genomes, whereas the MAF of rs2395029 that tags HLA-B*57:01, is 0.004. Conversely, HLA-B*57:03 that has been associated with the control
<table><caption>TABLE 2 List of non-synonymous variations linked to the haploblock.</caption><thead><tr><th>SNP</th><th>Position</th><th>MAF CTR</th><th>p value</th><th>Gene</th><th>Protein mutation</th><th>PolyPhen-2 mutation report</th></tr></thead><tbody><tr><td>rs7750641</td><td>31161533</td><td>11.4%</td><td>2.04E-08</td><td>TCF19</td><td>P109S</td><td>Benign (0%)</td></tr><tr><td>rs41558312</td><td>31411087</td><td>3.2%</td><td>2.08E-46</td><td>MICA</td><td>Q17R</td><td>Benign (9%)</td></tr><tr><td>rs41258944</td><td>32049465</td><td>5.9%</td><td>2.29E-14</td><td>TNXB</td><td>V3186I</td><td>Benign (11%)</td></tr><tr><td>rs41293883</td><td>31507043</td><td>3.6%</td><td>5.62E-38</td><td>MICB</td><td>T212I</td><td>Benign (6%)</td></tr><tr><td>rs12211410</td><td>32081646</td><td>9.2%</td><td>7.22E-13</td><td>TNXB</td><td>R1255H</td><td>Benign (0%)</td></tr><tr><td>rs2280801</td><td>31624287</td><td>5.5%</td><td>4.46E-19</td><td>PRRC2A</td><td>P106L</td><td>Benign (34%)</td></tr><tr><td>rs28732158</td><td>31662541</td><td>5.5%</td><td>4.48E-19</td><td>GPANK1</td><td>G266R</td><td>No prediction available</td></tr><tr><td>rs17207867</td><td>31970635</td><td>9.4%</td><td>1.27E-17</td><td>DXO</td><td>H261Q</td><td>Possibly damaging (66%)</td></tr><tr><td>rs8192591</td><td>32218019</td><td>3.3%</td><td>1.03E-17</td><td>NOTCH4</td><td>G534S</td><td>Possibly damaging (94%)</td></tr></tbody></table>

List of all the MHC SNPs in LD ($r^2>0.8$) with SNPs of the HLA-B*57:01 haploblock in Europeans that corresponds to a non-synonymous variation of the encoded protein. The last column presents the protein impact of the amino-acid variation predicted by Polyphen-2 (in%), and one can remark damaging variations for the proteins NOTCH4 and DXO, the other variations seeming less impacting.

of viral load (24, 25) in Africans is nearly absent in Europeans. Since we had access to African American populations in ICGH, we tried to replicate the observations made in the European population. We analyzed 2 independent groups of HIV-1 positive subjects of African descent available from ICGH, CS1 and CS2, consisting of respectively 380 and 376 subjects (Table 1A). These groups were enriched with elite controllers (ECs) (Table 1A). In both CS1 and CS2, the SNP allele rs1131446-T on chromosome 6 showed the most significant association with viral load at setpoint by regression analysis (data not shown). This SNP exists in Europeans and is in strong LD with rs2395029 ($R^2 = 0.88, D'=0.96$). In the two African American groups CS1 and CS2, rs1131446 had a MAF of 9%, whereas rs2395029 that tags HLA-B*57:01 was not observed. Rs1131446 represents a synonymous change at amino acid position 291 and serves as an excellent proxy for the HLA-B*57:03 allele, as LD assessment between rs1131446 and sequence-based HLA types in a subset of the ICGH sample of African descent (n=789) confirmed that the T allele is indeed linked to HLA-B*57:03 ($r^2 = 0.81$).

To identify SNPs in LD with rs1131446 in both Cs1 and Cs2 populations, we considered all SNPs with an $r^2$ value greater than 0.2 and a D' value greater than 0.8. Through this analysis, we found a total of 44 SNPs that were located between positions 31,274,985 and 31,403,294 on chromosome 6, spanning a distance of 128 kilobases. As shown in Figure 4A, the African American individuals who carry the minor allele rs1131446-T, marking HLA-B*57:03, carry simultaneously most of the 44 SNP minor alleles which are in strong LD (median is 37 SNPs carried out of 44). Reciprocally, each one of the 44 SNP minor alleles are carried by most HLA-B*57:03 carriers in both Cs1 and Cs2 groups (Figure 4B).

Figure 3B presents the mRNA expression pattern in PBMCs for the 44 SNPs, as determined by the GTEx database (see methods). Importantly, a very similar pattern to the one observed for the SNPs of the HLA-B*57:01 haploblock (Figure 3A) was found, with the same differentially transcribed proteins and non-coding RNAs (Figure 3B and Supplementary Table 2B), including MICB, HLA-B, HLA-S, ZBTB12, PSORS1C3, HCG27, and XXbacBPG181B23.7.

A discrepancy was nevertheless observed at the expression level of the 3 mRNA C4A, CCHCR1, and *Mir6891*. The analysis of the functional impact of all the genetic variants in LD with the SNPs of the haploblock revealed no damaging protein mutations. There were only 3 exonic variations observed that were localized in HLA-B and corresponded to the allele HLA-B*57:03, the other genetic variations annotated in the identified haploblock are described in Supplementary Tables 1D, E. The similarity of mRNA expression patterns between Europeans and African Americans suggests a potential additional mechanism for HIV control that supplements the HLA-B*57 alleles response by inducing an antiviral state that limits HIV replication.

Discussion

A previous genetic association study made by the ICGH consortium had shown there was no significant SNP difference between HIV-1 infected individuals and uninfected controls, except at the level of the CCR5 gene region (9). Here, we have conducted a comprehensive analysis focusing specifically on the elite controller (EC) phenotype that we compared with uninfected controls, and we found 2,626 positive signals ($p<10^{-8}$) out of 35,552 SNPs analyzed in the MHC region. By looking for the single nucleotide polymorphisms (SNPs) differentially enriched between stratified ECs: the EC3 group (viral load between 100 and 1,000 copies/ml) and the EC2 group (viral load below 100 copies/ml), we could sub-select 379 SNPs specifically associated with the control of viral load. Unsurprisingly, the strongest signal found was the famous SNP rs2395029, known to tag HLA-B*57:01. There is indeed an enrichment of carriers of this allele among ECs which could be up to 6 times compared to a control population (Table 1A). Further investigation of SNPs in LD with rs2395029 revealed us the presence of a previously unknown large haploblock spanning 1.9 megabases (MB), strongly associated with HLA-B*57:01 in individuals of European descent. The existence of such a substantial haploblock is surprising given the high level of
A

<figure><img src="image_6.png" /></figure>

B

<figure><img src="image_7.png" /></figure>

FIGURE 3
Transcriptional impact of the main SNPs of the haploblocks in the European and African American cohorts. (A) Representation of the transcriptional impact of the 54 representative SNPs of the *HLA-B*\*57:01 haploblock in PBMCs according to GTEx. For each SNP (represented by small blue circles for the enriched alleles and pink circles for the decreased alleles), the 3 genes exhibiting the best p values for mRNA expression were selected, and each SNP was thus connected to 3 genes (represented in small rectangular tags). The link is blue when the SNP minor allele induces an increased mRNA expression, the link is red when the SNP minor allele induces a decreased mRNA expression. Hence, according to GTEx, 21 genes are differentially expressed by the 54 SNPs representing the HLA-B\*57:01 haploblock. (B) Representation of the transcriptional impact of the 44 representative SNPs of the *HLA-B*\*57:03 haploblock in African Americans in PBMCs according to GTEx. For each SNP (represented by small blue circles for the haploblock alleles), the 3 genes exhibiting the best p values for mRNA expression were selected, and each SNP was thus connected to 3 genes (represented in small rectangular tags). The link is blue when the SNP minor allele (or a SNP with r²>0,8) induces an increased mRNA expression, the link is red when the SNP minor allele (or a SNP with r²>0,8) induces a decreased mRNA expression. We observe a highly similar transcriptional profile for the African American *HLA-B*\*57:03 SNPs and for the European *HLA-B*\*57:01 SNPs marked by a green rectangle. Two proteins HCG27 and PSORS1C3 are marked by a magenta star because the mRNA are also expressed at a similar level in European *HLA-B*\*57:01 individuals, but are not in the top 3. A few discrepancies between Europeans and African Americans are observed for the mRNA expression of *C4A*, *Mir6891*, and *CCHCR1* which are marked by a red star.

polymorphism in this region. This haploblock corresponds to the 57.1 ancestral haplotype previously described for classical *HLA* alleles (26). Our investigation has revealed that this ancestral haplotype involves in fact many more gene variants of the *MHC* region and likely contributes to the modification of biological processes involving numerous genes and long non-coding RNAs (lncRNAs), either through mutational effects (Table 2) or differential expression (Figure 3 and Supplementary Figure 1).

Significantly, we extended our study to include two independent cohorts of African Americans, a population known to benefit from *HLA-B*\*57:03-mediated protection. Intriguingly, we discovered that the primary SNP associated with viral load control, rs1131446, which tags *HLA-B*\*57:03, was also linked to a haploblock comprising 44 SNPs in strong linkage disequilibrium, spanning a genomic distance of 128 kilobases. Analysis of mRNA expression patterns associated with this *HLA-B*\*57:03 haploblock in African Americans, using data from GTEx, revealed a striking similarity to the mRNA expression pattern observed in the *HLA-B*\*57:01 haploblock among Europeans (Figures 3A, B). The small discrepancy observed for the mRNA expression of *C4A*, *CCHCR1*, and *Mir6891* is of interest since it suggests that these proteins may not be directly involved in the effect of the *HLA-B*\*57 haploblock on elite control.

The remarkable convergence of mRNA/protein-level impacts between these haploblocks in individuals of European and African descent suggests they could be a common mechanism contributing to the massive effect of *HLA-B*\*57 haplotypes on HIV-1 control. Due to its anti-retroviral activity (ZBTB12, MICB, HLA-B), and to its conservation across ancestries (shared by Europeans and Africans), one could foresee that this haploblock was perhaps involved in the control of retroviruses knowing that the human genome and particularly the MHC region contains an important proportion of endogenous retrovirus (27).

We have thus confirmed the massive effect of HLA-B*57:01 for EC and found a large haploblock that provides new clues to explain this effect. One of our initial queries was why only 1-2% of the *HLA-B*\*57:01 subjects ultimately become EC? To investigate this, we examined potential genetic variations between *HLA-B*\*57:01 CTR and *HLA-B*\*57:01 ECs, but we found no significant differences among the *HLA* common genetic variants between these two groups. The effect of *HLA-B*\*57:01 for the control of HIV-1 replication at viral setpoint (meaning in average 18 months after infection by HIV-1) has been demonstrated with numerous cohorts of seropositive subjects, often involving no more than 500 seropositive subjects (6, 24). The small proportion of elite controllers in such cohorts (0.2% to 0.5% of EC would correspond approximately to 2 or 3 subjects) cannot solely account for the robust association observed in these cohorts for *HLA-B*\*57:01 in relation to viral load control during the early stages of infection (viral setpoint). To explain this robust association, it is thus necessary to consider all carriers of the *HLA-B*\*57 haploblock (a few dozen subjects, given a minor allele frequency of 3% in Europeans and 8% in African Americans). Therefore, it is likely that
<figure><img src="image_8.png" /><figcaption>FIGURE 4</figcaption></figure>

Description of the representative SNPs of the African haploblock. (A) Box-plot presenting the numbers of minor alleles of the 44 SNPs in high LD with rs1131446 (tagging *HLA-B*57:03 in African Americans) carried for each individual in a given group: *HLA-B*57:03 individuals of the Cs1 cohort in blue, and *HLA-B*57:03 individuals of the Cs2 cohort in green. We see that most of the minor alleles of the 44 SNPs are carried simultaneously by the *HLA-B*57:03 individuals of the Cs1 and Cs2 groups (median of the box-plot = 37 for the Cs1 and Cs2 cohorts, out of 44 SNPs). (B) For each of the 44 SNPs of the *HLA-B*57:03 haploblock, histogram presenting the percentage of individuals carrying the minor allele among *HLA-B*57:03 African American subjects: in blue for the Cs1 group, in red for the Cs2 group. One can see that most of the SNP minor alleles are carried by a majority of *HLA-B*57:03 subjects, thus forming a haploblock. The SNPs are represented by their localization in chromosome 6 (GRCh38) on the x-axis.

the haploblock has antiviral effects against HIV-1 in all individuals carrying the *HLA-B*57:01 allele during the early stages of infection and it is not surprising to see a similar SNP distribution within the haploblock in both the control group and the elite controller group (Figure 2). This suggests a shift in our understanding, where this haploblock provides early protection to all carriers, followed by various disease progression patterns that may eventually lead to the elite controller phenotype. A majority of haploblock carriers will likely face additional factors, whether environmental or genetic, that will ultimately result in disease progression and prevent them from achieving EC status. Otherwise, all individuals carrying the *HLA-B*57:01 allele would become ECs. For example, the initial infectious dose could be a critical factor, as a low infectious dose could facilitate the antiviral activity of the haploblock, or on the contrary a higher infectious dose could facilitate viral dissemination despite the antiviral activity conferred by the haploblock. Another impacting factor could be the presence of rare *HLA* genetic variants -which were not considered in this study-yielding for some individuals additional protection alongside the haploblock. Further research examining genetic differences between *HLA-B*57:01 or *HLA-B*57:03 CTR and *HLA-B*57:01 or *HLA-B*57:03 ECs across the entire genome may also help answer this latter question in the future.

In summary, our study has uncovered an haploblock and several protein layers that contribute to the elite controller phenotype of *HLA-B*57 individuals. A key framework emerges, involving the interplay of direct antiviral factors such as the ZBTB12 protein, and the modulation of innate/adaptive immune responses through increased transcription of *HLA-B* and decreased

transcription of *MICB*. The identification of this haploblock within the major histocompatibility complex (MHC) region raises broader questions for future research, such as why the haploblock in Europeans is larger than that in Africans. Additionally, it prompts exploration of haploblocks in Asian populations and the potential existence of other multifunctional haploblocks within the MHC region.

## Materials and methods

For sake of clarity, a flowchart summarizing the study and its main steps is provided in Figure 5.

### Participant phenotypes and case/
control matching

In this study, we used genotyped data from the International Collaboration on HIV-1 Genomics (ICGH). The ICGH project, initiated in 2012 with the support of the National Institutes of Health, aimed to consolidate genomic datasets from HIV-1 infected individuals worldwide. The consortium comprised 26 cohorts of seropositive subjects genotyped on diverse platforms, representing four continents (US, Europe, Australia, Africa). Genotypes for uninfected control individuals were obtained from three participating centers, the Illumina genotype control database and the Myocardial Infarction Genetics Consortium (MIGen) (9, 28). Each dataset underwent preliminary quality control procedures
<figure><img src="image_9.png" /><figcaption>FIGURE 5</figcaption></figure>

Flowchart summarizing the study and its main steps. The left panel shows the workflow performed on the ECs of European descent from ICGH to identify and characterize the haploblock linked to HLA-B*57:01, and the right panel describes the replication study performed to identify a similar haploblock linked to HLA-B*57:03 in two African American cohorts from ICGH.

before centralizing all the data for combined analysis. However, to ensure consistency, additional quality control measures were implemented once the data were submitted, as described in the initial publication of ICGH (9).

All the individual cohorts contributing to the ICGH effort obtained ethical approval from their respective country institutions. The initial publication of ICGH in 2013 focused on susceptibility to HIV-1 infection (9) and included a comparison of 6,300 seropositive individuals with 7,200 uninfected controls of European descent. A subsequent publication in 2015 investigated the genetic association with viral load setpoint among the 6,300 seropositive subjects of European descent (9). In these studies, the cohorts were categorized into six groups of matched cases and controls based on genotyping platforms and geographic origin.

The present study focuses on two groups of Elite Controllers (EC) from the European cohorts, characterized by a viral load of less than 1,000 copies/mm3. Sufficient cases were available for our analysis in ill1 (418 ECs) and ill2 (125 ECs). Corresponding uninfected matched controls included 2,759 subjects in CTR1 and 513 subjects in CTR2. The matching process was described in the initial ICGH publication (9).

To replicate the findings related to the *HLA-B*\*57 haploblock found in the European population, we included two African American cohorts of seropositive subjects available within the ICGH consortium: cs1, consisting of 380 seropositive individuals genotyped with the Illumina chip 650, and cs2, consisting of 379 seropositive individuals genotyped with the Illumina chip 1M. Both cs1 and cs2 cohorts encompass Elite Controllers (ECs) with viral

load levels below 1,000 copies/ml, with 110 and 89 ECs, respectively, as well as regular seropositive individuals.

A summary of the various groups utilized in this study is provided in Table 1A.

## Proportion of elite controllers and proportion of HLA-B*57 carriers among ECs

The GRIV cohort, which includes Long-Term Non-Progressors (LTNP) and Elite Controllers (ECs), is part of the ICGH consortium (9). The GRIV cohort was collected in France based on the LTNP phenotype, defined as documented infection for more than 8 years, maintaining a T lymphocyte count above 500/mm3, and no clinical symptoms (29). Among the 250 GRIV LTNPs with documented viral load, 50 (20%) were identified as ECs with a viral load below 1,000/mm3. At the time of data collection of the GRIV cohort, it was estimated that LTNP represented 1 to 2% of the known HIV-1 infected patients. Consequently, the ECs in the GRIV cohort correspond approximately to 0.2 to 0.4% of the known seropositive patients. This observation aligns well with the findings reported by Okulicz et al. in 2011 (2).

Knowing that 30% of the GRIV ECs carry the *HLA-B*\*57:01 allele, while the prevalence of *HLA-B*\*57:01 in the general French population is 6%, this suggests that the proportion of ECs among *HLA-B*\*57:01 subjects falls within the range of 1 to 2%.
## Participant genotypes

The participants of European descent within the ICGH consortium were organized into matched case-control groups using a two-stage case/control matching strategy, as outlined in the initial publication of ICGH (9). This resulted in four clusters: Group 1 consisted of participants from the Netherlands genotyped on the Illumina platform, Group 2 included participants from France genotyped on the Illumina platform, Groups 3 and 4 encompassed participants from North America and non-Dutch/non-French European regions genotyped on the Illumina platform, and Groups 5 and 6 comprised participants from North America and non-Dutch/non-French European regions genotyped on the Affymetrix platform.

The ill1 EC group and its matching controls (CTR1) originated from Group 3, while the ill2 EC group and its matching controls (CTR2) were derived from Group 4. Genotypic data for these groups were generated using three genotyping arrays: Illumina 550, Illumina 650, and Illumina 1M. Prior to analysis, we conducted standard quality control (QC) procedures on each set (CTR1/ill1 and CTR2/ill2) to ensure the use of clean data.

Specifically, we filtered out monomorphic or rare variants with a minor allele frequency (MAF) less than 1%, as well as structural variations such as insertion-deletions. We also excluded sites with a missingness rate above 0.02 or a Hardy-Weinberg equilibrium (HWE) p-value lower than 10-6. Moreover, we examined the consistency between recorded sex information and sex inferred from the X chromosome, assessed heterozygosity and homozygosity rates, and verified the level of relatedness by calculating identity by descent (IBD). A similar QC process was separately conducted for the Cs1 and Cs2 groups of African American descent, who were genotyped on the Illumina 650 and 1M platforms.

## Imputation

Before imputation, additional quality control (QC) steps were applied to the two sets of European descent (ill1/CTR1 and ill2/CTR2). These QC measures utilized the checkbim steps of the McCarthy Group Tools (30), which involved removing single nucleotide polymorphisms (SNPs) with A/T and G/C alleles if the minor allele frequency (MAF) exceeded 40%. SNPs with discordant alleles, those with more than a 20% difference in allele frequency, and SNPs not present in the 1000 Genomes reference panel were also excluded. These QC steps ensured the utilization of high-quality data for the subsequent imputation process, enabling accurate results.

For ill1/CTR1, a total of 387,495 common SNPs were imputed, while for ill2/CTR2, 271,572 common SNPs were imputed. The imputation process followed a specific protocol: phasing was performed using Eagle2.4 (31), and imputation was carried out on the TOPMed Imputation Server (32) using Minimac4 (33) and the TOPMed reference panel (34).

The African American groups Cs1 and Cs2 underwent a similar imputation protocol, conducted separately from the European groups.

## Stratification

To assess patterns of population structure in the ill1/CTR1 and ill2/CTR2 datasets, a principal component analysis (PCA) was conducted. The analysis aimed to quantify the genetic ancestry of the participants. A set of 510,420 informative SNPs for ill1/CTR1 and 523,850 informative SNPs for ill2/CTR2 was selected for ancestral origin determination. To mitigate the influence of linkage disequilibrium, pruning was applied using an r² threshold of 0.3, a sliding window size of 50, and a step size of 5. Regions on chromosomes 6, 8, and 17 with high linkage disequilibrium were excluded from the analysis.

The results of the PCA confirmed the homogeneity of both the ill1/CTR1 and ill2/CTR2 datasets, consistent with the findings reported in the first ICGH publication (9).

## Association testing by meta-analysis

For the ill1/CTR1 and ill2/CTR2 datasets, a logistic regression analysis comparing Elite Controllers (ECs) and controls (Ill1 vs CTR1 and ill2 vs CTR2) was conducted using SNPtest (35). The analysis focused on 35,552 variant dosages from the Major Histocompatibility Complex (MHC) region under an additive model. To account for population structure and minimize potential confounding effects, the first five principal components (PCs) were included as covariates in both analyses.

Following the individual dataset analyses, a meta-analysis was performed using GWAMA software (36) to combine the p-values from each dataset. The meta-analysis aimed to identify significant associations that satisfy the following conditions: a combined p-value less than or equal to 5.10⁻⁸, both individual p-values less than 0.05, an infotest value greater than 0.75 (indicating a good imputation quality), and the effect sizes (OR) in the same direction across datasets.

## Selection of SNPs significantly associated with EC and EC2

As mentioned in the previous paragraph, after conducting the meta-analysis on the set of 35,552 SNPs imputed in the MHC region, we identified 2,626 significant SNPs. In our analysis, we focused on SNPs with a minor allele frequency (MAF) that was enriched among subjects with the lowest viral loads. To quantify this enrichment, we calculated an enrichment score (ES) between the EC3 group (viral load between 100 and 1,000 copies/ml) and the EC2 group (viral load below 100 copies/ml).

Equation (1) defines the enrichment score ($E_{score}$) as the percentage increase in the SNP's minor allele frequency (MAF).

$$ MAF_{ctr} > MAF_{cv3} > MAF_{cv2} \quad \text{or} \quad MAF_{ctr} < MAF_{cv3} < MAF_{cv2} \quad (1) $$

$$ E_{score} = \frac{MAF_{cv2} - MAF_{cv3}}{MAF_{cv3}} * 100 $$
$MAFcv2$ = MAF of the SNP in EC with viral load< 100 copies/ml.
$MAFcv3$ = MAF of the SNP in EC with 100 copies/ml< viral load< 1000 copies/ml.

$MAFctr$ = MAF of the SNP in controls (uninfected).

Following the identification of 2,626 significant SNPs from the meta-analysis, we further examined the subset of SNPs that exhibited an increase (or decrease) in minor allele frequency (MAF) of more than 30% in EC2 compared to EC3 and controls. Out of the 2,626 significant SNPs, 379 SNPs met this criterion, indicating a substantial increase (or decrease) in MAF in EC2.

To work with a more manageable set of representative SNPs, we performed a clumping analysis using the program PLINK (37) on these 379 SNPs. The clumping analysis grouped together SNPs with an $r^2$ correlation of 0.8 or higher and located within a physical distance of 200,000 base pairs. This clustering process resulted in a final set of 54 representative SNPs that could be used for further analysis.

Among these 54 representative SNPs, 44 exhibited a MAF increase greater than 30% between EC2 and EC3, corresponding to “protective” alleles. Additionally, 9 SNPs showed a MAF decrease greater than 30% between EC2 and EC3, corresponding to “risk” alleles. These representative SNPs provide a focused subset for subsequent investigations and analysis.

## Computation of the LD in various subgroups

To identify SNPs in linkage disequilibrium (LD) with SNP rs1131446, which serves as a marker for HLA-B*57:03 in African Americans, we employed the software PLINK (37) to conduct LD analysis on the genotypes of our cohorts. PLINK utilizes the maximum likelihood method to estimate LD coefficients, including the correlation coefficient ($r^2$) and normalized LD coefficient (D'), between pairs of SNPs (38).

To establish appropriate threshold values for LD analysis, we referred to values observed in the European haploblock with SNP rs2395029, the reference marker for HLA-B*57:01 in Europeans. Based on this, we selected a threshold of an $r^2$ correlation coefficient greater than 0.2 and a D' coefficient greater than 0.8 with rs1131446 in African Americans. We applied these criteria to both the Cs1 and Cs2 populations, and it resulted in the identification of 44 SNPs that exhibited significant LD with rs1131446 in both populations. These 44 SNPs provide valuable candidates for further investigation and analysis.

Consistent tools and methods were utilized for all $r^2$ calculations in our analysis, ensuring uniformity and comparability across the study.

## Score per individual and score per SNP allele, in each subgroup

For each individual (id), we computed 2 scores as the number of carried “protective” alleles, or alternatively, the number of carried “risk” alleles, from the 54 representative SNPs.

Knowing that 0 represents the alternate allele and 1, the reference allele, knowing that n is the number of protective or negative SNPs, knowing that $P_s(0/0)$ and $P_s(0/1)$ are the estimated posterior genotype probabilities for the SNP s imputed for individual id, we obtained a score per individual as presented in equation (2).

$$ Score_{id} = \sum_{s=1}^{s=n} P_s(0/0) + P_s(0/1) \quad (2) $$

We negate this score for under-enriched SNPs.

Each subgroup was represented by a boxplot in which each dot corresponded to the score of an individual (Figure 2). The Mann-Whitney U test was then employed to compare the distribution of the scores between subgroups.

For each SNP, it was also interesting to investigate the number of subjects carrying its minor allele in the various groups: the control group CTR, the elite control group EC (as well as EC2 and EC3), among the carriers of HLA-B*57:01 in Europeans or among the carriers of HLA-B*57:03 in African Americans (CTR_HLA-B*57:03, EC2_HLA-B*57:03, CS1_HLA-B*57:03, CS2_HLA-B*57:03). The HLA-B*57:01 carriers are those who have at least one G for SNP rs2395029, and the HLA-B*57:03 carriers are those who carry at least one T for SNP rs1131446 (since it is a good proxy of HLA-B*57:03 with $r^2>0,8$).

For each SNP, we thus computed a score in a given subgroup, corresponding to the global percentage of subjects wearing its minor allele, as presented in equation (3).

$$ Score_s = \frac{\sum_{id=1}^{id=N} P_{s,id}(0/0) + P_{s,id}(0/1)}{N} * 100 \quad (3) $$

## Biological exploration of the SNPs with genetic annotation (Annovar)

We annotated the 379 SNPs from the European haploblock and the 44 SNPs from the African haploblock using Annovar (39) on dbSNP (avSNP150) (40) and the RefGene gene database (41). In order to include all mutants of interest, we considered in our analysis all SNPs with a $r^2$ LD value greater than 0.8 with any SNP of the haploblock during the annotation process. We used Polyphen 2 (42) to predict the potential effect of amino acid substitutions on the structure and function of the protein for nonsynonymous SNPs identified.

## Transcriptional impact of the SNPs with GTEx

We used Ldexpress (43) from Ldlink (44), based on Genotype-Tissue Expression project data (GTEx, v8 release) (19), to investigate the genes significantly and differentially expressed in whole blood according to the alleles of the haploblock SNPs (as well as SNPs with $r^2 > 0.8$). GTEx uses the genetic LD found in specified populations (either of European descent or of African descent) to
select the SNPs of interest (with an r²>= 0.8 in our case), and then
computes the differential transcription of genes for these SNPs
among all samples (838 samples including 103 samples of African
American descent).

For each SNP, we picked up the 3 genes most differentially
expressed, i.e. that exhibited the lowest p-values.

The results obtained for the SNPs of the haploblock were
then visualized in a network created with Cytoscape tool (45). On
the representation of Figure 3, we use colours to present the
direction of the enrichment in EC2 vs EC3 for the minor SNP
allele ("protective" or "risk" allele), and the way by which it
influences the mRNA gene expression in GTEx (increased or
decreased expression).

Data availability statement

The data presented in the study are deposited in the Figshare
repository, 10.6084/m9 figshare.24615549.

Ethics statement

Ethical approval for this study involving human genetic data
was obtained for all cohorts by the local institutional review board
of each group/center/study participating to the ICGH collective
effort. All subjects provided written informed consent. The list of
the ICGH cohorts/participants is as follows : 1. The AIDS clinical
Trial Group (ACTG) in the USA 2. The AIDS Linked to the
IntraVenous Experience (ALIVE) Cohort in Baltimore, USA 3.
The Amsterdam Cohort Studies on HIV infection and AIDS (ACS)
in the Netherlands 4. The ANRS CO18 in France 5. The ANRS
PRIMO Cohort in France 6. The Center for HIV/AIDS Vaccine
Immunology (CHAVI) in the USA 7. The Danish HIV Cohort
Study in Denmark 8. The Genetic and Immunological Studies of
European and African HIV-1+ Long Term Non-Progressors
(GISHEAL) Study, in France and Italy 9. The GRIV Cohort in
France 10. The Hemophilia Growth and Development Study
(HGDS) in the USA 11. The Hospital Clinic-IDIBAPS Acute/
Recent HIV-1 Infection cohort in Barcelona, Spain 12. The Icona
Foundation Study in Italy 13. The International HIV Controllers
Study in Boston, USA 14. The IrsiCaixa Foundation Acute/Recent
HIV-1 Infection cohort in Barcelona, Spain 15. The Modena Cohort
in Modena, Italy 16. The Multicenter AIDS Cohort Study (MACS),
in Baltimore, Chicago, Pittsburgh and Los Angeles, USA 17. The
Multicenter Hemophilia Cohort Studies (MHCS) 18. The NCI
Laboratory of Genomic Diversity in Frederick, USA 19. The
Pumwani Sex Workers Cohort in Nairobi, Kenya, and Winnipeg,
Canada 20. The San Francisco City Clinic Cohort (SFCCC) in San
Francisco, USA 21. The Sanger RCC Study in Oxford, UK, and in
Uganda 22. The Swiss HIV Cohort Study (SHCS), in Switzerland
23. The US military HIV Natural History Study (NHS) 24. The
Wellcome Trust Case Control Consortium (WTCCC3) study of
the genetics of host control of HIV-1 infection in the Gambia 25.

The West Australian HIV cohort Study. The studies were
conducted in accordance with the local legislation and
institutional requirements.

Author contributions

MR: Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Data curation. LM: Data curation, Formal analysis, Visualization, Writing – original draft. J-LS: Formal analysis, Software, Validation, Writing – review & editing. MT: Data curation, Formal analysis, Visualization, Writing – review & editing. RM-S: Formal analysis, Methodology, Writing – review & editing. TL: Formal analysis, Data curation, Writing – review & editing. JN: Data curation, Writing – review & editing. RT: Formal analysis, Writing – review & editing. SL: Data curation, Validation, Writing – original draft. OD: Data curation, Writing – review & editing, Methodology. JF: Writing – review & editing, Funding acquisition, Project administration, Resources, Validation. AB: Writing – review & editing, Project administration, Supervision, Validation. SC: Data curation, Formal analysis, Writing – review & editing, Supervision. PM: Writing – original draft, Writing – review & editing, Data curation, Methodology, Resources, Validation, Visualization. J-FZ: Writing – original draft, Writing – review & editing, Investigation, Methodology, Validation, Visualization, Conceptualization, Funding acquisition, Project administration, Resources, Supervision.

Funding

The author(s) declare financial support was received for the
research, authorship, and/or publication of this article. Initial
funding for the International Collaboration for the Genomics of
HIV was provided by the NIH Office of AIDS Research. This project
has been funded in whole or in part with federal funds from the
National Cancer Institute, National Institutes of Health, under
contract HHSN26120080001E. MT is recipient of a fellowship
from the Foundation FundaMental. The Laboratory GBCM
acknowledges the program Mécénat-Santé of Mutuelles AXA for
funding its research. The funders had no role in study design, data
analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

The authors are grateful to all the patients and medical staff who
contributed to the collection of the cohorts. The authors thank
Stuart Z. Shapiro (Program Officer, Division of AIDS, National
Institute of Allergy and Infectious Diseases) and Stacy Carrington-
Lawrence (Chair of Etiology and Pathogenesis, NIH Office of AIDS
Research) who accompanied the creation of the International
Collaboration on the Genomics of HIV (ICGH). The authors are
also grateful to all the scientists who participated to the ICGH
international effort and helped constitute this exceptional collection of data. MT is recipient of a fellowship from the Foundation FundaMental. The Laboratory GBCM is grateful to the program Mécénat-Santé of Mutuelles AXA for funding its research.

## Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1305856/full#supplementary-material

## References

1. UNAIDS. *estimates* (2021). Available at: https://www.unaids.org/en.

2. Okulicz JF, Lambotte O. Epidemiology and clinical characteristics of elite controllers. *Curr Opin HIV AIDS* (2011) 6(3):1638. doi: 10.1097/COH.0b013e328344f35e

3. Berg MG, Olivo A, Harris BJ, Rodgers MA, James L, Mampunza S, et al. A high prevalence of potential HIV elite controllers identified over 30 years in Democratic Republic of Congo. *EBioMedicine. mars* (2021) 65:103258. doi: 10.1016/j.ebiom.2021.103258

4. Kaslow RA, Carrington M, Apple R, Park L, Muñoz A, Saah AJ, et al. Influence of combinations of human major histocompatibility complex genes on the course of HIV-1 infection. *Nat Med* (1996) 2(4):40511. doi: 10.1038/nm0496-405

5. Hendel H, Caillat-Zucman S, Lebuanec H, Carrington M, O'Brien S, Andrieu JM, et al. New class I and II HLA alleles strongly associated with opposite patterns of progression to AIDS. *J Immunol* (1999) 162(11):69426. doi: 10.4049/jimmunol.162.11.6942

6. Fellay J, Shianna KV, Ge D, Colombo S, Ledergerber B, Weale M, et al. A whole-genome association study of major determinants for host control of HIV-1. *Science* (2007) 317(5840):9447. doi: 10.1126/science.1143767

7. Limou S, Le Clerc S, Coulonges C, Carpentier W, Dina C, Delaneau O, et al. Genomewide association study of an AIDS-nonprogression cohort emphasizes the role played by HLA genes (ANRS Genomewide Association Study 02). *J Infect Dis* (2009) 199(3):41926. doi: 10.1086/596067

8. International HIV Controllers Study, Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PIW, et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. *Science* (2010) 330(6010):15517. doi: 10.1126/science.1195271

9. McLaren PJ, Coulonges C, Ripke S, van den Berg L, Buchbinder S, Carrington M, et al. Association study of common genetic variants and HIV-1 acquisition in 6,300 infected cases and 7,200 controls. *PloS Pathog* (2013) 9(7):e1003515. doi: 10.1371/journal.ppat.1003515

10. McLaren PJ, Coulonges C, Bartha I, Lenz TL, Deutsch AJ, Bashirova A, et al. Polymorphisms of large effect explain the majority of the host genetic contribution to variation of HIV-1 virus load. *Proc Natl Acad Sci U.S.A.* (2015) 112(47):1465863. doi: 10.1073/pnas.1514867112

11. Le Clerc S, Limou S, Zagury JF. Large-scale « OMICS » Studies to explore the physiopathology of HIV-1 infection. *Front Genet* (2019) 10:799. doi: 10.3389/fgene.2019.00799

12. McLaren PJ, Fellay J. HIV-1 and human genetic variation. *Nat Rev Genet* (2021) 22(10):64557. doi: 10.1038/s41576-021-00378-0

13. Naranbhai V, Carrington M. Host genetic variation and HIV disease: from mapping to mechanism. *Immunogenetics* (2017) 69(89):48998. doi: 10.1007/s00251-017-1000-z

14. Frankish A, Diekhans M, Jungreis I, Lagarde J, Loveland JE, Mudge JM, et al. Gencode 2021. *Nucleic Acids Res* (2021) 49(D1):D91623. doi: 10.1093/nar/gkaa1087

15. Le Clerc S, Delaneau O, Coulonges C, Spadoni JL, Labib T, Laville V, et al. Evidence after imputation for a role of MICA variants in nonprogression and elite control of HIV type 1 infection. *J Infect Dis* (2014) 210(12):194650. doi: 10.1093/infdis/jiu342

16. Naik S, MacFarlane M, Sarin A. Notch4 signaling confers susceptibility to TRAIL-induced apoptosis in breast cancer cells. *J Cell Biochem* (2015) 116(7):137180. doi: 10.1002/jcb.25094

17. Puri RV, Yerrathota S, Home T, Idowu JY, Chakravarthi VP, Ward CJ, et al. Notch4 activation aggravates NF-κB-mediated inflammation in HIV-1-associated nephropathy. *Dis Model Mech* (2019) 12(12):dmm040642. doi: 10.1242/dmm.040642

18. Benoni B, Benoni R, Trylcova J, Grab K, Pačes J, Weber J, et al. HIV-1 infection reduces NAD capping of host cell snRNA and snoRNA. *Biochemistry* (2022). doi: 10.1101/2022.11.10.515957

19. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. *Science* (2020) 369(6509):131830. doi: 10.1126/science.aaz1776

20. Markiewski MM, Lambris JD. The role of complement in inflammatory diseases from behind the scenes into the spotlight. *Am J Pathol* (2007) 171(3):71527. doi: 10.2353/ajpath.2007.070166

21. Noro F, Gianfagna F, Gialluisi A, De Curtis A, Di Castelnuovo A, Napoleone E, et al. ZBTB12 DNA methylation is associated with coagulation- and inflammation-related blood cell parameters: findings from the Moli-family cohort. *Clin Epigenet* (2019) 11(1):74. doi: 10.1186/s13148-019-0665-6

22. Han D, Liu G, Oh Y, Oh S, Yang S, Mandjikian L, et al. ZBTB12 is a molecular barrier to dedifferentiation in human pluripotent stem cells. *Nat Commun* (2023) 14(1):632. doi: 10.1038/s41467-023-36178-9

23. Plant K, Fairfax BP, Makino S, Vandiedonck C, Radhakrishnan J, Knight JC. Fine mapping genetic determinants of the highly variably expressed MHC gene ZFP57. *Eur J Hum Genet* (2014) 22(4):56871. doi: 10.1038/ejhg.2013.244

24. Pelak K, Goldstein DB, Walley NM, Fellay J, Ge D, Shianna KV, et al. Host determinants of HIV-1 control in African Americans. *J Infect Dis* (2010) 201(8):11419. doi: 10.1086/651382

25. McLaren PJ, Ripke S, Pelak K, Weintrobl AC, Patsopoulos NA, Jia X, et al. Fine-mapping classical HLA variation associated with durable host control of HIV-1 infection in African Americans. *Hum Mol Genet* (2012) 21(19):433447. doi: 10.1093/hmg/dds226

26. Degli-Esposti MA, Leaver AL, Christiansen FT, Witt CS, Abraham LJ, Dawkins RL. Ancestral haplotypes: conserved population MHC haplotypes. *Hum Immunol* (1992) 34(4):24252. doi: 10.1016/0198-8859(92)90023-G

27. Dawkins R, Leelayuwat C, Gaudieri S, Tay G, Hui J, Cattley S, et al. Genomics of the major histocompatibility complex: haplotypes, duplication, retroviruses and disease. *Immunol Rev* (1999) 167:275304. doi: 10.1111/j.1600-065X.1999.tb01399.x

28. Kathiresan SMyocardial Infarction Genetics Consortium. A PCSK9 missense variant associated with a reduced risk of early-onset myocardial infarction. *N Engl J Med* (2008) 358(21):2299300. doi: 10.1056/NEJMcr0707445

29. Winkler CA, Hendel H, Carrington M, Smith MW, Nelson GW, O'Brien SJ, et al. Dominant effects of CCR2-CCR5 haplotypes in HIV-1 disease progression. *J Acquir Immune Defic Syndr* (2004) 37(4):15348. doi: 10.1097/01.qai.0000127353.01578.63

30. McCarthy tools. Available at: https://www.well.ox.ac.uk/~wrayner/tools/.

31. Loh PR, Danecek P, Palamara PF, Fuchsberger C, A Reshef Y, K Finucane H, et al. Reference-based phasing using the Haplotype Reference Consortium panel. *Nat Genet* (2016) 48(11):14438. doi: 10.1038/ng.3679

32. Das S, Forer L, Schönherr S, Sidore C, Locke AE, Kwong A, et al. Next-generation genotype imputation service and methods. *Nat Genet* (2016) 48(10):12847. doi: 10.1038/ng.3656

33. Fuchsberger C, Abecasis GR, Hinds DA. minimac2: faster genotype imputation. *Bioinformatics* (2015) 31(5):7824. doi: 10.1093/bioinformatics/btu704
34. Taliun D, Harris DN, Kessler MD, Carlson J, Szpiech ZA, Torres R, et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. *Nature* (2021) 590(7845):2909. doi: 10.1038/s41586-021-03205-y

35. Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. *Nat Genet* (2007) 39(7):90613. doi: 10.1038/ng2088

36. Mägi R, Morris AP. GWAMA: software for genome-wide association meta-analysis. *BMC Bioinf* (2010) 11:288. doi: 10.1186/1471-2105-11-288

37. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. *Am J Hum Genet* (2007) 81(3):55975. doi: 10.1086/519795

38. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. *Gigascience* (2015) 4:7. doi: 10.1186/s13742-015-0047-8

39. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. *Nucleic Acids Res* (2010) 38(16):e164. doi: 10.1093/nar/gkq603

40. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. *Nucleic Acids Res* (2001) 29(1):30811. doi: 10.1093/nar/29.1.308

41. Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez Gene: gene-centered information at NCBI. *Nucleic Acids Res* (2011) 39 part 2: Database issue:D52–57. doi: 10.1093/nar/gkq1237

42. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. *Nat Methods* (2010) 7(4):2489. doi: 10.1038/nmeth0410-248

43. Lin SH, Thakur R, Machiela MJ. LDexpress: an online tool for integrating population-specific linkage disequilibrium patterns with tissue-specific expression data. *BMC Bioinf* (2021) 22(1):608. doi: 10.1186/s12859-021-04531-8

44. Machiela MJ, Chanock SJ. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. *Bioinformatics* (2015) 31(21):35557. doi: 10.1093/bioinformatics/btv402

45. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. *Genome Res* (2003) 13(11):2498504. doi: 10.1101/gr.1239303