Advances in the Evolutionary Understanding of MHC Polymorphism

Expression and Immune Characterization of Major Histocompatibility Complex in Paralichthys olivaceus after Antigen Stimulation

The Major histocompatibility complex (Mhc) is an important molecule for antigen presenting and binds to T cell receptors, activating T lymphocytes and triggering specific immune responses. To investigate the role of MhcII in adaptive immunity, in this study, mhcIIα and mhcIIβ of flounder (Paralichthys olivaceus) were cloned, polyclonal antibodies (Abs) against their extracellular regions were produced, respectively, and their distribution on cells and tissues and expression patterns, which varied by antigen stimulation or pathogen infection, were investigated. The results showed that the open reading frame (ORF) of mhcIIα is 708 bp, including 235 amino acids (aa); and the ORF of mhcIIβ is 741 bp, encoding 246aa. The mhcIIα and mhcIIβ were significantly expressed in gills, spleen, and peripheral blood leukocytes (PBLs). Their antibodies could specifically recognize eukaryotic expressed MhcIIα and MhcIIβ. MhcIIα+ and MhcIIβ+ cells were 30.2 ± 2.9% of the percentage in peripheral blood leukocytes. MhcII molecules were co-localized with CD83 and IgM on leukocytes, respectively, but not on CD4+ or CD8+ T lymphocyte subpopulations. The expression of both mhcIIα and mhcIIβ were significantly upregulated in flounder after bacteria and virus challenges. The percentages of MhcII+ cells, MhcII+/CD83+, and MhcII+/IgM+ double-positive cells increased significantly after PHA and ConA stimulation, respectively; they varied significantly in PBLs after polyI:C stimulation, and no variations were found after LPS treatment. In the meantime, variations in MhcII+ cells were consistent with that of CD4+ T lymphocytes. These results suggest that MhcII, mainly expressed in B cells and dendritic cells, play an essential role in antigen presentation, and respond significantly to exogenous antigens and T cell-dependent antigens. These results may provide an important reference for the study of cellular immunity in teleosts.

Bw4 ligand and direct T-cell receptor binding induced selection on HLA A and B alleles

Introduction

The HLA region is the hallmark of balancing selection, argued to be driven by the pressure to present a wide variety of viral epitopes. As such selection on the peptide-binding positions has been proposed to drive HLA population genetics. MHC molecules also directly binds to the T-Cell Receptor and killer cell immunoglobulin-like receptors (KIR).

Methods

We here combine the HLA allele frequencies in over six-million Hematopoietic Stem Cells (HSC) donors with a novel machine-learning-based method to predict allele frequency.

Results

We show for the first time that allele frequency can be predicted from their sequences. This prediction yields a natural measure for selection. The strongest selection is affecting KIR binding regions, followed by the peptide-binding cleft. The selection from the direct interaction with the KIR and TCR is centered on positively charged residues (mainly Arginine), and some positions in the peptide-binding cleft are not associated with the allele frequency, especially Tyrosine residues.

Discussion

These results suggest that the balancing selection for peptide presentation is combined with a positive selection for KIR and TCR binding.

Human HLA prolongs the host inflammatory response in Streptococcus suis serotype 2 infection compared to mouse H2 molecules

