The MHC, disease and selection

A novel framework for human leukocyte antigen (HLA) genotyping using probe capture-based targeted next-generation sequencing and computational analysis

Human leukocyte antigen (HLA) genes play pivotal roles in numerous immunological applications. Given the immense number of polymorphisms, achieving accurate high-throughput HLA typing remains challenging. This study aimed to harness the human pan-genome reference consortium (HPRC) resources as a potential benchmark for HLA reference materials. We meticulously annotated specific four field-resolution alleles for 11 HLA genes (HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3, -DRB4 and -DRB5) from 44 high-quality HPRC personal genome assemblies. For sequencing, we crafted HLA-specific probes and conducted capture-based targeted sequencing of the genomic DNA of the HPRC cohort, ensuring focused and comprehensive coverage of the HLA region of interest. We used publicly available short-read whole-genome sequencing (WGS) data from identical samples to offer a comparative perspective. To decipher the vast amount of sequencing data, we employed seven distinct software tools: OptiType, HLA-VBseq, HISAT genotype, SpecHLA, T1K, QzType, and DRAGEN. Each tool offers unique capabilities and algorithms for HLA genotyping, allowing comprehensive analysis and validation of the results. We then compared these results with benchmarks derived from personal genome assemblies. Our findings present a comprehensive four-field-resolution HLA allele annotation for 44 HPRC samples. Significantly, our innovative targeted next-generation sequencing (NGS) approach for HLA genes showed superior accuracy compared with conventional short-read WGS. An integrated analysis involving QzType, T1K, and DRAGEN was developed, achieving 100% accuracy for all 11 HLA genes. In conclusion, our study highlighted the combination of targeted short-read sequencing and astute computational analysis as a robust approach for HLA genotyping. Furthermore, the HPRC cohort has emerged as a valuable assembly-based reference in this realm.

Variations in the germinal center response revealed by genetically diverse mouse strains

The humoral response is complex and involves multiple cellular populations and signaling pathways. Bacterial and viral infections, as well as immunization regimens, can trigger this type of response, promoting the formation of microanatomical cellular structures called germinal centers (GCs). GCs formed in secondary lymphoid organs support the differentiation of high-affinity plasma cells and memory B cells. There is growing evidence that the quality of the humoral response is influenced by genetic variants. Using 12 genetically divergent mouse strains, we assessed the impact of genetics on GC cellular traits. At steady state, in the spleen, lymph nodes and Peyer's patches, we quantified GC B cells, plasma cells and follicular helper T cells. These traits were also quantified in the spleen of mice following immunization with a foreign antigen, namely, sheep red blood cells, in addition to the number and size of GCs. We observed both strain- and organ-specific variations in cell type abundance, as well as for GC number and size. Moreover, we find that some of these traits are highly heritable. Importantly, the results of this study inform on the impact of genetic diversity in shaping the GC response and identify the traits that are the most impacted by genetic background.

A synthetic review: natural history of amniote reproductive modes in light of comparative evolutionary genomics

There is a current lack of consensus on whether the ancestral parity mode was oviparity (egg-laying) or viviparity (live-birth) in amniotes and particularly in squamates (snakes, lizards, and amphisbaenids). How transitions between parity modes occur at the genomic level has primary importance for how science conceptualises the origin of amniotes, and highly variable parity modes in Squamata. Synthesising literature from medicine, poultry science, reproductive biology, and evolutionary biology, I review the genomics and physiology of five broad processes (here termed the 'Main Five') expected to change during transitions between parity modes: eggshell formation, embryonic retention, placentation, calcium transport, and maternal-fetal immune dynamics. Throughout, I offer alternative perspectives and testable hypotheses regarding proximate causes of parity mode evolution in amniotes and squamates. If viviparity did evolve early in the history of lepidosaurs, I offer the nucleation site hypothesis as a proximate explanation. The framework of this hypothesis can be extended to amniotes to infer their ancestral state. I also provide a mechanism and hypothesis on how squamates may transition from viviparity to oviparity and make predictions about the directionality of transitions in three species. After considering evidence for differing perspectives on amniote origins, I offer a framework that unifies (i) the extended embryonic retention model and (ii) the traditional model which describes the amniote egg as an adaptation to the terrestrial environment. Additionally, this review contextualises the origin of amniotes and parity mode evolution within Medawar's paradigm. Medawar posited that pregnancy could be supported by immunosuppression, inertness, evasion, or immunological barriers. I demonstrate that this does not support gestation or gravidity across most amniotes but may be an adequate paradigm to explain how the first amniote tolerated internal fertilization and delayed egg deposition. In this context, the eggshell can be thought of as an immunological barrier. If serving as a barrier underpins the origin of the amniote eggshell, there should be evidence that oviparous gravidity can be met with a lack of immunological responses in utero. Rare examples of two species that differentially express very few genes during gravidity, suggestive of an absent immunological reaction to oviparous gravidity, are two skinks Lampropholis guichenoti and Lerista bougainvillii. These species may serve as good models for the original amniote egg. Overall, this review grounds itself in the historical literature while offering a modern perspective on the origin of amniotes. I encourage the scientific community to utilise this review as a resource in evolutionary and comparative genomics studies, embrace the complexity of the system, and thoughtfully consider the frameworks proposed.

Human Leucocyte Antigen Class II Risk and Protective Alleles in Women with Cervical Intraepithelial Neoplasia

Background

Persistent human papillomavirus (HPV) infection is a necessary cause for development of cervical precancerous lesions and cervical cancer, however, only a small percentage of women progress to cervical cancer. The local immune response, determined, among other factors, by Human Leucocyte Antigen (HLA) genes, is thought to be significant. Still the results of genome studies are inconsistent and differ between ethnical populations. The aim of the study was to assess an association between HLA-DQA1*; DQB1*; DRB1* allele's genetic variants between women with cervical precancerous lesions and healthy controls in Latvia.

Materials and methods

From January until April 2017 we enrolled 84 consecutive patients referred for colposcopy to Riga East University Hospital (Latvia) due to abnormal cervical cytology results. 57 women who came for a regular check-up and had normal cytology smears were included in the control group. Material from the cervix was taken for subsequent HLA genotyping of 13 DRB1*, 8 DQA1*, and 12 DQB1* alleles. Colposcopy was performed on all participants. In case of visual suspicion for CIN cervical biopsy was done.

