Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

Exploring the HLA complex in autoimmunity: From the risk haplotypes to the modulation of expression

The genes mapping at the HLA region show high density, strong linkage disequilibrium and high polymorphism, which affect the association of HLA class I and class II genes with autoimmunity. We focused on the HLA haplotypes, genomic structures consisting of an array of specific alleles showing some degrees of genetic association with different autoimmune disorders. GWASs in many pathologies have identified variants in either the coding loci or the flanking regulatory regions, both in linkage disequilibrium in haplotypes, that are frequently associated with increased risk and may influence gene expression. We discuss the relevance of the HLA gene expression because the level of surface heterodimers determines the number of complexes presenting self-antigen and, thus, the strength of pathogenic autoreactive T cells immune response.

HLA-DQB1*05 subtypes and not DRB1*10:01 mediates risk in anti-IgLON5 disease

Anti-IgLON5 disease is a rare and likely underdiagnosed subtype of autoimmune encephalitis. The disease displays a heterogeneous phenotype that includes sleep, movement and bulbar-associated dysfunction. The presence of IgLON5-antibodies in CSF/serum, together with a strong association with HLA-DRB1*10:01∼DQB1*05:01, supports an autoimmune basis. In this study, a multicentric human leukocyte antigen (HLA) study of 87 anti-IgLON5 patients revealed a stronger association with HLA-DQ than HLA-DR. Specifically, we identified a predisposing rank-wise association with HLA-DQA1*01:05∼DQB1*05:01, HLA-DQA1*01:01∼DQB1*05:01 and HLA-DQA1*01:04∼DQB1*05:03 in 85% of patients. HLA sequences and binding cores for these three DQ heterodimers were similar, unlike those of linked DRB1 alleles, supporting a causal link to HLA-DQ. This association was further reflected in an increasingly later age of onset across each genotype group, with a delay of up to 11 years, while HLA-DQ-dosage dependent effects were also suggested by reduced risk in the presence of non-predisposing DQ1 alleles. The functional relevance of the observed HLA-DQ molecules was studied with competition binding assays. These proof-of-concept experiments revealed preferential binding of IgLON5 in a post-translationally modified, but not native, state to all three risk-associated HLA-DQ receptors. Further, a deamidated peptide from the Ig2-domain of IgLON5 activated T cells in two patients, compared with one control carrying HLA-DQA1*01:05∼DQB1*05:01. Taken together, these data support a HLA-DQ-mediated T-cell response to IgLON5 as a potentially key step in the initiation of autoimmunity in this disease.

The effect of type 1 diabetes protection and susceptibility associated HLA class II genotypes on DNA methylation in immune cells

The HLA region, especially HLA class I and II genes, which encode molecules for antigen presentation to T cells, plays a major role in the predisposition to autoimmune disorders. To clarify the mechanisms behind this association, we examined genome-wide DNA methylation by microarrays to cover over 850,000 CpG sites in the CD4+ T cells and CD19+ B cells of healthy subjects homozygous either for DRB1*15-DQA1*01-DQB1*06:02 (DR2-DQ6, n = 14), associated with a strongly decreased T1D risk, DRB1*03-DQA1*05-DQB1*02 (DR3-DQ2, n = 19), or DRB1*04:01-DQA1*03-DQB1*03:02 (DR4-DQ8, n = 17), associated with a moderately increased T1D risk. In total, we discovered 14 differentially methylated CpG probes, of which 10 were located in the HLA region and six in the HLA-DRB1 locus. The main differences were between the protective genotype DR2-DQ6 and the risk genotypes DR3-DQ2 and DR4-DQ8, where the DR2-DQ6 group was hypomethylated compared to the other groups in all but four of the differentially methylated probes. The differences between the risk genotypes DR3-DQ2 and DR4-DQ8 were small. Our results indicate that HLA variants have few systemic effects on methylation and that their effect on autoimmunity is conveyed directly by HLA molecules, possibly by differences in expression levels or function.

Nucleotide alterations in the HLA-C class I gene can cause aberrant splicing and marked changes in RNA levels in a polymorphic context-dependent manner

Polymorphisms of HLA genes, which play a crucial role in presenting peptides with diverse sequences in their peptide-binding pockets, are also thought to affect HLA gene expression, as many studies have reported associations between HLA gene polymorphisms and their expression levels. In this study, we devised an ectopic expression assay for the HLA class I genes in the context of the entire gene, and used the assay to show that the HLA-C*03:03:01 and C*04:01:01 polymorphic differences observed in association studies indeed cause different levels of RNA expression. Subsequently, we investigated the C*03:23N null allele, which was previously noted for its reduced expression, attributed to an alternate exon 3 3' splice site generated by G/A polymorphism at position 781 within the exon 3. We conducted a thorough analysis of the splicing patterns of C*03:23N, and revealed multiple aberrant splicing, including the exon 3 alternative splicing, which overshadowed its canonical counterpart. After confirming a significant reduction in RNA levels caused by the G781A alteration in our ectopic assay, we probed the function of the G-rich sequence preceding the canonical exon 3 3' splice site. Substituting the G-rich sequence with a typical pyrimidine-rich 3' splice site sequence on C*03:23N resulted in a marked elevation in RNA levels, likely due to the enhanced preference for the canonical exon 3 3' splice site over the alternate site. However, the same substitution led to a reduction in RNA levels for C*03:03:01. These findings suggested the dual roles of the G-rich sequence in RNA expression, and furthermore, underscore the importance of studying polymorphism effects within the framework of the entire gene, extending beyond conventional mini-gene reporter assays.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