Streptococcus suis (S. suis) is widely acknowledged as a significant zoonotic pathogen in Southeast Asia and China, which has led to a substantial number of fatalities in both swine and humans. Despite the prevalent use of mice as the primary animal model to study S. suis pathogenesis, the substantial differences in the major histocompatibility complex (MHC) between humans and mice underscore the ongoing exploration for a more suitable and effective animal model. In this study, humanized transgenic HLA-A11/DR1 genotypes mice were used to evaluate the differences between humanized HLA and murine H2 in S. suis infection. Following intravenous administration of S. suis suspensions, we investigated bacterial load, cytokine profiles, pathological alterations, and immune cell recruitment in both Wild-type (WT) and humanized mice across different post-infection time points. Relative to WT mice, humanized mice exhibited heightened pro-inflammatory cytokines, exacerbated tissue damage, increased granulocyte recruitment with impaired resolution, notably more pronounced during the late infection stage. Additionally, our examination of bacterial clearance rates suggests that HLA-A11/DR1 primarily influences cell recruitment and mitochondrial reactive oxygen species (ROS) production, which affects the bacterial killing capacity of macrophages in the late stage of infection. The reduced IL-10 production and lower levels of regulatory T cells in humanized mice could underlie their compromised resolution ability. Intervention with IL-10 promotes bacterial clearance and inflammatory regression in the late stages of infection in transgenic mice. Our findings underscore the heightened sensitivity of HLA-A11/DR1 mice with impaired resolution to S. suis infection, effectively mirroring the immune response seen in humans during infection. The humanized HLA-A11/DR1 mice could serve as an optimal animal model for investigating the pathogenic and therapeutic mechanisms associated with sepsis and other infectious diseases.

Assessing micro-macroparasite selective pressures and anthropogenic disturbance as drivers of immune gene diversity in a Neotropical wild cat

Anthropogenic environmental change is reducing available habitat for wild species, providing novel selection pressures such as infectious diseases and causing species to interact in new ways. The potential for emerging infectious diseases and zoonoses at the interface between humans, domestic animals, and wild species is a key global concern. In vertebrates, diversity at the major histocompatibility complex MHC is critical to disease resilience, and its study in wild populations provides insights into eco-evolutionary dynamics that human activities alter. In natural populations, variation at MHC loci is partly maintained by balancing selection, driven by pathogenic selective pressures. We hypothesize that MHC genetic diversity differs between guigna populations inhabiting human-dominated landscapes (higher pathogen pressures) versus more natural habitats (lower pathogen pressures). We predict that MHC diversity in guignas would be highest in human-dominated landscapes compared with continuous forest habitats. We also expected to find higher MHC diversity in guignas infected with micro and macro parasites (higher parasite load) versus non infected guignas. We characterized for the first time the genetic diversity at three MHC class I and II exons in 128 wild guignas (Leopardus guigna) across their distribution range in Chile (32-46° S) and Argentina, representing landscapes with varying levels of human disturbance. We integrated MHC sequence diversity with multiple measures of anthropogenic disturbance and both micro and macro parasite infection data. We also assessed signatures of positive selection acting on MHC genes. We found significantly higher MHC class I diversity in guignas inhabiting landscapes where houses were present, and with lower percentage of vegetation cover, and also in animals with more severe cardiorespiratory helminth infection (richness and intensity) and micro-macroparasite co-infection. This comprehensive, landscape-level assessment further enhances our knowledge on the evolutionary dynamics and adaptive potential of vertebrates in the face of emerging infectious disease threats and increasing anthropogenic impacts.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

MHC class Ib-restricted CD8+ T cells possess strong tumoricidal activities

The importance of classical CD8+ T cells in tumor eradication is well acknowledged. However, the anti-tumor activity of MHC (major histocompatibility complex) Ib-restricted CD8+ T (Ib-CD8+ T) cells remains obscure. Here, we show that CX3CR1-expressing Ib-CD8+ T cells (Ib-restricted CD8+ T cells) highly express cytotoxic factors, austerely resist exhaustion, and effectively eliminate various tumors. These Ib-CD8+ T cells can be primed by MHC Ia (MHC class Ia molecules) expressed on various cell types for optimal activation in a Tbet-dependent manner. Importantly, MHC Ia does not allogeneically activate Ib-CD8+ T cells, rather, sensitizes these cells for T cell receptor activation. Such effects were observed when MHC Ia+ cells were administered to tumor-bearing Kb-/-Db-/-mice. A similar population of tumoricidal CX3CR1+CD8+ T cells was identified in wild-type mice and melanoma patients. Adoptive transfer of Ib-CD8+ T cells to wild-type mice inhibited tumor progression without damaging normal tissues. Taken together, we demonstrate that MHC class Ia can prime Ib-CD8+ T cells for robust tumoricidal activities.

