Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II

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.

The 18th International HLA & Immunogenetics workshop project report: Creating fully representative MHC reference haplotypes

A fundamental endeavor of the International Histocompatibility and Immunogenetics Workshop (IHIW) was assembling a collection of DNA samples homozygous through the MHC genomic region. This collection proved invaluable for assay development in the histocompatibility and immunogenetics field, for generating the human reference genome, and furthered our understanding of MHC diversity. Defined by their HLA-A, -B, -C and -DRB1 alleles, the combined frequency of the haplotypes from these individuals is ~20% in Europe. Thus, a significant proportion of MHC haplotypes, both common and rare throughout the world, and including many associated with disease, are not yet represented. In this workshop component, we are collecting the next generation of MHC -homozygous samples, to expand, diversify and modernize this critical community resource that has been foundational to the field. We asked laboratories worldwide to identify samples homozygous through all HLA class I and/or HLA class II genes, or through whole-genome SNP genotyping or sequencing, to have extensive homozygosity tracts within the MHC region. The focus is non-Europeans or those having HLA haplotypes less common in Europeans. Through this effort, we have obtained samples from 537 individuals representing 294 distinct haplotypes, as determined by their HLA class I and II alleles, and an additional 50 haplotypes distinct in HLA class I or II alleles. Although we have expanded the diversity, many populations remain underrepresented, particularly from Africa, and we encourage further participation. The data will serve as a resource for investigators seeking to characterize variation across the MHC genomic region for disease and population studies.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Common risk alleles for schizophrenia within the major histocompatibility complex predict white matter microstructure

Recent research has highlighted the role of complement genes in shaping the microstructure of the brain during early development, and in contributing to common allele risk for Schizophrenia. We hypothesised that common risk variants for schizophrenia within complement genes will associate with structural changes in white matter microstructure within tracts innervating the frontal lobe. Results showed that risk alleles within the complement gene set, but also intergenic alleles, significantly predict axonal density in white matter tracts connecting frontal cortex with parietal, temporal and occipital cortices. Specifically, risk alleles within the Major Histocompatibility Complex region in chromosome 6 appeared to drive these associations. No significant associations were found for the orientation dispersion index. These results suggest that changes in axonal packing - but not in axonal coherence - determined by common risk alleles within the MHC genomic region - including variants related to the Complement system - appear as a potential neurobiological mechanism for schizophrenia.

Investigating the genetic makeup of the major histocompatibility complex (MHC) in the United Arab Emirates population through next-generation sequencing

The Human leukocyte antigen (HLA) molecules are central to immune response and have associations with the phenotypes of various diseases and induced drug toxicity. Further, the role of HLA molecules in presenting antigens significantly affects the transplantation outcome. The objective of this study was to examine the extent of the diversity of HLA alleles in the population of the United Arab Emirates (UAE) using Next-Generation Sequencing methodologies and encompassing a larger cohort of individuals. A cohort of 570 unrelated healthy citizens of the UAE volunteered to provide samples for Whole Genome Sequencing and Whole Exome Sequencing. The definition of the HLA alleles was achieved through the application of the bioinformatics tools, HLA-LA and xHLA. Subsequently, the findings from this study were compared with other local and international datasets. A broad range of HLA alleles in the UAE population, of which some were previously unreported, was identified. A comparison with other populations confirmed the current population's unique intertwined genetic heritage while highlighting similarities with populations from the Middle East region. Some disease-associated HLA alleles were detected at a frequency of > 5%, such as HLA-B*51:01, HLA-DRB1*03:01, HLA-DRB1*15:01, and HLA-DQB1*02:01. The increase in allele homozygosity, especially for HLA class I genes, was identified in samples with a higher level of genome-wide homozygosity. This highlights a possible effect of consanguinity on the HLA homozygosity. The HLA allele distribution in the UAE population showcases a unique profile, underscoring the need for tailored databases for traditional activities such as unrelated transplant matching and for newer initiatives in precision medicine based on specific populations. This research is part of a concerted effort to improve the knowledge base, particularly in the fields of transplant medicine and investigating disease associations as well as in understanding human migration patterns within the Arabian Peninsula and surrounding regions.

Full resolution HLA and KIR genes annotation for human genome assemblies

The HLA (Human Leukocyte Antigen) genes and the KIR (Killer cell Immunoglobulin-like Receptor) genes are critical to immune responses and are associated with many immune-related diseases. Located in highly polymorphic regions, they are hard to be studied with traditional short-read alignment-based methods. Although modern long-read assemblers can often assemble these genes, using existing tools to annotate HLA and KIR genes in these assemblies remains a non-trivial task. Here, we describe Immuannot, a new computation tool to annotate the gene structures of HLA and KIR genes and to type the allele of each gene. Applying Immuannot to 56 regional and 212 whole-genome assemblies from previous studies, we annotated 9,931 HLA and KIR genes and found that almost half of these genes, 4,068, had novel sequences compared to the current Immuno Polymorphism Database (IPD). These novel gene sequences were represented by 2,664 distinct alleles, some of which contained non-synonymous variations resulting in 92 novel protein sequences. We demonstrated the complex haplotype structures at the two loci and reported the linkage between HLA/KIR haplotypes and gene alleles. We anticipate that Immuannot will speed up the discovery of new HLA/KIR alleles and enable the association of HLA/KIR haplotype structures with clinical outcomes in the future.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

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.

HLA Genes: A Hallmark of Functional Genetic Variation and Complex Evolution

The major histocompatibility complex (MHC) with its highly polymorphic HLA genes represents one of the most intensely studied genomic regions in the genome. MHC proteins play a key role in antigen-specific immunity and are associated with a wide range of complex diseases. Despite decades of research and many advances in the field, the characterization and interpretation of its genetic and genomic variability remain challenging. Here an overview is provided of the MHC, the nature of its exceptional variability, and the complex evolutionary processes assumed to drive this variability. Highlighted are also recent advances in the field that promise to improve our understanding of the variability in the MHC and in antigen-specific immunity more generally.

High-throughput complement component 4 genomic sequence analysis with C4Investigator

