Distinct transcriptome architectures underlying lupus establishment and exacerbation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease involving multiple immune cells. To elucidate SLE pathogenesis, it is essential to understand the dysregulated gene expression pattern linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 136 SLE and 89 healthy donors. We profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We then identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and therapeutic responses. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that current genetic studies may not well capture clinically vital biology. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic and genetic studies.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases.

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

Analysis of antibody binding specificities in twin and SNP-genotyped cohorts reveals that antiviral antibody epitope selection is a heritable trait.

Human immune responses to viral infections are highly variable, but the genetic factors that contribute to this variability are not well characterized. We used VirScan, a high-throughput epitope scanning technology, to analyze pan-viral antibody reactivity profiles of twins and SNP-genotyped individuals. Using these data, we determined the heritability and genomic loci associated with antibody epitope selection, response breadth, and control of Epstein-Barr virus (EBV) viral load. 107 EBV peptide reactivities were heritable and at least two Epstein-Barr nuclear antigen 2 (EBNA-2) reactivities were associated with variants in the MHC class II locus. We identified an EBV serosignature that predicted viral load in peripheral blood mononuclear cells and was associated with variants in the MHC class I locus. Our study illustrates the utility of epitope profiling to investigate the genetics of pathogen immunity, reports heritable features of the antibody response to viruses, and identifies specific HLA loci important for EBV epitope selection.

Phage display of environmental protein toxins and virulence factors reveals the prevalence, persistence, and genetics of antibody responses.

Microbial exposures are crucial environmental factors that impact healthspan by sculpting the immune system and microbiota. Antibody profiling via Phage ImmunoPrecipitation Sequencing (PhIP-Seq) provides a high-throughput, cost-effective approach for detecting exposure and response to microbial protein products. We designed and constructed a library of 95,601 56-amino acid peptide tiles spanning 14,430 proteins with "toxin" or "virulence factor" keyword annotations. We used PhIP-Seq to profile the antibodies of ∼1,000 individuals against this "ToxScan" library. In addition to enumerating immunodominant antibody epitopes, we studied the age-dependent stability of the ToxScan profile and used a genome-wide association study to find that the MHC-II locus modulates bacterial epitope selection. We detected previously described anti-flagellin antibody responses in a Crohn's disease cohort and identified an association between anti-flagellin antibodies and juvenile dermatomyositis. PhIP-Seq with the ToxScan library is thus an effective tool for studying the environmental determinants of health and disease at cohort scale.

Joint host-pathogen genomic analysis identifies hepatitis B virus mutations associated with human NTCP and HLA class I variation.

Evolutionary changes in the hepatitis B virus (HBV) genome could reflect its adaptation to host-induced selective pressure. Leveraging paired human exome and ultra-deep HBV genome-sequencing data from 567 affected individuals with chronic hepatitis B, we comprehensively searched for the signatures of this evolutionary process by conducting "genome-to-genome" association tests between all human genetic variants and viral mutations. We identified significant associations between an East Asian-specific missense variant in the gene encoding the HBV entry receptor NTCP (rs2296651, NTCP S267F) and mutations within the receptor-binding region of HBV preS1. Through in silico modeling and in vitro preS1-NTCP binding assays, we observed that the associated HBV mutations are in proximity to the NTCP variant when bound and together partially increase binding affinity to NTCP S267F. Furthermore, we identified significant associations between HLA-A variation and viral mutations in HLA-A-restricted T cell epitopes. We used in silico binding prediction tools to evaluate the impact of the associated HBV mutations on HLA presentation and observed that mutations that result in weaker binding affinities to their cognate HLA alleles were enriched. Overall, our results suggest the emergence of HBV escape mutations that might alter the interaction between HBV PreS1 and its cellular receptor NTCP during viral entry into hepatocytes and confirm the role of HLA class I restriction in inducing HBV epitope variations.

Homozygous variant in TKFC abolishing triokinase activities is associated with isolated immunodeficiency.