An integrated approach to the characterization of immune repertoires using AIMS: An Automated Immune Molecule Separator

The adaptive immune system employs an array of receptors designed to respond with high specificity to pathogens or molecular aberrations faced by the host organism. Binding of these receptors to molecular fragments-collectively referred to as antigens-initiates immune responses. These antigenic targets are recognized in their native state on the surfaces of pathogens by antibodies, whereas T cell receptors (TCR) recognize processed antigens as short peptides, presented on major histocompatibility complex (MHC) molecules. Recent research has led to a wealth of immune repertoire data that are key to interrogating the nature of these molecular interactions. However, existing tools for the analysis of these large datasets typically focus on molecular sets of a single type, forcing researchers to separately analyze strongly coupled sequences of interacting molecules. Here, we introduce a software package for the integrated analysis of immune repertoire data, capable of identifying distinct biophysical differences in isolated TCR, MHC, peptide, antibody, and antigen sequence data. This integrated analytical approach allows for direct comparisons across immune repertoire subsets and provides a starting point for the identification of key interaction hotspots in complementary receptor-antigen pairs. The software (AIMS-Automated Immune Molecule Separator) is freely available as an open access package in GUI or command-line form.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

Non-MHC immunity genes do not affect parasite load in European invasive populations of common raccoon

Understanding the evolutionary mechanisms behind invasion success enables predicting which alien species and populations are the most predisposed to become invasive. Parasites may mediate the success of biological invasions through their effect on host fitness. The evolution of increased competitive ability (EICA) hypothesis assumes that escape from parasites during the invasion process allows introduced species to decrease investment in immunity and allocate resources to dispersal and reproduction. Consequently, the selective pressure of parasites on host species in the invasive range should be relaxed. We used the case of the raccoon Procyon lotor invasion in Europe to investigate the effect of gastrointestinal pathogen pressure on non-MHC immune genetic diversity of newly established invasive populations. Despite distinct differences in parasite prevalence between analysed populations, we detected only marginal associations between two analysed SNPs and infection intensity. We argue that the differences in parasite prevalence are better explained by detected earlier associations with specific MHC-DRB alleles. While the escape from native parasites seems to allow decreased investment in overall immunity, which relaxes selective pressure imposed on immune genes, a wide range of MHC variants maintained in the invasive range may protect from newly encountered parasites.

Evolutionary trade-off between innate and acquired immune defences in birds

BACKGROUND

The development, maintenance, and use of immune defences are costly. Therefore, animals face trade-offs in terms of resource allocation within their immune system and between their immune system and other physiological processes. To maximize fitness, evolution may favour investment in one immunological defence or subsystem over another in a way that matches a species broader life history strategy. Here, we used phylogenetically-informed comparative analyses to test for relationships between two immunological components. Natural antibodies and complement were used as proxies for the innate branch; structural complexity of the major histocompatibility complex (MHC) region was used for the acquired branch.

RESULTS

We found a negative association between the levels of natural antibodies (i.e., haemagglutination titre) and the total MHC gene copy number across the avian phylogeny, both at the species and family level. The family-level analysis indicated that this association was apparent for both MHC-I and MHC-II, when copy numbers within these two MHC regions were analysed separately. The association remained significant after controlling for basic life history components and for ecological traits commonly linked to pathogen exposure.

CONCLUSION

Our results provide the first phylogenetically robust evidence for an evolutionary trade-off within the avian immune system, with a more developed acquired immune system (i.e., more complex MHC architecture) in more derived bird lineages (e.g., passerines) being accompanied by an apparent downregulation of the innate immune system.

Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90β1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIβ1 sequences. The G90β1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

Tutorial: a statistical genetics guide to identifying HLA alleles driving complex disease

The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.

The quest for good genes: Epigamic traits, fitness, MHC and multilocus heterozygosity in the guppy