The complement component 4 gene loci, composed of the C4A and C4B genes and located on chromosome 6, encodes for complement component 4 (C4) proteins, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. C4A and C4B gene loci exhibit copy number variation, with each composite gene varying between 0 and 5 copies per haplotype. C4A and C4B genes also vary in size depending on the presence of the human endogenous retrovirus (HERV) in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which affects expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4A and C4B copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4A and C4B variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4A and C4B sequences from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines the overall gene copy numbers, as well as C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S). Additionally, C4Ivestigator reports the full overall C4A and C4B aligned sequence, enabling nucleotide level analysis. To demonstrate the utility of this workflow we have analyzed C4A and C4B variation in the 1000 Genomes Project Data set, showing that these genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

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

Background

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

Methods

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

Results

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

Results

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

High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

HLA-F transcriptional and protein differential expression according to its genetic polymorphisms

Many specificities single out HLA-F: its structure, expression regulation at cell membrane and function. HLA-F mRNA is detected in the most cell types and the protein is localized in the ER and Golgi apparatus. When expressed at cell surface, HLA-F may be associated to β2-microglobulin and peptide or expressed as an open-conformer molecule. HLA-F reaches the membrane upon activation of different primary cell types and cell-lines. HLA-F has its highest affinity for the KIR3DS1-activating NK receptor, but also binds inhibitory immune receptors. Some studies reported that HLA-F expression is associated with its genotype. Higher HLA-F mRNA expression associated with F*01:01:02, and 3 noncoding SNPs, rs1362126, rs2523405, and rs2523393, located in HLA-F-AS1 or upstream the HLA-F sequence were associated with HLA-F mRNA expression. Given the implication of HLA-F in many clinical setting, and the undisclosed process of its expression regulation, we aim to confirm the effect of the aforementioned SNPs with HLA-F transcriptional and protein expression. We analyzed the distribution, frequency and linkage disequilibrium of these SNPs at worldwide scale in the 1000 Genomes Project samples. Influence on the genotype of each SNP on HLA-F expression was explored using RNAseq data from the 1000 Genomes Project, and using Q-PCR and intracellular cytometry in PBMC from healthy individuals. Our results show that the SNPs under studied displayed remarkably different allelic proportion according to geography and confirm that rs1362126, rs2523405, and rs2523393 displayed the most concordant results, with the highest effect size and a double-dose effect.

HLA-OLI: A new MHC class I pseudogene and HLA-Y are located on a 60 kb indel in the human MHC between HLA-W and HLA-J

Analysis of publicly available whole-genome sequence data from the Human Pangenome Project and the 1000 Genomes Project has identified a DNA segment of approximately 60 kb in the major histocompatibility complex (MHC) between HLA-W and HLA-J that is present in some MHC haplotypes but not others. This DNA segment is largely repeat element-rich but includes the pseudogene HLA-Y, thus pinpointing the location of this pseudogene, and a new HLA class I sequence we have called HLA-OLI. HLA-OLI clusters phylogenetically with the HLA class I pseudogenes, HLA-P and HLA-W, and appears to have a similar genetic structure. The availability of whole-genome sequence data from diverse populations enables a detailed characterization of the MHC at the population level and will have implications for understanding MHC disease associations and the non-HLA MHC factors that impact unrelated hematopoietic cell transplant outcomes.

Evaluation of Density-Based Spatial Clustering for Identifying Genomic Loci Associated with Ischemic Stroke in Genome-Wide Data

The genetic architecture of ischemic stroke (IS), which is one of the leading causes of death worldwide, is complex and underexplored. The traditional approach for associative gene mapping is genome-wide association studies (GWASs), testing individual single-nucleotide polymorphisms (SNPs) across the genomes of case and control groups. The purpose of this research is to develop an alternative approach in which groups of SNPs are examined rather than individual ones. We proposed, validated and applied to real data a new workflow consisting of three key stages: grouping SNPs in clusters, inferring the haplotypes in the clusters and testing haplotypes for the association with phenotype. To group SNPs, we applied the clustering algorithms DBSCAN and HDBSCAN to linkage disequilibrium (LD) matrices, representing pairwise r2 values between all genotyped SNPs. These clustering algorithms have never before been applied to genotype data as part of the workflow of associative studies. In total, 883,908 SNPs and insertion/deletion polymorphisms from people of European ancestry (4929 cases and 652 controls) were processed. The subsequent testing for frequencies of haplotypes restored in the clusters of SNPs revealed dozens of genes associated with IS and suggested the complex role that protocadherin molecules play in IS. The developed workflow was validated with the use of a simulated dataset of similar ancestry and the same sample sizes. The results of classic GWASs are also provided and discussed. The considered clustering algorithms can be applied to genotypic data to identify the genomic loci associated with different qualitative traits, using the workflow presented in this research.

SpecHLA enables full-resolution HLA typing from sequencing data

Reconstructing diploid sequences of human leukocyte antigen (HLA) genes, i.e., full-resolution HLA typing, from sequencing data is challenging. The high homogeneity across HLA genes and the high heterogeneity within HLA alleles complicate the identification of genomic source loci for sequencing reads. Here, we present SpecHLA, which utilizes fine-tuned reads binning and local assembly to achieve accurate full-resolution HLA typing. SpecHLA accepts sequencing data from paired-end, 10×-linked-reads, high-throughput chromosome conformation capture (Hi-C), Pacific Biosciences (PacBio), and Oxford Nanopore Technology (ONT). It can also incorporate pedigree data and genotype frequency to refine typing. In 32 Human Genome Structural Variation Consortium, Phase 2 (HGSVC2) samples, SpecHLA achieved 98.6% accuracy for G-group-resolution HLA typing, inferring entire HLA alleles with an average of three mismatches fewer, ten gaps fewer, and 590 bp less edit distance than HISAT-genotype per allele. Additionally, SpecHLA exhibited a 2-field typing accuracy of 98.6% in 875 real samples. Finally, SpecHLA detected HLA loss of heterozygosity with 99.7% specificity and 96.8% sensitivity in simulated samples of cancer cell lines.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

Complete sequencing of a cynomolgus macaque major histocompatibility complex haplotype

