Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations

Genome-wide association analysis reveals the associations of NPHP4, TYW1-AUTS2 and SEMA6D for Behçet's disease and HLA-B*46:01 for its intestinal involvement

BACKGROUND

Intestinal involvement in Behçet's disease (BD) is associated with poor prognosis and is more prevalent in East Asian than in Mediterranean populations. Identifying the genetic causes of intestinal BD is important for understanding the pathogenesis and for appropriate treatment of BD patients.

METHODS

We performed genome-wide association studies (GWAS) and imputation/replication genotyping of human leukocyte antigen (HLA) alleles for 1,689 Korean and Turkish patients with BD (including 379 patients with intestinal BD) and 2,327 healthy controls, followed by replication using 593 Japanese patients with BD (101 patients with intestinal BD) and 737 healthy controls. Stratified cross-phenotype analyses were performed for 1) overall BD, 2) intestinal BD, and 3) intestinal BD without association of overall BD.

RESULTS

We identified three novel genome-wide significant susceptibility loci including NPHP4 (rs74566205; P=1.36 × 10-8), TYW1-AUTS2 (rs60021986; P=1.14 × 10-9), and SEMA6D (rs4143322; P=5.54 × 10-9) for overall BD, and a new association with HLA-B*46:01 for intestinal BD (P=1.67 × 10-8) but not for BD without intestinal involvement. HLA peptide binding analysis revealed that Mycobacterial peptides, have a stronger binding affinity to HLA-B*46:01 compared to the known risk allele HLA-B*51:01.

CONCLUSIONS

HLA-B*46:01 is associated with the development of intestinal BD; NPHP4, TYW1-AUTS2, and SEMA6D are susceptibility loci for overall BD.

Regulatory SVA retrotransposons and classical HLA genotyped-transcripts associated with Parkinson's disease

Introduction

Parkinson's disease (PD) is a neurodegenerative and polygenic disorder characterised by the progressive loss of neural dopamine and onset of movement disorders. We previously described eight SINE-VNTR-Alu (SVA) retrotransposon-insertion-polymorphisms (RIPs) located and expressed within the Human Leucocyte Antigen (HLA) genomic region of chromosome 6 that modulate the differential co-expression of 71 different genes including the HLA classical class I and class II genes in a Parkinson's Progression Markers Initiative (PPMI) cohort.

Aims and methods

In the present study, we (1) reanalysed the PPMI genomic and transcriptomic sequencing data obtained from whole blood of 1521 individuals (867 cases and 654 controls) to infer the genotypes of the transcripts expressed by eight classical HLA class I and class II genes as well as DRA and the DRB3/4/5 haplotypes, and (2) examined the statistical differences between three different PD subgroups (cases) and healthy controls (HC) for the HLA and SVA transcribed genotypes and inferred haplotypes.

Results

Significant differences for 57 expressed HLA alleles (21 HLA class I and 36 HLA class II alleles) up to the three-field resolution and four of eight expressed SVA were detected at p<0.05 by the Fisher's exact test within one or other of three different PD subgroups (750 individuals with PD, 57 prodromes, 60 individuals who had scans without evidence of dopamine deficits [SWEDD]), when compared against a group of 654 HCs within the PPMI cohort and when not corrected by the Bonferroni test for multiple comparisons. Fourteen of 20 significant alleles were unique to the PD-HC comparison, whereas 31 of the 57 alleles overlapped between two or more different subgroup comparisons. Only the expressed HLA-DRA*01:01:01 and -DQA1*03:01:01 protective alleles (PD v HC), the -DQA1*03:03:01 risk (HC v Prodrome) or protective allele (PD v Prodrome), the -DRA*01:01:02 and -DRB4*01:03:02 risk alleles (SWEDD v HC), and the NR_SVA_381 present genotype (PD v HC) at a 5% homozygous insertion frequency near HLA-DPA1, were significant (Pc<0.1) after Bonferroni corrections. The homologous NR_SVA_381 insertion significantly decreased the transcription levels of HLA-DPA1 and HLA-DPB1 in the PPMI cohort and its presence as a homozygous genotype is a risk factor (Pc=0.012) for PD. The most frequent NR_SVA_381 insertion haplotype in the PPMI cohort was NR_SVA_381/DPA1*02/DPB1*01 (3.7%). Although HLA C*07/B*07/DRB5*01/DRB1*15/DQB1*06 was the most frequent HLA 5-loci phased-haplotype (n, 76) in the PPMI cohort, the NR_SVA_381 insertion was present in only six of them (8%).