Crossmatch assays in transplantation: Physical or virtual?: A review

The value of the crossmatch test in assessing pretransplant immunological risk is vital for clinical decisions, ranging from the indication of the transplant to the guidance of induction protocols and treatment with immunosuppressants. The crossmatch tests in transplantation can be physical or virtual, each with its advantages and limitations. Currently, the virtual crossmatch stands out for its sensitivity and specificity compared to the physical tests. Additionally, the virtual crossmatch can be performed in less time, allowing for a reduction in cold ischemia time. It shows a good correlation with the results of physical tests and does not negatively impact graft survival. Proper communication between clinicians and the transplant immunology laboratory will lead to a deeper understanding of each patient's immunological profile, better donor-recipient selection, and improved graft survival.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Killer instincts: natural killer cells as multifactorial cancer immunotherapy

Natural killer (NK) cells integrate heterogeneous signals for activation and inhibition using germline-encoded receptors. These receptors are stochastically co-expressed, and their concurrent engagement and signaling can adjust the sensitivity of individual cells to putative targets. Against cancers, which mutate and evolve under therapeutic and immunologic pressure, the diversity for recognition provided by NK cells may be key to comprehensive cancer control. NK cells are already being trialled as adoptive cell therapy and targets for immunotherapeutic agents. However, strategies to leverage their naturally occurring diversity and agility have not yet been developed. In this review, we discuss the receptors and signaling pathways through which signals for activation or inhibition are generated in NK cells, focusing on their roles in cancer and potential as targets for immunotherapies. Finally, we consider the impacts of receptor co-expression and the potential to engage multiple pathways of NK cell reactivity to maximize the scope and strength of antitumor activities.

Development of a semi-automated MHC-associated peptide proteomics (MAPPs) method using streptavidin bead-based immunoaffinity capture and nano LC-MS/MS to support immunogenicity risk assessment in drug development

Major histocompatibility complex (MHC)-Associated Peptide Proteomics (MAPPs) is an ex vivo method used to assess the immunogenicity risk of biotherapeutics. MAPPs can identify potential T-cell epitopes within the biotherapeutic molecule. Using adalimumab treated human monocyte derived dendritic cells (DCs) and a pan anti-HLA-DR antibody (Ab), we systematically automated and optimized biotin/streptavidin (SA)-capture antibody coupling, lysate incubation with capture antibody, as well as the washing and elution steps of a MAPPs method using functionalized magnetic beads and a KingFisher Magnetic Particle processor. Automation of these steps, combined with capturing using biotinylated-Ab/SA magnetic beads rather than covalently bound antibody, improved reproducibility as measured by minimal inter-and intra-day variability, as well as minimal analyst-to-analyst variability. The semi-automated MAPPs workflow improved sensitivity, allowing for a lower number of cells per analysis. The method was assessed using five different biotherapeutics with varying immunogenicity rates ranging from 0.1 to 48% ADA incidence in the clinic. Biotherapeutics with ≥10%immunogenicity incidence consistently presented more peptides (1.8-28 fold) and clusters (10-21 fold) compared to those with <10% immunogenicity incidence. Our semi-automated MAPPs method provided two main advantages over a manual workflow- the robustness and reproducibility affords confidence in the epitopes identified from as few as 5 to 10 donors and the method workflow can be readily adapted to incorporate different capture Abs in addition to anti-HLA-DR. The incorporation of semi-automated MAPPs with biotinylated-Ab/SA bead-based capture in immunogenicity screening strategies allows the generation of more consistent and reliable data, helping to improve immunogenicity prediction capabilities in drug development. MHC associated peptide proteomics (MAPPs), Immunogenicity risk assessment, in vitro/ex vivo, biotherapeutics, Major Histocompatibility Complex Class II (MHC II), LC-MS, Immunoaffinity Capture, streptavidin magnetic beads.

Pilot study of genome-wide DNA methylation and gene expression for treatment response to escitalopram in panic disorder

BACKGROUND

Antidepressants, particularly selective serotonin reuptake inhibitors, are currently considered the first-line treatment for panic disorder (PD). However, little is known about the relationship between the biomarkers that may predict better treatment.

AIM

To compare genome-wide methylation and gene expression patterns between responsive and non-responsive patients with PD after 4 wk of escitalopram treatment.

METHODS

Thirty patients with PD were enrolled in this study (responders = 13; non-responders = 17). All patients were assessed using the PD Severity Scale-Chinese version before and after treatment. The Illumina Infinium MethylationEPIC (850k) BeadChip for genome-wide methylation screening and mRNA sequencing was used in all patients with PD.