Macaques provide the most widely used nonhuman primate models for studying the immunology and pathogenesis of human diseases. Although the macaque major histocompatibility complex (MHC) region shares most features with the human leukocyte antigen (HLA) region, macaques have an expanded repertoire of MHC class I genes. Although a chimera of two rhesus macaque MHC haplotypes was first published in 2004, the structural diversity of MHC genomic organization in macaques remains poorly understood owing to a lack of adequate genomic reference sequences. We used ultralong Oxford Nanopore and high-accuracy Pacific Biosciences (PacBio) HiFi sequences to fully assemble the ∼5.2-Mb M3 haplotype of an MHC-homozygous, Mauritian-origin cynomolgus macaque (Macaca fascicularis). The MHC homozygosity allowed us to assemble a single MHC haplotype unambiguously and avoid chimeric assemblies that hampered previous efforts to characterize this exceptionally complex genomic region in macaques. The high quality of this new assembly is exemplified by the identification of an extended cluster of six Mafa-AG genes that contains a recent duplication with a highly similar ∼48.5-kb block of sequence. The MHC class II region of this M3 haplotype is similar to the previously sequenced rhesus macaque haplotype and HLA class II haplotypes. The MHC class I region, in contrast, contains 13 MHC-B genes, four MHC-A genes, and three MHC-E genes (vs. 19 MHC-B, two MHC-A, and one MHC-E in the previously sequenced haplotype). These results provide an unambiguously assembled single contiguous cynomolgus macaque MHC haplotype with fully curated gene annotations that will inform infectious disease and transplantation research.

View from the Biological Property: Insight into the Functional Diversity and Complexity of the Gut Mucus

Due to mucin's important protective effect on epithelial tissue, it has garnered extensive attention. The role played by mucus in the digestive tract is undeniable. On the one hand, mucus forms "biofilm" structures that insulate harmful substances from direct contact with epithelial cells. On the other hand, a variety of immune molecules in mucus play a crucial role in the immune regulation of the digestive tract. Due to the enormous number of microorganisms in the gut, the biological properties of mucus and its protective actions are more complicated. Numerous pieces of research have hinted that the aberrant expression of intestinal mucus is closely related to impaired intestinal function. Therefore, this purposeful review aims to provide the highlights of the biological characteristics and functional categorization of mucus synthesis and secretion. In addition, we highlight a variety of the regulatory factors for mucus. Most importantly, we also summarize some of the changes and possible molecular mechanisms of mucus during certain disease processes. All these are beneficial to clinical practice, diagnosis, and treatment and can provide some potential theoretical bases. Admittedly, there are still some deficiencies or contradictory results in the current research on mucus, but none of this diminishes the importance of mucus in protective impacts.

Decoupling peptide binding from T cell receptor recognition with engineered chimeric MHC-I molecules

Major Histocompatibility Complex class I (MHC-I) molecules display self, viral or aberrant epitopic peptides to T cell receptors (TCRs), which employ interactions between complementarity-determining regions with both peptide and MHC-I heavy chain 'framework' residues to recognize specific Human Leucocyte Antigens (HLAs). The highly polymorphic nature of the HLA peptide-binding groove suggests a malleability of interactions within a common structural scaffold. Here, using structural data from peptide:MHC-I and pMHC:TCR structures, we first identify residues important for peptide and/or TCR binding. We then outline a fixed-backbone computational design approach for engineering synthetic molecules that combine peptide binding and TCR recognition surfaces from existing HLA allotypes. X-ray crystallography demonstrates that chimeric molecules bridging divergent HLA alleles can bind selected peptide antigens in a specified backbone conformation. Finally, in vitro tetramer staining and biophysical binding experiments using chimeric pMHC-I molecules presenting established antigens further demonstrate the requirement of TCR recognition on interactions with HLA framework residues, as opposed to interactions with peptide-centric Chimeric Antigen Receptors (CARs). Our results underscore a novel, structure-guided platform for developing synthetic HLA molecules with desired properties as screening probes for peptide-centric interactions with TCRs and other therapeutic modalities.

CRISPR-based targeted haplotype-resolved assembly of a megabase region

Constructing high-quality haplotype-resolved genome assemblies has substantially improved the ability to detect and characterize genetic variants. A targeted approach providing readily access to the rich information from haplotype-resolved genome assemblies will be appealing to groups of basic researchers and medical scientists focused on specific genomic regions. Here, using the 4.5 megabase, notoriously difficult-to-assemble major histocompatibility complex (MHC) region as an example, we demonstrated an approach to construct haplotype-resolved assembly of the targeted genomic region with the CRISPR-based enrichment. Compared to the results from haplotype-resolved genome assembly, our targeted approach achieved comparable completeness and accuracy with reduced computing complexity, sequencing cost, as well as the amount of starting materials. Moreover, using the targeted assembled personal MHC haplotypes as the reference both improves the quantification accuracy for sequencing data and enables allele-specific functional genomics analyses of the MHC region. Given its highly efficient use of resources, our approach can greatly facilitate population genetic studies of targeted regions, and may pave a new way to elucidate the molecular mechanisms in disease etiology.

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.

Quantitative PCR from human genomic DNA: The determination of gene copy numbers for congenital adrenal hyperplasia and RCCX copy number variation