Results

There were 57 "no CIN" patients, 23 histologically proven CIN 1 and 61 CIN2+ cases in the study population. CIN2+ was more often associated with DQA1*0401 (OR 6.68, 95% CI 1.47-30.29, p=0.014), DRB*15 (OR 2.99, 95% CI 1.22-7.39, p=0.017), DQB1*0401 (OR 2.91, 95%CI 1.11-7.68, p=0.03), DQA1*0103 (OR 2.72, 95% CI 1.02-7.21, p=0.045), DRB1*11 (OR 2.42, 95% CI 1.10-5.33, p=0.029) and DQB1*0301 (OR 1.94, 95% CI 1.12-3.38, p=0.018). Women with "no CIN" more often had DQB1*0501 (OR 0.17, 95% CI 0.04-0.81, p=0.026), DRB1*16 (OR 0.21, 95% CI 0.06-0.78, p=0.019), DQA1*0301 (OR 0.35, 95% CI 0.14-0.87, p=0.024) and DRB1*14 (OR 0.59, 95% CI 0.01-0.46, p=0.007).

Conclusions

In the current study we have demonstrated a strong association with risk and protective HLA class II alleles that are determined by the HLA-DRB1*; DQA1*; DQB1*.

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.

Specific T-cell receptor beta-rearrangements of gluten-triggered CD8+ T-cells are enriched in celiac disease patients' duodenal mucosa

Celiac disease (CeD) is an autoimmune disorder affecting the small intestine with gluten as disease trigger. Infections including Influenza A, increase the CeD risk. While gluten-specific CD4+ T-cells, recognizing HLA-DQ2/DQ8 presented gluten-peptides, initiate and sustain the celiac immune response, CD8+ α/β intraepithelial T-cells elicit mucosal damage. Here, we subjected TCRs from a cohort of 56 CeD patients and 22 controls to an analysis employing 749 published CeD-related TCRβ-rearrangements derived from gluten-specific CD4+ T-cells and gluten-triggered peripheral blood CD8+ T-cells. We show, that in addition to TCRs from gluten-specific CD4+ T-cells, TCRs of gluten-triggered CD8+ T-cells are significantly enriched in CeD duodenal tissue samples. TCRβ-rearrangements of gluten-triggered CD8+ T-cells were even more expanded in patients than TCRs from gluten-specific CD4+ T-cells (p < 0.0002) and highest in refractory CeD. Sequence alignments with TCR-antigen databases suggest that a subgroup of these most likely indirectly gluten-triggered TCRs recognize microbial, viral, and autoantigens.

Epstein-Barr virus-acquired immunodeficiency in myalgic encephalomyelitis-Is it present in long COVID?

Both myalgic encephalomyelitis or chronic fatigue syndrome (ME/CFS) and long COVID (LC) are characterized by similar immunological alterations, persistence of chronic viral infection, autoimmunity, chronic inflammatory state, viral reactivation, hypocortisolism, and microclot formation. They also present with similar symptoms such as asthenia, exercise intolerance, sleep disorders, cognitive dysfunction, and neurological and gastrointestinal complaints. In addition, both pathologies present Epstein-Barr virus (EBV) reactivation, indicating the possibility of this virus being the link between both pathologies. Therefore, we propose that latency and recurrent EBV reactivation could generate an acquired immunodeficiency syndrome in three steps: first, an acquired EBV immunodeficiency develops in individuals with "weak" EBV HLA-II haplotypes, which prevents the control of latency I cells. Second, ectopic lymphoid structures with EBV latency form in different tissues (including the CNS), promoting inflammatory responses and further impairment of cell-mediated immunity. Finally, immune exhaustion occurs due to chronic exposure to viral antigens, with consolidation of the disease. In the case of LC, prior to the first step, there is the possibility of previous SARS-CoV-2 infection in individuals with "weak" HLA-II haplotypes against this virus and/or EBV.

Identification of long non-coding RNA MSTRG.5748.1 and MSTRG.7894.1 from Megalobrama amblycephala and their potential roles in innate immunity

Megalobrama amblycephala is one of the most economically important freshwater fish in China, and the bacterial septicemia caused by Aeromonas hydrophila is a serious threat to the breeding industry of M. amblycephala. Unfortunately, the characterization of long noncoding RNA (lncRNA) in response to A. hydrophila infection has not been performed in M. amblycephala. To better understand the biological significance of lncRNA in the immune system, we identified two lncRNA, named MSTRG.5748.1 and MSTRG.7894.1, as playing critical roles in the antibacterial response of M. amblycephala. After separating the nucleus and cytoplasm of the hepatocytes from M. amblycephala, cellular localization of MSTRG.5748.1 and MSTRG.7894.1 was performed to predict their functions. The results showed that MSTRG.5748.1 was mainly expressed in the nucleus, suggesting that its functions are mostly to regulate the expression of downstream genes through epistasis and transcription. MSTRG.7894.1 existed in both the nucleus and cytoplasm, which indicated that it has many regulatory modes. qPCR analysis showed that MSTRG.5748.1 and MSTRG.7894.1 were expressed in the immune-related organs of M. amblycephala, and significantly changed in the liver after A. hydrophila infection. RNA-seq analysis revealed that differentially expressed genes (DEGs) were mainly enriched in antigen processing and presentation via MHC class I, RIG-I-like receptor (RLR) signaling pathway, and IFN-related pathway, and a large number of pathway-related genes were significantly regulated after lncRNA overexpression in muscle cell of M. amblycephala. Overexpression of MSTRG.5748.1 and MSTRG.7894.1 significantly inhibited the expression of STING and IFN, significantly upregulated muscle cell viability, and promoted cell proliferation by targeting STING and IFN.

Multi-functional adaptor SKAP1: regulator of integrin activation, the stop-signal, and the proliferation of T cells

T-cell activation is a complex process involving a network of kinases and downstream molecular scaffolds or adaptors that integrate surface signals with effector functions. One key immune-specific adaptor is Src kinase-associated phosphoprotein 1 (SKAP1), which is also known as src kinase-associated protein of 55 kDa (SKAP55). This mini-review explains how SKAP1 plays multiple roles in regulating integrin activation, the "stop-signal", and the optimization of the cell cycling of proliferating T cells through interactions with various mediators, including the Polo-like kinase 1 (PLK1). Ongoing research on SKAP1 and its binding partners will likely provide important insights into the regulation of immune function and have implications for the development of new treatments for disease states such as cancer and autoimmunity.