RESULTS

A total of 701 differentially methylated positions (DMPs) were found between responders and non-responders (|Δβ| ≥ 0.06, q < 0.05), and the hyper- and hypomethylated CpG sites were 511 (72.9%) and 190 (27.1%), respectively. Relative to non-responders, there were 59 differential transcripts, of which 20 were downregulated and 39 were upregulated (q < 0.05). However, no differentially expressed genes were identified by mRNA sequencing after correcting for multiple testing (|log2(FC)| > 1, q > 0.05).

CONCLUSION

This preliminary study showed that DMPs might be associated with the treatment response to escitalopram in PD; however, these DMPs need to be verified in large samples.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Maternal carriage of H-Y restricting HLA class II alleles is a negative prognostic factor for women with recurrent pregnancy loss after birth of a boy

Immune system aberrations are suggested to be an important factor in the pathophysiology of unexplained secondary recurrent pregnancy loss (sRPL). The objective was to investigate if the sex ratio of the firstborn child in sRPL patients differs from the background population and whether the sex of the firstborn child has a negative impact on the pregnancy prognosis alone and/or in combination with carriage of male-specific minor histocompatibility (H-Y) restricting HLA class II alleles. From January 2016 to October 2022, 582 patients with unexplained RPL were admitted to the RPL Center of Western Denmark and continuously followed-up. HLA-DRB1 and -DQB1 typing was performed as part of the routine diagnostic work-up. In sRPL patients, a history of a firstborn boy was significantly more frequent than in the Danish background population and was associated with significantly lower odds of a successful reproductive outcome in the first pregnancy after admission compared to a firstborn girl (OR=0.41, 95% CI: 0.20-0.83, p = 0.014). The odds of a successful reproductive outcome were enhanced in patients carrying ≥ 1 H-Y-restricting HLA class II alleles with a first-born girl compared to a firstborn boy (OR=3.33, 95% CI: 1.40-7.88, p = 0.005), while no difference in successful reproductive outcome was seen in sRPL patients not carrying these alleles (OR=1.20, 95% CI: 0.33-4.43, p = 0.781). The sex ratio of children born after RPL was similar to the Danish background population. These findings confirm previous findings and suggests that a harmful immune response triggered by H-Y-antigen exposure during a previous pregnancy in preconditioned women may cause sRPL.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes

The human Major Histocompatibility Complex (MHC) or Human Leukocyte Antigen (HLA) super-locus is a highly polymorphic genomic region that encodes more than 140 coding genes including the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. In recent years, MHC genotyping and haplotyping using Sanger sequencing and next-generation sequencing (NGS) methods have produced many hundreds of genomic sequences of the HLA super-locus for comparative studies of the genetic architecture and diversity between the same and different haplotypes. In this special issue on 'The Current Landscape of HLA Genomics and Genetics', we provide a short review of some of the recent analytical developments used to investigate the SNP polymorphisms, structural variants (indels), transcription and haplotypes of the HLA super-locus. This review highlights the importance of using reference cell-lines, population studies, and NGS methods to improve and update our understanding of the mechanisms, architectural structures and combinations of human MHC genomic alleles (SNPs and indels) that better define and characterise haplotypes and their association with various phenotypes and diseases.

MUC22, HLA-A, and HLA-DOB variants and COVID-19 in resilient super-agers from Brazil

BACKGROUND

Although aging correlates with a worse prognosis for Covid-19, super elderly still unvaccinated individuals presenting mild or no symptoms have been reported worldwide. Most of the reported genetic variants responsible for increased disease susceptibility are associated with immune response, involving type I IFN immunity and modulation; HLA cluster genes; inflammasome activation; genes of interleukins; and chemokines receptors. On the other hand, little is known about the resistance mechanisms against SARS-CoV-2 infection. Here, we addressed polymorphisms in the MHC region associated with Covid-19 outcome in super elderly resilient patients as compared to younger patients with a severe outcome.

METHODS

SARS-CoV-2 infection was confirmed by RT-PCR test. Aiming to identify candidate genes associated with host resistance, we investigated 87 individuals older than 90 years who recovered from Covid-19 with mild symptoms or who remained asymptomatic following positive test for SARS-CoV-2 as compared to 55 individuals younger than 60 years who had a severe disease or died due to Covid-19, as well as to the general elderly population from the same city. Whole-exome sequencing and an in-depth analysis of the MHC region was performed. All samples were collected in early 2020 and before the local vaccination programs started.

RESULTS