Quantitative PCR (qPCR) is used for the determination of gene copy number (GCN). GCNs contribute to human disorders, and characterize copy number variation (CNV). The single laboratory method validations of duplex qPCR assays with hydrolysis probes on CYP21A1P and CYP21A2 genes, residing a CNV (RCCX CNV) and related to congenital adrenal hyperplasia, were performed using 46 human genomic DNA samples. We also performed the verifications on 5 qPCR assays for the genetic elements of RCCX CNV; C4A, C4B, CNV breakpoint, HERV-K(C4) CNV deletion and insertion alleles. Precision of each qPCR assay was under 1.01 CV%. Accuracy (relative error) ranged from 4.96±4.08% to 9.91±8.93%. Accuracy was not tightly linked to precision, but was significantly correlated with the efficiency of normalization using the RPPH1 internal reference gene (Spearman's ρ: 0.793-0.940, p>0.0001), ambiguity (ρ = 0.671, p = 0.029) and misclassification (ρ = 0.769, p = 0.009). A strong genomic matrix effect was observed, and target-singleplex (one target gene in one assay) qPCR was able to appropriately differentiate 2 GCN from 3 GCN at best. The analysis of all GCNs from the 7 qPCR assays using a multiplex approach increased the resolution of differentiation, and produced 98% of GCNs unambiguously, and all of which were in 100% concordance with GCNs measured by Southern blot, MLPA and aCGH. We conclude that the use of an internal (in one assay with the target gene) reference gene, the use of allele-specific primers or probes, and the multiplex approach (in one assay or different assays) are crucial for GCN determination using qPCR or other methods.

Mobilisation of HLA-F on the surface of bronchial epithelial cells and platelets in asthmatic patients

Uncontrolled inflammation of the airways in chronic obstructive lung diseases leads to exacerbation, accelerated lung dysfunction and respiratory insufficiency. Among these diseases, asthma affects 358 million people worldwide. Human bronchial epithelium cells (HBEC) express both anti-inflammatory and activating molecules, and their deregulated expression contribute to immune cell recruitment and activation, especially platelets (PLT) particularly involved in lung tissue inflammation in asthma context. Previous results supported that HLA-G dysregulation in lung tissue is associated with immune cell activation. We investigated here HLA-F expression, reported to be mobilised on immune cell surface upon activation and displaying its highest affinity for the KIR3DS1-activating NK receptor. We explored HLA-F transcriptional expression in HBEC; HLA-F total expression in PBMC and HBEC collected from healthy individuals at rest and upon chemical activation and HLA-F membrane expression in PBMC, HBEC and PLT collected from healthy individuals at rest and upon chemical activation. We compared HLA-F transcriptional expression in HBEC from healthy individuals and asthmatic patients and its surface expression in HBEC and PLT from healthy individuals and asthmatic patients. Our results support that HLA-F is expressed by HBEC and PLT under healthy physiological conditions and is retained in cytoplasm, barely expressed on the surface, as previously reported in immune cells. In both cell types, HLA-F reaches the surface in the inflammatory asthma context whereas no effect is observed at the transcriptional level. Our study suggests that HLA-F surface expression is a ubiquitous post-transcriptional process in activated cells. It may be of therapeutic interest in controlling lung inflammation.

Molecular Characterization of MHC Class I Alpha 1 and 2 Domains in Asian Seabass (Lates calcarifer)

The Asian seabass is of importance both as a farmed and wild animal. With the emergence of infectious diseases, there is a need to understand and characterize the immune system. In humans, the highly polymorphic MHC class I (MHC-I) molecules play an important role in antigen presentation for the adaptive immune system. In the present study, we characterized a single MHC-I gene in Asian seabass (Lates calcarifer) by amplifying and sequencing the MHC-I alpha 1 and alpha 2 domains, followed by multi-sequence alignment analyses. The results indicated that the Asian seabass MHC-I α1 and α2 domain sequences showed an overall similarity within Asian seabass and retained the majority of the conserved binding residues of human leukocyte antigen-A2 (HLA-A2). Phylogenetic tree analysis revealed that the sequences belonged to the U lineage. Mapping the conserved binding residue positions on human HLA-A2 and grass carp crystal structure showed a high degree of similarity. In conclusion, the availability of MHC-I α1 and α2 sequences enhances the quality of MHC class I genetic information in Asian seabass, providing new tools to analyze fish immune responses to pathogen infections, and will be applicable in the study of the phylogeny and the evolution of antigen-specific receptors.

Sequence Variations Within HLA-G and HLA-F Genomic Segments at the Human Leukocyte Antigen Telomeric End Associated With Acute Graft-Versus-Host Disease in Unrelated Bone Marrow Transplantation

Acute graft-versus-host disease (aGVHD) is defined as a syndrome of an immunological response of graft to the host that occurs early after allogeneic hematopoietic stem cell transplantation (HCT). This disease is frequently observed even in HCT matched for human leukocyte antigen (HLA) alleles at multiple gene loci. Although the HLA region represents complex and diverse genomic characteristics, detailed association analysis is required for the identification of uncharacterized variants that are strongly associated with aGVHD. We genotyped three loci, OR2H2, HLA-F-AS1, and HLA-G, that are located in the 460 kb of HLA telomeric region and statistically analyzed the genotypes including HLA-DPB1 with clinical and transplantation outcomes using 338 unrelated bone marrow transplantation (UR-BMT) patient–donor pairs who were matched for HLA-A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 (HLA-10/10). Multivariate analyses demonstrated that HLA-F-AS1 and HLA-DPB1 mismatches were associated with grade II–IV aGVHD (hazard ratio (HR), 1.76; 95% CI, 1.07–2.88; p = 0.026; and HR, 1.59; CI, 1.02–2.49; p = 0.042, respectively). There was no confounding between HLA-F-AS1 and HLA-DPB1 (p = 0.512), suggesting that the HLA-F-AS1 mismatch has a strong effect on aGVHD independently of HLA-DPB1. Moreover, a stratified analysis suggested possible associations of HLA-F-AS1, HLA-DPB1, and/or HLA-G mismatches with grade II–IV aGVHD and the more severe grade III–IV aGVHD. These findings provide new insights into understanding the molecular mechanism of aGVHD caused by HLA-matched UR-BMT.

The Role and Function of Mucins and Its Relationship to Inflammatory Bowel Disease

Mucus is present throughout the gastrointestinal tract and is essential for regulating gut microbiota homeostasis and preventing disease by protecting the gastrointestinal barrier from microorganisms, pathogens and toxins or other irritants. Mucin (MUC)-2 is a secreted protein produced by epithelial goblet cells as the main component of mucus. Defects in the gastrointestinal tract, such as inflammation and ulcers, cause damage to the mucus barrier, which can worsen mucus quality and reduce mucus production. Therefore, we would like to review the characteristics of MUC2 and its role in intestinal disorders and highlight the importance of further studies. We also investigated whether the role of MUC2 differs between children and adults, ulcerative colitis (UC) and Crohn's disease (CD).

Modulation of Tumor Immune Microenvironment and Prognostic Value of Ferroptosis-Related Genes, and Candidate Target Drugs in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most common type of malignant brain tumor, among which IDH1-wild type GBM has a poor prognosis. Recent studies have shown that ferroptosis-related genes (FRGs) are correlated with the development and progression of cancer. In GBM, the role of FRGs associated with IDH1 status as biological indicators and therapeutic targets remains to be clarified. Ten of FRGs (STEAP3, HSPB1, MAP1LC3A, SOCS1, LOX, CAPG, CP, GDF15, CDKN1A, and CD44) associated with IDH1 status in GBM were identified as key genes through screening by survival analysis and Random Forest using The Cancer Genome Atlas (TCGA) datasets, and the protein expressions of key genes were verified. Transwell and qPCR results showed that ferroptosis promoted the migration of glioblastoma cells and affected the expression of key genes. Our study established the ferroptosis-related prognostic model for GBM patients based on ten key genes by a different modeling method from previous study, the GSVA algorithm. Further, we took the methods of functional enrichment analysis, clinical characteristics, immune cell infiltration, immunomodulator, ESTIMATE and single nucleotide variant (SNV) analysis to study the molecular mechanisms of prognostic model and key genes. The results showed that ten key genes were strongly associated with immune-related factors and were significantly involved in the p53 signaling pathway, senescence and autophagy in cancer, and in the negative regulation of protein kinase activity. Moreover, potential therapeutic drugs were identified by Virtual Screening and Molecular Docking. Our study indicated that the novel ferrotosis-related prognostic model for GBM patients and key genes possessed the prognostic and therapeutic values.

Immunotherapeutic Value of MAP1LC3C and Its Candidate FDA-Approved Drugs Identified by Pan-Cancer Analysis, Virtual Screening and Sensitivity Analysis

Background: The autophagy pathway within the tumour microenvironment can be regulated to inhibit or promote tumour development. In the fight against tumour growth, immunotherapy induces an anti-tumour immune response, whereas autophagy modulates this immune response. A key protein in the autophagy pathway, microtubule-associated protein 1 light chain 3 (MAP1LC3), has recently become a hotspot for tumour research. As a relatively novel member, the function of MAP1LC3C in tumours still need to be investigated. Therefore, the goal of this study was to look into the possible link between MAP1LC3C and immunotherapy for 33 kinds of human malignancies by using pan-cancer analysis. Methods: High-throughput sequencing data from The Cancer Genome Atlas, Genotype-Tissue Expression Project and Cancer Cell Line Encyclopedia databases, combined with clinical data, were used to analyze the expression of MAP1LC3C in 33 types of cancer, as well as patient prognosis and neoplasm staging. Activity scores were calculated using ssGSEA to assess the MAP1LC3C activity in pan-cancer. Associations between MAP1LC3C and the tumour microenvironment, including immune cell infiltration and immunomodulators, were analyzed. Moreover, tumour tissue ImmuneScores and StromalScores were analyzed using the ESTIMATE algorithm. Additionally, associations between MAP1LC3C and tumour mutational burden/microsatellite instability, were investigated. Finally, based on the expression and structure of MAP1LC3C, the United States Food and Drug Administration (FDA)-approved drugs, were screened by virtual screening, molecular docking and NCI-60 drug sensitivity analysis. Results: Our study found that MAP1LC3C was differentially expressed in tumour and normal tissues in 23 of 33 human cancer types, among which MAP1LC3C had prognostic effects in 12 cancer types, and MAP1LC3C expression was significantly correlated with tumour stage in four cancer types. In addition, MAP1LC3C activity in 14 cancer types was consistent with changes in transcription levels. Moreover, MAP1LC3C strongly correlated with immune infiltration, immune modulators and immune markers. Finally, a number of FDA-approved drugs were identified via virtual screening and drug sensitivity analysis. Conclusion: Our study investigated the prognostic and immunotherapeutic value of MAP1LC3C in 33 types of cancer, and several FDA-approved drugs were identified to be highly related to MAP1LC3C and can be potential cancer therapeutic candidates.

Development of an Immortalized Porcine Fibroblast Cell Panel With Different Swine Leukocyte Antigen Genotypes

Immortalized cell lines are valuable resources to expand the molecular characterization of major histocompatibility complex genes and their presented antigens. We generated a panel of immortalized cell lines by transfecting human telomerase reverse transcriptase (hTERT) into primary fibroblast cells prepared from ear, fetal, and lung tissues of 10 pigs from five breeds and successfully cultured them for 30-45 passages. The cell growth characteristic of the immortalized fibroblasts was similar to that of primary fibroblast, which was unable to form colonies on soft agar. The genotypes of major swine leukocyte antigen (SLA) genes, including three classical class I (SLA-1, -2, and -3) and three class II genes (DQB1, DRB1, and DQA), were determined using high-resolution typing. A total of 58 alleles, including a novel allele for SLA-2, were identified. Each cell line was unique. A cell line derived from a National Institutes of Health miniature pig was homozygous across the six major SLA genes. The expression levels of SLA classical class I genes varied among the cell lines and were slightly upregulated in the immortalized compared to the primary cells based on semiquantitative reverse transcription polymerase chain reaction. The immortalized porcine fibroblast cell lines with diverse SLA haplotypes that were developed in this study have potential to be applied in studies regarding the molecular characteristics and genetic structure of SLA genes and epitope-major histocompatibility complex interactions in pigs.

Genetic diversity and disease: The past, present, and future of an old idea

Why do infectious diseases erupt in some host populations and not others? This question has spawned independent fields of research in evolution, ecology, public health, agriculture, and conservation. In the search for environmental and genetic factors that predict variation in parasitism, one hypothesis stands out for its generality and longevity: genetically homogeneous host populations are more likely to experience severe parasitism than genetically diverse populations. In this perspective piece, I draw on overlapping ideas from evolutionary biology, agriculture, and conservation to capture the far-reaching implications of the link between genetic diversity and disease. I first summarize the development of this hypothesis and the results of experimental tests. Given the convincing support for the protective effect of genetic diversity, I then address the following questions: (1) Where has this idea been put to use, in a basic and applied sense, and how can we better use genetic diversity to limit disease spread? (2) What new hypotheses does the established disease-diversity relationship compel us to test? I conclude that monitoring, preserving, and augmenting genetic diversity is one of our most promising evolutionarily informed strategies for buffering wild, domesticated, and human populations against future outbreaks.

Functional analysis of a novel MHC-Iα genotype in orange-spotted grouper: Effects on Singapore grouper iridovirus (SGIV) replication and apoptosis

The classical major histocompatibility complex class I (MHC-Ⅰ) molecule plays a key role in vertebrate immune response for its important functions in antigen presentation and immune regulation. MHC pathway is closely related to many diseases involving autoimmunity, antigen intrusion and inflammation. However, rare literatures about the effect of MHC-I on fish cells apoptosis were reported. In this study, a novel type of MHC-Ⅰα genotype from orange-spotted grouper (named EcMHC-ⅠA*01) were cloned and characterized. It shared a 77% identity to its Epinephelus coioides MHC-Iα homology that has been uploaded to NCBI (ACZ97571.1). Molecular characterization analysis showed that EcMHC-ⅠA*01 encodes a 357-amino-acid protein, containing a signal peptide，α1，α2，α3, Cytoplasmic (Cyt) and Transmembrane (TM) domains. Tissue expression pattern showed that EcMHC-ⅠA*01 was extensively distributed in twelve selected tissues, with higher expression in the gill, intestine and skin. The expression of EcMHC-ⅠA*01 in grouper liver and spleen tissues were significantly induced by different stimuli (Zymosan A, LPS, Ploy I:C, RGNNV and SGIV). Comparing with the EcMHC-ⅠA*01 expression levels induced by Zymosan A, Ploy I:C and RGNNV, the effects induced by SGIV and LPS were more significant. Subcellular localization analysis showed that EcMHC-ⅠA*01 localizes throughout the cytoplasm appeared both diffuse and focal intracellular expression pattern. Overexpression of EcMHC-ⅠA*01 inhibited the CPE progression, the mRNA expression of the SGIV related genes (MCP, LITAF, ICP-18 and VP19) and the protein expression of MCP. Meanwhile, qRT-PCR result showed that EcMHC-ⅠA*01 overexpression upregulated the expression of interferon signaling molecules (IFN-γ, ISG56, MDA5 and MXI) and inflammatory cytokines (IL-1β, IL-6, TNF-α and TRAF6). In addition, our results showed that overexpression of EcMHC-ⅠA*01 promoted the apoptosis of normal fathead minnow (FHM) cells as well as the apoptosis of FHM cells induced by SGIV. However, there was no significant change in the activity of caspase 3 between control group and EcMHC-ⅠA*01 overexpression group, suggesting that EcMHC-ⅠA*01-induced apoptosis may not depend on the caspase 3 pathway. Taken together, these data in our study provide new insights into the role of MHC-I in antiviral immune response and apoptosis in fish.

HLA-H*02:07 Is a Membrane-Bound Ligand of Denisovan Origin That Protects against Lysis by Activated Immune Effectors

The biological relevance of genes initially categorized as “pseudogenes” is slowly emerging, notably in innate immunity. In the HLA region on chromosome 6, HLA-H is one such pseudogene; yet, it is transcribed, and its variation is associated with immune properties. Furthermore, two HLA-H alleles, H*02:07 and H*02:14, putatively encode a complete, membrane-bound HLA protein. Here we thus hypothesized that HLA-H contributes to immune homeostasis similarly to tolerogenic molecules HLA-G, -E, and -F. We tested if HLA-H*02:07 encodes a membrane-bound protein that can inhibit the cytotoxicity of effector cells. We used an HLA-null human erythroblast cell line transduced with HLA-H*02:07 cDNA to demonstrate that HLA-H*02:07 encodes a membrane-bound protein. Additionally, using a cytotoxicity assay, our results support that K562 HLA-H*02:07 inhibits human effector IL-2–activated PBMCs and human IL-2–independent NK92-MI cell line activity. Finally, through in silico genotyping of the Denisovan genome and haplotypic association with Denisovan-derived HLA-A*11, we also show that H*02:07 is of archaic origin. Hence, admixture with archaic humans brought a functional HLA-H allele into modern European and Asian populations.

Haplotypic Associations and Differentiation of MHC Class II Polymorphic Alu Insertions at Five Loci With HLA-DRB1 Alleles in 12 Minority Ethnic Populations in China

The analysis of polymorphic variations in the human major histocompatibility complex (MHC) class II genomic region on the short-arm of chromosome 6 is a scientific enquiry to better understand the diversity in population structure and the effects of evolutionary processes such as recombination, mutation, genetic drift, demographic history, and natural selection. In order to investigate associations between the polymorphisms of HLA-DRB1 gene and recent Alu insertions (POALINs) in the HLA class II region, we genotyped HLA-DRB1 and five Alu loci (AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10), and determined their allele frequencies and haplotypic associations in 12 minority ethnic populations in China. There were 42 different HLA-DRB1 alleles for ethnic Chinese ranging from 12 alleles in the Jinuo to 28 in the Yugur with only DRB1^∗08:03, DRB1^∗09:01, DRB1^∗12:02, DRB1^∗14:01, DRB1^∗15:01, and DRB1^∗15:02 present in all ethnic groups. The POALINs varied in frequency between 0.279 and 0.514 for AluDPB2, 0 and 0.127 for AluDQA2, 0.777 and 0.995 for AluDQA1, 0.1 and 0.455 for AluDRB1 and 0.084 and 0.368 for AluORF10. By comparing the data of the five-loci POALIN in 13 Chinese ethnic populations (including Han-Yunnan published data) against Japanese and Caucasian published data, marked differences were observed between the populations at the allelic or haplotypic levels. Five POALIN loci were in significant linkage disequilibrium with HLA-DRB1 in different populations and AluDQA1 had the highest percentage association with most of the HLA-DRB1 alleles, whereas the nearby AluDRB1 indel was strongly haplotypic for only DRB1^∗01, DRB1^∗10, DRB1^∗15 and DRB1^∗16. There were 30 five-locus POALIN haplotypes inferred in all populations with H5 (no Alu insertions except for AluDQA1) and H21 (only AluDPB2 and AluDQA1 insertions) as the two predominant haplotypes. Neighbor joining trees and principal component analyses of the Alu and HLA-DRB1 polymorphisms showed that genetic diversity of these genomic markers is associated strongly with the population characteristics of language family, migration and sociality. This comparative study of HLA-DRB1 alleles and multilocus, lineage POALIN frequencies of Chinese ethnic populations confirmed that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic and evolutionary markers for the identification of allele and haplotype lineages and genetic variations within the same and/or different populations.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggests that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

HLA Variants and Inhibitor Development in Hemophilia A: A Retrospective Case-Controlled Study Using the ATHNdataset

In hemophilia A (HA) patients, F8 gene-defects as genetic risk-factors for developing inhibitors to Factor VIII have been extensively studied. Here we provide estimates of inhibitor-risk associated with the patient's Human Leukocyte Antigen (HLA). We used next generation sequencing for high-resolution HLA Class II typing of 997 HA patients. Using inhibitor prevalence reports from the My Life Our Future (MLOF) research repository, we calculated Odds Ratios (OR) for inhibitor development in a multivariate model considering HLA-DRB1/3/4/5, HLA-DPB1, HLA-DQB1, race, F8 pathogenic variant type, and age. Participants with 1 HLA variant (DPB1*02:02) had developed inhibitors at a higher rate while participants with 2 HLA variants (DRB1*04:07; DRB1*11:04) had developed inhibitors at a lower rate. Additionally, patients with missense variants had developed inhibitors at a lower rate and participants with large structural changes (>50 bp) had developed inhibitors at a higher rate (both compared to Intron 22 inversion). Using a cohort of participants with a distribution of HLA-DRB1 alleles comparable to that in the North American population we show that the HLA repertoire of a HA patient can be a risk-factor for inhibitor development.

Role for Mucin-5AC in Upper and Lower Airway Pathogenesis in Mice

Mucin-5AC (MUC5AC) is a major secreted mucin in pathogenic airways. To determine its role in mucus-related airway disorders, Muc5ac-deficient (Muc5ac^-/-) and wild-type (Muc5ac^+/+) mice were compared in bleomycin-induced pulmonary fibrosis, respiratory syncytial virus (RSV) disease, and ozone toxicity. Significantly greater inflammation and fibrosis by bleomycin were developed in Muc5ac^-/- lungs compared to Muc5ac^+/+ lungs. More severe mucous cell metaplasia in fibrotic Muc5ac^-/- lungs coincided with bronchial Muc2, Muc4, and Muc5b overexpression. Airway RSV replication was higher in Muc5ac^-/- than in Muc5ac^+/+ during early infection. RSV-caused pulmonary epithelial death, bronchial smooth muscle thickening, and syncytia formation were more severe in Muc5ac^-/- compared to Muc5ac^+/+. Nasal septal damage and subepithelial mucoserous gland enrichment by RSV were greater in Muc5ac^-/- than in Muc5ac^+/+. Ozone exposure developed more severe nasal airway injury accompanying submucosal gland hyperplasia and pulmonary proliferation in Muc5ac^-/- than in Muc5ac^+/+. Ozone caused periodic acid-Schiff-positive secretion only in Muc5ac^-/- nasal airways. Lung E-cadherin level was relatively lower in Muc5ac^-/- than in Muc5ac^+/+ basally and after bleomycin, RSV, and ozone exposure. Results indicate that MUC5AC is an essential mucosal component in acute phase airway injury protection. Subepithelial gland hyperplasia and adaptive increase of other epithelial mucins may compensate airway defense in Muc5ac^-/- mice.

Challenges in the application of NGS in the clinical laboratory

Next-generation sequencing (NGS), also known as massively parallel sequencing, has revolutionized genomic research. The current advances in NGS technology make it possible to provide high resolution, high throughput HLA typing in clinical laboratories. The focus of this review is on the recent development and implementation of NGS in clinical laboratories. Here, we examine the critical role of NGS technologies in clinical immunology for HLA genotyping. Two major NGS platforms (Illumina and Ion Torrent) are characterized including NGS library preparation, data analysis, and validation. Challenges of NGS implementation in the clinical laboratory are also discussed, including sequencing error rate, bioinformatics, result interpretation, analytic sensitivity, as well as large data storage. This review aims to promote the broader applications of NGS technology in clinical laboratories and advocate for the novel applications of NGS to drive future research.

Reduced PCR-generated errors from a hybrid capture-based NGS assay for HLA typing

Next generation sequencing (NGS) assays are state of the art for HLA genotyping. To sequence on an Illumina sequencer, the DNA of interest must be enriched, fragmented, and bookended with known oligonucleotide sequences, a process known as library construction. Many HLA genotyping assays enrich the target loci by long-range PCR (LR-PCR), prior to fragmentation. This PCR step has been reported to introduce errors in the DNA to be sequenced, including inaccurate replication of repeated sequences, and the in vitro recombination of alleles encoded on separate chromosomes. An alternative library construction method involves fragmentation of genomic DNA, followed by hybrid-capture (HC) enrichment of target HLA loci. This HC-based method involves PCR, but with far fewer cycles. Consequently, the HC method had significantly fewer PCR-induced errors, including more faithful replication of repeated sequences, and the near elimination of recombinant sequences. These improvements likely produce more accurate NGS sequencing data of HLA loci.

State-of-the-art genome inference in the human MHC

The Major Histocompatibility Complex (MHC) on the short arm of chromosome 6 is associated with more diseases than any other region of the genome; it encodes the antigen-presenting Human Leukocyte Antigen (HLA) proteins and is one of the key immunogenetic regions of the genome. Accurate genome inference and interpretation of MHC association signals have traditionally been hampered by the region's uniquely complex features, such as high levels of polymorphism; inter-gene sequence homologies; structural variation; and long-range haplotype structures. Recent algorithmic and technological advances have, however, significantly increased the accessibility of genetic variation in the MHC; these developments include (i) accurate SNP-based HLA type imputation; (ii) genome graph approaches for variation-aware genome inference from next-generation sequencing data; (iii) long-read-based diploid de novo assembly of the MHC; (iv) cost-effective targeted MHC sequencing methods. Applied to hundreds of thousands of samples over the last years, these technologies have already enabled significant biological discoveries, for example in the field of autoimmune disease genetics. Remaining challenges concern the development of integrated methods that leverage haplotype-resolved de novo assembly of the MHC for the development of improved MHC genotyping methods for short reads and the construction of improved reference panels for SNP-based imputation. Improved genome inference in the MHC can crucially contribute to an improved genetic and functional understanding of many immune-related phenotypes and diseases.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Mucin signature as a potential tool to predict susceptibility to COVID-19

The Corona Virus Infectious Disease-19 (COVID-19) pandemic has played havoc on both the global health and economy. It is necessary to find a molecular signature that differentiates between low-risk and high-risk individuals. Pathogens, including viruses of the upper respiratory tract, utilize mucin proteins to enter into host cells. Mucins are critical components of innate immunity and also play important roles in infectious disease progression. Their expression is regulated by different cytokines during infection and inflammation. A comparison of mucin signatures between an asymptomatic versus symptomatic and between patients with mild versus severe symptoms could help identify other important proteins involved in the pathology of this new virus. Recent studies on the pathogenicity of the SARS-CoV-2 have found receptors that help its entry into the cells. In this review, we present an overview of how mucins are connected to the pathogenicity of the virus and propose that studying the glycome and mucin signature may lead to the development of a biomarker in predicting the susceptibility, progression, and response to therapy in COVID-19 patients.

Association of HLA polymorphisms and acetaminophen-related Steven-Johnson syndrome with severe ocular complications in Thai population

Background/aims

To investigate the association of genetic polymorphisms of human leucocyte antigens (HLA) class I and II genes with acetaminophen-related Steven-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN) who developed severe ocular complications (SOC) in the Thai population.