The novel HLA-DQB1*03:02:01:14 allele was possibly generated by a recombination event

The novel HLA-DQB1*03:02:01:14 was likely generated by a recombination event between DQB1*03:02:01:01 and DQB1*03:03:02:01.

Insights into the pleiotropic relationships between chronic back pain and inflammation-related musculoskeletal conditions: rheumatoid arthritis and osteoporotic abnormalities

ABSTRACT

The ageing process includes the development of debilitating musculoskeletal (MSK) conditions, including chronic back pain (CBP), rheumatoid arthritis (RA), and osteoporosis (OP). The mechanisms involved in the genetic-epidemiological relationships between these MSK phenotypes are controversial and limited and thus require clarification, in particular, between CBP and the other MSK phenotypes. A cross-sectional statistical analysis was conducted using Europeans from the UK Biobank data collection, including 73,794 CBP, 4883 RA, and 7153 OP cases as well as 242,216 calcaneus bone mineral density scores. C-reactive protein (CRP) was measured for 402,165 subjects in this sample. Genetic correlations were assessed to evaluate shared genetic background between traits. Mendelian randomization was performed to assess a causal relationship between CBP and RA and OP along with other risk factors, such as CRP. Colocalization analysis was conducted to identify shared pleiotropic regions between the examined traits. Bayesian modelling was performed to determine a potential pathway that may explain the interrelationships among these traits. Mendelian randomization analyses revealed that CRP causally predicts CBP only (β = 0.183, 95% CI = 0.077-0.290, P -value = 0.001). Horizontally pleiotropy appeared to explain the relationship between CBP and RA and OP. Through colocalization analysis, several genomic regions emerged describing common genetic influences between CBP and its proposed risk factors, including HLA-DQA1/HLA-DQB1, APOE , SOX5, and MYH7B as well as Histone 1 genes. We speculate that among other factors, CBP and its MSK comorbidities may arise from common inflammatory mechanisms. Colocalized identified genes may aid in advancing or improving the mode of treatment in patients with CBP.

Malaria-driven adaptation of MHC class I in wild bonobo populations

The malaria parasite Plasmodium falciparum causes substantial human mortality, primarily in equatorial Africa. Enriched in affected African populations, the B*53 variant of HLA-B, a cell surface protein that presents peptide antigens to cytotoxic lymphocytes, confers protection against severe malaria. Gorilla, chimpanzee, and bonobo are humans' closest living relatives. These African apes have HLA-B orthologs and are infected by parasites in the same subgenus (Laverania) as P. falciparum, but the consequences of these infections are unclear. Laverania parasites infect bonobos (Pan paniscus) at only one (TL2) of many sites sampled across their range. TL2 spans the Lomami River and has genetically divergent subpopulations of bonobos on each side. Papa-B, the bonobo ortholog of HLA-B, includes variants having a B*53-like (B07) peptide-binding supertype profile. Here we show that B07 Papa-B occur at high frequency in TL2 bonobos and that malaria appears to have independently selected for different B07 alleles in the two subpopulations.

A Genetic Risk Score Distinguishes Different Types of Autoantibody-Mediated Membranous Nephropathy

Introduction

Membranous nephropathy (MN) is the leading cause of nephrotic syndrome in adults and is characterized by detectable autoantibodies against glomerular antigens, most commonly phospholipase A2 receptor 1 (PLA2R1) and thrombospondin type-1 domain containing 7A (THSD7A). In Europeans, genetic variation in at least five loci, PLA2R1, HLA-DRB1, HLA-DQA1, IRF4, and NFKB1, affects the risk of disease. Here, we investigated the genetic risk differences between different autoantibody states.

Methods

1,409 MN individuals were genotyped genome-wide with a dense SNV array. The genetic risk score (GRS) was calculated utilizing the previously identified European MN loci, and results were compared with 4,929 healthy controls and 422 individuals with steroid-sensitive nephrotic syndrome.

Results

GRS was calculated in the 759 MN individuals in whom antibody status was known. The GRS for MN was elevated in the anti-PLA2R1 antibody-positive (N = 372) compared with both the unaffected control (N = 4,929) and anti-THSD7A-positive (N = 31) groups (p < 0.0001 for both comparisons), suggesting that this GRS reflects anti-PLA2R1 MN. Among PLA2R1-positive patients, GRS was inversely correlated with age of disease onset (p = 0.009). Further, the GRS in the dual antibody-negative group (N = 355) was intermediate between controls and the PLA2R1-positive group (p < 0.0001).

Conclusion

We demonstrate that the genetic risk factors for PLA2R1- and THSD7A-antibody-associated MN are different. A higher GRS is associated with younger age of onset of disease. Further, a proportion of antibody-negative MN cases have an elevated GRS similar to PLA2R1-positive disease. This suggests that in some individuals with negative serology the disease is driven by autoimmunity against PLA2R1.

Origin, distribution, and function of three frequent coding polymorphisms in the gene for the human P2X7 ion channel

P2X7, an ion channel gated by extracellular ATP, is widely expressed on the plasma membrane of immune cells and plays important roles in inflammation and apoptosis. Several single nucleotide polymorphisms have been identified in the human P2RX7 gene. In contrast to other members of the P2X family, non-synonymous polymorphisms in P2X7 are common. Three of these occur at overall frequencies of more than 25% and affect residues in the extracellular “head”-domain of P2X7 (155 Y/H), its “lower body” (270 R/H), and its “tail” in the second transmembrane domain (348 T/A). Comparison of the P2X7 orthologues of human and other great apes indicates that the ancestral allele is Y—R—T (at 155–270–348). Interestingly, each single amino acid variant displays lower ATP-sensitivity than the ancestral allele. The originally published reference sequence of human P2X7, often referred to as “wildtype,” differs from the ancestral allele at all three positions, i.e. H—H—A. The 1,000 Genome Project determined the sequences of both alleles of 2,500 human individuals, including roughly 500 persons from each of the five major continental regions. This rich resource shows that the ancestral alleles Y155, R270, and T348 occur in all analyzed human populations, albeit at strikingly different frequencies in various subpopulations (e.g., 25%–59% for Y155, 59%–77% for R270, and 13%–47% for T348). BLAST analyses of ancient human genome sequences uncovered several homozygous carriers of variant P2X7 alleles, possibly reflecting a high degree of inbreeding, e.g., H—R—T for a 50.000 year old Neanderthal, H—R—A for a 24.000 year old Siberian, and Y—R—A for a 7,000 year old mesolithic European. In contrast, most present-day individuals co-express two copies of P2X7 that differ in one or more amino acids at positions 155, 270, and 348. Our results improve the understanding of how P2X7 structure affects its function and suggest the importance of considering P2X7 variants of participants when designing clinical trials targeting P2X7.