We found that the resilient super elderly group displayed a higher frequency of some missense variants in the MUC22 gene (a member of the mucins' family) as one of the strongest signals in the MHC region as compared to the severe Covid-19 group and the general elderly control population. For example, the missense variant rs62399430 at MUC22 is two times more frequent among the resilient super elderly (p = 0.00002, OR = 2.24).

CONCLUSION

Since the pro-inflammatory basal state in the elderly may enhance the susceptibility to severe Covid-19, we hypothesized that MUC22 might play an important protective role against severe Covid-19, by reducing overactive immune responses in the senior population.

HLA allele-specific expression: Methods, disease associations, and relevance in hematopoietic stem cell transplantation

Varying HLA allele-specific expression levels are associated with human diseases, such as graft versus host disease (GvHD) in hematopoietic stem cell transplantation (HSCT), cytotoxic T cell response and viral load in HIV infection, and the risk of Crohn’s disease. Only recently, RNA-based next generation sequencing (NGS) methodologies with accompanying bioinformatics tools have emerged to quantify HLA allele-specific expression replacing the quantitative PCR (qPCR) -based methods. These novel NGS approaches enable the systematic analysis of the HLA allele-specific expression changes between individuals and between normal and disease phenotypes. Additionally, analyzing HLA allele-specific expression and allele-specific expression loss provide important information for predicting efficacies of novel immune cell therapies. Here, we review available RNA sequencing-based approaches and computational tools for NGS to quantify HLA allele-specific expression. Moreover, we explore recent studies reporting disease associations with differential HLA expression. Finally, we discuss the role of allele-specific expression in HSCT and how considering the expression quantification in recipient-donor matching could improve the outcome of HSCT.

Distribution of HLA-A, -B, -C, -DRB1, -DQB1, -DPB1 allele frequencies in patients with COVID-19 bilateral pneumonia in Russians, living in the Chelyabinsk region (Russia)

In this population-based case-control study conducted in the Chelyabinsk region of Russia, we examined the distribution of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1, in a group of 100 patients with confirmed COVID-19 bilateral pneumonia. Typing was performed by NGS and statistical calculations were carried out with the Arlequin program. HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 alleles were compared between patients with COVID-19 and 99 healthy controls. We identified that COVID-19 susceptibility is associated with alleles and genotypes rs9277534A (disequilibrium with HLA-DPB1*02:01, -02:02, -04:01, -04:02, -17:01 alleles) with low expression of protein products HLA-DPB1 (pc < 0.028) and homozygosity at HLA-C*04 (p = 0.024, pc = 0.312). Allele HLA-A*01:01 was decreased in a group of patients with severe forms of bilateral pneumonia, and therefore it may be considered as a protective factor for the development of severe symptoms of COVID-19 (p = 0.009, pc = 0.225). Our studies provide further evidence for the functional association between HLA genes and COVID-19.

How genetic risk contributes to autoimmune liver disease

Genome-wide association studies (GWAS) for autoimmune hepatitis (AIH) and GWAS/genome-wide meta-analyses (GWMA) for primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC) have been successful over the past decade, identifying about 100 susceptibility loci in the human genome, with strong associations with the HLA locus and many susceptibility variants outside the HLA locus with relatively low risk. However, identifying causative variants and genes and determining their effects on liver cells and their immunological microenvironment is far from trivial. Polygenic risk scores (PRSs) based on current genome-wide data have limited potential to predict individual disease risk. Interestingly, results of mediated expression score regression analysis provide evidence that a substantial portion of gene expression at susceptibility loci is mediated by genetic risk variants, in contrast to many other complex diseases. Genome- and transcriptome-wide comparisons between AIH, PBC, and PSC could help to better delineate the shared inherited component of autoimmune liver diseases (AILDs), and statistical fine-mapping, chromosome X-wide association testing, and genome-wide in silico drug screening approaches recently applied to GWMA data from PBC could potentially be successfully applied to AIH and PSC. Initial successes through single-cell RNA sequencing (scRNA-seq) experiments in PBC and PSC now raise high hopes for understanding the impact of genetic risk variants in the context of liver-resident immune cells and liver cell subpopulations, and for bridging the gap between genetics and disease.

Regulation of HLA class I expression by non-coding gene variations

The Human Leukocyte Antigen (HLA) is a critical genetic system for different outcomes after solid organ and hematopoietic cell transplantation. Its polymorphism is usually determined by molecular technologies at the DNA level. A potential role of HLA allelic expression remains under investigation in the context of the allogenic immune response between donors and recipients. In this study, we quantified the allelic expression of all three HLA class I loci (HLA-A, B and C) by RNA sequencing and conducted an analysis of expression quantitative traits loci (eQTL) to investigate whether HLA expression regulation could be associated with non-coding gene variations. HLA-B alleles exhibited the highest expression levels followed by HLA-C and HLA-A alleles. The max fold expression variation was observed for HLA-C alleles. The expression of HLA class I loci of distinct individuals demonstrated a coordinated and paired expression of both alleles of the same locus. Expression of conserved HLA-A~B~C haplotypes differed in distinct PBMC's suggesting an individual regulated expression of both HLA class I alleles and haplotypes. Cytokines TNFα /IFNβ, which induced a very similar upregulation of HLA class I RNA and cell surface expression across alleles did not modify the individually coordinated expression at the three HLA class I loci. By identifying cis eQTLs for the HLA class I genes, we show that the non-coding eQTLs explain 29%, 13%, and 31% of the respective HLA-A, B, C expression variance in unstimulated cells, and 9%, 23%, and 50% of the variance in cytokine-stimulated cells. The eQTLs have significantly higher effect sizes in stimulated cells compared to unstimulated cells for HLA-B and HLA-C genes expression. Our data also suggest that the identified eQTLs are independent from the coding variation which defines HLA alleles and thus may be influential on intra-allele expression variability although they might not represent the causal eQTLs.