BACKGROUND: Triokinase and FMN cyclase (TKFC) is a bifunctional enzyme involved in fructose metabolism. Triokinase catalyses the phosphorylation of fructose-derived glyceraldehyde (GA) and exogenous dihydroxyacetone (DHA), while FMN cyclase generates cyclic FMN. TKFC regulates the antiviral immune response by interacting with IFIH1 (MDA5). Previously reported pathogenic variants in TKFC are associated with either a multisystemic disease or isolated hypotrichosis with loose anagen hairs. METHODS: Whole-exome sequencing identified a homozygous novel variant in TKFC (c.1624G>A; p.Gly542Arg) in an individual with a complex primary immunodeficiency disorder. The variant was characterised using enzymatic assays and yeast studies of mutant recombinant proteins. RESULTS: The individual presented with chronic active Epstein-Barr virus disease and multiple bacterial and viral infections. Clinical investigations revealed hypogammaglobulinaemia, near absent natural killer cells and decreased memory B cells. Enzymatic assays showed that this variant displayed defective DHA and GA kinase activity while maintaining FMN cyclase activity. An allogenic bone marrow transplantation corrected the patient's immunodeficiency. CONCLUSION: Our report suggests that TKFC may have a role in the immunological system. The pathological features associated with this variant are possibly linked with DHA/GA kinase inactivation through a yet an unknown mechanism. This report thus adds a possible new pathway of immunometabolism to explore further.

GenIA, the Genetic Immunology Advisor database for inborn errors of immunity.

BACKGROUND: To date, no publicly accessible platform has captured and synthesized all of the layered dimensions of genotypic, phenotypic, and mechanistic information published in the field of inborn errors of immunity (IEIs). Such a platform would represent the extensive and complex landscape of IEIs and could increase the rate of diagnosis in patients with a suspected IEI, which remains unacceptably low. OBJECTIVE: Our aim was to create an expertly curated, patient-centered, multidimensional IEI database that enables aggregation and sophisticated data interrogation and promotes involvement from diverse stakeholders across the community. METHODS: The database structure was designed following a subject-centered model and written in Structured Query Language (SQL). The web application is written in Hypertext Preprocessor (PHP), Hypertext Markup Language (HTML), Cascading Style Sheets (CSS), and JavaScript. All data stored in the Genetic Immunology Advisor (GenIA) are extracted by manually reviewing published research articles. RESULTS: We completed data collection and curation for 24 pilot genes. Using these data, we have exemplified how GenIA can provide quick access to structured, longitudinal, more thorough, comprehensive, and up-to-date IEI knowledge than do currently existing databases, such as ClinGen, Human Phenotype Ontology (HPO), ClinVar, or Online Mendelian Inheritance in Man (OMIM), with which GenIA intends to dovetail. CONCLUSIONS: GenIA strives to accurately capture the extensive genetic, mechanistic, and phenotypic heterogeneity found across IEIs, as well as genetic paradigms and diagnostic pitfalls associated with individual genes and conditions. The IEI community's involvement will help promote GenIA as an enduring resource that supports and improves knowledge sharing, research, diagnosis, and care for patients with genetic immune disease.

HLA-DPA1*02:01~B1*01:01 is a risk haplotype for primary sclerosing cholangitis mediating activation of NKp44+ NK cells.