EBV protection- and susceptibility-related HLA alleles and EBV status in the Chinese population: A single-center study

BACKGROUND

Although most adults are infected by Epstein-Barr virus (EBV), some patients develop highly lethal diseases associated with EBV infection, including EBV-hemophagocytic lymphohistiocytosis (EBV-HLH), chronic active EBV infections (CAEBV), and lymphoma, the pathogeneses of which remain to be investigated. The human leukocyte antigen (HLA) complex may be associated with the viral infection pathway, and, therefore, HLA alleles may be associated with EBV-related diseases and subpopulations of infected cells, studies related to EBV-associated diseases, and subpopulations of infected cells that were conducted in China are scarce.

METHODS

In this study, we analyzed the high-resolution HLA genotypes of 269 patients with EBV-associated diseases and 213 EBV-seronegative hematopoietic stem cell donors using PCR-SBT assay and elucidated the associations of HLA-A, -B, -C, -DRB1, and -DQB1 alleles with EBV-associated diseases in the Chinese population, Benjamini-Hochberg correction to adjust for multiple testing. HLA genotypes were also analyzed in patients with EBV-associated diseases showing EBV-infected lymphocyte subpopulations.

RESULTS

We found that individuals carrying the following alleles showed the following levels of risks: HLA-DRB1*11 allele, reduced risk of EBV-related disease (OR [odds ratio]: 0.56; 95% confidence interval [95% CI]: 0.32-0.99; p < .05; Adjust p = .71); HLA-DQB1*06:02 allele, reduced risk (OR: 0.5699; 95% CI: 0.3486-0.9317; p < .05; Adjust p = .57); and HLA-B*15:01 allele, increased risk (OR: 1.763; 95% CI: 0.3486-0.9317; p < .05; Adjust p = .57). Patients with EBV-associated diseases showing the B*15:01 genotype had a higher risk of T-cell, NK-cell, and multicell infections than those with other genotype subgroups.

CONCLUSIONS

These findings highlight the importance of HLA genotype for assessing EBV infectivity.

Some Peculiarities of Anthrax Epidemiology in Herbivorous and Carnivorous Animals

Anthrax is an especially dangerous zooanthroponosis caused by the Gram-positive spore-forming bacterium Bacillus anthracis. A notable feature of this disease is the difference in susceptibility to it among different groups of animals. Anthrax primarily affects herbivorous ungulate mammals; they are easily infected, and their disease often leads to rapid, even sudden, death. However, predators and scavengers are extremely resistant to anthrax, and if they become infected, they usually become mildly ill. As the result of the increased sensitivity of ungulates to anthrax and the possibility of disease transmission from them to humans, most studies of anthrax have focused on the diagnosis, prevention, and treatment of infection in farm animals and humans. The issues of anthrax in other animals, such as predators, and the peculiarities of anthrax epidemiology in wild ungulates have not been sufficiently detailed in the literature. In this article, we provide a review of literature sources that describe the differential susceptibility to infection of various groups of animals to anthrax and some epidemiological features of anthrax in animals that are not the main hosts of B. anthracis.

Genomic deciphering of sex determination and unique immune system of a potential model species rare minnow (Gobiocypris rarus)

Gobiocypris rarus is sensitive to environmental pollution, especially to heavy metal and grass carp reovirus (GCRV). Hence, it has potential utility as a biological monitor. Genetic deciphering of its unique immune system will advance our understanding of its unique adaptive strategies, which provide cues for its better application. A de novo genome of rare minnow was obtained, and its sex determination mechanism is ZZ/ZW. We identified several specific mutation genes and specific lost genes of rare minnow, and these might be related to the sensitivity of rare minnow to environmental stimuli. We also analyzed the gene expression level of different organs/tissues and found that several IFIT genes may play key roles in GCRV resistance. In addition, knockout of the gene PCDH10L indicates that PCDH10L affects Pb²⁺-induced mortality in rare minnow. Rare minnow is ready for genetic manipulation and shows potential as an emerging experimental model.

The Immunogenetics of Vitiligo: An Approach Toward Revealing the Secret of Depigmentation

Vitiligo is a hypomelanotic skin disease and considered to be of autoimmune origin due to breaching of immunological self-tolerance, resulting in inappropriate immune responses against melanocytes. The development of vitiligo includes a strong heritable component. Different strategies ranging from linkage studies to genome-wide association studies are used to explore the genetic factors responsible for the disease. Several vitiligo loci containing the respective genes have been identified which contribute to vitiligo and genetic variants for some of the genes are still unknown. These genes include mainly the proteins that play a role in immune regulation and a few other genes important for apoptosis and regulation of melanocyte functions. Despite the available data on genetic variants and risk alleles which influence the biological processes, only few immunological pathways have been found responsible for all ranges of severity and clinical manifestations of vitiligo. However, studies have concluded that vitiligo is of autoimmune origin and manifests due to complex interactions in immune components and their inappropriate response toward melanocytes. The genes involved in the immune regulation and processing the melanocytes antigen and its presentation can serve as effective immune-therapeutics that can target specific immunological pathways involved in vitiligo. This chapter highlights those immune-regulatory genes involved in vitiligo susceptibility and loci identified to date and their implications in vitiligo pathogenesis.

Improved high-throughput MHC typing for non-model species using long-read sequencing