Conclusions

These data suggest that expressed SVA and HLA gene alleles in circulating white blood cells are coordinated differentially in the regulation of immune responses and the long-term onset and progression of PD, the mechanisms of which have yet to be elucidated.

18th International HLA and Immunogenetics Workshop: Report on the SNP-HLA Reference Consortium (SHLARC) component

The SNP-HLA Reference Consortium (SHLARC), a component of the 18th International HLA and Immunogenetics Workshop, is aimed at collecting diverse and extensive human leukocyte antigen (HLA) data to create custom reference panels and enhance HLA imputation techniques. Genome-wide association studies (GWAS) have significantly contributed to identifying genetic associations with various diseases. The HLA genomic region has emerged as the top locus in GWAS, particularly in immune-related disorders. However, the limited information provided by single nucleotide polymorphisms (SNPs), the hallmark of GWAS, poses challenges, especially in the HLA region, where strong linkage disequilibrium (LD) spans several megabases. HLA imputation techniques have been developed using statistical inference in response to these challenges. These techniques enable the prediction of HLA alleles from genotyped GWAS SNPs. Here we present the SHLARC activities, a collaborative effort to create extensive, and multi-ethnic reference panels to enhance HLA imputation accuracy.

Mass Spectrometry-Based Immunopeptidomics of Peptides Presented on Human Leukocyte Antigen Proteins

Human leukocyte antigen (HLA) proteins are a group of glycoproteins that are expressed at the cell surface, where they present peptides to T cells through physical interactions with T-cell receptors (TCRs). Hence, characterizing the set of peptides presented by HLA proteins, referred to hereafter as the immunopeptidome, is fundamental for neoantigen identification, immunotherapy, and vaccine development. As a result, different methods have been used over the years to identify peptides presented by HLA proteins, including competition assays, peptide microarrays, and yeast display systems. Nonetheless, over the last decade, mass spectrometry-based immunopeptidomics (MS-immunopeptidomics) has emerged as the gold-standard method for identifying peptides presented by HLA proteins. MS-immunopeptidomics enables the direct identification of the immunopeptidome in different tissues and cell types in different physiological and pathological states, for example, solid tumors or virally infected cells. Despite its advantages, it is still an experimentally and computationally challenging technique with different aspects that need to be considered before planning an MS-immunopeptidomics experiment, while conducting the experiment and with analyzing and interpreting the results. Hence, we aim in this chapter to provide an overview of this method and discuss different practical considerations at different stages starting from sample collection until data analysis. These points should aid different groups aiming at utilizing MS-immunopeptidomics, as well as, identifying future research directions to improve the method.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

Sex-Dimorphic Analyses Identify Novel and Sex-Specific Genetic Associations in Inflammatory Bowel Disease

BACKGROUND

Sex is an integral variable often overlooked in complex disease genetics. Differences between sexes have been reported in natural history, disease complications, and age of onset in inflammatory bowel disease (IBD). While association studies have identified >230 IBD loci, there have been a limited number of studies investigating sex differences underlying these genetic associations.

METHODS

We report the first investigation of sex-dimorphic associations via meta-analysis of a sex-stratified association study (34 579 IBD cases, 39 125 controls). In addition, we performed chromosome (chr) X-specific analyses, considering models of X inactivation (XCI) and XCI escape. Demographic and clinical characteristics were also compared between sexes.

RESULTS

We identified significant differences between sexes for disease location and perianal complication in Crohn's disease and disease extent in ulcerative colitis. We observed genome-wide-significant sex-dimorphic associations (P < 5 × 10-8) at loci not previously reported in large-scale IBD genetic studies, including at chr9q22, CARMIL1, and UBASH3A. We identified variants in known IBD loci, including in chr2p15 and within the major histocompatibility complex on chr6, exhibiting sex-specific patterns of association (P < 5 × 10-7 in one sex only). We identified 3 chrX associations with IBD, including a novel Crohn's disease susceptibility locus at Xp22.

CONCLUSIONS

