Genetic diversity of the LILRB1 and LILRB2 coding regions in an admixed Brazilian population sample

Abstract

Leukocyte Immunoglobulin (Ig)-like Receptors (LILR) LILRB1 and LILRB2 play a pivotal role in maintaining self-tolerance and modulating the immune response through interaction with classical and non-classical Human Leukocyte Antigen (HLA) molecules. Although both diversity and natural selection patterns over HLA genes have been extensively evaluated, little information is available concerning the genetic diversity and selection signatures on the LIRB1/2 regions. Therefore, we identified the LILRB1/2 genetic diversity using next-generation sequencing in a population sample comprising 528 healthy control individuals from Sao Paulo State, Brazil. We identified 58 LILRB1 Single Nucleotide Variants (SNVs), which gave rise to 13 haplotypes with at least 1% of frequency. For LILRB2, we identified 41 SNVs arranged into 11 haplotypes with frequencies above 1%. We found evidence of either positive or purifying selection on LILRB1/2 coding regions. Some residues in both proteins showed to be under the effect of positive selection, suggesting that amino acid replacements in these proteins resulted in beneficial functional changes. Finally, we have shown that allelic variation (six and five amino acid exchanges in LILRB1 and LILRB2, respectively) affects the structure and/or stability of both molecules. Nonetheless, LILRB2 has shown higher average stability, with no D1/D2 residue affecting protein structure. Taken together, our findings demonstrate that LILRB1 and LILRB2 are highly polymorphic and provide strong evidence supporting the directional selection regime hypothesis.