The 'good genes' hypothesis for the evolution of male secondary sexual traits poses that female preferences for such traits are driven by indirect genetic benefits. However, support for the hypothesis remains ambiguous, and, in particular, the genetic basis for the benefits has rarely been investigated. Here, we use seminatural populations of Trinidadian guppies to investigate whether sexually selected traits (orange, black and iridescent colouration, gonopodium length and body size) predict fitness measured as the number of grandoffspring, a metric that integrates across fitness components and sexes. Furthermore, we tested whether two potential sources of genetic benefits-major histocompatibility complex (MHC) genotypes and multilocus heterozygosity (MLH)-are significant predictors of fitness and of the size of sexually selected traits. We found a significant, nonlinear effect of the area of black pigmentation and male body size on the number of grandoffspring, suggesting stabilizing selection on black area, and nonlinear selection favouring small body size. MLH was heritable (h2  = 0.14) and significantly predicted the number of grandoffspring, indicating the potential for genetic benefits based on heterozygosity. We also found support for local heterozygosity effects, which may reflect a noneven distribution of genetic load across the genome. MHC genotype was not significantly associated with any tested fitness component, or with the load of Gyrodactylus parasites. Neither MHC nor MLH was significant predictor of sexually selected traits. Overall, our results highlight the role of heterozygosity in determining fitness, but do not provide support for male sexually selected traits being indicators of genetic quality.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Adaptive genetic variations are key for understanding evolutionary processes influencing invasions. However, we have limited knowledge on how adaptive genetic diversity in invasive species responds to new pathogenic environments. Here, we compared variations in immune major histocompatibility complex (MHC) class-II β gene and neutral loci in relation to pathogenic chytrid fungus (Batrachochytrium dendrobatidis, Bd) infection across invasive and native populations of American bullfrog between China and United States (US). Chinese invasive populations show a 60% reduction in neutral cytb variations relative to US native populations, and there were similar MHC variation and functional diversity between them. One MHC allele private to China was under recent positive selection and associated with decreased Bd infection, partly explaining the lower Bd prevalence for Chinese populations than for native US populations. These results suggest that pathogen-mediated selection favors adaptive MHC variations and functional diversity maintenance against serious bottlenecks during the early invasions (within 15 generations) of bullfrogs.

Blue light promotes vitamin C-mediated ferroptosis of melanoma through specifically upregulating transporter SVCT2 and generating Fe2

Vitamin C (VC)-based cancer therapy is a promising therapeutic approach for a variety of cancers due to its profound effects on redox reactions and metabolic pathways. However, high administration dosage of VC for necessary therapeutic efficacy for cancers increases the risk of overt side effects and limits its clinical use. Here, we show cutaneous blue light irradiation can specifically upregulate the sodium-dependent vitamin C transporter 2 (SVCT2) of the tumor and increase effectively the VC concentration at the tumor sites by an overall low dosage administration. In the mouse melanoma model, blue light stimulates the SVCT2 expression through the nuclear factor-kappa B (NF-κB) signaling pathway both in vitro and in vivo. The increased cellular VC together with Fe2+ generated by blue light simultaneously elevate cellular oxidative stress and trigger the ferroptosis of melanoma. With this revealed mechanism, the synergistic actions of blue light on the VC transporter and Fe2+ generation lead to a ca. 20-fold reduction in the administration dosage of VC with an effective melanoma elimination and prolonged survival. The work defines the killing mechanism of blue light on VC-based cancer therapy and provides a practical approach for promoting VC uptake. This light-assisted VC therapy is not only highly efficient for melanoma but also considerable for a broad clinical utility.

Associations between HLA class II alleles and IgE sensitization to allergens in the Qatar Biobank cohort

Background

Allergic disorders are the consequence of IgE sensitization to allergens. Population studies have shown that certain human leukocyte antigen (HLA) alleles are associated with increased or decreased risk of developing allergy.

Objective

We aimed to characterize the relationship between HLA class II allelic diversity and IgE sensitization in an understudied Arab population.

Methods