OBJECTIVE: Primary sclerosing cholangitis (PSC) is characterised by bile duct strictures and progressive liver disease, eventually requiring liver transplantation. Although the pathogenesis of PSC remains incompletely understood, strong associations with HLA-class II haplotypes have been described. As specific HLA-DP molecules can bind the activating NK-cell receptor NKp44, we investigated the role of HLA-DP/NKp44-interactions in PSC. DESIGN: Liver tissue, intrahepatic and peripheral blood lymphocytes of individuals with PSC and control individuals were characterised using flow cytometry, immunohistochemical and immunofluorescence analyses. HLA-DPA1 and HLA-DPB1 imputation and association analyses were performed in 3408 individuals with PSC and 34 213 controls. NK cell activation on NKp44/HLA-DP interactions was assessed in vitro using plate-bound HLA-DP molecules and HLA-DPB wildtype versus knock-out human cholangiocyte organoids. RESULTS: NKp44+NK cells were enriched in livers, and intrahepatic bile ducts of individuals with PSC showed higher expression of HLA-DP. HLA-DP haplotype analysis revealed a highly elevated PSC risk for HLA-DPA1*02:01~B1*01:01 (OR 1.99, p=6.7×10-50). Primary NKp44+NK cells exhibited significantly higher degranulation in response to plate-bound HLA-DPA1*02:01-DPB1*01:01 compared with control HLA-DP molecules, which were inhibited by anti-NKp44-blocking. Human cholangiocyte organoids expressing HLA-DPA1*02:01-DPB1*01:01 after IFN-γ-exposure demonstrated significantly increased binding to NKp44-Fc constructs compared with unstimulated controls. Importantly, HLA-DPA1*02:01-DPB1*01:01-expressing organoids increased degranulation of NKp44+NK cells compared with HLA-DPB1-KO organoids. CONCLUSION: Our studies identify a novel PSC risk haplotype HLA-DP A1*02:01~DPB1*01:01 and provide clinical and functional data implicating NKp44+NK cells that recognise HLA-DPA1*02:01-DPB1*01:01 expressed on cholangiocytes in PSC pathogenesis.

Progression to type 1 diabetes in the DPT-1 and TN07 clinical trials is critically associated with specific residues in HLA-DQA1-B1 heterodimers.