HLA-DQ heterodimers in hematopoietic cell transplantation

HLA-DQ heterodimers increase the susceptibility to autoimmune diseases, but their role in hematopoietic cell transplantation is unknown. We tested the hypothesis that outcome after HLA-matched and HLA-DQ-mismatched hematopoietic cell transplantation is influenced by HLA-DQ heterodimers. Heterodimers were defined in 5164 HLA-matched and 520 HLA-DQ-mismatched patients and their transplant donors according to well-established crystallographic criteria. Group 1 (G1) heterodimers are any DQA1*02/03/04/05/06α paired with any DQB1*02/03/04β. Group 2 (G2) heterodimers are DQA1*01α paired with any DQB1*05/06β. Multivariable models identified significantly higher relapse risk in G1G2 and G2G2 compared with G1G1 HLA-matched patients with malignant disease; risk increased with an increasing number of G2 molecules. In HLA-DQ-mismatched transplantation for malignant diseases, matching or mismatching for G2 increased relapse risk. G2 lowered disease-free survival after both HLA-matched and HLA-DQ-mismatched transplantation. A paradigm based on HLA-DQ heterodimers provides a functional definition of the hematopoietic cell transplantation barrier and a means to lower risks for future patients.

HLA informs risk predictions after haploidentical stem cell transplantation with posttransplantation cyclophosphamide

Hematopoietic cell transplantation from HLA-haploidentical related donors is increasingly used to treat hematologic cancers; however, characteristics of the optimal haploidentical donor have not been established. We studied the role of donor HLA mismatching in graft-versus-host disease (GVHD), disease recurrence, and survival after haploidentical donor transplantation with posttransplantation cyclophosphamide (PTCy) for 1434 acute leukemia or myelodysplastic syndrome patients reported to the Center for International Blood and Marrow Transplant Research. The impact of mismatching in the graft-versus-host vector for HLA-A, -B, -C, -DRB1, and -DQB1 alleles, the HLA-B leader, and HLA-DPB1 T-cell epitope (TCE) were studied using multivariable regression methods. Outcome was associated with HLA (mis)matches at individual loci rather than the total number of HLA mismatches. HLA-DRB1 mismatches were associated with lower risk of disease recurrence. HLA-DRB1 mismatching with HLA-DQB1 matching correlated with improved disease-free survival. HLA-B leader matching and HLA-DPB1 TCE-nonpermissive mismatching were each associated with improved overall survival. HLA-C matching lowered chronic GVHD risk, and the level of HLA-C expression correlated with transplant-related mortality. Matching status at the HLA-B leader and HLA-DRB1, -DQB1, and -DPB1 predicted disease-free survival, as did patient and donor cytomegalovirus serostatus, patient age, and comorbidity index. A web-based tool was developed to facilitate selection of the best haploidentical-related donor by calculating disease-free survival based on these characteristics. In conclusion, HLA factors influence the success of haploidentical transplantation with PTCy. HLA-DRB1 and -DPB1 mismatching and HLA-C, -B leader, and -DQB1 matching are favorable. Consideration of HLA factors may help to optimize the selection of haploidentical related donors.

Downregulation of HLA class II is associated with relapse after allogeneic stem cell transplantation and alters recognition by antigen-specific T cells

Genomic deletion of donor-patient-mismatched HLA alleles in leukemic cells is a major cause of relapse after allogeneic hematopoietic stem cell transplantation (HSCT). Mismatched HLA is frequently lost as an individual allele or a whole region in HLA-class I, however, it is downregulated in HLA-class II. We hypothesized that there might be a difference in T cell recognition capacity against epitopes associated with HLA-class I and HLA-class II and consequently such allogeneic immune pressure induced HLA alterations in leukemic cells. To investigate this, we conducted in vitro experiments with T cell receptor-transduced T (TCR-T) cells. The cytotoxic activity of NY-ESO-1-specific TCR-T cells exhibited similarly against K562 cells with low HLA-A*02:01 expression. However, we demonstrated that the cytokine production against low HLA-DPB1*05:01 expression line decreased gradually from the HLA expression level approximately 2-log lower than normal expressors. Using sort-purified leukemia cells before and after HSCT, we applied the next-generation sequencing, and revealed that there were several marked downregulations of HLA-class II alleles which demonstrated consistently low expression from pre-transplantation. The marked downregulation of HLA-class II may lead to decreased antigen recognition ability of antigen-specific T cells and may be one of immune evasion mechanism associated with HLA-class II downregulation.

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.

MHC Variants Associated With Symptomatic Versus Asymptomatic SARS-CoV-2 Infection in Highly Exposed Individuals