These analyses identified novel IBD loci, in addition to characterizing sex-specific patterns of associations underlying sex-dimorphic associations. By elucidating the role of sex in IBD genetics, our study will help enhance our understanding of the differences between the sexes in IBD biology and underscores a need to move beyond conventional sex-combined analyses to appreciate the genetic architecture of IBD more comprehensively.

HLA-DP on Epithelial Cells Enables Tissue Damage by NKp44+ Natural Killer Cells in Ulcerative Colitis

BACKGROUND & AIMS

Ulcerative colitis (UC) is characterized by severe inflammation and destruction of the intestinal epithelium, and is associated with specific risk single nucleotide polymorphisms in HLA class II. Given the recently discovered interactions between subsets of HLA-DP molecules and the activating natural killer (NK) cell receptor NKp44, genetic associations of UC and HLA-DP haplotypes and their functional implications were investigated.

METHODS

HLA-DP haplotype and UC risk association analyses were performed (UC: n = 13,927; control: n = 26,764). Expression levels of HLA-DP on intestinal epithelial cells (IECs) in individuals with and without UC were quantified. Human intestinal 3-dimensional (3D) organoid cocultures with human NK cells were used to determine functional consequences of interactions between HLA-DP and NKp44.

RESULTS

These studies identified HLA-DPA1∗01:03-DPB1∗04:01 (HLA-DP401) as a risk haplotype and HLA-DPA1∗01:03-DPB1∗03:01 (HLA-DP301) as a protective haplotype for UC in European populations. HLA-DP expression was significantly higher on IECs of individuals with UC compared with controls. IECs in human intestinal 3D organoids derived from HLA-DP401pos individuals showed significantly stronger binding of NKp44 compared with HLA-DP301pos IECs. HLA-DP401pos IECs in organoids triggered increased degranulation and tumor necrosis factor production by NKp44+ NK cells in cocultures, resulting in enhanced epithelial cell death compared with HLA-DP301pos organoids. Blocking of HLA-DP401-NKp44 interactions (anti-NKp44) abrogated NK cell activity in cocultures.

CONCLUSIONS

We identified an UC risk HLA-DP haplotype that engages NKp44 and activates NKp44+ NK cells, mediating damage to intestinal epithelial cells in an HLA-DP haplotype-dependent manner. The molecular interaction between NKp44 and HLA-DP401 in UC can be targeted by therapeutic interventions to reduce NKp44+ NK cell-mediated destruction of the intestinal epithelium in UC.

Pathological mechanism and targeted drugs of ulcerative colitis: A review

Ulcerative colitis (UC) is a chronic inflammatory disease of the colon with abdominal pain, diarrhea, and mucopurulent stools as the main symptoms. Its incidence is increasing worldwide, and traditional treatments have problems such as immunosuppression and metabolic disorders. In this article, the etiology and pathogenesis of ulcerative colitis are reviewed to clarify the targeted drugs of UC in the latest research. Our aim is to provide more ideas for the clinical treatment and new drug development of UC, mainly by analyzing and sorting out the relevant literature on PubMed, summarizing and finding that it is related to the main genetic, environmental, immune and other factors, and explaining its pathogenesis from the NF-κB pathway, PI3K/Akt signaling pathway, and JAK/STAT signaling pathway, and obtaining anti-TNF-α monoclonal antibodies, integrin antagonists, IL-12/IL-23 antagonists, novel UC-targeted drugs such as JAK inhibitors and SIP receptor agonists. We believe that rational selection of targeted drugs and formulation of the best dosing strategy under the comprehensive consideration of clinical evaluation is the best way to treat UC.

An Extremes of Phenotype Approach Confirms Significant Genetic Heterogeneity in Patients with Ulcerative Colitis

BACKGROUND AND AIMS

Ulcerative colitis [UC] is a major form of inflammatory bowel disease globally. Phenotypic heterogeneity is defined by several variables including age of onset and disease extent. The genetics of disease severity remains poorly understood. To further investigate this, we performed a genome wide association [GWA] study using an extremes of phenotype strategy.

METHODS

We conducted GWA analyses in 311 patients with medically refractory UC [MRUC], 287 with non-medically refractory UC [non-MRUC] and 583 controls. Odds ratios [ORs] were calculated for known risk variants comparing MRUC and non-MRUC, and controls.