The major histocompatibility complex (MHC) plays a critical role in the vertebrate immune system. Accurate MHC typing is critical to understanding not only host fitness and disease susceptibility, but also the mechanisms underlying host‐pathogen co‐evolution. However, due to the high degree of gene duplication and diversification of MHC genes, it is often technically challenging to accurately characterise MHC genetic diversity in non‐model species. Here we conducted a systematic review to identify common issues associated with current widely used MHC typing approaches. Then to overcome these challenges, we developed a long‐read based MHC typing method along with a new analysis pipeline. Our approach enables the sequencing of fully phased MHC alleles spanning all key functional domains and the separation of highly similar alleles as well as the removal of technical artefacts such as PCR heteroduplexes and chimeras. Using this approach, we performed population‐scale MHC typing in the Tasmanian devil (Sarcophilus harrisii), revealing previously undiscovered MHC functional diversity in this endangered species. Our new method provides a better solution for addressing research questions that require high MHC typing accuracy. Since the method is not limited by species or the number of genes analysed, it will be applicable for studying not only the MHC but also other complex gene families.

Negative trade-off between neoantigen repertoire breadth and the specificity of HLA-I molecules shapes antitumor immunity

Human leukocyte antigen class I (HLA-I) genes shape our immune response against pathogens and cancer. Certain HLA-I variants can bind a wider range of peptides than others, a feature that could be favorable against a range of viral diseases. However, the implications of this phenomenon on cancer immune response are unknown. Here we quantified peptide repertoire breadth (or promiscuity) of a representative set of HLA-I alleles and found that patients with cancer who were carrying HLA-I alleles with high peptide-binding promiscuity have significantly worse prognosis after immune checkpoint inhibition. This can be explained by a reduced capacity of the immune system to discriminate tumor neopeptides from self-peptides when patients carry highly promiscuous HLA-I variants, shifting the regulation of tumor-infiltrating T cells from activation to tolerance. In summary, HLA-I peptide-binding specificity shapes neopeptide immunogenicity and the self-immunopeptidome repertoire in an antagonistic manner, and could underlie a negative trade-off between antitumor immunity and genetic susceptibility to viral infections.

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

"HLA-10-SNP": A new qPCR-based system for rapid, effective, and accurate detection of 10 important SNPs in the HLA region

Mastering the SNP content in the HLA region can be based on it to judiciously select unrelated donor stem cells with preferable MHC matching to lower postoperative complications. Herein, quantitative PCR-based primers and probes were designed for 10 transplant outcome-associated SNP loci with two-allelic polymorphism, and then a new detection system ("HLA-10-SNP") was established. Compared with Sanger sequencing, its accuracy has been proven to reach 100%. Additionally, fluorescent PCR typing of 10 important SNPs via this system expressed excellent repeatability (sensitivity, 20 ng). Overall, the new system achieves single-sample classification precision and easily distinguishable results, equipped with the advantages of simple, rapid, accurate and effective, promising to acquire widespread popularization and application in clinical settings. This article is protected by copyright. All rights reserved.

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.

Unique pathogen peptidomes facilitate pathogen-specific selection and specialization of MHC alleles

A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes—the pool of short peptides derived from the pathogen’s proteome—so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.

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.

Major histocompatibility complex (MHC) associations with diseases in ethnic groups of the Arabian Peninsula

Since the discovery of human leukocyte antigens (HLAs), the function of major histocompatibility complex (MHC) gene families in a wide range of diseases have been the subject of research for decades. In particular, the associations of autoimmune disorders to allelic variants and candidate genes encoding the MHC are well documented. However, despite decades of research, the knowledge of MHC associations with human disease susceptibility have been predominantly studied in European origin, with limited understanding in different populations and ethnic groups. This is particularly evident in countries and ethnic populations of the Arabian Peninsula. Human MHC haplotypes, and its association with diseases, of the variable ethnic groups of this region are poorly studied. This review compiled published manuscripts that have reported a list of autoimmune diseases (insulin-dependent diabetes mellitus, systemic lupus erythematosus, myasthenia gravis, rheumatoid arthritis, psoriasis vulgaris, and multiple sclerosis) associated with MHC class I and class II in the populations of the Arabian Peninsula, specifically Bahrain, Kuwait, Oman, Qatar, Saudi Arabia, the United Arab Emirates, and Yemen. Data available was compared with other three ethnic groups, namely Caucasians, Asians, and Africans. The limited data available in the public domain on the association between MHC gene and autoimmune diseases highlight the challenges in the Middle Eastern region.

Genetic Bias, Diversity Indices, Physiochemical Properties and CDR3 Motifs Divide Auto-Reactive from Allo-Reactive T-Cell Repertoires

The distinct properties of allo-reactive T-cell repertoires are not well understood. To investigate whether auto-reactive and allo-reactive T-cell repertoires encoded distinct properties, we used dextramer enumeration, enrichment, single-cell T-cell receptor (TCR) sequencing and multiparameter analysis. We found auto-reactive and allo-reactive T-cells differed in mean ex vivo frequency which was antigen dependent. Allo-reactive T-cells showed clear differences in TCR architecture, with enriched usage of specific T-cell receptor variable (TRBJ) genes and broader use of T-cell receptor variable joining (TRBJ) genes. Auto-reactive T-cell repertoires exhibited complementary determining regions three (CDR3) lengths using a Gaussian distribution whereas allo-reactive T-cell repertoires exhibited distorted patterns in CDR3 length. CDR3 loops from allo-reactive T-cells showed distinct physical-chemical properties, tending to encode loops that were more acidic in charge. Allo-reactive T-cell repertoires differed in diversity metrics, tending to show increased overall diversity and increased homogeneity between repertoires. Motif analysis of CDR3 loops showed allo-reactive T-cell repertoires differed in motif preference which included broader motif use. Collectively, these data conclude that allo-reactive T-cell repertoires are indeed different to auto-reactive repertoires and provide tangible metrics for further investigations and validation. Given that the antigens studied here are overexpressed on multiple cancers and that allo-reactive TCRs often show increased ligand affinity, this new TCR bank also has translational potential for adoptive cell therapy, soluble TCR-based therapy and rational TCR design.

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.

The adverse events of haematopoietic stem cell transplantation are associated with gene polymorphism within human leukocyte antigen region