Despite the high number of individuals infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) who develop coronavirus disease 2019 (COVID-19) symptoms worldwide, many exposed individuals remain asymptomatic and/or uninfected and seronegative. This could be explained by a combination of environmental (exposure), immunological (previous infection), epigenetic, and genetic factors. Aiming to identify genetic factors involved in immune response in symptomatic COVID-19 as compared to asymptomatic exposed individuals, we analyzed 83 Brazilian couples where one individual was infected and symptomatic while the partner remained asymptomatic and serum-negative for at least 6 months despite sharing the same bedroom during the infection. We refer to these as "discordant couples". We performed whole-exome sequencing followed by a state-of-the-art method to call genotypes and haplotypes across the highly polymorphic major histocompatibility complex (MHC) region. The discordant partners had comparable ages and genetic ancestry, but women were overrepresented (65%) in the asymptomatic group. In the antigen-presentation pathway, we observed an association between HLA-DRB1 alleles encoding Lys at residue 71 (mostly DRB1*03:01 and DRB1*04:01) and DOB*01:02 with symptomatic infections and HLA-A alleles encoding 144Q/151R with asymptomatic seronegative women. Among the genes related to immune modulation, we detected variants in MICA and MICB associated with symptomatic infections. These variants are related to higher expression of soluble MICA and low expression of MICB. Thus, quantitative differences in these molecules that modulate natural killer (NK) activity could contribute to susceptibility to COVID-19 by downregulating NK cell cytotoxic activity in infected individuals but not in the asymptomatic partners.

Identification of three novel HLA-DRA alleles by next-generation sequencing

Three novel HLA-DRA alleles, DRA*01:03, DRA*01:04, and DRA*01:05 alleles with unique amino acid sequences.

Haplotype-Specific Expression Analysis of MHC Class II Genes in Healthy Individuals and Rheumatoid Arthritis Patients

HLA-DRB1 alleles have been associated with several autoimmune diseases. For anti-citrullinated protein antibody positive rheumatoid arthritis (RA), HLA-DRB1 shared epitope (SE) alleles are the major genetic risk factors. In order to study the genetic regulation of major histocompatibility complex (MHC) Class II gene expression in immune cells, we investigated transcriptomic profiles of a variety of immune cells from healthy individuals carrying different HLA-DRB1 alleles. Sequencing libraries from peripheral blood mononuclear cells, CD4+ T cells, CD8+ T cells, and CD14+ monocytes of 32 genetically pre-selected healthy female individuals were generated, sequenced and reads were aligned to the standard reference. For the MHC region, reads were mapped to available MHC reference haplotypes and AltHapAlignR was used to estimate gene expression. Using this method, HLA-DRB and HLA-DQ were found to be differentially expressed in different immune cells of healthy individuals as well as in whole blood samples of RA patients carrying HLA-DRB1 SE-positive versus SE-negative alleles. In contrast, no genes outside the MHC region were differentially expressed between individuals carrying HLA-DRB1 SE-positive and SE-negative alleles, thus HLA-DRB1 SE alleles have a strong cis effect on gene expression. Altogether, our findings suggest that immune effects associated with different allelic forms of HLA-DR and HLA-DQ may be associated not only with differences in the structure of these proteins, but also with differences in their expression levels.

Plasma cell biology: Foundations for targeted therapeutic development in transplantation

Solid organ transplantation is a life-saving procedure for patients with end-stage organ disease. Over the past 70 years, tremendous progress has been made in solid organ transplantation, particularly in T-cell-targeted immunosuppression and organ allocation systems. However, humoral alloimmune responses remain a major challenge to progress. Patients with preexisting antibodies to human leukocyte antigen (HLA) are at significant disadvantages in regard to receiving a well-matched organ, moreover, those who develop anti-HLA antibodies after transplantation face a significant foreshortening of renal allograft survival. Historical therapies to desensitize patients prior to transplantation or to treat posttransplant AMR have had limited effectiveness, likely because they do not significantly reduce antibody levels, as plasma cells, the source of antibody production, remain largely unaffected. Herein, we will discuss the significance of plasma cells in transplantation, aspects of their biology as potential therapeutic targets, clinical challenges in developing strategies to target plasma cells in transplantation, and lastly, novel approaches that have potential to advance the field.

Mapping the SARS-CoV-2 spike glycoprotein-derived peptidome presented by HLA class II on dendritic cells

Understanding and eliciting protective immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an urgent priority. To facilitate these objectives, we profile the repertoire of human leukocyte antigen class II (HLA-II)-bound peptides presented by HLA-DR diverse monocyte-derived dendritic cells pulsed with SARS-CoV-2 spike (S) protein. We identify 209 unique HLA-II-bound peptide sequences, many forming nested sets, which map to sites throughout S including glycosylated regions. Comparison of the glycosylation profile of the S protein to that of the HLA-II-bound S peptides reveals substantial trimming of glycan residues on the latter, likely induced during antigen processing. Our data also highlight the receptor-binding motif in S1 as a HLA-DR-binding peptide-rich region and identify S2-derived peptides with potential for targeting by cross-protective vaccine-elicited responses. Results from this study will aid analysis of CD4⁺ T cell responses in infected individuals and vaccine recipients and have application in next-generation vaccine design.