Methods

A prospective case–control study including 20 unrelated Thai acetaminophen-related SJS/TEN patients with SOC and 60 Thai healthy volunteers, recruited at three university hospitals in Bangkok, Thailand, from September 2014 to August 2019. HLA genes were analysed using PCR amplification followed by hybridisation with sequence-specific oligonucleotide (SSO) probes with bead-based typing kits. The carrier and gene frequencies of individual HLA alleles in patients were compared with those in control volunteers based on dominant assumption using Fisher’s exact test.

Results

Among HLA class I polymorphisms, HLA-A*33:03, HLA-B*44:03 and HLA-C*07:01 were significantly associated with acetaminophen-related SJS/TEN and SOC with high ORs (95% CI, corrected p value; Pc) in carrier frequency of 5.4 (1.8 to 16.3, Pc=0.0274), 9.0 (95% CI 2.7 to 30.4, Pc=0.0034), and 9.3 (2.8 to 30.2, Pc=0.0022), respectively. There were no significant HLA class II associations with the disease after corrected for a total number of alleles tested.

Conclusion

HLA-B*44:03 was strongly associated with acetaminophen-related SJS/TEN patients who developed SOC in Thai population. In addition, we also found moderate to strong associations with HLA-A*33:03 and HLA-C*07:01 suggesting their potential roles in the pathogenesis of SOC in acetaminophen-related SJS/TEN.