Adverse reactions may still occur in some patients after receiving haematopoietic stem cell transplantation (HSCT), even when choosing a human leukocyte antigen (HLA)-matched donor. The adverse reactions of transplantation include disease relapse, graft-versus-host disease (GVHD), mortality and CMV infection. However, only the relapse was discussed in our previous study. Therefore, in this study, we investigated the correlation between the gene polymorphisms within the HLA region and the adverse reactions of post-HSCT in patients with acute leukaemia (n = 176), where 72 patients were diagnosed with acute lymphocytic leukaemia (ALL) and 104 were acute myeloid leukaemia (AML). The candidate single nucleotide polymorphisms were divided into three models: donor, recipient, and donor-recipient pairs and the data of ALL and AML were analysed individually. Based on the results, we found 16 SNPs associated with the survival rates, the risk of CMV infection, or the grade of GVHD in either donor, recipient, or donor-recipient matching models. In the ALL group, the rs209132 of TRIM27 in the donor group was related to CMV infection (p = 0.021), the rs213210 of RING1 in the recipient group was associated with serious GVHD (p = 0.003), and the rs2227956 of HSPA1L in the recipient group correlated with CMV infection (p = 0.001). In the AML group, the rs3130048 of BAG6 in the donor-recipient pairs group was associated with serious GVHD (p = 0.048). Moreover, these SNPs were further associated with the duration time of survival after transplantation. These results could be applied to select the best donor in HSCT.

Dog leucocyte antigen (DLA) class II haplotypes and risk of canine diabetes mellitus in specific dog breeds

Background

Canine diabetes mellitus (DM) is a common endocrine disease in domestic dogs. A number of pathological mechanisms are thought to contribute to the aetiopathogenesis of relative or absolute insulin deficiency, including immune-mediated destruction of pancreatic beta cells. DM risk varies considerably between different dog breeds, suggesting that genetic factors are involved and contribute susceptibility or protection. Associations of particular dog leucocyte antigen (DLA) class II haplotypes with DM have been identified, but investigations to date have only considered all breeds pooled together. The aim of this study was to analyse an expanded data set so as to identify breed-specific diabetes-associated DLA haplotypes.

Methods

The 12 most highly represented breeds in the UK Canine Diabetes Register were selected for study. DLA-typing data from 646 diabetic dogs and 912 breed-matched non-diabetic controls were analysed to enable breed-specific analysis of the DLA. Dogs were genotyped for allelic variation at DLA-DRB1, -DQA1, -DQB1 loci using DNA sequence-based typing. Genotypes from all three loci were combined to reveal three-locus DLA class II haplotypes, which were evaluated for statistical associations with DM. This was performed for each breed individually and for all breeds pooled together.

Results

Five dog breeds were identified as having one or more DLA haplotype associated with DM susceptibility or protection. Four DM-associated haplotypes were identified in the Cocker Spaniel breed, of which one haplotype was shared with Border Terriers. In the three breeds known to be at highest risk of DM included in the study (Samoyed, Tibetan Terrier and Cairn Terrier), no DLA haplotypes were found to be associated with DM.

Conclusions

Novel DLA associations with DM in specific dog breeds provide further evidence that immune response genes contribute susceptibility to this disease in some cases. It is also apparent that DLA may not be contributing obvious or strong risk for DM in some breeds, including the seven breeds analysed for which no associations were identified.

In silico tools for accurate HLA and KIR inference from clinical sequencing data empower immunogenetics on individual-patient and population scales

Immunogenetic variation in humans is important in research, clinical diagnosis and increasingly a target for therapeutic intervention. Two highly polymorphic loci play critical roles, namely the human leukocyte antigen (HLA) system, which is the human version of the major histocompatibility complex (MHC), and the Killer-cell immunoglobulin-like receptors (KIR) that are relevant for responses of natural killer (NK) and some subsets of T cells. Their accurate classification has typically required the use of dedicated biological specimens and a combination of in vitro and in silico efforts. Increased availability of next generation sequencing data has led to the development of ancillary computational solutions. Here, we report an evaluation of recently published algorithms to computationally infer complex immunogenetic variation in the form of HLA alleles and KIR haplotypes from whole-genome or whole-exome sequencing data. For both HLA allele and KIR gene typing, we identified tools that yielded >97% overall accuracy for four-digit HLA types, and >99% overall accuracy for KIR gene presence, suggesting the readiness of in silico solutions for use in clinical and high-throughput research settings.

Next-generation sequencing of 11 HLA loci in a large dengue vaccine cohort from the Philippines

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4-6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls.

Amino Acid Polymorphisms in Hla Class II Differentiate Between Thyroid and Polyglandular Autoimmunity

CONTEXT

The structure of the human leucocyte antigen (HLA) peptide-binding clefts strongly contributes to monoglandular and polyglandular autoimmunity (AP).

OBJECTIVE

To investigate the impact of amino acid polymorphisms on the peptide-binding interactions within HLA class II and its association with AP.

DESIGN

Immunogenetic study.

SETTING

Tertiary referral center for autoimmune endocrine diseases.

SUBJECTS

587 subjects with AP, autoimmune thyroid disease (AITD), type 1 diabetes (T1D), and healthy unrelated controls were typed for HLA class II.

METHODS

Amino acids within the peptide binding cleft that are encoded by HLA class II exon 2 were listed for all codon positions in all subjects. Overall comparisons between disease and control groups with respect to allele distribution at a given locus were performed by assembling rare alleles applying an exact Freeman Halton contingency table test with Monte-Carlo P values based on 150 000 samples.

RESULTS