HLA Alleles B*53:01 and C*06:02 Are Associated With Higher Risk of P. falciparum Parasitemia in a Cohort in Uganda

Variation within the HLA locus been shown to play an important role in the susceptibility to and outcomes of numerous infections, but its influence on immunity to P. falciparum malaria is unclear. Increasing evidence indicates that acquired immunity to P. falciparum is mediated in part by the cellular immune response, including NK cells, CD4 and CD8 T cells, and semi-invariant γδ T cells. HLA molecules expressed by these lymphocytes influence the epitopes recognized by P. falciparum-specific T cells, and class I HLA molecules also serve as ligands for inhibitory receptors including KIR. Here we assessed the relationship of HLA class I and II alleles to the risk of P. falciparum infection and symptomatic malaria in a cohort of 892 Ugandan children and adults followed prospectively via both active and passive surveillance. We identified two HLA class I alleles, HLA-B*53:01 and HLA-C*06:02, that were associated with a higher prevalence of P. falciparum infection. Notably, no class I or II HLA alleles were found to be associated with protection from P. falciparum parasitemia or symptomatic malaria. These findings suggest that class I HLA plays a role in the ability to restrict parasitemia, supporting an essential role for the cellular immune response in P. falciparum immunity. Our findings underscore the need for better tools to enable mechanistic studies of the T cell response to P. falciparum at the epitope level and suggest that further study of the role of HLA in regulating pre-erythrocytic stages of the P. falciparum life cycle is warranted.

Next-generation sequencing technologies: An overview

Since the days of Sanger sequencing, next-generation sequencing technologies have significantly evolved to provide increased data output, efficiencies, and applications. These next generations of technologies can be categorized based on read length. This review provides an overview of these technologies as two paradigms: short-read, or "second-generation," technologies, and long-read, or "third-generation," technologies. Herein, short-read sequencing approaches are represented by the most prevalent technologies, Illumina and Ion Torrent, and long-read sequencing approaches are represented by Pacific Biosciences and Oxford Nanopore technologies. All technologies are reviewed along with reported advantages and disadvantages. Until recently, short-read sequencing was thought to provide high accuracy limited by read-length, while long-read technologies afforded much longer read-lengths at the expense of accuracy. Emerging developments for third-generation technologies hold promise for the next wave of sequencing evolution, with the co-existence of longer read lengths and high accuracy.

HLA RNA Sequencing With Unique Molecular Identifiers Reveals High Allele-Specific Variability in mRNA Expression

The HLA gene complex is the most important single genetic factor in susceptibility to most diseases with autoimmune or autoinflammatory origin and in transplantation matching. Most studies have focused on the vast allelic variation in these genes; only a few studies have explored differences in the expression levels of HLA alleles. In this study, we quantified mRNA expression levels of HLA class I and II genes from peripheral blood samples of 50 healthy individuals. The gene- and allele-specific mRNA expression was assessed using unique molecular identifiers, which enabled PCR bias removal and calculation of the number of original mRNA transcripts. We identified differences in mRNA expression between different HLA genes and alleles. Our results suggest that HLA alleles are differentially expressed and these differences in expression levels are quantifiable using RNA sequencing technology. Our method provides novel insights into HLA research, and it can be applied to quantify expression differences of HLA alleles in various tissues and to evaluate the role of this type of variation in transplantation matching and susceptibility to autoimmune diseases.

Single-Cell RNA Sequencing of Tocilizumab-Treated Peripheral Blood Mononuclear Cells as an in vitro Model of Inflammation

COVID-19 has posed a significant threat to global health. Early data has revealed that IL-6, a key regulatory cytokine, plays an important role in the cytokine storm of COVID-19. Multiple trials are therefore looking at the effects of Tocilizumab, an IL-6 receptor antibody that inhibits IL-6 activity, on treatment of COVID-19, with promising findings. As part of a clinical trial looking at the effects of Tocilizumab treatment on kidney transplant recipients with subclinical rejection, we performed single-cell RNA sequencing of comparing stimulated PBMCs before and after Tocilizumab treatment. We leveraged this data to create an in vitro cytokine storm model, to better understand the effects of Tocilizumab in the presence of inflammation. Tocilizumab-treated cells had reduced expression of inflammatory-mediated genes and biologic pathways, particularly amongst monocytes. These results support the hypothesis that Tocilizumab may hinder the cytokine storm of COVID-19, through a demonstration of biologic impact at the single-cell level.

Class II HLA mismatch improves outcomes following haploidentical transplantation with posttransplant cyclophosphamide

