key: cord-0715584-hdaa45ah
authors: Liu, Fuguo; Cocker, Alexander T. H.; Pugh, Jason L.; Djaoud, Zakia; Parham, Peter; Guethlein, Lisbeth A.
title: Natural LILRB1 D1-D2 variants show frequency differences in populations and bind to HLA class I with various avidities
date: 2022-05-13
journal: Immunogenetics
DOI: 10.1007/s00251-022-01264-7
sha: 986d46ca9bd041dfb96c9081ba83a48fda4bc9d4
doc_id: 715584
cord_uid: hdaa45ah

Leukocyte immunoglobulin-like receptor B1 (LILRB1) is widely expressed on various immune cells and the engagement of LILRB1 to HLA class I and pathogen-derived proteins can modulate the immune response. In the current study, 108 LILRB1 alleles were identified by screening the LILRB1 locus from the 1000 Genomes Phase 3 database. Forty-six alleles that occurred in three or more individuals encode 28 LILRB1 allotypes, and the inferred LILRB1 allotypes were then grouped into 9 LILRB1 D1-D2 variants for further analysis. We found that variants 1, 2, and 3 represent the three most frequent LILRB1 D1-D2 variants and the nine variants show frequency differences in populations. The binding assay demonstrated that variant 1 bound to HLA class I with the highest avidity, and all tested LILRB1 D1-D2 variants bound to HLA-C with lower avidity than to HLA-A and -B. Locus-specific polymorphisms at positions 183, 189, and 268 in HLA class I and dimorphisms in HLA-A (positions 207 and 253) and in HLA-B (position 194) affect their binding to LILRB1. Notably, the electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I as revealed by electrostatic analysis and by comparison of different binding avidities caused by polymorphisms at positions 72 and 103 of LILRB1. In this paper, we present a comprehensive study of the population genetics and binding abilities of LILRB1. The data will help us better understand the LILRB1-related diversity of the immune system and lay a foundation for functional studies. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00251-022-01264-7.

Leukocyte immunoglobulin-like receptor B1 (LILRB1), also known as LIR-1, ILT-2, MIR-7, and CD85j, is an inhibitory receptor that was first identified by Cosman et al. in 1997 . LILRB1 is encoded within the leukocyte receptor complex, which harbors the LILR and killer cell immunoglobulin-like receptor (KIR) families, and occupies about 1 Mb of human chromosome 19q13.4 (Wende et al. 1999; Canavez et al. 2001) . LILRB1 is widely expressed on various immune cells that include macrophages, monocytes, dendritic cells, B cells, and subsets of NK cells and T cells Burshtyn and Morcos 2016) . It is also expressed on decidual immune cells (McIntire et al. 2008) . The expression pattern of LILRB1 suggests that it plays crucial roles in both immune responses and in supporting fetal tolerance in pregnancy.

LILRB1 is a type I transmembrane glycoprotein with four extracellular immunoglobulin-like (Ig-like) domains and four immunoreceptor tyrosine-based inhibitory motifs (ITIMs) within the cytoplasmic region for signal transduction . The two membrane distal Ig domains (D1-D2) are responsible for interacting with HLA class I, with the D1 domain engaging with the α3 region of the HLA class I heavy chain, and the D2 domain interacting with beta 2 microglobulin (β2m) . No major binding sites have been identified within the D3 and D4 domains of LILRB1 (Nam et al. 2013) . Since the α3 region of HLA class I antigens and β2m are relatively conserved, LILRB1 binds to a broad range of HLA class I. The LILRB1-HLA class I binding pattern is different from that of KIRs and T cell receptors (TCR), which recognize the polymorphic α1 and α2 domains of subsets of HLA class I (Djaoud and Parham 2020) . Jones et al. found that LILRB1 bound to a broad range of HLA class I by testing the binding of one LILRB1 variant (unspecified, but most likely variant 1) to a panel of 97 HLA class I antigens in a single antigen bead assay (Jones et al. 2011) . The effect of amino acid changes in the D1-D2 domains of LILRB1 on binding to select HLA class I has been investigated. The four most-studied amino acid changes in this region are L68P, A93T, I142T, and S155I, with the numbers corresponding to the full-length protein sequence. Yu et al. found that LILRB1 bound to HLA-B*58, with the LILRB1 LAIS variant binding to a lesser extent than the PTTI variant (Yu et al. 2018) , and Kuroki et al. demonstrated that the LILRB1 LAIS and PTTI variants bind to HLA-A*11, HLA-C*04, HLA-B*35, and HLA-G1 to a similar degree (Kuroki et al. 2005) . These studies focused on the two common LILRB1 variants and a small number of HLA ligands. As the differences in binding avidities between LILRB1 and HLA class I may substantially be related to the diversity of immune response among individuals, we determined the range of natural variation in LILRB1 and tested binding against a comprehensive panel of HLA allotypes.

