The MHC in the era of next-generation sequencing: Implications for bridging structure with function

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

HLA-DPA1*02:01:25, a novel allele with a synonymous transition in exon 2 identified by next-generation sequencing

HLA-DPA1*02:01:25 differs from DPA1*02:01:01:02 by a synonymous transition in exon 2.

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.

De novo genome assembly depicts the immune genomic characteristics of cattle

Immunogenomic loci remain poorly understood because of their genetic complexity and size. Here, we report the de novo assembly of a cattle genome and provide a detailed annotation of the immunogenomic loci. The assembled genome contains 143 contigs (N50 ~ 74.0 Mb). In contrast to the current reference genome (ARS-UCD1.2), 156 gaps are closed and 467 scaffolds are located in our assembly. Importantly, the immunogenomic regions, including three immunoglobulin (IG) loci, four T-cell receptor (TR) loci, and the major histocompatibility complex (MHC) locus, are seamlessly assembled and precisely annotated. With the characterization of 258 IG genes and 657 TR genes distributed across seven genomic loci, we present a detailed depiction of immune gene diversity in cattle. Moreover, the MHC gene structures are integrally revealed with properly phased haplotypes. Together, our work describes a more complete cattle genome, and provides a comprehensive view of its complex immune-genome.

Simultaneous detection of DNA variation and methylation at HLA class II locus and immune gene promoters using targeted SureSelect Methyl-Sequencing

The Human Leukocyte Antigen (HLA) locus associates with a variety of complex diseases, particularly autoimmune and inflammatory conditions. The HLA-DR15 haplotype, for example, confers the major risk for developing Multiple Sclerosis in Caucasians, pinpointing an important role in the etiology of this chronic inflammatory disease of the central nervous system. In addition to the protein-coding variants that shape the functional HLA-antigen-T cell interaction, recent studies suggest that the levels of HLA molecule expression, that are epigenetically controlled, also play a role in disease development. However, deciphering the exact molecular mechanisms of the HLA association has been hampered by the tremendous genetic complexity of the locus and a lack of robust approaches to investigate it. Here, we developed a method to specifically enrich the genomic DNA from the HLA class II locus (chr6:32,426,802-34,167,129) and proximal promoters of 2,157 immune-relevant genes, utilizing the Agilent RNA-based SureSelect Methyl-Seq Capture related method, followed by sequencing to detect genetic and epigenetic variation. We demonstrated successful simultaneous detection of the genetic variation and quantification of DNA methylation levels in HLA locus. Moreover, by the detection of differentially methylated positions in promoters of immune-related genes, we identified relevant pathways following stimulation of cells. Taken together, we present a method that can be utilized to study the interplay between genetic variance and epigenetic regulation in the HLA class II region, potentially, in a wide disease context.

Detection and characterization of the novel HLA-DPA1*02:66:02N allele, with a premature stop codon in exon 2

The failure to identify HLA null alleles in bone marrow transplantation could be life-threatening because this could result in an HLA mismatch with the ability to trigger the graft-vs-host disease (GVHD) and to reduce patient's survival. In this report we describe the identification and characterization of the novel HLA-DPA1*02:66:02N allele with a non-sense codon in exon 2. This new allele was discovered in two unrelated bone marrow donors during routine HLA-typing using next-generation sequencing (NGS). DPA1*02:66:02N is homologous to DPA1*02:01:01:03 with a single nucleotide difference in exon 2, codon 50, where the replacement of C located at genomic position 3825 by T, causes the formation of a premature stop codon (TGA), resulting in a null allele. This description illustrates the benefits of HLA typing by NGS since it permits to reduce ambiguities, identify new alleles, analyze multiple HLA loci and improve transplantation outcome.

Identification of the novel HLA-DQA1*02:01:09:01 allele by two different next-generation sequencing platforms

A synonymous nucleotide substitution in exon 3 results in the novel HLA-DQA1*02:01:09:01 allele.

Genetics of neurosarcoidosis

Sarcoidosis is a systemic, inflammatory, granulomatous disease characterized by great variability in organ involvement, clinical course, and severity. While pulmonary manifestations are almost universal, the central and peripheral nervous systems can also be affected. Neurosarcoidosis occurs in ∼5-15% of cases and is among the manifestations with the highest morbidity and mortality. It is known that sarcoidosis has genetic underpinnings and while multiple studies aimed at identifying associations to sarcoidosis susceptibility and prognosis, very few studies have focused on neurosarcoidosis. This review summarizes the genetic studies to date, compares and contrasts those findings with other genetic effects in sarcoidosis, and offers ideas for moving the field forward.