AIMS/HYPOTHESIS: The aim of this work was to explore molecular amino acids (AAs) and related structures of HLA-DQA1-DQB1 that underlie its contribution to the progression from stages 1 or 2 to stage 3 type 1 diabetes. METHODS: Using high-resolution DQA1 and DQB1 genotypes from 1216 participants in the Diabetes Prevention Trial-Type 1 and the Diabetes Prevention Trial, we applied hierarchically organised haplotype association analysis (HOH) to decipher which AAs contributed to the associations of DQ with disease and their structural properties. HOH relied on the Cox regression to quantify the association of DQ with time-to-onset of type 1 diabetes. RESULTS: By numerating all possible DQ heterodimers of α- and ß-chains, we showed that the heterodimerisation increases genetic diversity at the cellular level from 43 empirically observed haplotypes to 186 possible heterodimers. Heterodimerisation turned several neutral haplotypes (DQ2.2, DQ2.3 and DQ4.4) to risk haplotypes (DQ2.2/2.3-DQ4.4 and DQ4.4-DQ2.2). HOH uncovered eight AAs on the α-chain (-16α, -13α, -6α, α22, α23, α44, α72, α157) and six AAs on the ß-chain (-18ß, ß9, ß13, ß26, ß57, ß135) that contributed to the association of DQ with progression of type 1 diabetes. The specific AAs concerned the signal peptide (minus sign, possible linkage to expression levels), pockets 1, 4 and 9 in the antigen-binding groove of the α1ß1 domain, and the putative homodimerisation of the αß heterodimers. CONCLUSIONS/INTERPRETATION: These results unveil the contribution made by DQ to type 1 diabetes progression at individual residues and related protein structures, shedding light on its immunological mechanisms and providing new leads for developing treatment strategies. DATA AVAILABILITY: Clinical trial data and biospecimen samples are available through the National Institute of Diabetes and Digestive and Kidney Diseases Central Repository portal ( https://repository.niddk.nih.gov/studies ).

A pan-cancer analysis of prognostic significance and immunological role of lysosomal-associated membrane protein 3.

Lysosomal dysfunction can drive carcinogenesis. Lysosomal-associated membrane protein 3 (LAMP3), is a member of the Lysosome Associated Membrane Proteins and is involved in the malignant phenotype such as tumour metastasis and drug resistance, while the mechanisms that regulate the malignant progression of tumour remain vague. Our study aims to provide a more systematic and comprehensive understanding of the role of LAMP3 in the progression of various cancers by various databases.We explored the role of LAMP3 in pan-cancer using The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database. Multiple online web platforms and software were used for data analysis, including HPA, TIMER, TISIDB, GEPIA, UALCAN, Kaplan-Meier plotter, DAVID and TIGER. The immunohistochemistry was used to quantify the LAMP3 and PD-L1 expression levels in cancer.High LAMP3 expression was found in most cancers and differentially expressed across molecular and immune subtypes. The expression of LAMP3 was involved in the immune-associated processes of Antigen processing and presentation, Th17 cell differentiation, Th1 and Th2 cell differentiation, and the immune-associated pathways of T cell receptor and B cell receptor signalling pathways in most cancers. It also correlated with genetic markers of immunomodulators in various cancers. LAMP3 and PD-L1 expression in BRCA and HNSC tissues was higher than that in corresponding adjacent normal tissues by immunohistochemistry. There is a significant correlation between the expression of LAMP3 and PD-L1.Our study elucidates that LAMP3 has different expression patterns and genetic alteration patterns in different tumours. It is a potential biomarker for immune-related cancer diagnosis, prognosis and efficacy prediction.

Selection and Interpretation of Molecular Diagnostics in Heart Transplantation.

The number of heart transplants performed annually in the United States and worldwide continues to increase, but there has been little change in graft longevity and patient survival over the past 2 decades. The reference standard for diagnosis of acute cellular and antibody-mediated rejection includes histologic and immunofluorescence evaluation of endomyocardial biopsy samples, despite invasiveness and high interrater variability for grading histologic rejection. Circulating biomarkers and molecular diagnostics have shown substantial predictive value in rejection monitoring, and emerging data support their use in diagnosing other posttransplant complications. The use of genomic (cell-free DNA), transcriptomic (mRNA and microRNA profiling), and proteomic (protein expression quantitation) methodologies in diagnosis of these posttransplant outcomes has been evaluated with varying levels of evidence. In parallel, growing knowledge about the genetically mediated immune response leading to rejection (immunogenetics) has enhanced understanding of antibody-mediated rejection, associated graft dysfunction, and death. Antibodies to donor human leukocyte antigens and the technology available to evaluate these antibodies continues to evolve. This review aims to provide an overview of biomarker and immunologic tests used to diagnose posttransplant complications. This includes a discussion of pediatric heart transplantation and the disparate rates of rejection and death experienced by Black patients receiving a heart transplant. This review describes diagnostic modalities that are available and used after transplant and the landscape of future investigations needed to enhance patient outcomes after heart transplantation.

Systematic characterization of immunoglobulin loci and deep sequencing of the expressed repertoire in the Atlantic cod (Gadus morhua).

BACKGROUND: The Atlantic cod is a prolific species in the Atlantic, despite its inconsistent specific antibody response. It presents a peculiar case within vertebrate immunology due to its distinct immune system, characterized by the absence of MHCII antigen presentation pathway, required for T cell-dependent antibody responses. Thorough characterisation of immunoglobulin loci and analysis of the antibody repertoire is necessary to further our understanding of the Atlantic cod's immune response on a molecular level. RESULTS: A comprehensive search of the cod genome (gadmor3.0) identified the complete set of IgH genes organized into three sequential translocons on chromosome 2, while IgL genes were located on chromosomes 2 and 5. The Atlantic cod displayed a moderate germline V gene diversity, comprising four V gene families for both IgH and IgL, each with distinct chromosomal locations and organizational structures. 5'RACE sequencing revealed a diverse range of heavy chain CDR3 sequences and relatively limited CDR3 diversity in light chains. The analysis highlighted a differential impact of V-gene germline CDR3 length on receptor CDR3 length between heavy and light chains, underlining different recombination processes. CONCLUSIONS: This study reveals that the Atlantic cod, despite its inconsistent antibody response, maintains a level of immunoglobulin diversity comparable to other fish species. The findings suggest that the extensive recent duplications of kappa light chain genes do not result in increased repertoire diversity. This research provides a comprehensive view of the Atlantic cod's immunoglobulin gene organization and repertoire, necessary for future studies of antibody responses at the molecular level.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

BACKGROUND: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). RESULTS: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages. CONCLUSIONS: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