RESULTS

MRUC-control analysis had the greatest yield of genome-wide significant single nucleotide polymorphisms [SNPs] [2018], including lead SNP = rs111838972 [OR = 1.82, p = 6.28 × 10-9] near MMEL1 and a locus in the human leukocyte antigen [HLA] region [lead SNP = rs144717024, OR = 12.23, p = 1.7 × 10-19]. ORs for the lead SNPs were significantly higher in MRUC compared to non-MRUC [p < 9.0 × 10-6]. No SNPs reached significance in the non-MRUC-control analysis (top SNP, rs7680780 [OR 2.70, p = 5.56 × 10-8). We replicate findings for rs4151651 in the Complement Factor B [CFB] gene and demonstrate significant changes in CFB gene expression in active UC. Detailed HLA analyses support the strong associations with MHC II genes, particularly HLA-DQA1, HLA-DQB1 and HLA-DRB1 in MRUC.

CONCLUSIONS

Our MRUC subgroup replicates multiple known UC risk variants in contrast to non-MRUC and demonstrates significant differences in effect sizes compared to those published. Non-MRUC cases demonstrate lower ORs similar to those published. Additional risk and prognostic loci may be identified by targeted recruitment of individuals with severe disease.

A novel unconventional T cell population enriched in Crohn's disease

OBJECTIVE

One of the current hypotheses to explain the proinflammatory immune response in IBD is a dysregulated T cell reaction to yet unknown intestinal antigens. As such, it may be possible to identify disease-associated T cell clonotypes by analysing the peripheral and intestinal T-cell receptor (TCR) repertoire of patients with IBD and controls.

DESIGN

We performed bulk TCR repertoire profiling of both the TCR alpha and beta chains using high-throughput sequencing in peripheral blood samples of a total of 244 patients with IBD and healthy controls as well as from matched blood and intestinal tissue of 59 patients with IBD and disease controls. We further characterised specific T cell clonotypes via single-cell RNAseq.

RESULTS

We identified a group of clonotypes, characterised by semi-invariant TCR alpha chains, to be significantly enriched in the blood of patients with Crohn's disease (CD) and particularly expanded in the CD8⁺ T cell population. Single-cell RNAseq data showed an innate-like phenotype of these cells, with a comparable gene expression to unconventional T cells such as mucosal associated invariant T and natural killer T (NKT) cells, but with distinct TCRs.

CONCLUSIONS

We identified and characterised a subpopulation of unconventional Crohn-associated invariant T (CAIT) cells. Multiple evidence suggests these cells to be part of the NKT type II population. The potential implications of this population for CD or a subset thereof remain to be elucidated, and the immunophenotype and antigen reactivity of CAIT cells need further investigations in future studies.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

HLA-DRB1∗16 and -DQB1∗05 alleles are strongly associated with autoimmune pancreatitis in a cohort of hundred patients

BACKGROUND/OBJECTIVES

Autoimmune diseases are often associated with human leukocyte antigen (HLA) haplotypes, indicating that changes in major histocompatibility complex (MHC)-dependent self-peptide or antigen presentation contribute to autoimmunity. In our study, we aimed to investigate HLA alleles in a large European cohort of autoimmune pancreatitis (AIP) patients.

METHODS

Hundred patients with AIP, diagnosed and classified according to the International Consensus Diagnostic Criteria (ICDC), were prospectively enrolled in the study. Forty-four patients with chronic pancreatitis (CP) and 254 healthy subjects served as control groups. DNA was isolated from blood samples and two-digit HLA typing was performed with sequence-specific primer (SSP-) PCR. HLA allele association strength to AIP was calculated as odds ratio.

RESULTS

We uncovered a strong enrichment of HLA-DQB1 homozygosity in type 1 and type 2 AIP patients. Moreover, a significantly increased incidence of the HLA-DRB1∗16 and HLA-DQB1∗05 alleles and a concomitant lack of the HLA-DRB1∗13 allele was detected in AIP type 1 and type 2 patients. In contrast, the HLA-DQB1∗02 allele was underrepresented in the 'not otherwise specified' (NOS) AIP subtype. We detected no significant difference in the HLA-DRB3, HLA-DRB4 and HLA-DRB5 allele frequency in our cohort.

CONCLUSIONS

Although AIP type 1 and type 2 are characterized by distinct histopathological characteristics, both subtypes are associated with the same HLA alleles, indicating that the disease might rely on similar immunogenic mechanisms. However, AIP NOS represented another subclass of AIP.

A genome-wide association study in autoimmune neurological syndromes with anti-GAD65 autoantibodies

Autoimmune neurological syndromes (AINS) with autoantibodies against the 65 kDa isoform of the glutamic acid decarboxylase (GAD65) present with limbic encephalitis, including temporal lobe seizures or epilepsy, cerebellitis with ataxia, and stiff-person-syndrome or overlap forms. Anti-GAD65 autoantibodies are also detected in autoimmune diabetes mellitus, which has a strong genetic susceptibility conferred by human leukocyte antigen (HLA) and non-HLA genomic regions. We investigated the genetic predisposition in patients with anti-GAD65 AINS. We performed a genome-wide association study (GWAS) and an association analysis of the HLA region in a large German cohort of 1214 individuals. These included 167 patients with anti-GAD65 AINS, recruited by the German Network for Research on Autoimmune Encephalitis (GENERATE), and 1047 individuals without neurological or endocrine disease as population-based controls. Predictions of protein expression changes based on GWAS findings were further explored and validated in the CSF proteome of a virtually independent cohort of 10 patients with GAD65-AINS and 10 controls. Our GWAS identified 16 genome-wide significant (P < 5 × 10-8) loci for the susceptibility to anti-GAD65 AINS. The top variant, rs2535288 [P = 4.42 × 10-16, odds ratio (OR) = 0.26, 95% confidence interval (CI) = 0.187-0.358], localized to an intergenic segment in the middle of the HLA class I region. The great majority of variants in these loci (>90%) mapped to non-coding regions of the genome. Over 40% of the variants have known regulatory functions on the expression of 48 genes in disease relevant cells and tissues, mainly CD4+ T cells and the cerebral cortex. The annotation of epigenomic marks suggested specificity for neural and immune cells. A network analysis of the implicated protein-coding genes highlighted the role of protein kinase C beta (PRKCB) and identified an enrichment of numerous biological pathways participating in immunity and neural function. Analysis of the classical HLA alleles and haplotypes showed no genome-wide significant associations. The strongest associations were found for the DQA1*03:01-DQB1*03:02-DRB1*04:01HLA haplotype (P = 4.39 × 10-4, OR = 2.5, 95%CI = 1.499-4.157) and DRB1*04:01 allele (P = 8.3 × 10-5, OR = 2.4, 95%CI = 1.548-3.682) identified in our cohort. As predicted, the CSF proteome showed differential levels of five proteins (HLA-A/B, C4A, ATG4D and NEO1) of expression quantitative trait loci genes from our GWAS in the CSF proteome of anti-GAD65 AINS. These findings suggest a strong genetic predisposition with direct functional implications for immunity and neural function in anti-GAD65 AINS, mainly conferred by genomic regions outside the classical HLA alleles.

HLA imputation and its application to genetic and molecular fine-mapping of the MHC region in autoimmune diseases

Variations of human leukocyte antigen (HLA) genes in the major histocompatibility complex region (MHC) significantly affect the risk of various diseases, especially autoimmune diseases. Fine-mapping of causal variants in this region was challenging due to the difficulty in sequencing and its inapplicability to large cohorts. Thus, HLA imputation, a method to infer HLA types from regional single nucleotide polymorphisms, has been developed and has successfully contributed to MHC fine-mapping of various diseases. Different HLA imputation methods have been developed, each with its own advantages, and recent methods have been improved in terms of accuracy and computational performance. Additionally, advances in HLA reference panels by next-generation sequencing technologies have enabled higher resolution and a more reliable imputation, allowing a finer-grained evaluation of the association between sequence variations and disease risk. Risk-associated variants in the MHC region would affect disease susceptibility through complicated mechanisms including alterations in peripheral responses and central thymic selection of T cells. The cooperation of reliable HLA imputation methods, informative fine-mapping, and experimental validation of the functional significance of MHC variations would be essential for further understanding of the role of the MHC in the immunopathology of autoimmune diseases.

Diversity of HLA class I and class II alleles in Iran populations: Systematic review and Meta-Analaysis

The human leukocyte antigen (HLA) system plays an essential role in the peptides antigen presentation and more regulation of immune responses. Regarding all HLA molecules' associations with various diseases and their clinical utilities in understanding drug reactions or prediction of transplantation outcome, there is a need for much more extensive HLA data generated from Asian countries.

METHOD

A comprehensive search was conducted in electronic databases between 1990 and 2021 to identify relevant articles to HLA frequency in the normal Iranian population. Two independent reviewers screened and selected the eligible studies. After data extraction, the meta-analysis was performed using STATA version 14. The overall frequencies and their 95% confidence intervals (CIs) were obtained using the random-effects model.

RESULTS

Among 1141 studies 78 were eligible for this study and the sample sizes varied from 14 to 15,600. The most frequent alleles of HLA class I were HLA-A*02 (22%; 95%CI: 20-24%; I² = 88.63%), -B*35 (18%; 95%CI: 16-21%; I² = 90.95%), -C*12 (18%; 95%CI: 13-22%; I² = 89.51%). HLA-DQA1*01 (42%; 95%CI: 40-44%; I² = 56.80%), -DQB1*03 (38%; 95%CI: 35-42%; I² = 92.38%), and -DRB1*11 (24%; 95%CI: 22-26%; I² = 90.72%) were the most frequent alleles of HLA class II in Iran.

DISCUSSION

Our meta-analysis results point out that the comprehensive report of HLA allele frequency in the Iranian population could be helpful as reference data for planning and managing transplantation and immune disease treatment in Iran.

Immunopeptidomics toolkit library (IPTK): a python-based modular toolbox for analyzing immunopeptidomics data

Background

The human leukocyte antigen (HLA) proteins play a fundamental role in the adaptive immune system as they present peptides to T cells. Mass-spectrometry-based immunopeptidomics is a promising and powerful tool for characterizing the immunopeptidomic landscape of HLA proteins, that is the peptides presented on HLA proteins. Despite the growing interest in the technology, and the recent rise of immunopeptidomics-specific identification pipelines, there is still a gap in data-analysis and software tools that are specialized in analyzing and visualizing immunopeptidomics data.

Results

We present the IPTK library which is an open-source Python-based library for analyzing, visualizing, comparing, and integrating different omics layers with the identified peptides for an in-depth characterization of the immunopeptidome. Using different datasets, we illustrate the ability of the library to enrich the result of the identified peptidomes. Also, we demonstrate the utility of the library in developing other software and tools by developing an easy-to-use dashboard that can be used for the interactive analysis of the results.

Conclusion

IPTK provides a modular and extendable framework for analyzing and integrating immunopeptidomes with different omics layers. The library is deployed into PyPI at https://pypi.org/project/IPTKL/ and into Bioconda at https://anaconda.org/bioconda/iptkl, while the source code of the library and the dashboard, along with the online tutorials are available at https://github.com/ikmb/iptoolkit.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04315-0.

HLA class II immunopeptidomics reveals that co-inherited HLA-allotypes within an extended haplotype can improve proteome coverage for immunosurveillance

Human leucocyte antigen (HLA) class II molecules in humans are encoded by three different loci, HLA-DR, -DQ, and -DP. These molecules share approximately 70% sequence similarity and all present peptide ligands to circulating T cells. While the peptide repertoires of numerous HLA-DR, -DQ, and -DP allotypes have been examined, there have been few reports on the combined repertoire of these co-inherited molecules expressed in a single cell as an extended HLA haplotype. Here we describe the endogenous peptide repertoire of a human B lymphoblastoid cell line (C1R) expressing the class II haplotype HLA-DR12/DQ7/DP4. We have identified 71350 unique naturally processed peptides presented collectively by HLA-DR12, HLA-DQ7, or HLA-DP4. The resulting "haplodome" is complemented by the cellular proteome defined by standard LC-MS/MS approaches. This large dataset has shed light on properties of these class II ligands especially the preference for membrane and extracellular source proteins. Our data also provides insights into the co-evolution of these conserved haplotypes of closely linked and co-inherited HLA molecules; which together increase sequence coverage of cellular proteins for immune surveillance with minimal overlap between each co-inherited HLA-class II allomorph.