Exceptional diversity of KIR and HLA class I in Egypt

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis

Background

The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods.

Methods

We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls.

Results

We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk.

Results

The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles

OBJECTIVES

Drug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still's disease with atypical lung disease. We sought to characterise features of patients with Still's disease with DRESS compared with drug-tolerant Still's controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort.

METHODS

In a case/control study, we collected a multicentre series of patients with Still's disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still's controls (n=65). We retrospectively analysed clinical data from all Still's subjects and typed 94/131 for HLA. European Still's-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still's cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still's-DRESS cases (n=64) compared with drug-tolerant Still's controls (n=30). KD subjects (n=19) were similarly studied.

RESULTS

Still's-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still's-DRESS (64%) versus drug-tolerant Still's (3%; p=1.2×10-14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still's-DRESS cases versus INCHARGE Still's controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still's controls (p=6.3×10-10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions.

CONCLUSIONS

DRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.

High-resolution HLA allele and haplotype frequencies in several unrelated populations determined by next generation sequencing: 17th International HLA and Immunogenetics Workshop joint report

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity

The KIR (killer-cell immunoglobulin-like receptor) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies.

Behçet disease, new insights in disease associations and manifestations: a next-generation sequencing study

Behçet disease is a multi-system disease associated with human leukocyte antigen (HLA) class I polymorphism. High-resolution next-generation sequencing (NGS) with haplotype analysis has not been performed previously for this disease. Sixty Egyptian patients diagnosed according to the International Study Group (ISG) criteria for Behçet disease and 160 healthy geographic and ethnic-matched controls were genotyped for HLA class I loci (HLA-A, B, C). For HLA class II loci (DRB1, DRB3/4/5, DQA1, DQB1, DPA1, DPB1), 40 control samples were genotyped. High-resolution HLA genotyping was performed using NGS and the results were analyzed. Clinical manifestations were oral ulcers (100%), genital ulcers (100%), eye (55%) and neurological (28%) and vascular involvement (35%). HLA-B*51:08 [odds ratio (OR) = 19·75, 95% confidence interval (CI) = 6·5-79; P < 0·0001], HLA-B*15:03 (OR = 12·15, 95% CI = 3·7-50·7; P < 0·0001), HLA-C*16:02 (OR = 6·53, 95% CI = 3-14; P < 0·0001), HLA-A*68:02 (OR = 3·14, 95% CI = 1·1-8·9; P < 0·01) were found to be associated with Behçet disease, as were HLA-DRB1*13:01 and HLA-DQB1*06:03 (OR = 3·39, 95% CI = 0·9-18·9; P = 0·04 for both). By contrast, HLA-A*03:01 (OR = 0·13, 95% CI = 0-0·8; P = 0·01) and HLA-DPB1*17:01 were found to be protective (OR = 0·27, 95% CI = 0·06-1·03; P = 0·02). We identified strong linkage disequilibrium between HLA-B*51:08 and C*16:02 and A*02:01 in a haplotype associated with Behçet disease. HLA-B*51:08 was significantly associated with legal blindness (OR = 2·98, 95% CI = 1·06-8·3; P = 0·01). In Egyptian Behçet patients, HLA-B*51:08 is the most common susceptibility allele and holds poor prognosis for eye involvement.

High-throughput SARS-CoV-2 and host genome sequencing from single nasopharyngeal swabs

During COVID19 and other viral pandemics, rapid generation of host and pathogen genomic data is critical to tracking infection and informing therapies. There is an urgent need for efficient approaches to this data generation at scale. We have developed a scalable, high throughput approach to generate high fidelity low pass whole genome and HLA sequencing, viral genomes, and representation of human transcriptome from single nasopharyngeal swabs of COVID19 patients.

Multiplex family with GAD65-Abs neurologic syndromes

Objective:

Neurologic autoimmune syndromes associated with anti–glutamate acid decarboxylase 65 antibodies (GAD65-Abs) are rare and mostly sporadic.

Methods:

We describe a niece and her aunt with GAD65-Abs neurologic syndromes. High-resolution HLA typing of Class I and Class II alleles was performed using next-generation sequencing.

Results:

The proband had cerebellar ataxia and probable limbic encephalitis features, whereas her niece had stiff-person syndrome. Both had a high titer of GAD65-Abs in serum and CSF and showed signs of inflammation in CSF. Both affected members carried the same rare recombinant DRB1*15:01:01∼DQA1*01:02:01∼DQB1*05:02:01 haplotype, which may or may not be involved in disease susceptibility. Of interest, other unaffected members of the family either had the same HLA haplotype but normal serum GAD65-Abs or had different HLA types but a high titer of serum GAD65-Abs without neurologic symptoms, suggesting cumulative effects.

Conclusions:

This unique association strengthens the concept that hereditary factors, possibly including specific HLA haplotypes, play a role in neurologic syndromes associated with GAD65-Abs.