It has also been shown that LILRB1 is related to pathogen infection and several diseases. LILRB1 binds to the HLA class I mimic UL18 encoded by human cytomegalovirus (Yang and Bjorkman 2008) with a binding affinity > 1000-fold higher than that for HLA class I (Chapman et al. 1999) . LILRB1 also binds to a subset of repetitive interspersed families of polypeptides (RIFINs) of the P. falciparum parasite, which can inhibit LILRB1-expressing B cells and NK cells and help to evade the host immune system (Saito et al. 2017) . LILRB1 is also associated with several autoimmune diseases, such as systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) (Monsiváis-Urenda et al. 2007; Kuroki et al. 2005) . LILRB1 plays dual roles in cancer biology, acting both as an immune checkpoint molecule and as a tumor-sustaining factor. These observations of the roles of LILRB1 in disease have highlighted the importance of analyzing the LILRB1 diversity across populations to lay a foundation for functional studies to determine disease mechanism. However, the natural LILRB1 alleles and their frequencies in each human population largely remain unknown.

In the present study, we analyzed the LILRB1 locus from the 1000 Genomes Phase 3 database. LILRB1 alleles and allotypes were inferred, and their frequencies in populations were calculated. Eight natural LILRB1 D1-D2 variants were produced and their binding avidities for a panel of 97 HLA class I antigens by single antigen bead assay were tested. The data will help us better understand the LILRB1-related diversity of the immune system and lay a foundation for functional studies.