A snapshot of human leukocyte antigen (HLA) diversity using data from the Allele Frequency Net Database

The extensive allelic variability observed in several genes related to the immune response and its significance in different areas including transplantation, disease association studies, diversity in human populations, among many others, has led the scientific community to analyse these variants among individuals. Serving as an electronic data warehouse, the Allele Frequency Net Database (AFND, http://www.allelefrequencies.net) contains data on the frequency of immune related genes and their corresponding alleles from more than 1700 worldwide population samples covering more than ten million unrelated individuals. The collection of population data sets available in AFND encompasses different polymorphic regions including the highly-polymorphic human leukocyte antigen (HLA) system for which more than 1200 populations are available. In this article, we provide an insight of the high diversity found in the HLA region by examining population data sets stored in AFND, as well as a description of the available data sets for further analyses.

Review of Genetic Variation as a Predictive Biomarker for Chronic Graft-Versus-Host-Disease After Allogeneic Stem Cell Transplantation

Chronic graft-versus-host disease (cGvHD) is one of the major complications of allogeneic stem cell transplantation (HSCT). cGvHD is an autoimmune-like disorder affecting multiple organs and involves a dermatological rash, tissue inflammation and fibrosis. The incidence of cGvHD has been reported to be as high as 30% to 60% and there are currently no reliable tools for predicting the occurrence of cGvHD. There is therefore an important unmet clinical need for predictive biomarkers. The present review summarizes the state of the art for genetic variation as a predictive biomarker for cGvHD. We discuss three different modes of action for genetic variation in transplantation: genetic associations, genetic matching, and pharmacogenetics. The results indicate that currently, there are no genetic polymorphisms or genetic tools that can be reliably used as validated biomarkers for predicting cGvHD. A number of recommendations for future studies can be drawn. The majority of studies to date have been under-powered and included too few patients and genetic markers. Like in all complex multifactorial diseases, large collaborative genome-level studies are now needed to achieve reliable and unbiased results. Some of the candidate genes, in particular, CTLA4, HSPE, IL1R1, CCR6, FGFR1OP, and IL10, and some non-HLA variants in the HLA gene region have been replicated to be associated with cGvHD risk in independent studies. These associations should now be confirmed in large well-characterized cohorts with fine mapping. Some patients develop cGvHD despite very extensive immunosuppression and other treatments, indicating that the current therapeutic regimens may not always be effective enough. Hence, more studies on pharmacogenetics are also required. Moreover, all of these studies should be adjusted for diagnostic and clinical features of cGvHD. We conclude that future studies should focus on modern genome-level tools, such as machine learning, polygenic risk scores and genome-wide association study-transcription meta-analyses, instead of focusing on just single variants. The risk of cGvHD may be related to the summary level of immunogenetic differences, or whole genome histocompatibility between each donor-recipient pair. As the number of genome-wide analyses in HSCT is increasing, we are approaching an era where there will be sufficient data to incorporate these approaches in the near future.

Role of Human Leukocyte Antigen System as A Predictive Biomarker for Checkpoint-Based Immunotherapy in Cancer Patients

Recent advances in cancer immunotherapy have clearly shown that checkpoint-based immunotherapy is effective in a small subgroup of cancer patients. However, no effective predictive biomarker has been identified so far. The major histocompatibility complex, better known in humans as human leukocyte antigen (HLA), is a very polymorphic gene complex consisting of more than 200 genes. It has a crucial role in activating an appropriate host immune response against pathogens and tumor cells by discriminating self and non-self peptides. Several lines of evidence have shown that down-regulation of expression of HLA class I antigen derived peptide complexes by cancer cells is a mechanism of tumor immune escape and is often associated to poor prognosis in cancer patients. In addition, it has also been shown that HLA class I and II antigen expression, as well as defects in the antigen processing machinery complex, may predict tumor responses in cancer immunotherapy. Nevertheless, the role of HLA in predicting tumor responses to checkpoint-based immunotherapy is still debated. In this review, firstly, we will describe the structure and function of the HLA system. Secondly, we will summarize the HLA defects and their clinical significance in cancer patients. Thirdly, we will review the potential role of the HLA as a predictive biomarker for checkpoint-based immunotherapy in cancer patients. Lastly, we will discuss the potential strategies that may restore HLA function to implement novel therapeutic strategies in cancer patients.

A diploid assembly-based benchmark for variants in the major histocompatibility complex

Most human genomes are characterized by aligning individual reads to the reference genome, but accurate long reads and linked reads now enable us to construct accurate, phased de novo assemblies. We focus on a medically important, highly variable, 5 million base-pair (bp) region where diploid assembly is particularly useful - the Major Histocompatibility Complex (MHC). Here, we develop a human genome benchmark derived from a diploid assembly for the openly-consented Genome in a Bottle sample HG002. We assemble a single contig for each haplotype, align them to the reference, call phased small and structural variants, and define a small variant benchmark for the MHC, covering 94% of the MHC and 22368 variants smaller than 50 bp, 49% more variants than a mapping-based benchmark. This benchmark reliably identifies errors in mapping-based callsets, and enables performance assessment in regions with much denser, complex variation than regions covered by previous benchmarks.