We explored associations between IgE sensitization to 7 allergen mixes and mesquite (comprising 41 food or aeroallergens) and 45 common classical HLA class II alleles in a well-defined cohort of 797 individuals representing the general adult population of Qatari nationals and long-term residents. To do so, we performed HLA calling from whole genome sequencing data at 2-field resolution using 2 independent algorithms. We then applied 3 different regression models to assess either each allergen mix independently, in the context of IgE sensitization to other allergens tested, or polysensitization.

Results

More than half (n = 447) of the study participants showed IgE sensitization to at least 1 allergen, most of them (n = 400) to aeroallergens (Phadiatop). We identified statistically significant negative and positive associations with 24 HLA class II alleles. These have been reported to confer risk or protection from variety of diseases; however, only a few have previously been associated with allergy in other populations.

Conclusions

Our study reveals several new risk and protective genetic markers for allergen-specific IgE sensitization. This is a first and essential step toward a better understanding of the origins of allergic diseases in this understudied population.

Neutral Forces and Balancing Selection Interplay to Shape the Major Histocompatibility Complex Spatial Patterns in the Striped Hamster in Inner Mongolia: Suggestive of Broad-Scale Local Adaptation

BACKGROUND

The major histocompatibility complex (MHC) plays a key role in the adaptive immune response to pathogens due to its extraordinary polymorphism. However, the spatial patterns of MHC variation in the striped hamster remain unclear, particularly regarding the relative contribution of the balancing selection in shaping MHC spatial variation and diversity compared to neutral forces.

METHODS

In this study, we investigated the immunogenic variation of the striped hamster in four wild populations in Inner Mongolia which experience a heterogeneous parasitic burden. Our goal was to identify local adaptation by comparing the genetic structure at the MHC with that at seven microsatellite loci, taking into account neutral processes.

RESULTS

We observed significant variation in parasite pressure among sites, with parasite burden showing a correlation with temperature and precipitation. Molecular analysis revealed a similar co-structure between MHC and microsatellite loci. We observed lower genetic differentiation at MHC loci compared to microsatellite loci, and no correlation was found between the two.

CONCLUSIONS

Overall, these results suggest a complex interplay between neutral evolutionary forces and balancing selection in shaping the spatial patterns of MHC variation. Local adaptation was not detected on a small scale but may be applicable on a larger scale.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Antigen presentation in vertebrates: Structural and functional aspects

Antigen presentation is a key process of the immune system and is responsible for the activation of T cells. The main characters are the major histocompatibility complex class I (MHC-I) and class II (MHC-II) molecules, and accessory proteins that act as chaperones for these glycoproteins. Current knowledge of this process and also the elucidation of the structural features of these proteins, has been extensively reviewed in humans. Unfortunately, this is not the case for non-human species, wherein the function and structural characteristic of the antigen presentation proteins is far from being understood. The majority of previous studies in non-human species, especially in teleost fish and lower vertebrates, are limited to the transcriptomic level, which leads to gaps in the knowledge about the functional process of antigen presentation in these species. This review summarizes what is known so far about antigen presentation pathways in vertebrates from a structural and functional perspective. The focus is not only on the MHC receptors, but also, on the forgotten characters of these pathways such as the proteins of the peptide loading complex, and the MHC-II chaperone invariant chain.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

Balancing selection on the complement system of a wild rodent

BACKGROUND

Selection pressure exerted by pathogens can influence patterns of genetic diversity in the host. In the immune system especially, numerous genes encode proteins involved in antagonistic interactions with pathogens, paving the way for coevolution that results in increased genetic diversity as a consequence of balancing selection. The complement system is a key component of innate immunity. Many complement proteins interact directly with pathogens, either by recognising pathogen molecules for complement activation, or by serving as targets of pathogen immune evasion mechanisms. Complement genes can therefore be expected to be important targets of pathogen-mediated balancing selection, but analyses of such selection on this part of the immune system have been limited.

RESULTS