The Monte Carlo exact Fisher test demonstrated marked differences in all 3 loci, DQA1, DQB1, and DRB1 (P < .0001) between AP and both AITD and controls, as well as between AP type II (Addison's disease as a major endocrine component) and AP type III (T1D + AITD). Differences were also noted between AP and T1D pertaining to the DRB1 allele (P < .041). Seven amino acid positions, DRB1-13, DRB1-26, DRB1-71, DRB1-74, DQA1-47, DQA1-56, and DQB1-57, significantly contributed to AP. Five positions in DQA1 (11, 47, 50, 56, and 69) completely correlated (P < .0001).

CONCLUSION

Amino acid polymorphisms within HLA class II exon 2 mediate the AP risk and differentiate between thyroid and polyglandular autoimmunity.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.

Divergent Allele Advantage at Human MHC Genes: Signatures of Past and Ongoing Selection

The highly polymorphic genes of the major histocompatibility complex (MHC) play a key role in adaptive immunity. Divergent allele advantage, a mechanism of balancing selection, is proposed to contribute to their exceptional polymorphism. It assumes that MHC genotypes with more divergent alleles allow for broader antigen-presentation to immune effector cells, by that increasing immunocompetence. However, the direct correlation between pairwise sequence divergence and the corresponding repertoire of bound peptides has not been studied systematically across different MHC genes. Here, we investigated this relationship for five key classical human MHC genes (human leukocyte antigen; HLA-A, -B, -C, -DRB1, and -DQB1), using allele-specific computational binding prediction to 118,097 peptides derived from a broad range of human pathogens. For all five human MHC genes, the genetic distance between two alleles of a heterozygous genotype was positively correlated with the total number of peptides bound by these two alleles. In accordance with the major antigen-presentation pathway of MHC class I molecules, HLA-B and HLA-C alleles showed particularly strong correlations for peptides derived from intracellular pathogens. Intriguingly, this bias coincides with distinct protein compositions between intra- and extracellular pathogens, possibly suggesting adaptation of MHC I molecules to present specifically intracellular peptides. Eventually, we observed significant positive correlations between an allele’s average divergence and its population frequency. Overall, our results support the divergent allele advantage as a meaningful quantitative mechanism through which pathogen-mediated selection leads to the evolution of MHC diversity.

Inter- and intrabreed diversity of the major histocompatibility complex (MHC) in primitive and draft horse breeds

Background

Polymorphism of major histocompatibility complex (MHC) genes ensures effective immune responses against a wide array of pathogens. However, artificial selection, as performed in the case of domestic animals, may influence MHC diversity. Here, we investigate and compare the MHC diversity of three populations of horses, for which different breeding policies were applied, to evaluate the impact of artificial selection and the environment on MHC polymorphism.

Methods

Samples of DNA were taken from 100 Polish draft horses, 38 stabled Konik Polski horses and 32 semiferal Konik Polski horses. MHC alleles and haplotype diversity within and between these populations of horses was estimated from 11 MHC microsatellite loci.

Results

MHC diversity measured based on allelic richness, observed heterozygosity, expected heterozygosity and polymorphism content was similar across the MHC microsatellite loci in all three populations. The highest expected heterozygosity was detected in semiferal primitive horses (He = 0.74), while the lowest was calculated for draft horses (He = 0.65). In total, 203 haplotypes were determined (111 in Polish draft horses, 43 in semiferal Konik Polski horses and 49 in stabled Konik Polski horses), and four haplotypes were shared between the two populations of Koniks. None of these haplotypes were present in any of the previously investigated horse breeds. Intra-MHC recombination events were detected in all three populations. However, the population of semiferal Konik horses showed the highest recombination frequency among the three horse populations. In addition, three recombination events were detected.

Conclusions

These results showed that despite the different breeding policies, the MHC allele and haplotype diversity was similarly high in all three horse populations. Nevertheless, the proportion of new haplotypes in the offspring was the highest in semiferal Konik Polski horses, which indicates the influence of the environment on MHC diversity in horses. Thus, we speculate that the genetic makeup of the domestic horse MHC might be more strongly influenced by the environment than by artificial selection. Moreover, intra-MHC conversion, insertion, and deletion and intra-MHC recombination may be proposed as mechanisms underlying the generation of new MHC haplotypes in horses.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6

BACKGROUND

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.

RESULTS

Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.

CONCLUSIONS

Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Association between single nucleotide polymorphisms within HLA region and disease relapse for patients with hematopoietic stem cell transplantation

Disease relapse occurs in patients with leukemia even hematopoietic stem cell transplantation (HSCT) was performed with human leukocyte antigen (HLA)-matched donors. As revealed previously by Petersdorf et al., there are nine single nucleotide polymorphisms (SNPs) located in the HLA region that potentially modulate the efficacy of HSCT. In this study, we investigated whether or not the genomic variants 500 base pairs flanking the nine transplantation-related SNPs were related to the risk of post-HSCT relapse for patients with leukemia (n = 141). The genomic DNAs collected from 85 patients with acute myeloid leukemia (AML), 56 patients with acute lymphocytic leukemia (ALL), and their respective HLA-matched donors were subject to SNPs analysis, conferred by the mode of mismatch between donor-recipient pair or by recipient or donor genotype analysis. Seven SNPs were revealed to associate with the risk of relapse post-HSCT. For patients with AML, the increased risk of post-HSCT relapse was associated with the donor SNP of rs111394117 in the intron of NOTCH4 gene, and the recipient SNPs of rs213210 in the ring finger protein 1 (RING1) gene promoter, and rs17220087 and rs17213693 in the intron of HLA-DOB gene. For patients with ALL, the increased risk of post-HSCT relapse was associated with the donor SNP of rs213210 in the RING1 gene promoter, and the recipient SNPs of rs79327197 in the HLA-DOA gene promoter, rs2009658 in the telomeric end of lymphotoxin-alpha (LTA) gene, rs17220087 and rs17213693 in the intron of HLA-DOB gene, and rs2070120 in the 3′-UTR of HLA-DOB gene. This study sheds new insight into selecting better candidate donors for performing HSCT in patients with AML and ALL.

Thinking About the Evolution of Complex Traits in the Era of Genome-Wide Association Studies

Many traits of interest are highly heritable and genetically complex, meaning that much of the variation they exhibit arises from differences at numerous loci in the genome. Complex traits and their evolution have been studied for more than a century, but only in the last decade have genome-wide association studies (GWASs) in humans begun to reveal their genetic basis. Here, we bring these threads of research together to ask how findings from GWASs can further our understanding of the processes that give rise to heritable variation in complex traits and of the genetic basis of complex trait evolution in response to changing selection pressures (i.e., of polygenic adaptation). Conversely, we ask how evolutionary thinking helps us to interpret findings from GWASs and informs related efforts of practical importance.

Using next-generation sequencing for characterising HLA-DRB1 and DQB1 loci in a cohort of Colombian women

The Colombian population is characterised by a high genetic diversity, secondary to the ethnic mixture arising from colonisation. Unfortunately, few reports are available regarding HLA-DRB1 and DQB1 diversity in Colombia to date. HLA-DRB1 and DQB1 diversity was identified in this study using next-generating sequencing (NGS) on a cohort of Colombian women. Cervical samples taken from 276 women were used for typing DRB1 and DQB1 loci by Illumina MiSeq. Allele and haplotype frequencies were calculated using an expectation-maximisation algorithm. Hardy-Weinberg Equilibrium and linkage disequilibrium (LD) between loci were evaluated. Forty-seven DRB1 alleles and 14 DQB1 alleles were identified. DRB1*04:07:01G and DQB1*03:02:01G alleles occurred most frequently in the target population. Significant LD was found in 44 out of the 144 identified haplotypes, within which DRB1*04:07:01G-DQB1*03:02:01G occurred most frequently (6.56%). The alleles and haplotypes found with NGS agreed with that found in previous reports involving lower resolution for the Colombian population, and greater genetic variability was found, especially concerning DRB1. Comparing allele and haplotype frequency distribution in the target population to that of other populations denoted HLA system intra- and inter-population diversity.

Epstein-Barr Virus and the Human Leukocyte Antigen Complex

Purpose

While most adults are infected Epstein-Barr virus (EBV), 3-5% remain uninfected. The human leukocyte antigen (HLA) complex, which controls many pathogens, may influence infection and disease associated with EBV.

Recent Findings

Numerous EBV proteins and miRNAs down-regulate HLA class I and II expression on the cell surface. HLA class II functions as a receptor for EBV entry into B cells. Specific HLA class II alleles correlate with the susceptibility of B cells to EBV infection in vitro and with EBV seropositivity or seronegativity of humans. HLA class I polymorphisms correlate with development and severity of EBV infectious mononucleosis and with the risk of several virus-associated malignancies including nasopharyngeal carcinoma, Hodgkin lymphoma, and post-transplant lymphoproliferative disease.

Significance

These findings indicate that while EBV has evolved to use MHC class II as a receptor for virus entry, polymorphisms in MHC class II and class I influence virus infection and disease.

Genetic Identification of Two Novel Loci Associated with Steroid-Sensitive Nephrotic Syndrome

BACKGROUND

Steroid-sensitive nephrotic syndrome (SSNS), the most common form of nephrotic syndrome in childhood, is considered an autoimmune disease with an established classic HLA association. However, the precise etiology of the disease is unclear. In other autoimmune diseases, the identification of loci outside the classic HLA region by genome-wide association studies (GWAS) has provided critical insights into disease pathogenesis. Previously conducted GWAS of SSNS have not identified non-HLA loci achieving genome-wide significance.

METHODS

In an attempt to identify additional loci associated with SSNS, we conducted a GWAS of a large cohort of European ancestry comprising 422 ethnically homogeneous pediatric patients and 5642 ethnically matched controls.

RESULTS

The GWAS found three loci that achieved genome-wide significance, which explain approximately 14% of the genetic risk for SSNS. It confirmed the previously reported association with the HLA-DR/DQ region (lead single-nucleotide polymorphism [SNP] rs9273542, P=1.59×10^-43; odds ratio [OR], 3.39; 95% confidence interval [95% CI], 2.86 to 4.03) and identified two additional loci outside the HLA region on chromosomes 4q13.3 and 6q22.1. The latter contains the calcium homeostasis modulator family member 6 gene CALHM6 (previously called FAM26F). CALHM6 is implicated in immune response modulation; the lead SNP (rs2637678, P=1.27×10^-17; OR, 0.51; 95% CI, 0.44 to 0.60) exhibits strong expression quantitative trait loci effects, the risk allele being associated with lower lymphocytic expression of CALHM6.

CONCLUSIONS

Because CALHM6 is implicated in regulating the immune response to infection, this may provide an explanation for the typical triggering of SSNS onset by infections. Our results suggest that a genetically conferred risk of immune dysregulation may be a key component in the pathogenesis of SSNS.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

An exome wide association study of pulmonary tuberculosis patients and their asymptomatic household contacts

Tuberculosis is a leading cause of death in India. To identify genetic variants associated with susceptibility or resistance to Mycobacterium tuberculosis infection, we have performed an exome-wide association study with 0.2 million exonic variants among 119 pairs of tuberculosis patients and their clinically asymptomatic household contacts. The strongest association was identified for rs61104666[A], a synonymous variant (p.E292E) of exon 5 of the gene SIGLEC15 (OR = 2.4, p = 1.49 × 10^-5). We also found association of non-coding variants in the 3'UTR region of a gene encoding the class II human leukocyte antigens (HLAs), HLA-DRA. rs13209234[A] (minor allele frequency (MAF) = 13.8%) (OR = 0.35, P = 2.5 × 10^-4) and rs3177928[A] (minor allele frequency (MAF) = 13.7%) (OR = 0.35, P = 3.3 × 10^-4) were associated with protection from tuberculosis. These two SNPs, rs13209234 and rs3177928, are in complete linkage disequilibrium. These associations remained valid when additional data on freshly recruited individuals were jointly analyzed on 250 patient-control pairs. The identified gene, HLA-DRA, suggest involvement of immune regulation, indicating pathways associated with antigen presentation in tuberculosis infection.

Functional variation at an expressed MHC class IIβ locus associates with Ranavirus infection intensity in larval anuran populations

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIβ locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIβ alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIβ functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Whereas specialist major histocompatibility complex (MHC) molecules initiate immune response against only relatively few pathogens, generalists provide protection against a broad range. Accordingly, this study shows that the geographical distribution of generalist MHC alleles in human populations reflects exposure to diverse infectious diseases.

Variation in the human genome influences our susceptibility to infectious diseases, but the causal link between disease and underlying mutation often remains enigmatic. Major histocompatibility complex II (MHC class II) molecules shape both our immune response against pathogens and our tolerance of self-peptides. The genomic region that encodes MHC molecules is of particular interest, as it is home to more genetic disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Here, we propose that MHC class II molecules can be categorized into two major types; specialists initiate effective immune response against only relatively few pathogens, while generalists provide protection against a broad range of pathogens. As support, we demonstrate that generalist MHC class II variants are more prevalent in human populations residing in pathogen-rich areas.