Genotyping of canine MHC gene DLA-88 by next-generation sequencing reveals high frequencies of new allele discovery and gene duplication

Dogs have served as one of the most reliable preclinical models for a variety of diseases and treatments, including stem/progenitor cell transplantation. At the genetic epicenter of dog transplantation models, polymorphic major histocompatibility complex (MHC) genes are most impactful on transplantation success. Among the canine class I and class II genes, DLA-88 has been best studied in transplantation matching and outcomes, with 129 DLA-88 alleles identified. In this study we developed and tested a next generation (NGS) sequencing protocol for rapid identification of DLA-88 genotypes in dogs and compared the workflow and data generated with an established DLA-88 Sanger sequencing protocol that has been in common prior use for clinical studies. By testing the NGS protocol on a random population of 382 dogs, it was possible to demonstrate superior efficacy based on laboratory execution and overall cost. In addition, NGS proved far more effective at discovering new alleles and detecting multiple alleles associated with gene duplication. A total of 51 new DLA-88 alleles are reported here. This rate of new allele discovery indicates that a large pool of yet un-discovered DLA-88 alleles exists in the domestic dog population. In addition, more than 46% of dogs carried three or more copies of DLA-88, further emphasizing the need for more sensitive and cost-effective DLA typing methodology for the dog clinical model.

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.

The novel HLA-DQB1*03:493 allele, the first with glutamic acid at position-10 in the leader peptide

A nonsynonymous nucleotide substitution in exon 1 results in the novel HLA-DQB1*03:493 allele.

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.

Race and Survival in Unrelated Hematopoietic Cell Transplantation

Survival after hematopoietic cell transplantation depends on race/ethnicity and histocompatibility (HLA) between the patient and transplant donor. HLA sequence variation is a genetic construct of continental populations, but its role in accounting for racial disparities of transplant outcome is unknown. To determine disparities in transplant survivorship among patients of diverse race while accounting for patient and donor HLA variation. A total of 26,945 self-described Japanese, U.S. Asian, White, Hispanic, and Black patients received an unrelated donor transplant for the treatment of a life-threatening blood disorder. The risk of mortality with and without adjustment for known HLA risk factors (number and location of donor mismatches; patient HLA-B leader genotype and HLA-DRβ peptide-binding motif) was studied using multivariable models. Survival after HLA-matched transplantation for patients with optimal leader and peptide-binding features was estimated for each race, as was the improvement in survival over calendar-year time by considering year of transplantation as a continuous linear variable. The number, location, and nature of donor HLA mismatches and the frequency of patient HLA-B and HLA-DRB1 sequence motifs differed by race. Japanese patients had superior survival compared to other races without consideration of HLA. After HLA adjustment, three mortality risk strata were identified: Japanese and U.S. Asian (low-risk); White and Hispanic (intermediate-risk), and Black patients (high-risk). Survival for patients with optimal donor and HLA characteristics was superior for Japanese, intermediate for U.S. Asian, White, and Hispanic, and lowest for Black patients. Five-year increments of transplant year were associated with greater decreases in mortality hazards for Black and Hispanic patients than for Japanese, U.S. Asian and White patients. Transplant survivorship disparities are influenced by HLA as a genetic construct of race. A more complete understanding of the factors that influence transplant outcomes provides opportunities to narrow disparities for future patients.

Comparative gene expression profile in circulating PBMCs of Bos indicus and crossbred cattle to understand disease tolerance mechanism

The present investigation was performed to compare the global gene expression profile in peripheral blood mononuclear cells (PBMCs) of Bos indicus and crossbred (Bos taurus × B. indicus) cattle. Previously, several studies revealed the disease tolerance potential of B. indicus cattle but underlying genetic mechanism is still not fully explored. The PBMCs model was used for this investigation as it plays crucial role in the immune system regulation. Transcriptomic analysis revealed total 6767 significantly differentially expressed transcripts (fold change (absolute) >2.0, p < .05). In addition, 4149 transcripts were upregulated, 2618 transcripts were downregulated and fold change (absolute) of differentially expressed transcript varied from -223.32 to 213.63. Functional annotation analysis of differentially expressed genes confirmed their role in various molecular pathways viz. innate immune response, antigen processing and presentation, MHC protein complex, defense response to bacterium, regulation of immune response, positive regulation of JAK-STAT cascade, cytoskeletal protein binding, etc. Protein-protein interaction network analysis provided understanding of inter-relationship of immune genes with differentially expressed genes. In conclusion, this study could provide comprehensive information about the dysregulated genes and biological pathways in PBMCs which might be responsible for disease tolerance in B. indicus cattle.

Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines

Introduction

Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs.

Methods

We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments.

Results

We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes.

Discussion

Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Validation of Single Nucleotide Variant Assays for Human Leukocyte Antigen Haplotypes HLA-B*15:02 and HLA-A*31:01 Across Diverse Ancestral Backgrounds

The human leukocyte antigen haplotypes HLA-B*15:02 and HLA-A*31:01 have been linked to life-threatening adverse drug reactions to the anticonvulsants carbamazepine and oxcarbazepine. Identification of these haplotypes via pharmacogenetic techniques facilitates implementation of precision medicine to prevent such reactions. Using reference samples from diverse ancestral origins, we investigated the test analytical validity (i.e., ability to detect whether or not the haplotypes were present or absent) of TaqMan assays for single nucleotide variants previously identified as potentially being able to "tag" these haplotypes. A TaqMan custom assay for rs10484555 and an inventoried assay for rs17179220 and were able to identify with 100% sensitivity and 100% specificity HLA-B*15:02 and HLA-A*31:01 respectively. A custom assay for rs144012689 that takes into account a neighboring single nucleotide variant with manual calling was also able to identify HLA-B*15:02 with 100% sensitivity and 100% specificity. A custom assay for rs106235 identified HLA-A*31:01 with 100% sensitivity and 95% specificity. The slight reduction in specificity for the latter was owing to another haplotype (HLA-A*33:03) also being detected. While any positive call using the rs106235 assay could therefore be further investigated, as the presence of the HLA-A*31:01 haplotype confers adverse drug reaction risk, the absence of false negatives (indexed by sensitivity) is more important than false positives. In summary, we present validated TaqMan assay methodology for efficient detection of HLA haplotypes HLA-B*15:02 and HLA-A*31:01. Our data are relevant for other genotyping technologies that identify, or have the potential to identify, these haplotypes using single nucleotide variants.

Exonic SNP in MHC-DMB2 is associated with gene expression and humoral immunity in Japanese quails

The DMB2 gene is widely expressed at high levels in avian. This gene plays an important role in humoral immunity. The aim of this study was to investigate the effects of 361 G > C Single nucleotide polymorphism (SNP) on DMB2 protein structure and gene expression to determine how the 361 G > C SNP affects humoral immune response in Japanese quails. 0.2 mL of 5% sheep red blood cell (SRBC) was injected into breast muscle of 130 Japanese quails on 28 days. After DNA extraction, PCR was carried out to amplify a 333-base pair DNA fragment from the exon 2 of DMB2 gene. The pattern of all samples was determined through RFLP technique. PCR-RFLP results identified two alleles segregating (C, G) as three genotypes (CC, CG and GG) in Japanese Quails. The antibody response to SRBC with CC genotype was significantly higher than the CG and GG genotypes (P < 0.01). In silico analysis showed that the 361 G > C SNP has no effect on the physicochemical properties and 3D structure. The results of RT-qPCR indicated that the effect of genotype on gene expression is significant, so that the expression of CC genotype is more than CG and GG genotype. It can be inferred that the 361 G > C SNP in the exon 2 of MHC-DMB2 gene is not desirable. This mutation decreases humoral immune response by reducing DMB2 gene expression.

Impact of Previously Unrecognized HLA Mismatches Using Ultrahigh Resolution Typing in Unrelated Donor Hematopoietic Cell Transplantation

PURPOSE

Ultrahigh resolution (UHR) HLA matching is reported to result in better outcomes following unrelated donor hematopoietic cell transplantation, improving survival and reducing post-transplant complications. However, most studies included relatively small numbers of patients. Here we report the findings from a large, multicenter validation study.

METHODS

UHR HLA typing was available on 5,140 conventionally 10 out of 10 HLA-matched patients with malignant disease transplanted between 2008 and 2017.

RESULTS

After UHR HLA typing, 82% of pairs remained 10 out of 10 UHR-matched; 12.3% of patients were 12 out of 12 UHR HLA-matched. Compared with 12 out of 12 UHR-matched patients, probabilities of grade 2-4 acute graft-versus-host disease (aGVHD) were significantly increased with UHR mismatches (overall P = .0019) and in those patients who were HLA-DPB1 T-cell epitope permissively mismatched or nonpermissively mismatched (overall P = .0011). In the T-cell-depleted subset, the degree of UHR HLA mismatch was only associated with increased transplant-related mortality (TRM) (overall P = .0068). In the T-cell-replete subset, UHR HLA matching was associated with a lower probability of aGVHD (overall P = .0020); 12 out of 12 UHR matching was associated with reduced TRM risk when compared with HLA-DPB1 T-cell epitope permissively mismatched patients, whereas nonpermissive mismatching resulted in a greater risk (overall P = .0003).