HLA disparity is the major predictor of outcome following unrelated donor (UD) transplantation, where a single mismatch (mm) at the HLA-A, HLA-B, HLA-C, or HLA-DRB1 locus leads to increased mortality, and mismatching at multiple loci compounds this effect. In contrast, HLA disparity has not been shown to increase mortality in the context of haploidentical transplant using posttransplant cyclophosphamide (PTCy). To better define the consequences of loci-specific HLA mm, we analyzed 208 consecutive patients undergoing haploidentical transplantation for hematologic malignancy using PTCy at our institution (median age, 52 years [range, 19-75 years]; peripheral blood stem cell, 66%; reduced-intensity conditioning, 59%). Median follow-up was 65.4 months (range, 34.3-157.2 months). In univariate analysis, a single class II HLA mm at HLA-DR, HLA-DQ or a nonpermissive (np) HLA-DP mm had a protective effect on disease-free and overall survival (OS), primarily a result of reduced relapse risk. Furthermore, this survival effect was cumulative, so that patients with 3 class II mm (HLA-DR, HLA-DQ, and np HLA-DP) had the best OS. In multivariate analysis, HLA-DR mm and np HLA-DP mm were both independently associated with improved OS (hazard ratio [HR], 0.43; P =.001; and HR, 0.47; P =.011, respectively). In contrast, single or multiple mm at HLA-A, HLA-B, or HLA-C loci had no effect on acute graft-versus-host disease (GVHD), nonrelapse mortality (NRM), relapse, or survival, although the presence of an HLA-A mm was associated with increased chronic GVHD incidence. The association of class II mm with lower relapse occurred without a corresponding increase in NRM or acute or chronic GVHD. These findings will require validation in larger registry studies.

Single-cell RNA sequencing of Tocilizumab-treated peripheral blood mononuclear cells as an in vitro model of inflammation

COVID-19 has posed a significant threat to global health. Early data has revealed that IL-6, a key regulatory cytokine, plays an important role in the cytokine storm of COVID-19. Multiple trials are therefore looking at the effects of Tocilizumab, an IL-6 receptor antibody that inhibits IL-6 activity, on treatment of COVID-19, with promising findings. As part of a clinical trial looking at the effects of Tocilizumab treatment on kidney transplant recipients with subclinical rejection, we performed single-cell RNA sequencing of comparing stimulated PBMCs before and after Tocilizumab treatment. We leveraged this data to create an in vitro cytokine storm model, to better understand the effects of Tocilizumab in the presence of inflammation. Tocilizumab-treated cells had reduced expression of inflammatory-mediated genes and biologic pathways, particularly amongst monocytes. These results support the hypothesis that Tocilizumab may hinder the cytokine storm of COVID-19, through a demonstration of biologic impact at the single-cell level.

Mapping the SARS-CoV-2 spike glycoprotein-derived peptidome presented by HLA class II on dendritic cells

Understanding and eliciting protective immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an urgent priority. To facilitate these objectives, we have profiled the repertoire of human leukocyte antigen class II (HLA-II)-bound peptides presented by HLA-DR diverse monocyte-derived dendritic cells pulsed with SARS-CoV-2 spike (S) protein. We identify 209 unique HLA-II-bound peptide sequences, many forming nested sets, which map to sites throughout S including glycosylated regions. Comparison of the glycosylation profile of the S protein to that of the HLA-II-bound S peptides revealed substantial trimming of glycan residues on the latter, likely introduced during antigen processing. Our data also highlight the receptor-binding motif in S1 as a HLA-DR-binding peptide-rich region. Results from this study have application in vaccine design, and will aid analysis of CD4+ T cell responses in infected individuals and vaccine recipients.

Identification of Novel Alleles and Structural Haplotypes of Major Histocompatibility Complex Class I and DRB Genes in Domestic Cat (Felis catus) by a Newly Developed NGS-Based Genotyping Method

The major histocompatibility complex (MHC) is a highly polymorphic and duplicated genomic region that encodes transplantation and immune regulatory molecules. Although it is well-known that particular MHC allelic polymorphisms and haplotypes are genetically relate to immune-mediated diseases detailed information of the cat MHC (Feline Leukocyte Antigen; FLA) genetic and haplotypic structure and diversity is limited in comparison to humans and many other species. In this study, to better understand the degree and types of allele and allelic haplotype diversity of FLA-class I (FLA-I) and FLA-DRB loci in domestic cats, we identified six expressible FLA-I loci in peripheral white blood cells by in silico estimation of the coding exons and NGS-based amplicon sequencing using five unrelated cats. We then used a newly developed NGS-based genotyping method to genotype and annotate 32 FLA-I and 16 FLA-DRB sequences in two families of 20 domestic cats. A total of 14 FLA-I and seven FLA-DRB were identified as novel polymorphic sequences. Phylogenetic analyses grouped the sequences into six FLA-I (FLA-E/H/K, FLA-A, FLA-J, FLA-L, FLA-O and a tentatively named FLA-E/H/K_Rec) and four FLA-DRB (FLA-DRB1, FLA-DRB3, FLA-DRB4, and FLA-DRB5) lineages. Pedigree analysis of two cat families revealed eight distinct FLA structural haplotypes (Class I - DRB) with five to eight FLA-I and two to three FLA-DRB transcribed loci per haplotype. It is evident that the eight FLA haplotypes were generated by gene duplications and deletions, and rearrangements by genetic recombination with the accumulation and/or inheritance of novel polymorphisms. These findings are useful for further genetic diversity analysis and disease association studies among cat breeds and in veterinary medicine.