The genetics of Graves' disease.

Graves' disease (GD) is the commonest cause of hyperthyroidism and has a strong female preponderance. Everyday clinical practice suggests strong aggregation within families and twin studies demonstrate that genetic factors account for 60-80% of risk of developing GD. In this review, we collate numerous genetic studies and outline the discoveries over the years, starting with historic candidate gene studies and then exploring more recent genome-wide linkage and association studies, which have involved substantial cohorts of East Asian patients as well as those of European descent. Variants in genes including HLA, CTLA4, and PTPN22 have been shown to have substantial individual effects on disease susceptibility. In addition, we examine emerging evidence concerning the possibility that genetic variants may correlate with relevant clinical phenotypes including age of onset of GD, severity of thyrotoxicosis, goitre size and relapse of hyperthyroidism following antithyroid drug therapy, as well as thyroid eye disease. This review supports the inheritance of GD as a complex genetic trait, with a growing number of more than 80 susceptibility loci identified so far. Future implementation of more targeted clinical therapies requires larger studies investigating the influence of these genetic variants on the various phenotypes and different outcomes of conventional treatments.

Predicting flow cytometry crossmatch results from single-antigen bead testing.

The aim of this study was to devise an algorithm that would predict flow cytometry crossmatch (FCXM) results using single-antigen bead (SAB) mean fluorescent intensity (MFI) levels using samples received through the National External Quality Assurance Scheme (NEQAS) 2B external proficiency testing scheme between 2019 and 2023. A total of 159 serum samples were retrospectively screened using LABScreen Single Antigen Class I and II (SAB), and 40 peripheral blood samples were human leucocyte antigen (HLA) typed with LABType SSO. Donor-specific antibodies were identified for each cell-serum combination tested, and cumulative MFI values were calculated for each test before correlating the screening result with the consensus crossmatch results for this scheme. HLA Class I MFIs were combined to predict the T cell crossmatch. For the B cell crossmatch prediction, two options were considered: (i) HLA Class II MFI values alone and (ii) HLA Class I + Class II MFIs. Receiver operating characteristic analysis was carried out to identify the combined MFI threshold that predicted NEQAS consensus results with the greatest sensitivity and specificity. HLA Class I combined MFI >5000 predicted T cell crossmatch results with 96% sensitivity, 100% specificity, 100% positive predictive value (PPV) and 92% negative predictive value (NPV). For B cell results, HLA Class I + Class II combined MFIs >11,000 gave the best model, showing 97% sensitivity, 82% specificity, 96% PPV and 85% NPV. However, for samples with only HLA Class II sensitization, combined MFIs >13,000 improved the B cell crossmatch predictions: 92% sensitivity, 95% specificity, 96% PPV and 91% NPV. Using this model, combined MFI can be used to predict the immunological risk posed by donor-specific antibodies when it is not possible to carry out an FCXM.

Forming a new perspective: Post-structural approaches to determination of donor compatibility and post-transplant assessment of allograft health.

The purpose of this review is to encourage a new perspective on the question of donor-recipient compatibility and post-transplant assessment of graft health based on functional measures. The premise is that we should be better sighted on what (and how) the immune system responds toward rather than what is merely there. Continuance of the pursuit of further and better definition of antigens and antibodies is not however discouraged but seen as necessary to improved understanding of the structural correlates of functional immunity. There currently exists, in the opinion of the authors, an opportunity for histocompatibility and immunogenetics laboratories to develop and widen their scope of involvement into these new areas of laboratory activity in support and to the benefit of the transplant programmes they serve.