CONCLUSION

This study did not confirm that UHR 12 out of 12 HLA matching increases the probability of overall survival but does demonstrate that aGVHD risk, and in certain settings TRM, is lowest in UHR HLA-matched pairs and thus warrants consideration when multiple 10 out of 10 HLA-matched donors of equivalent age are available.

High throughput analysis of MHC-I and MHC-DR diversity of Brazilian cattle populations

The major histocompatibility complex (MHC) contains many genes that play key roles in initiating and regulating immune responses. This includes the polymorphic MHCI and MHCII genes that present epitopes to CD8+ and CD4+ T-cells, respectively. Consequently, the characterisation of the repertoire of MHC genes is an important component of improving our understanding of the genetic variation that determines the outcomes of immune responses. In cattle, MHC (BoLA) research has predominantly focused on Holstein-Friesian animals (as the most economically important breed globally), although the development of high-throughput approaches has allowed the BoLA-DRB3 repertoire to be studied in a greater variety of breeds. In a previous study we reported on the development of a MiSeq-based method to enable high-throughput and high-resolution analysis of bovine MHCI repertoires. Herein, we report on the expansion of this methodology to incorporate analysis of the BoLA-DRB3 and its application to analyse MHC diversity in a large cohort of cattle from Brazil (>500 animals), including representatives from the three major Bos indicus breeds present in Brazil - Guzerat, Gir and Nelore. This large-scale description of paired MHCI-DRB3 repertoires in Bos indicus cattle has identified a small number of novel DRB3 alleles, a large number of novel MHCI alleles and haplotypes, and provided novel insights into MHCI-MHCII association - further expanding our knowledge of bovine MHC diversity.

The changing landscape of HLA typing: Understanding how and when HLA typing data can be used with confidence from bench to bedside

Human leukocyte antigen (HLA) genes are extraordinary for their extreme diversity and widespread impact on human health and disease. More than 30,000 HLA alleles have been officially named and more alleles continue to be discovered at a rapid pace. HLA typing systems which have been developed to detect HLA diversity have advanced rapidly and are revolutionizing our understanding of HLA's clinical importance. However, continuous improvements in knowledge and technology have created challenges for clinicians and scientists. This review explains how differences in HLA typing systems can impact the HLA types that are assigned. The consequences of differences in laboratory testing methods and reference databases are described. The challenges of using HLA types that are not equivalent are illustrated. A fundamental understanding of the continual expansion of our understanding of HLA diversity and limitations in some of the typing data is essential for using typing data appropriately in clinical and research settings.

BSHI guideline: HLA matching and donor selection for haematopoietic progenitor cell transplantation

A review of the British Society for Histocompatibility and Immunogenetics (BSHI) Guideline 'HLA matching and donor selection for haematopoietic progenitor cell transplantation' published in 2016 was undertaken by a BSHI appointed writing committee. Literature searches were performed and the data extracted were presented as recommendations according to the GRADE nomenclature.

Genome-Wide Association Study of Root System Development at Seedling Stage in Rice

Root network structure plays a crucial role in growth and development processes in rice. Longer, more branched root structures help plants to assimilate water and nutrition from soil, support robust plant growth, and improve resilience to stresses such as disease. Understanding the molecular basis of root development through screening of root-related traits in rice germplasms is critical to future rice breeding programs. This study used a small germplasm collection of 137 rice varieties chosen from the Korean rice core set (KRICE_CORE) to identify loci linked to root development. Two million high-quality single nucleotide polymorphisms (SNPs) were used as the genotype, with maximum root length (MRL) and total root weight (TRW) in seedlings used as the phenotype. Genome-wide association study (GWAS) combined with Principal Components Analysis (PCA) and Kinship matrix analysis identified four quantitative trait loci (QTLs) on chromosomes 3, 6, and 8. Two QTLs were linked to MRL and two were related to TRW. Analysis of Linkage Disequilibrium (LD) decay identified a 230 kb exploratory range for detection of candidate root-related genes. Candidates were filtered using RNA-seq data, gene annotations, and quantitative real-time PCR (qRT-PCR), and five previously characterized genes related to root development were identified, as well as four novel candidate genes. Promoter analysis of candidate genes showed that LOC_Os03g08880 and LOC_Os06g13060 contained SNPs with the potential to impact gene expression in root-related promoter motifs. Haplotype analysis of candidate genes revealed diverse haplotypes that were significantly associated with phenotypic variation. Taken together, these results indicate that LOC_Os03g08880 and LOC_Os06g13060 are strong candidate genes for root development functions. The significant haplotypes identified in this study will be beneficial in future breeding programs for root improvement.