Using a population sample of whole-genome resequencing data from wild bank voles (n = 31), we estimated the extent of genetic diversity and tested for signatures of balancing selection in multiple complement genes (n = 44). Complement genes showed higher values of standardised β (a statistic expected to be high under balancing selection) than the genome-wide average of protein coding genes. One complement gene, FCNA, a pattern recognition molecule that interacts directly with pathogens, was found to have a signature of balancing selection, as indicated by the Hudson-Kreitman-Aguadé test (HKA) test. Scans for localised signatures of balancing selection in this gene indicated that the target of balancing selection was found in exonic regions involved in ligand binding.

CONCLUSION

The present study adds to the growing evidence that balancing selection may be an important evolutionary force on components of the innate immune system. The identified target in the complement system typifies the expectation that balancing selection acts on genes encoding proteins involved in direct interactions with pathogens.

The acute inflammatory response of teleost fish

Acute inflammation is crucial to the immune responses of fish. The process protects the host from infection and is central to induction of subsequent tissue repair programs. Activation of proinflammatory signals reshapes the microenvironment within an injury/infection site, initiates leukocyte recruitment, promotes antimicrobial mechanisms and contributes to the resolution of inflammation. Inflammatory cytokines and lipid mediators are primary contributors to these processes. Uncontrolled or persistent induction results in delayed tissue healing. The kinetics by which inducers and regulators of acute inflammation exert their actions is essential for understanding the pathogenesis of fish diseases and identifying potential treatments. Although, a number of these are well-conserved across, others are not, reflecting the unique physiologies and life histories of members of this unique animal group.

The biogenesis of the immunopeptidome

The immunopeptidome is the repertoire of peptides bound and presented by the MHC class I, class II, and non-classical molecules. The peptides are produced by the degradation of most cellular proteins, and in some cases, peptides are produced from extracellular proteins taken up by the cells. This review attempts to first describe some of its known and well-accepted concepts, and next, raise some questions about a few of the established dogmas in this field: The production of novel peptides by splicing is questioned, suggesting here that spliced peptides are extremely rare, if existent at all. The degree of the contribution to the immunopeptidome by degradation of cellular protein by the proteasome is doubted, therefore this review attempts to explain why it is likely that this contribution to the immunopeptidome is possibly overstated. The contribution of defective ribosome products (DRiPs) and non-canonical peptides to the immunopeptidome is noted and methods are suggested to quantify them. In addition, the common misconception that the MHC class II peptidome is mostly derived from extracellular proteins is noted, and corrected. It is stressed that the confirmation of sequence assignments of non-canonical and spliced peptides should rely on targeted mass spectrometry using spiking-in of heavy isotope-labeled peptides. Finally, the new methodologies and modern instrumentation currently available for high throughput kinetics and quantitative immunopeptidomics are described. These advanced methods open up new possibilities for utilizing the big data generated and taking a fresh look at the established dogmas and reevaluating them critically.

Crystal Structure of a Classical MHC Class I Molecule in Dogs; Comparison of DLA-88*0 and DLA-88*5 Category Molecules

DLA-88 is a classical major histocompatibility complex (MHC) class I gene in dogs, and allelic DLA-88 molecules have been divided into two categories named "DLA-88*0" and "DLA-88*5." The defining difference between the two categories concerns an LQW motif in the α2 domain helical region of the DLA-88*5 molecules that includes the insertion of an extra amino acid compared to MHC class I consensus length. We here show that this motif has been exchanged by recombination between different DLA-88 evolutionary lineages. Previously, with pDLA-88*508:01, the structure of a molecule of the DLA-88*5 category was elucidated. The present study is the first to elucidate a structure, using X-ray crystallography, of the DLA-88*0 category, namely DLA-88*001:04 complexed with β2m and a nonamer peptide derived from canine distemper virus (CDV). The LQW motif that distinguishes DLA-88*5 from DLA-88*0 causes a shallower peptide binding groove (PBG) and a leucine exposed at the top of the α2 domain helix expected to affect T cell selection. Peptide ligand amino acid substitution and pMHC-I complex formation and stability analyses revealed that P2 and P3 are the major anchor residue positions for binding to DLA-88*001:04. We speculate that the distribution pattern of the LQW motif among canine classical MHC class I alleles represents a strategy to enhance allogeneic rejection by T cells of transmissible cancers such as canine transmissible venereal tumor (CTVT).