In order to assess global genetic variation in LILRB1 locus, we selected individuals from the 1000 Genomes Project database (http:// www. inter natio nalge nome. org) that had read depth sufficient for robust SNP determination (Clarke et al. 2017; Auton et al. 2015) . We correlated read depth to mapped.bam file size by selecting 15 individuals with file sizes ranging from 15 to 40 Gb and determining the read depth for each coding exon (E3-E16) of LILRB1 using MOSDEPTH (Pedersen and Quinlan 2018) . Sufficient read depth for each coding exon of LILRB1 was obtained when individuals had mapped.bam file sizes of 30 Gb and larger. Using this cutoff, 429 individuals (out of 2504 assessed) from 26 sub-populations were selected from the 1000 Genomes Project Phase 3 database.

We obtained the SNPs for the complete LILRB1 region on human chromosome 19 (55, 085, 148, 979) using Data Slicer (http:// grch37. ensem bl. org/ Homo_ sapie ns/ Tools/ DataS licer). The focus of this analysis was variation in the coding region so we excluded SNPs in the introns, 5ʹ-UTR, and 3ʹ-UTR. In the coding region, 149 SNPs were found, 66 of which only occur once and were excluded in the present study. The remaining 83 SNPs were analyzed by PHASE2.1 (Stephens et al. 2001 ) using the default parameters to infer the alleles. The inferred alleles were then translated into corresponding allotypes for further analysis. This analysis provided all common naturally occurring variants and their population frequencies.

The sequence encoding the D1-D2 domains of LILRB1 variant 1 and the hinge-CH2-CH3 region of human IgG1 was synthesized by GenScript, with XbaI and NotI restriction sites at the 5ʹ-end and 3ʹ-end respectively. The synthesized sequence was first cloned into the NotI/XbaI sites of the pcDNA3.1( +) vector and the sequence confirmed by Sanger sequencing (MCLAB, CA, USA). The sequencing primers are listed in Table 1 . Site-directed mutagenesis was conducted using the plasmid pcDNA3.1-LILRB1-D1D2-Fc to generate the sequences encoding the remaining eight LILRB1-D1D2-Fc proteins. The primers used for site-directed mutagenesis are listed in Table 1 . Sequences of all variants were confirmed by Sanger sequencing. The pcDNA3.1-LILRB1-D1D2-Fc constructs were digested with XbaI and NotI (New England Biolabs, UK) and the target sequences were cloned into XbaI and NotI digested pAcGP67A vector to be used for protein production. The sequences in the pAcGP67A-LILRB1-D1D2-Fc constructs were again confirmed by Sanger sequencing. Despite repeated attempts at mutagenesis, we were unable to successfully produce the variant 8 clone and continued the analysis with the remaining eight variants.

Protein production and purification was conducted as described (Hilton et al. 2015) . Briefly, the pAcGP67A-LILRB1-D1D2-Fc constructs were co-transferred into Sf9 insect cells with linearized baculovirus (Expression System, Cat. #91-002, USA) using Cellfectin (Invitrogen, CA, USA). SDS-PAGE was performed to measure the protein yield after each round of viral amplification. For large-scale protein production, 1 L of Sf9 cells at a cell density of 2 M/mL was infected with 1 mL of high viral titer P2 supernatant for 96 h. Cell supernatant was collected by centrifugation (2500 g, 10 min, 4 °C) and filtered through a 0.2-µm filter (Nalgene, Cat. #450-0020). After incubation with protein A sepharose beads overnight at 4 °C, the beads were washed with 500 mL of PBS. The LILRB1-D1D2-Fc proteins were eluted from the beads with 0.1 M glycine (pH 2.7), with the eluted protein being immediately neutralized with 0.2 M Tris-base (pH 9). The proteins were adjusted to 1 mg/mL and stored at −20 °C for further use.

The integrity of the purified proteins was assessed by flow cytometry using a Cytek Aurora spectral cytometer. Briefly, 20 µL of anti-human IgG-coated beads was added to 500 µL of LILRB1-D1D2-Fc fusion protein at a concentration of 100 µg/mL. After incubation at 4 °C for 30 min, the beads were collected by centrifugation (250 g for 3 min) and washed 3 times with FACS buffer (FCB). Following the last wash, the beads were resuspended in 50 µL of FCB with 2 µL of APC conjugated anti-human LILRB1 antibody (Thermo Fisher, Cat. #17-5129-42). After incubation at 4 °C for 30 min, the beads were collected by centrifugation (250 g for 3 min) and washed 3 times with FCB. The beads were then resuspended in 150 µL of PBS and analyzed by flow cytometry.

LILRB1-D1D2-Fc proteins were tested for binding to 97 classical HLA class I antigens (31 HLA-A, 50 HLA-B, and 16 HLA-C) using the single antigen bead assay developed by One Lambda (Thermo Fisher Scientific, Cat. #LS1A04) as described (Hilton et al. 2017) . Briefly, LILRB1-D1D2-Fc proteins were diluted to 100 µg/mL with PBS and incubated with the LABScreen microbeads (60 min at 4 °C, gently shaking). As a positive control for antigen density on individual beads, beads were incubated with W6/32, which is a mouse mAb that binds a common epitope to all HLA-A, -B, and -C (W6/32 (Barnstable et al. 1978) ). Beads incubated with same volume of PBS were the negative control. After 4 washes, LILRB1-D1D2-Fc binding beads and negative control beads were resuspended in 100 µL PBS with 1% PE-conjugated goat anti-human IgG-Fc antibody and W6/32 binding beads were resuspended in 100 µL PBS -T7  TAA TAC GAC TCA CTA TAG G  pcDNA3.1-BGH  TAG AAG GCA CAG TCG AGG  pAcGP67A-F  AAA TGA TAA CCA TCT CGC  pAcGP67A-R  GTC CAA GTT TCC CTG  Internal-F  AAC CAG GCT CTG TGA TCA with 1% PE-conjugated goat anti-mouse IgG-Fc antibody. They were then incubated for 60 min at 4 °C, gently shaking. After 4 washes, the beads were resuspended in 100 µL PBS and amount of bound antibody was quantitated using a BioPlex 200. The binding of LILRB1-D1D2-Fc to HLA class I was normalized to W6/32 binding, calculated using the following formula: (LILRB1-D1D2-Fc binding −negative control binding) / (W6/32 binding −negative control binding).

The data were analyzed statistically using the SPSS statistics package 25 (SPSS Inc. Chicago, USA). One-way analysis of variance (ANOVA) and LSD post hoc test were used to analyze the binding data, with P ≤ 0.05 between groups considered significant.

The overall goals of our study were to determine frequency differences for LILRB1 polymorphisms in worldwide populations, and to assess the effect of these polymorphism on LILRB1 binding to ligands. We analyzed individuals from the 1000 Genomes Project database to determine the number and frequency of SNPs in the coding region of LILRB1.

The results of this study were used to determine the common LILRB1 D1-D2 allotypes which were used in binding experiments.

The 1000 Genomes Project was designed to provide a comprehensive description of human genetic variation by sequencing multiple individuals from worldwide populations. Over three phases of the 1000 Genomes Project, the genomes of 2504 individuals from 26 sub-populations were sequenced and reconstructed. The average read depth for whole-genome sequencing of the 2504 individuals is 7.4X (Auton et al. 2015) . While Auton et al. have already qualified the imputation of accuracy for the whole dataset, we still want to know the sequencing quality of our targeted LILRB1 region. We expected that read depth in the coding region of LILRB1 would correlate with the file size of available mapped reads. To determine a cutoff for file size that would be sufficient for robust SNP calling, fifteen individuals with mapped.bam files ranging from 15 to 40 Gb were selected from the 1000 Genomes Phase 3 database:

HG01912, HG03432, HG02502, HG02580, HG02477, HG01956, HG02558, HG03950, HG02508, HG03085, HG02974, NA19984, HG02420, NA19900, and HG02284 were analyzed. The read depth for each coding exon (E3-E16) of LILRB1 was determined using MOSDEPTH. We found that an average read depth of about 8X, sufficient for robust SNP calling, for the complete coding region of LILRB1 was obtained when individuals had mapped.bam file sizes of 30 Gb and larger (Sup. Fig. 1A) . A representative read depth analysis is shown in Sup. Fig. 1A . Using this file size cutoff, 429 individuals were selected from 26 sub-populations (Fig. 1 ).

Analysis of the coding region of LILRB1 in the 429 selected individuals identified 149 SNPs. Of these, 66 SNPs occurred only once and were excluded from further analysis. The remaining 83 SNPs were found to be 62.7% nonsynonymous and 37.3% synonymous (Sup. Fig. 1B) . Analysis of the 83 SNPs by PHASE revealed 108 alleles, of which 46 alleles occur once, 16 alleles occur twice, and 46 alleles occur three or more times (Sup. Figs. 1C and 2). These 46 LILRB1 alleles that occur three or more times encode 28 LILRB1 allotypes and are listed in Fig. 2 . As our goal was to assess the effect of LILRB1 polymorphism on binding to HLA class I, we then grouped the 28 LILRB1 allotypes into nine LILRB1 D1-D2 variants based on the amino acid sequences of the D1-D2 region (Fig. 3A, B) for assessment by binding assay, as crystal structure analysis of LILRB1-HLA class I showed these domains are responsible for ligand interactions ).

We also calculated the frequency of these nine natural LILRB1 D1-D2 variants in populations (Fig. 3C) 

All LILRB1 D1-D2 variants except variant 8 were successfully produced and purified from Sf9 cells (Fig. 4A) . The integrity of the purified proteins was assessed by flow cytometry using APC conjugated anti-human LILRB1 antibody. The data demonstrated that all produced proteins bind to anti-human LILRB1 antibody efficiently, indicating that the D1-D2 portion is functional (Fig. 4B) . The binding of LILRB1 D1-D2 variants to HLA class I was tested by single antigen bead assay (Fig. 4C ). Variant 1 binds to HLA class I with the highest avidity. Variant 5 binds to HLA class I with higher avidity than other variants except variant 1. Variants 2, 3, 6, and 9 have similar binding avidities for HLA class I, with the binding avidity of variant 2 for HLA class I slightly lower than that of variant 9. Variants 4 and 7 bind to HLA class I with the lowest binding avidity, and there is no significant difference between their binding avidities. The detailed binding data of LILRB1 D1-D2 variants for HLA class I are shown in Sup. Fig. 3 . Although their binding avidities differ to some extent, they have a similar binding specificity for HLA class I.

We analyzed the binding to HLA-A, -B, and -C separately to assess binding differences for the LILRB1 D1-D2 variants (Fig. 5) . Our data demonstrated that all tested LILRB1 D1-D2 variants bound to HLA-C with significantly lower avidity compared to binding to HLA-A or -B. Variants 3, 4, 5, 6, 7, and 9 bound to HLA-A with slightly higher avidity than to HLA-B. However, variants 1 and 2 binding to HLA-A and HLA-B had no significant difference.

Since the α3 domain of HLA class I is responsible for interacting with the D1 domain of LILRB1 as shown by the crystal structure , we then exam- ined whether the amino acid changes in the α3 domain of HLA class I affect binding to natural LILRB1 D1-D2 variants. In total, there are 21 polymorphisms in the α3 domain for the 97 tested HLA class I antigens ( Fig. 6A-C) . Four of the polymorphisms were not dimorphic within at least one of the three HLA class I types. These polymorphisms distinguished the HLA-A, -B, and -C from each other and may be responsible for the avidity differences seen between them. All HLA-A and HLA-B have 183D, and they bind to all LILRB1 D1-D2 variants with significantly higher avidity than HLA-C with 183E (Fig. 7A ). All HLA-A have 189M, and all HLA-B and HLA-C have 189V, and the binding avidity of HLA-A with 189M is significantly higher than that of HLA-B and HLA-C with 189V (Fig. 7B ). All HLA-A and HLA-C have 239G, and HLA-B (except HLA-B*73:01) have 239R, and there is no significant difference in the binding avidity for LILRB1 D1-D2 variants between the two groups (Sup. Fig. 4A ). HLA-A (except HLA-A*80:01) and HLA-B (except HLA-B*73:01) have 268K, and all HLA-C have 268E, and the binding avidity of HLA-C with 268E for LILRB1 D1-D2 variants is significantly lower than of HLA-A and HLA-B with 268K (Fig. 7C) . For the remaining 17 polymorphisms, 11 of them occur in less than two HLA class I allotypes and were not used for statistical comparison. Finally, the remaining six polymorphisms are dimorphic within at least one of the three HLA class I types. Our data demonstrated that 194I in HLA-B is associated with higher binding to all tested LILRB1 D1-D2 variants (Fig. 7D ). 207G in HLA-A is associated with higher binding to LILRB1 variants 2, 3, 4, 6, 7, and 9 (Fig. 7E) , and 253E in HLA-A is associated with higher binding to LILRB1 variants 3, 4, 6, 7, and 9 (Fig. 7F) . However, dimorphisms at positions 184 (184P vs 184A), 193 (193A vs 193P), 194 (194V vs 194I) , and 246 (246A vs 246S) in 

The structure of LILRB1 in complex with HLA-A*02 (PDB acc. no.: 1P7Q) was retrieved and visualized in ChimeraX (Pettersen et al. 2021) . Electrostatic analysis showed that the contact sites between the LILRB1 D1-D2 domains and α3 domain of HLA-A*02 and β2m are oppositely charged (Fig. 8A-D) . The contact sites of the D1 domain of LILRB1 are positively charged, binding to the negatively charged contact sites of the α3 domain of HLA class I. However, the contact sites of the D2 domain of LILRB1 are negatively charged, binding to the positively charged contact sites of the β2m domain. It is worth noting that there are two pairs of LILRB1 D1-D2 variants that differ at only one position and at that position the amino acids have different charges (Fig. 3B) . Variant 1 and variant 7 have arginine (R) and glutamine (Q) at position 72 respectively in the D1 domain, with arginine being positively charged and glutamine being neutral at physiological pH. Similarly, variant 3 and variant 5 have aspartate (D) and histidine (H) at position 103 respectively in the D1 domain, with aspartate being negatively charged and histidine being positively charged at physiological pH. The binding data demonstrated that variant 1 with 72R binds to HLA class I with significantly higher avidity compare to variant 7 with 72Q (Fig. 8E) . Variant 3 with 103D binds to HLA class I with significantly lower avidity than variant 5 with 103H (Fig. 8F) . These data suggest that electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I.

It is unfortunate that we were unable to produce variant 8. It has two substitutions at positions 107 and 110 differen- tiating it from the more common variant 3. Variant 8 has a cysteine at position 107 compared to the positively charged arginine present in all other variants. We would predict this to reduce the binding affinity of variant 8. In addition, at position 110, variant 8 has a proline instead of the serine that is present in all other variants. Introduction of the proline could cause an alteration in the structure that may affect binding. Future comparison of the binding avidity of variant 8 to other variants will give us more information on the interaction mechanism between LILRB1 and HLA class I.

This study investigated the genetic diversity of LILRB1 in worldwide populations and assessed the binding of eight natural LILRB1 D1-D2 variants to a panel of 97 classical HLA class I antigens. We showed that LILRB1 is moderately polymorphic and the frequency of LILRB1 D1-D2 variants varies in populations. Variant 1 binds to HLA class I with the highest avidity, and variants 2, 3, 5, 6, and 9 bind to HLA class I with moderate avidity, whereas variants 4 and 7 bind to HLA class I with the lowest avidity. All tested LILRB1 D1-D2 variants bind to HLA-C with lower avidity compared with binding to HLA-A and HLA-B. Additionally, polymorphisms in the α3 domain of HLA class I influence the binding to LILRB1, partially accounting for the variation in binding avidity of LILRB1 for HLA class I. Lastly, we found that the electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I, which is evidenced by electrostatic analysis and by comparison of different binding avidities caused by polymorphisms at positions 72 and 103 of LILRB1.

We comprehensively investigated the genetic diversity of human LILRB1. Although there are a few studies genotyping the LILRB1 locus, they are limited to a small number of cohorts or specific to a certain population. Kuroki et al. screened polymorphisms in the LILRB1 locus of 18 unrelated Japanese healthy individuals and revealed 17 SNPs, including 4 nonsynonymous SNPs (L68P, A93T, I142T, and S155I) in the D1-D2 domains (Kuroki et al. 2005 ). Davidson et al. genotyped a group of 11 healthy donors and found 10 LILRB1 alleles and 85 SNPs, which is greater than that reported in Japanese populations (Davidson et al. 2010 ). More recently, Yu et al. deduced the natural LILRB1 variants by analyzing the polymorphisms L68P, A93T, I142T, and S155I in the 1000 Genomes Project and calculated their proportions in worldwide populations (Yu et al. 2018) . The limitations in that study are as follows: (1) detailed LILRB1 polymorphisms in the entire coding region and the inferred alleles were not reported; (2) only 4 nonsynonymous SNPs in the D1-D2 domain were included to infer the D1-D2 variants; (3) the haplotype frequency was only calculated for the worldwide populations; and (4) only LAIS and PTTI variants were produced and functionally examined in the further analysis. We found 149 SNPs in the coding region of LILRB1, and 108 alleles were inferred by analyzing 83 SNPs that occurred twice or more. Nine LILRB1 D1-D2 variants were further revealed by analyzing all polymorphisms in the D1-D2 domain. We also found that the frequency of LILRB1 D1-D2 variants shows population difference and variants 1-3 represent the three most frequent LILRB1 D1-D2 variants in worldwide populations. Additionally, polymorphisms in the introns and UTRs were also examined in several studies (Kuroki et al. 2005; Davidson et al. 2010; Wiśniewski et al. 2015; Yu et al. 2020 ). However, we did not include the polymorphisms in the non-coding region in our study. Our findings will enable the investigation of LILRB1 variant association to diseases within global populations and lay a foundation for functional studies.

Several studies have demonstrated that LILRB1 binds to a broad range of HLA class I Jones et al. 2011; Burshtyn and Morcos 2016) . However, whether LILRB1 binds to HLA class I with various avidities and how natural LILRB1 variants differently bind to HLA class I largely remain unknown. Previously, the binding of LILRB1 to HLA-A*01:01, HLA-A*03:01, HLA-B*07:02, HLA-B*08:01, HLA-B*15:01, HLA-B*27:02, HLA-B*27:05, HLA-B*35, HLA-B*51, and HLA-B*58 has been determined by flow cytometry analysis (Fanger et al. 1998; Cosman et al. 1997; Baía et al. 2016; Yu et al. 2018) , measuring the binding of LILRB1 protein to HLA class I transfectants. The binding of HLA-A*11, HLA-C*04, HLA-C*06:02, and HLA-C*07:02 has been determined by surface plasmon resonance (Kuroki et al. 2005; Chapman et al. 1999; Shiroishi et al. 2003) . The binding of HLA-A*02:01, HLA-G1, and HLA-F was revealed by structural analysis Wang et al. 2020; Dulberger et al. 2017) . Nevertheless, only one LILRB1 variant was tested in their studies and the number of tested HLA class I antigens is limited. Our results demonstrated that LILRB1 D1-D2 variant 1 bound to HLA class I with the highest avidity, and variants 2, 3, 5, 6, and 9 bind to HLA class I with moderate avidity, whereas variants 4 and 7 bind to HLA class I with the lowest avidity. This comprehensive comparison of the binding avidity for HLA class I between LILRB1 D1-D2 variants will greatly benefit the future functional studies on LILRB1. Jones et al. tested the interaction of one LILRB1 variant (unspecified, but most likely variant 1) to a panel of 97 HLA class I antigens, and found that LILRB1 bound to HLA class I with low variations in their binding avidity (Jones et al. 2011) . Indeed, we found that variant 1 bound to HLA class I with low variation, which is consistent with their results to some extent. However, other tested LILRB1 D1-D2 variants bound to HLA class I with higher binding variations. What's more, we found that all tested LILRB1 D1-D2 variants bound to HLA-C with the lowest avidity, and LILRB1 D1-D2 vari-ants bound to HLA-A with slightly higher avidity than to binding to HLA-B. We identified that 183E and 268E only occur in the α3 domain of HLA-C, whereas HLA-A and HLA-B have 183D and 268 K instead. These locus-specific polymorphisms may be responsible for the low binding of HLA-C to LILRB1, although they are not in direct contact (Sup. Fig. 5) . Notably, three dimorphisms in HLA-A and HLA-B are also associated with binding to in HLA-B are associated with higher binding to LILRB1. Taken together, LILRB1 D1-D2 variants bind to HLA class I with various avidities and the polymorphisms in the α3 domain of HLA class I have significantly contributed to the binding variations for LILRB1.

Variations in the binding avidity of individual HLA class I antigens for LILRB1 D1-D2 variants were also observed in the binding assay, which may alter the balance between inhibitory and activating signals within immune cells and subsequently influence the progression of pathogen infection and development of disease. We found that LILRB1 D1-D2 variants displayed higher binding to HLA-A*11:01 and A*11:02. HLA-A*11 is part of a haplotype A11-Cw4-B35-DR1-DQ1 that is associated with rapid progression of HIV (Roger 1998) . It is possible that the strong interaction of HLA-A*11:01 and A*11:02 with LILRB1 on immune cells, such as macrophages, CD8 T cells, and NK cells, results in the inhibition of these cells' anti-viral responses. In contrast, several HLA class I antigens bind to LILRB1 D1-D2 variants with weaker avidity, and reduced inhibition of immune cells may enhance the production of cytokines and subsequently influence the disease outcome. HLA-C*04:01 and HLA-C*17:01 have been reported to be associated with severe outcomes in COVID-19 (Weiner et al. 2021; Bonaccorsi et al. 2021) . We found that both HLA-C*04:01 and HLA-C*17:01 bind to LILRB1 D1-D2 variants with extremely low avidity (Sup. Fig. 3 ), which could lead to hyper-responsive APCs, CD8 T cells, and NK cells, a subsequent cytokine storm, and poorer outcomes in COVID-19. In addition, we found that LILRB1 D1-D2 variants bind to HLA-A*26:01, HLA-B*13:01, HLA-B*15:02, HLA-B*15:11, HLA-B*15:13, HLA-B*58:01, HLA-C*06:02, and HLA-C*07:02 with weaker avidity (Sup. Fig. 3 ). HLA-A*26:01 has been reported to play a predominant role in causing eye inflammations, particularly in Northeast Asian populations (Nakamura et al. 2019; Kaburaki et al. 2010) . Several studies have demonstrated that HLA-B*13:01, HLA-B*15:02, HLA-B*15:11, HLA-B*15:13, and HLA-B*58:01 are strongly associated with severe skin conditions in East Asian populations (Satapornpong et al. 2021; Chang et al. 2017; Wong et al. 2022; Dean and Kane 2013) . It is worth noting that LILRB1 D1-D2 variant 2 is enriched in Eastern Asian populations with a frequency of 46.2%, and it binds to HLA class I with weaker avidities. Thus, it is very possible that low interaction of these HLA class I antigens with LILRB1 D1-D2 variants, particularly variant 2, plays critical roles in causing severe skin conditions (or other inflammations) in East Asian populations. Lastly, both HLA-C*06:02 and HLA-C*07:02 are risk alleles associated with psoriasis (Dand et al. 2019; Li et al. 2020) ; the low binding of them to LILRB1 D1-D2 variants may partially account for it.

We found that the electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I. Electrostatic analysis revealed that the contact sites of D1 domain of LILRB1 are mainly positively charged, while the contact sites of the α3 domain of HLA class I are negatively charged. What's more, the comparison of LILRB1 D1-D2 variants that vary at position 72 (72R vs 72Q) or position 103 (103D vs 103H) revealed that positively charged residues at these two positions significantly enhance the binding to HLA class I. In addition, several studies have demonstrated that electrostatic interaction can be important to protein structures because they are potent attractions (Sheinerman et al. 2000; Froloff et al. 1997) . One good example is that the ionic charge of the loaded peptide significantly influences the stability of HLA class I antigens (Shiroishi et al. 2006; Grifoni et al. 2015) . Thus, we concluded here that electrostatic interaction is an important mechanism that used by interactions between LILRB1 and HLA class I, and this interaction can be precisely modified by mutation of charged residues at the contact sites.

A global reference for human genetic variation

Interaction of the LILRB1 inhibitory receptor with HLA class Ia dimers

Production of monoclonal antibodies to group A erythrocytes, HLA and other human cell surface antigensnew tools for genetic analysis

HLA-C*17 in COVID-19 patients: hints for associations with severe clinical outcome and cardiovascular risk

A family of human lymphoid and myeloid Ig-like receptors, some of which bind to MHC class I molecules

The expanding spectrum of ligands for leukocyte Ig-like receptors

Comparison of chimpanzee and human leukocyte Ig-like receptor genes reveals framework and rapidly evolving genes

Association of HLA-B*15:13 and HLA-B*15:02 with phenytoin-induced severe cutaneous adverse reactions in a Malay population

The inhibitory receptor LIR-1 uses a common binding interaction to recognize class I MHC molecules and the viral homolog UL18

The international Genome sample resource (IGSR): a worldwide collection of genome variation incorporating the 1000 Genomes Project data

A novel immunoglobulin superfamily receptor for cellular and viral MHC class I molecules

HLA-C*06:02 genotype is a predictive biomarker of biologic treatment response in psoriasis

polymorphism and surface phenotypes of natural killer cells

Allopurinol therapy and HLA-B*58:01 Genotype

National Center for Biotechnology Information (US)

HLAs, TCRs, and KIRs, a triumvirate of human cell-mediated immunity

Human leukocyte antigen f presents peptides and regulates immunity through interactions with NK cell receptors

The MHC class I binding proteins LIR-1 and LIR-2 inhibit Fc receptor-mediated signaling in monocytes

On the calculation of binding free energies using continuum methods: application to MHC class I protein-peptide interactions

Structural Differences in KIR3DL1 and LILRB1 Interaction with HLA-B and the loading peptide polymorphisms: in silico evidences

Resurrecting KIR2DP1: a key intermediate in the evolution of human inhibitory NK cell receptors that recognize HLA-C

The production of KIR-Fc fusion proteins and their use in a multiplex HLA class I binding assay

HLA class I allelic sequence and conformation regulate leukocyte Ig-like receptor binding

Genetic association of HLA-A*2601 with ocular Behçet's disease in Japanese patients

Extensive polymorphisms of LILRB1 (ILT2, LIR1) and their association with HLA-DRB1 shared epitope negative rheumatoid arthritis

Cross-sectional study reveals that HLA-C*07:02 is a potential biomarker of early onset/lesion severity of psoriasis

Novel HLA-G-binding leukocyte immunoglobulin-like receptor (LILR) expression patterns in human placentas and umbilical cords

Analysis of expression and function of the inhibitory receptor ILT2 (CD85j/LILRB1/LIR-1) in peripheral blood mononuclear cells from patients with systemic lupus erythematosus (SLE)

The association analysis between HLA-A*26 and Behçet's disease

Crystal structures of the two membrane-proximal Ig-like domains (D3D4) of LILRB1/B2: alternative models for their involvement in peptide-HLA binding

Mosdepth: quick coverage calculation for genomes and exomes

UCSF ChimeraX: structure visualization for researchers, educators, and developers

Influence of host genes on HIV-1 disease progression

Immune evasion of Plasmodium falciparum by RIFIN via inhibitory receptors

HLA-B*13:01 is a predictive marker of dapsone-induced severe cutaneous adverse reactions in Thai patients

Electrostatic aspects of protein-protein interactions

Efficient leukocyte Ig-like receptor signaling and crystal structure of disulfide-linked HLA-G dimer

Human inhibitory receptors Ig-like transcript 2 (ILT2) and ILT4 compete with CD8 for MHC class I binding and bind preferentially to HLA-G

A new statistical method for haplotype reconstruction from population data

Structures of the four Ig-like domain LILRB2 and the four-domain LILRB1 and HLA-G1 complex

Increased risk of severe clinical course of COVID-19 in carriers of HLA-C*04:01

Organization of the leukocyte receptor cluster (LRC) on human chromosome 19q13

Crystal structure of HLA-A2 bound to LIR-1, a host and viral major histocompatibility complex receptor

Genetic polymorphisms and expression of HLA-G and its receptors, KIR2DL4 and LILRB1, in non-small cell lung cancer

HLA-B*15:11 status and carbamazepineinduced severe cutaneous adverse drug reactions in HLA-B*15:02 negative Chinese

Structure of UL18, a peptide-binding viral MHC mimic, bound to a host inhibitory receptor

LILRB1 intron 1 has a polymorphic regulatory region that enhances transcription in NK cells and recruits YY1

polymorphisms influence posttransplant HCMV susceptibility and ligand interactions

We would like to thank Dr. Rahul D Pawar and Dr. Jeanette Baker for advice on performing the single antigen beads assay, Dr. Jinfan Wang for critical discussion, and Niassan Beyzaie for advice on data extraction and analysis.Funding This work was supported by an NIH grant (R01 AI136952).

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00251-022-01264-7.