Saliva direct PCR protocol for HLA-DQB1*02 genotyping.

Celiac disease (CD) is an immune disorder, that is triggered by gluten ingestion in genetically predisposed individuals. The HLA-DQB1*02 allele is the main predisposing genetic factor and a candidate for first-line genotyping screening. We designed and validated a simple, DNA purification-free PCR protocol directly from crude saliva, enabling the detection of the DQB1*02 allele. This assay also distinguishes homozygous from heterozygous carriers. We propose this method for use in mass screening and/or epidemiological studies.

UK NEQAS and BSHI guideline: Laboratory testing and clinical interpretation of HLA genotyping results supporting the diagnosis of coeliac disease.

Coeliac disease is a common immune-mediated inflammatory disorder caused by dietary gluten in genetically susceptible individuals. While the diagnosis of coeliac disease is based on serological and histological criteria, HLA-DQ genotyping can be useful, especially in excluding the diagnosis in patients who do not carry the relevant DQ heterodimers: DQA1*05 DQB1*02, DQB1*03:02 or DQA1*02 DQB1*02 (commonly referred to as DQ2.5, DQ8 and DQ2.2, respectively). External quality assessment results for HLA genotyping in coeliac disease have revealed concerning errors in HLA genotyping, reporting and clinical interpretation. In response, these guidelines have been developed as an evidence-based approach to guide laboratories undertaking HLA genotyping for coeliac disease and provide recommendations for reports to standardise and improve the communication of results.

Homozygous HLA-DQB1*06:02 combined with T-cell receptor alpha polymorphism results in narcolepsy onset - A familial case report.

Narcolepsy is a life-long neurological disorder with well-established genetic risk factors. Human leukocyte antigen-DQB1*06:02 remains the strongest genetic predeterminant; however, polymorphisms in genes encoding the T-cell receptor alpha chain are also strongly linked. This case report shows the inheritance pathway of these genetic markers contributing to narcolepsy onset in a 17-year-old female.

Leukocyte immunoglobulin-like receptor A3 gene deletion in five Chinese populations and protective association with nasopharyngeal carcinoma.

Among the thirteen leukocyte Ig-like receptor (LILR) loci located at 19q13.4, LILRA3 is unique in that it encodes a soluble protein lacking the transmembrane and cytoplasmic domains, and a 6.7 kb deletion spanning the first seven exons has been detected in some human individuals. Presently, there is a lack of data about the distribution of LILRA3 gene deletion in more diverse ethnic groups. Also, no previous studies have investigated the correlation between copy number variation (CNV) of LILRA3 and nasopharyngeal carcinoma (NPC). In this study, five populations from China mainland: two Southern Han populations, Hunan (N = 1478) and Guandong (N = 107); one Southeastern Han population, Fujian (N = 439); and two Northern populations, Inner Mongolia Han (N = 104) and Mongol population from Inner Mongolia (N = 158) were investigated for CNV of LILRA3 using polymerase chain reaction-sequence-specific priming (PCR-SSP) method. LILRA3 variants were also examined in a cohort of NPC cases (N = 1142) in Hunan Han population. The five Chinese populations demonstrated northward increase in frequency of the deleted form of LILRA3 gene (LILRA3*Del) (all corrected p values < 0.05). Inter-population comparison also uncovered significant differentiation in the distribution of CNV of LILRA3 among modern human populations. LILRA3*Del was found to confer significantly reduced risk to NPC in Hunan Han population (at allelic level: OR = 0.79, 95% CI = 0.71-0.89, p < 0.0001; at genotype level: OR = 0.63, 95% CI = 0.51-0.79, p < 0.0001). No interaction was found between LILRA3 variants and HLA-A*02:07, HLA-A*11:01, HLA-B*13 and HLA-B*46:01 alleles in susceptibility to NPC. Our study constitutes the first demonstration of LILRA3 gene as a locus linked to NPC susceptibility in a southern Chinese population. Future independent studies in other populations are warranted to confirm the findings reported in this study.