Identification of the hybrid gene LILRB5-3 by long-read sequencing and implication of its novel signaling function.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) on human chromosome 19q13.4 encode 11 immunoglobulin superfamily receptors, exhibiting genetic diversity within and between human populations. Among the LILR genes, the genomic region surrounding LILRB3 and LILRA6 has yet to be fully characterized due to their significant sequence homology, which makes it difficult to differentiate between them. To examine the LILRB3 and LILRA6 genomic region, a tool named JoGo-LILR CN Caller, which can call copy number from short-read whole genome sequencing (srWGS) data, was applied to an extensive international srWGS dataset comprising 2,504 samples. During this process, a previously unreported loss of both LILRB3 and LILRA6 was detected in three samples. Using long-read sequencing of these samples, we have discovered a novel large deletion (33,692 bp) in the LILRB3 and LILRA6 genomic regions in the Japanese population. This deletion spanned three genes, LILRB3, LILRA6, and LILRB5, resulting in LILRB3 exons 12-13 being located immediately downstream of LILRB5 exons 1-12 with the loss of LILRA6, suggesting the potential expression of a hybrid gene between LILRB5 and LILRB3 (LILRB5-3). Transcription and subsequent translation of the LILRB5-3 hybrid gene were also verified. The hybrid junction was located within the intracellular domain, resulting in an LILRB5 extracellular domain fused to a partial LILRB3 intracellular domain with three immunoreceptor tyrosine-based inhibitory motifs (ITIMs), suggesting that LILRB5-3 acquired a novel signaling function. Further application of the JoGo-LILR tool to srWGS samples suggested the presence of the LILRB5-3 hybrid gene in the CEU population. Our findings provide insight into the genetic and functional diversity of the LILR family.

Characterization of LILRB3 and LILRA6 allelic variants in the Japanese population.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) are encoded by members of a human multigene family, comprising 11 protein-coding genes and two pseudogenes. Among the LILRs, LILRB3 and LILRA6 show the highest homology with each other, along with high allelic and copy number variations. Therefore, it has been difficult to discriminate between them, both genetically and functionally, precluding disease association studies of LILRB3 and LILRA6. In this study, we carefully performed variant screening of LILRB3 and LILRA6 by cDNA cloning from Japanese individuals and identified four allelic lineages showing significantly high non-synonymous-to-synonymous ratios in pairwise comparisons. Furthermore, the extracellular domains of the LILRB3*JP6 and LILRA6*JP1 alleles were identical at the DNA level, suggesting that gene conversion-like events diversified LILRB3 and LILRA6. To determine the four allelic lineages from genomic DNA, we established a lineage typing method that accurately estimated the four allelic lineages in addition to specific common alleles from genomic DNA. Analysis of LILRA6 copy number variation revealed one, two, and three copies of LILRA6 in the Japanese-in-Tokyo (JPT) population. Flow cytometric analysis showed that an anti-LILRB3 antibody did not recognize the second most common lineage in the Japanese population, indicating significant amino acid differences across the allelic lineages. Taken together, our findings indicate that our lineage typing is useful for classifying the lineage-specific functions of LILRB3 and LILRA6, serving as the basis for disease association studies.