Critical Review of Existing MHC I Immunopeptidome Isolation Methods

Major histocompatibility complex class I (MHC I) plays a crucial role in the development of adaptive immune response in vertebrates. MHC molecules are cell surface protein complexes loaded with short peptides and recognized by the T-cell receptors (TCR). Peptides associated with MHC are named immunopeptidome. The MHC I immunopeptidome is produced by the proteasome degradation of intracellular proteins. The knowledge of the immunopeptidome repertoire facilitates the creation of personalized antitumor or antiviral vaccines. A huge number of publications on the immunopeptidome diversity of different human and mouse biological samples-plasma, peripheral blood mononuclear cells (PBMCs), and solid tissues, including tumors-appeared in the scientific journals in the last decade. Significant immunopeptidome identification efficiency was achieved by advances in technology: the immunoprecipitation of MHC and mass spectrometry-based approaches. Researchers optimized common strategies to isolate MHC-associated peptides for individual tasks. They published many protocols with differences in the amount and type of biological sample, amount of antibodies, type and amount of insoluble support, methods of post-fractionation and purification, and approaches to LC-MS/MS identification of immunopeptidome. These parameters have a large impact on the final repertoire of isolated immunopeptidome. In this review, we summarize and compare immunopeptidome isolation techniques with an emphasis on the results obtained.

Identification of a single nucleotide deletion in the novel HLA-DQB1*06:379N allele, detected by Polymerase Chain Reaction-Sequence Based Typing but not by Next Generation Sequencing

In this report we describe the characterization of a novel null HLA-DQB1 allele, detected by polymerase chain reaction-sequence-based typing but not by next-generation sequencing (NGS). The new allele, HLA-DQB1*06:379N, was discovered in an Italian patient with acute myeloid leukemia and also in one of her healthy siblings. The new allele is largely homologous to DQB1*06:02:01:01 with a T deletion in exon 2, in codon 11, which causes a frameshift and the formation of an unusual and premature TGA STOP in codon 27. This case report underlines the importance of not removing Sanger method sequencing from routine use for high-resolution HLA typing, but to maintain it together with NGS technology.

Cell Communication-mediated Nonself-Recognition and -Intolerance in Representative Species of the Animal Kingdom

Why has histo-incompatibility arisen in evolution and can cause self-intolerance? Compatible/incompatible reactions following natural contacts between genetically-different (allogeneic) colonies of marine organisms have inspired the conception that self-nonself discrimination has developed to reduce invasion threats by migratory foreign germ/somatic stem cells, in extreme cases resulting in conquest of the whole body by a foreign genome. Two prominent model species for allogeneic discrimination are the marine invertebrates Hydractinia (Cnidaria) and Botryllus (Ascidiacea). In Hydractinia, self-nonself recognition is based on polymorphic surface markers encoded by two genes (alr1, alr2), with self recognition enabled by homophilic binding of identical ALR molecules. Variable expression patterns of alr alleles presumably account for the first paradigm of autoaggression in an invertebrate. In Botryllus, self-nonself recognition is controlled by a single polymorphic gene locus (BHF) with hundreds of codominantly expressed alleles. Fusion occurs when both partners share at least one BHF allele while rejection develops when no allele is shared. Molecules involved in allorecognition frequently contain immunoglobulin or Ig-like motifs, case-by-case supplemented by additional molecules enabling homophilic interaction, while the mechanisms applied to destroy allogeneic grafts or neighbors include taxon-specific tools besides common facilities of natural immunity. The review encompasses comparison with allorecognition in mammals based on MHC-polymorphism in transplantation and following feto-maternal cell trafficking.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.

Importance of the Major Histocompatibility Complex (Swine Leukocyte Antigen) in Swine Health and Biomedical Research