Immune genotypes, immune responses, and survival in a wild bird population

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The F_ST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

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

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

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

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.

Cytokine gene polymorphism and parasite susceptibility in free-living rodents: Importance of non-coding variants

Associations between genetic variants and susceptibility to infections have long been studied in free-living hosts so as to infer the contemporary evolutionary forces that shape the genetic polymorphisms of immunity genes. Despite extensive studies of proteins interacting with pathogen-derived ligands, such as MHC (major histocompatilbility complex) or TLR (Toll-like receptors), little is known about the efferent arm of the immune system. Cytokines are signalling molecules that trigger and modulate the immune response, acting as a crucial link between innate and adaptive immunity. In the present study we investigated how genetic variation in cytokines in bank voles Myodes glareolus affects their susceptibility to infection by parasites (nematodes: Aspiculuris tianjensis, Heligmosomum mixtum, Heligmosomoides glareoli) and microparasites (Cryptosporidium sp, Babesia microti, Bartonella sp.). We focused on three cytokines: tumour necrosis factor (TNF), lymphotoxin alpha (LTα), and interferon beta (IFNβ1). Overall, we identified four single nucleotide polymorphisms (SNPs) associated with susceptibility to nematodes: two located in LTα and two in IFNβ1. One of those variants was synonymous, another located in an intron. Each SNP associated with parasite load was located in or next to a codon under selection, three codons displayed signatures of positive selection, and one of purifying selection. Our results indicate that cytokines are prone to parasite-driven selection and that non-coding variants, although commonly disregarded in studies of the genetic background of host-parasite co-evolution, may play a role in susceptibility to infections in wild systems.

Complex patterns shape immune genes diversity during invasion of common raccoon in Europe - Selection in action despite genetic drift

Rapid adaptation is common in invasive populations and is crucial to their long-term success. The primary target of selection in the invasive species' new range is standing genetic variation. Therefore, genetic drift and natural selection acting on existing variation are key evolutionary processes through which invaders will evolve over a short timescale. In this study, we used the case of the raccoon Procyon lotor invasion in Europe to identify the forces shaping the diversity of immune genes during invasion. The genes involved in the defence against infection should be under intense selection pressure in the invasive range where novel pathogens are expected to occur. To disentangle the selective and demographic processes shaping the adaptive immune diversity of its invasive and expanding populations, we have developed species-specific single-nucleotide polymorphism markers located in the coding regions of targeted immune-related genes. We characterised the genetic diversity of 110 functionally important immune genes in two invasive and one native raccoon genetic clusters, each presenting a different demographic history. Despite the strong effect of demographic processes in the invasive clusters, we detected a subset of genes exhibiting the diversity pattern suggestive of selection. The most likely process shaping the variation in those genes was balancing selection. The selected genes belong to toll-like receptors and cytokine-related genes. Our results suggest that the prevalence of selection depends on the level of diversity, that is - less genetically diverse invasive population from the Czech Republic displayed fewer signs of selection. Our results highlight the role of standing genetic variation in adapting to new environment. Understanding the evolutionary mechanisms behind invasion success would enable predicting how populations may respond to environmental change.

Assessing the causes and consequences of gut mycobiome variation in a wild population of the Seychelles warbler

BACKGROUND

Considerable research has focussed on the importance of bacterial communities within the vertebrate gut microbiome (GM). However, studies investigating the significance of other microbial kingdoms, such as fungi, are notably lacking, despite their potential to influence host processes. Here, we characterise the fungal GM of individuals living in a natural population of Seychelles warblers (Acrocephalus sechellensis). We evaluate the extent to which fungal GM structure is shaped by environment and host factors, including genome-wide heterozygosity and variation at key immune genes (major histocompatibility complex (MHC) and Toll-like receptor (TLR)). Importantly, we also explore the relationship between fungal GM differences and subsequent host survival. To our knowledge, this is the first time that the genetic drivers and fitness consequences of fungal GM variation have been characterised for a wild vertebrate population.

RESULTS

Environmental factors, including season and territory quality, explain the largest proportion of variance in the fungal GM. In contrast, neither host age, sex, genome-wide heterozygosity, nor TLR3 genotype was associated with fungal GM differences in Seychelles warblers. However, the presence of four MHC-I alleles and one MHC-II allele was associated with changes in fungal GM alpha diversity. Changes in fungal richness ranged from between 1 and 10 sequencing variants lost or gained; in some cases, this accounted for 20% of the fungal variants carried by an individual. In addition to this, overall MHC-I allelic diversity was associated with small, but potentially important, changes in fungal GM composition. This is evidenced by the fact that fungal GM composition differed between individuals that survived or died within 7 months of being sampled.

CONCLUSIONS

Our results suggest that environmental factors play a primary role in shaping the fungal GM, but that components of the host immune system-specifically the MHC-may also contribute to the variation in fungal communities across individuals within wild populations. Furthermore, variation in the fungal GM can be associated with differential survival in the wild. Further work is needed to establish the causality of such relationships and, thus, the extent to which components of the GM may impact host evolution. Video Abstract.

Identification of comutation in signaling pathways to predict the clinical outcomes of immunotherapy

BACKGROUND

Immune checkpoint blockades (ICBs) have emerged as a promising treatment for cancer. Recently, tumour mutational burden (TMB) and neoantigen load (NAL) have been proposed to be potential biomarkers to predict the efficacy of ICB; however, they were limited by difficulties in defining the cut-off values and inconsistent detection platforms. Therefore, it is critical to identify more effective predictive biomarkers for screening patients who will potentially benefit from immunotherapy. In this study, we aimed to identify comutated signaling pathways to predict the clinical outcomes of immunotherapy.

METHODS

Here, we comprehensively analysed the signaling pathway mutation status of 9763 samples across 33 different cancer types from The Cancer Genome Atlas (TCGA) by mapping the somatic mutations to the pathways. We then explored the comutated pathways that were associated with increased TMB and NAL by using receiver operating characteristic (ROC) curve analysis and multiple linear regressions.

RESULTS

Our results revealed that comutation of the Spliceosome (Sp) pathway and Hedgehog (He) signaling pathway (defined as SpHe-comut⁺) could be used as a predictor of increased TMB and NAL and was associated with increased levels of immune-related signatures. In seven independent immunotherapy cohorts, we validated that SpHe-comut⁺ patients exhibited a longer overall survival (OS) or progression-free survival (PFS) and a higher objective response rate (ORR) than SpHe-comut^- patients. Moreover, a combination of SpHe-comut status with PD-L1 expression further improved the predictive value for ICB therapy.

CONCLUSION

Overall, SpHe-comut⁺ was demonstrated to be an effective predictor of immunotherapeutic benefit in seven independent immunotherapy cohorts and may serve as a potential and convenient biomarker for the clinical application of ICB therapy.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

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

Background

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

Methods

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

Results

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

Conclusion

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

Extending the Mass Spectrometry-Detectable Landscape of MHC Peptides by Use of Restricted Access Material

Mass spectrometry-based immunopeptidomics enables the comprehensive identification of major histocompatibility complex (MHC) peptides from a cell culture as well as from tissue or tumor samples and is applied for the identification of tumor-specific and viral T-cell epitopes. Although mass spectrometry is generally considered an "unbiased" method for MHC peptide identification, the physicochemical properties of MHC peptides can greatly influence their detectability. Here, we demonstrate that highly hydrophobic peptides are lost during sample preparation when C18 solid-phase extraction (SPE) is used for separating MHC peptides from proteins. To overcome this limitation, we established an optimized protocol involving restricted access material (RAM). Compared to C18-SPE, RAM-SPE improved the overall MHC peptide recovery and extended the landscape of mass spectrometry-detectable MHC peptides toward more hydrophobic peptides.

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.