In pigs, the major histocompatibility complex (MHC), or swine leukocyte antigen (SLA) complex, maps to Sus scrofa chromosome 7. It consists of three regions, the class I and class III regions mapping to 7p1.1 and the class II region mapping to 7q1.1. The swine MHC is divided by the centromere, which is unique among mammals studied to date. The SLA complexspans between 2.4 and 2.7 Mb, depending on haplotype, and encodes approximately 150 loci, with at least 120 genes predicted to be functional. Here we update the whole SLA complex based on the Sscrofa11.1 build and annotate the organization for all recognized SLA genes and their allelic sequences. We present SLA nomenclature and typing methods and discuss the expression of SLA proteins, as well as their role in antigen presentation and immune, disease, and vaccine responses. Finally, we explore the role of SLA genes in transplantation and xenotransplantation and their importance in swine biomedical models.

Human leucocyte antigen diversity: A biological gift to escape infections, no longer a barrier for haploidentical Hemopoietic Stem Cell Transplantation

Since the beginning of life, every multicellular organism appeared to have a complex innate immune system although the adaptive immune system, centred on lymphocytes bearing antigen receptors generated by somatic recombination, arose in jawed fish approximately 500 million years ago. The major histocompatibility complex MHC, named the Human leucocyte antigen (HLA) system in humans, represents a vital function structure in the organism by presenting pathogen-derived peptides to T cells as the main initial step of the adaptive immune response. The huge level of polymorphism observed in HLA genes definitely reflects selection, favouring heterozygosity at the individual or population level, in a pathogen-rich environment, although many are located in introns or in exons that do not code for the antigen-biding site of the HLA. Over the past three decades, the extent of allelic diversity at HLA loci has been well characterized using high-resolution HLA-DNA typing and the number of new HLA alleles, produced through next-generation sequencing methods, is even more rapidly increasing. The level of the HLA system polymorphism represents an obstacle to the search of potential compatible donors for patients affected by haematological disease proposed for a hematopoietic stem cell transplant (HSCT). Data reported in literature clearly show that antigenic and/or allelic mismatches between related or unrelated donors and patients influences the successful HSCT outcome. However, the recent development of the new transplant strategy based on the choice of haploidentical donors for HSCT is questioning the role of HLA compatibility, since the great HLA disparities present do not worsen the overall clinical outcome. Nowadays, NGS has contributed to define at allelic levels the HLA polymorphism and solve potential ambiguities. However, HLA functions and tissue typing probably need to be further investigated in the next future, to understand the reasons why in haploidentical transplants the presence of a whole mismatch haplotype between donors and recipients, both the survival rate and the incidence of acute GvHD or graft rejection are similar to those reported for unrelated HSCTs.

MHC Genomics and Disease: Looking Back to Go Forward

Ancestral haplotypes are conserved but extremely polymorphic kilobase sequences, which have been faithfully inherited over at least hundreds of generations in spite of migration and admixture. They carry susceptibility and resistance to diverse diseases, including deficiencies of CYP21 hydroxylase (47.1) and complement components (18.1), as well as numerous autoimmune diseases (8.1). The haplotypes are detected by segregation within ethnic groups rather than by SNPs and GWAS. Susceptibility to some other diseases is carried by specific alleles shared by multiple ancestral haplotypes, e.g., ankylosing spondylitis and narcolepsy. The difference between these two types of association may explain the disappointment with many GWAS. Here we propose a pathway for combining the two different approaches. SNP typing is most useful after the conserved ancestral haplotypes have been defined by other methods.

Factors affecting HLA expression: A review

The detection and semiquantitative measurement of circulating human leucocyte antigen (HLA)-specific antibodies is essential for the management of patients before and after transplantation. In addition, the pretransplant cross-match to assess the reactivity of recipient HLA antibody against donor lymphocytes has long been the gold standard to prevent hyperacute rejection. Whilst both of these tests assume that recipient HLA-specific antibody is the only variable in the assessment of transplant risk, this is not the case. Transplant immunologists recognize that some HLA antigens are expressed at levels a magnitude lower than others (e.g., HLA-C, HLA-DQ), but within loci, and between different cell types there are many factors that influence HLA expression in both resting and activated cells. HLA is not usually expressed without the specific promoter proteins NLRC5, for HLA class I, and CIITA, for class II. The quantity of HLA protein production is then affected by factors including promoter region polymorphisms, alternative exon splice sites, methylation and microRNA-directed degradation. Different loci are influenced by multiple combinations of these control mechanisms making prediction of HLA regulation difficult, but an ability to measure the cellular expression of each HLA antigen, in conjunction with knowledge of circulating HLA-specific antibody, would lead to a more informed algorithm to assess transplant risk.