Generalists and Specialists: A New View of How MHC Class I Molecules Fight Infectious Pathogens

A personal view on developmental and comparative immunology: What, how and why?

What are the future directions of the fields of developmental and comparative immunology? In thinking through this question as I write, I find myself marvelling at the very long ways that we have come since I began as a PhD student some 50 years ago. I think that we cannot know what technical and theoretical advances will emerge in the future, nor will our initial aims survive the realities of what appears in our sights, often from unexpected directions. I feel that we should not allow what we already know about some well-studied systems to blind us to the wide range of possibilities, and that remaining a humble seeker helps the uptake of new realities. Finally, it would be good to try answering the whole range of questions about developmental and comparative immunology, from what to how to why.

Early human migration determines the risk of being attacked by wolves: ethnic gene diversity on the development of systemic lupus erythematosus

The prevalence of systemic lupus erythematosus (SLE) varies significantly based on ethnicity rather than geographic distribution; thus, the prevalence is higher in Asian, Hispanic, and Black African populations than in European populations. The risk of developing lupus nephritis (LN) is the highest among Asian populations. Therefore, we hypothesize that human genetic diversity between races has occurred through the early human migration and human genetic adaptation to various environments, with a particular focus on pathogens. Additionally, we compile the currently available evidence on the ethnic gene diversity of SLE and how it relates to disease severity. The human leukocyte antigen (HLA) locus is well established as associated with susceptibility to SLE; specific allele distributions have been observed across diverse populations. Notably, specific amino acid residues within these HLA loci demonstrate significant associations with SLE risk. The non-HLA genetic loci associated with SLE risk also varies across diverse ancestries, implicating distinct immunological pathways, such as the type-I interferon and janus kinase-signal transducers and activators of transcription (JAK-STAT) pathways in Asians, the type-II interferon signaling pathway in Europeans, and B cell activation pathway in Africans. Furthermore, assessing individual genetic susceptibility using genetic risk scores (GRS) for SLE helps to reveal the diverse prevalence, age of onset, and clinical phenotypes across different ethnicities. A higher GRS increases the risk of LN and the severity of SLE. Therefore, understanding ethnic gene diversity is crucial for elucidating disease mechanisms and SLE severity, which could enable the development of novel drugs specific to each race.

Polymorphism in F pocket affects peptide selection and stability of type 1 diabetes-associated HLA-B39 allotypes

HLA-B*39:06, HLA-B*39:01, and HLA-B*38:01 are closely related HLA allotypes differentially associated with type 1 diabetes (T1D) risk and progression. B*39:06 is highly predisposing, while B*39:01 and B*38:01 are weakly predisposing and protective allotypes, respectively. Here, we aimed to decipher molecular mechanisms underlying the differential association of these allotypes with T1D pathogenesis. We addressed peptide binding and conformational stability of HLA-B allotypes using computational and experimental approaches. Computationally, we found that B*39:06 and B*39:01 allotypes had more rigid peptide-binding grooves and were more promiscuous in binding peptides than B*38:01. Peptidomes of B*39:06 and B*39:01 contained fewer strong binders and were of lower affinity than that of B*38:01. Experimentally, we demonstrated that B*39:06 and B*39:01 had a higher capacity to bind peptides and exit to the cell surface but lower surface levels and were degraded faster than B*38:01. In summary, we propose that promiscuous B*39:06 and B*39:01 may bind suboptimal peptides and transport them the cell surface, where such unstable complexes may contribute to the pathogenesis of T1D.

Large scale transcriptional analysis of MHC class I haplotype diversity in sheep

Domestic sheep (Ovis aries) have been an important component of livestock agricultural production for thousands of years. Preserving genetic diversity within livestock populations maintains a capacity to respond to changing environments and rapidly evolving pathogens. MHC genetic diversity can influence immune functionality at individual and population levels. Here, we focus on defining functional MHC class I haplotype diversity in a large cohort of Scottish Blackface sheep pre-selected for high levels of MHC class II DRB1 diversity. Using high-throughput amplicon sequencing with three independent sets of barcoded primers we identified 134 MHC class I transcripts within 38 haplotypes. Haplotypes were identified with between two and six MHC class I genes, plus variable numbers of conserved sequences with very low read frequencies. One or two highly transcribed transcripts dominate each haplotype indicative of two highly polymorphic, classical MHC class I genes. Additional clusters of medium, low, and very low expressed transcripts are described, indicative of lower transcribed classical, non-classical and genes whose function remains to be determined.

MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs

West African Mastomys rodents are the primary reservoir of the zoonotic Lassa virus (LASV). The virus causes haemorrhagic Lassa fever and considerable mortality in humans. To date, the role of Mastomys immunogenetics in resistance to, and persistence of, LASV infections is largely unknown. Here, we investigated the role of Major Histocompatibility Complex class I (MHC-I) on LASV infection status (i.e., active vs. cleared infection, determined via PCR and an immunofluorescence assay on IgG antibodies, respectively) in Mastomys natalensis and M. erythroleucus sampled within southwestern Nigeria. We identified more than 190 and 90 MHC-I alleles by Illumina high throughput-sequencing in M. natalensis and M. erythroleucus, respectively, with different MHC allele compositions and frequencies between LASV endemic and non-endemic sites. In M. natalensis, the MHC allele ManaMHC-I*006 was negatively associated with active infections (PCR-positive) and positively associated with cleared infections (IgG-positive) simultaneously, suggesting efficient immune responses that facilitate LASV clearance in animals carrying this allele. Contrarily, alleles ManaMHC-I*008 and ManaMHC-I*021 in M. natalensis, and MaerMHC-I*008 in M. erythroleucus, were positively associated with active infection, implying susceptibility. Alleles associated with susceptibility shared a glutamic acid at the positively selected codon 57, while ManaMHC-I*006 featured an arginine. There was no link between number of MHC alleles per Mastomys individual and LASV prevalence. Thus, specific alleles, but not MHC diversity per se, seem to mediate antibody responses to viremia. We conclude that co-evolution with LASV likely shaped the MHC-I diversity of the main LASV reservoirs in southwestern Nigeria, and that information on reservoir immunogenetics may hold insights into transmission dynamics and zoonotic spillover risks.

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.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

HLA variants and TCR diversity against SARS-CoV-2 in the pre-COVID-19 era

HLA antigen presentation and T-cell mediated immunity are critical to control acute viral infection such as COVID-19 caused by SARS-CoV-2. Recent data suggest that both the depth of peptide presentation and the breadth of the T-cell repertoire are associated with disease outcome. It has also been shown that unexposed subjects can develop strong T-cell responses against SARS-CoV-2 due to heterologous immunity. In this study, we explored the anti-SARS-CoV-2 T-cell repertoire by analyzing previously published T-cell receptor (TCR) CDR3β immunosequencing data in a cohort of 116 healthy donors and in the context of immune reconstitution after allogeneic hematopoietic stem cell transplantation in 116 recipients collected during the pre-COVID-19 era. For this, 143,310 publicly available SARS-CoV-2 specific T-cell sequences were investigated among the 3.5 million clonotypes in the cohort. We also performed HLA class I peptide binding predictions using the reference proteome of the virus and high resolution genotyping data in these patients. We could demonstrate that individuals are fully equipped at the genetic level to recognize SARS-CoV-2. This is evidenced by the 5% median cumulative frequency of clonotypes having their sequence matched to a SARS-CoV-2 specific T-cell. In addition, any combination of HLA class I variants in this cohort is associated with a broad capacity of presenting hundreds of SARS-CoV-2 derived peptides. These results could be explained by heterologous immunity and random somatic TCR recombination. We speculate that these observations could explain the efficacy of the specific immune response against SARS-CoV-2 in individuals without risk factors of immunodeficiency and infected prior to vaccination.

The two-stage molecular scenery of SARS-CoV-2 infection with implications to disease severity: An in-silico quest

Introduction

The two-stage molecular profile of the progression of SARS-CoV-2 (SCOV2) infection is explored in terms of five key biological/clinical questions: (a) does SCOV2 exhibits a two-stage infection profile? (b) SARS-CoV-1 (SCOV1) vs. SCOV2: do they differ? (c) does and how SCOV2 differs from Influenza/INFL infection? (d) does low viral-load and (e) does COVID-19 early host response relate to the two-stage SCOV2 infection profile? We provide positive answers to the above questions by analyzing the time-series gene-expression profiles of preserved cell-lines infected with SCOV1/2 or, the gene-expression profiles of infected individuals with different viral-loads levels and different host-response phenotypes.

Methods

Our analytical methodology follows an in-silico quest organized around an elaborate multi-step analysis pipeline including: (a) utilization of fifteen gene-expression datasets from NCBI's gene expression omnibus/GEO repository; (b) thorough designation of SCOV1/2 and INFL progression stages and COVID-19 phenotypes; (c) identification of differentially expressed genes (DEGs) and enriched biological processes and pathways that contrast and differentiate between different infection stages and phenotypes; (d) employment of a graph-based clustering process for the induction of coherent groups of networked genes as the representative core molecular fingerprints that characterize the different SCOV2 progression stages and the different COVID-19 phenotypes. In addition, relying on a sensibly selected set of induced fingerprint genes and following a Machine Learning approach, we devised and assessed the performance of different classifier models for the differentiation of acute respiratory illness/ARI caused by SCOV2 or other infections (diagnostic classifiers), as well as for the prediction of COVID-19 disease severity (prognostic classifiers), with quite encouraging results.

Results

The central finding of our experiments demonstrates the down-regulation of type-I interferon genes (IFN-1), interferon induced genes (ISGs) and fundamental innate immune and defense biological processes and molecular pathways during the early SCOV2 infection stages, with the inverse to hold during the later ones. It is highlighted that upregulation of these genes and pathways early after infection may prove beneficial in preventing subsequent uncontrolled hyperinflammatory and potentially lethal events.

Discussion

The basic aim of our study was to utilize in an intuitive, efficient and productive way the most relevant and state-of-the-art bioinformatics methods to reveal the core molecular mechanisms which govern the progression of SCOV2 infection and the different COVID-19 phenotypes.

Discovery of novel MHC-B haplotypes in Chantecler chickens

The major histocompatibility complex (MHC) is a highly polymorphic cluster of genes which contribute to immune response. Located on chromosome 16, the chicken MHC has great influence over disease resistance and susceptibility. Through the use of a high-density SNP panel which encompasses the MHC-B region, haplotypes can be easily identified. This study aims to use an MHC-B SNP panel to evaluate the MHC-B variability in the Chantecler breed. This breed is native to Quebec, Canada, and is a dual-purpose breed known for its strong resistance to extreme cold temperatures. The Chantecler breed faced a near extinction event in the 1970s, which most likely resulted in a genetic bottleneck and loss of diversity. Despite this, SNP haplotype diversity was observed among 4 Chantecler populations. A total of 8 haplotypes were observed. Of these haplotypes, 6 were previously defined in other breeds, and the other 2 were unique to the Chantecler. Within the populations, the number of haplotypes ranged from 4 to 7, with 3 haplotypes, including the novel BSNP-Chant01, being present in all the groups. This study shows existence of reasonable diversity in the MHC-B region of the Chantecler breed and our results further contribute to understanding the variability of this region in chickens.

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.

Visualising tapasin- and TAPBPR-assisted editing of major histocompatibility complex class-I immunopeptidomes

Which peptides are selected for presentation by major histocompatibility complex class-I (MHC-I) molecules is a key determinant of successful immune responses. Peptide selection is co-ordinated by the tapasin and TAP Binding PRotein (TAPBPR) proteins, which ensure MHC-I molecules preferentially acquire high-affinity-binding peptides. New structural analyses have offered insight into how tapasin achieves this function within the peptide-loading complex (PLC) (comprising the Transporter associated with Antigen Presentation (TAP) peptide transporter, tapasin-ERp57, MHC-I and calreticulin), and how TAPBPR performs a peptide editing function independently of other molecules. The new structures reveal nuances in how tapasin and TAPBPR interact with MHC-I, and how calreticulin and ERp57 complement tapasin to exploit the plasticity of MHC-I molecules to achieve peptide editing.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

Conserved Evolution of MHC Supertypes among Japanese Frogs Suggests Selection for Bd Resistance

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is a major threat to amphibians, yet there are no reports of major disease impacts in East Asian frogs. Genetic variation of the major histocompatibility complex (MHC) has been associated with resistance to Bd in frogs from East Asia and worldwide. Using transcriptomic data collated from 11 Japanese frog species (one individual per species), we isolated MHC class I and IIb sequences and validated using molecular cloning. We then compared MHC from Japanese frogs and other species worldwide, with varying Bd susceptibility. Supertyping analysis, which groups MHC alleles based on physicochemical properties of peptide binding sites, identified that all examined East Asian frogs contained at least one MHC-IIb allele belonging to supertype ST-1. This indicates that, despite the large divergence times between some Japanese frogs (up to 145 million years), particular functional properties in the peptide binding sites of MHC-II are conserved among East Asian frogs. Furthermore, preliminary analysis using NetMHCIIpan-4.0, which predicts potential Bd-peptide binding ability, suggests that MHC-IIb ST-1 and ST-2 have higher overall peptide binding ability than other supertypes, irrespective of whether the peptides are derived from Bd, other fungi, or bacteria. Our findings suggest that MHC-IIb among East Asian frogs may have co-evolved under the same selective pressure. Given that Bd originated in this region, it may be a major driver of MHC evolution in East Asian frogs.

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.

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.

A Wider and Deeper Peptide-Binding Groove for the Class I Molecules from B15 Compared with B19 Chickens Correlates with Relative Resistance to Marek's Disease

The chicken MHC is known to confer decisive resistance or susceptibility to various economically important pathogens, including the iconic oncogenic herpesvirus that causes Marek's disease (MD). Only one classical class I gene, BF2, is expressed at a high level in chickens, so it was relatively easy to discern a hierarchy from well-expressed thermostable fastidious specialist alleles to promiscuous generalist alleles that are less stable and expressed less on the cell surface. The class I molecule BF2*1901 is better expressed and more thermostable than the closely related BF2*1501, but the peptide motif was not simpler as expected. In this study, we confirm for newly developed chicken lines that the chicken MHC haplotype B15 confers resistance to MD compared with B19. Using gas phase sequencing and immunopeptidomics, we find that BF2*1901 binds a greater variety of amino acids in some anchor positions than does BF2*1501. However, by x-ray crystallography, we find that the peptide-binding groove of BF2*1901 is narrower and shallower. Although the self-peptides that bound to BF2*1901 may appear more various than those of BF2*1501, the structures show that the wider and deeper peptide-binding groove of BF2*1501 allows stronger binding and thus more peptides overall, correlating with the expected hierarchies for expression level, thermostability, and MD resistance. Our study provides a reasonable explanation for greater promiscuity for BF2*1501 compared with BF2*1901, corresponding to the difference in resistance to MD.

The mode of action of tapasin on major histocompatibility class I (MHC-I) molecules

Tapasin (Tsn) plays a critical role in antigen processing and presentation by major histocompatibility complex class I (MHC-I) molecules. The mechanism of Tsn-mediated peptide loading and exchange hinges on the conformational dynamics governing the interaction of Tsn and MHC-I with recent structural and functional studies pinpointing the critical sites of direct or allosteric regulation. In this review, we highlight these recent findings and relate them to the extensive molecular and cellular data that are available for these evolutionary interdependent proteins. Furthermore, allotypic differences of MHC-I with regard to the editing and chaperoning function of Tsn are reviewed and related to the mechanistic observations. Finally, evolutionary aspects of the mode of action of Tsn will be discussed, a short comparison with the Tsn-related molecule TAPBPR (Tsn-related protein) will be given, and the impact of Tsn on noncanonical MHC-I molecules will be described.

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

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

The influences of microbial colonisation and germ-free status on the chicken TCRβ repertoire

Microbial colonisation is paramount to the normal development of the immune system, particularly at mucosal sites. However, the relationships between the microbiome and the adaptive immune repertoire have mostly been explored in rodents and humans. Here, we report a high-throughput sequencing analysis of the chicken TCRβ repertoire and the influences of microbial colonisation on tissue-resident TCRβ+ cells. The results reveal that the microbiome is an important driver of TCRβ diversity in both intestinal tissues and the bursa of Fabricius, but not in the spleen. Of note, public TCRβ sequences (shared across individuals) make a substantial contribution to the repertoire. Additionally, different tissues exhibit biases in terms of their V family and J gene usage, and these effects were influenced by the gut-associated microbiome. TCRβ clonal expansions were identified in both colonised and germ-free birds, but differences between the groups were indicative of an influence of the microbiota. Together, these findings provide an insight into the avian adaptive immune system and the influence of the microbiota on the TCRβ repertoire.

Genetic variation in chicken interferon signalling pathway genes in research lines showing differential viral resistance

Avian viruses of economic interest are a significant burden on the poultry industry, affecting production traits and resulting in mortality. Furthermore, the zoonosis of avian viruses risks pandemics developing in humans. Vaccination is the most common method of controlling viruses; however current vaccines often lack cross-protection against multiple strains of each virus. The mutagenicity of these viruses has also led to virulent strains emerging that can overcome the protection offered by vaccines. Breeding chickens with a more robust innate immune response may help in tackling current and emerging viruses. Understanding the genetic evolution of different lines will thus provide a useful tool in helping the host in the fight against pathogens. This study focuses on the interferon genes and their receptors in different chicken lines that are known to be more resistant or susceptible to particular avian viruses. Comparing genotypic differences in these core immune genes between the chicken lines may explain the phenotypic differences observed and aid the identification of causative variations. The relative resistance/susceptibility of each line to viruses of interest (Marek's disease virus, infectious bursal disease, infectious bronchitis virus and avian influenza virus) has previously been determined. Here we identify single nucleotide polymorphisms in interferons and downstream genes. Functional prediction tools were used to identify variants that may be affecting protein structure, mRNA secondary structure or transcription factor and micro-RNA binding sites. These variants were then considered in the context of the research lines and their distribution between phenotypes. We highlight 60 variants of interest in the interferon pathway genes that may account for susceptibility/resistance to viral pathogens.

Some thoughts about what non-mammalian jawed vertebrates are telling us about antigen processing and peptide loading of MHC molecules

The major histocompatibility complex (MHC) of mammals encodes highly polymorphic classical class I and class II molecules with crucial roles in immune responses, as well as various nonclassical molecules encoded by the MHC and elsewhere in the genome that have a variety of functions. These MHC molecules are supported by antigen processing and peptide loading pathways which are well-understood in mammals. This review considers what has been learned about the MHC, MHC molecules and the supporting pathways in non-mammalian jawed vertebrates. From the initial understanding from work with the chicken MHC, a great deal of diversity in the structure and function has been found. Are there underlying principles?

•

The genomic organisation of the MHC varies enormously across jawed vertebrates.

•

Total numbers of MHC genes vary among vertebrates, with only a few classical MHC genes.

•

Some nonclassical MHC and classical pathway genes appear earlier than others.

•

Obvious co-evolution within MHC pathways occurs in some species, but not others.

•

The promiscuity of interactions may correlate with differences in genomic organisation.

Peptide Centric Vβ Specific Germline Contacts Shape a Specialist T Cell Response

Certain CD8 T cell responses are particularly effective at controlling infection, as exemplified by elite control of HIV in individuals harboring HLA-B57. To understand the structural features that contribute to CD8 T cell elite control, we focused on a strongly protective CD8 T cell response directed against a parasite-derived peptide (HF10) presented by an atypical MHC-I molecule, H-2Ld. This response exhibits a focused TCR repertoire dominated by Vβ2, and a representative TCR (TG6) in complex with Ld-HF10 reveals an unusual structure in which both MHC and TCR contribute extensively to peptide specificity, along with a parallel footprint of TCR on its pMHC ligand. The parallel footprint is a common feature of Vβ2-containing TCRs and correlates with an unusual Vα-Vβ interface, CDR loop conformations, and Vβ2-specific germline contacts with peptides. Vβ2 and Ld may represent "specialist" components for antigen recognition that allows for particularly strong and focused T cell responses.

New vistas unfold: Chicken MHC molecules reveal unexpected ways to present peptides to the immune system

The functions of a wide variety of molecules with structures similar to the classical class I and class II molecules encoded by the major histocompatibility complex (MHC) have been studied by biochemical and structural studies over decades, with many aspects for humans and mice now enshrined in textbooks as dogma. However, there is much variation of the MHC and MHC molecules among the other jawed vertebrates, understood in the most detail for the domestic chicken. Among the many unexpected features in chickens is the co-evolution between polymorphic TAP and tapasin genes with a dominantly-expressed class I gene based on a different genomic arrangement compared to typical mammals. Another important discovery was the hierarchy of class I alleles for a suite of properties including size of peptide repertoire, stability and cell surface expression level, which is also found in humans although not as extreme, and which led to the concept of generalists and specialists in response to infectious pathogens. Structural studies of chicken class I molecules have provided molecular explanations for the differences in peptide binding compared to typical mammals. These unexpected phenomena include the stringent binding with three anchor residues and acidic residues at the peptide C-terminus for fastidious alleles, and the remodelling binding sites, relaxed binding of anchor residues in broad hydrophobic pockets and extension at the peptide C-terminus for promiscuous alleles. The first few studies for chicken class II molecules have already uncovered unanticipated structural features, including an allele that binds peptides by a decamer core. It seems likely that the understanding of how MHC molecules bind and present peptides to lymphocytes will broaden considerably with further unexpected discoveries through biochemical and structural studies for chickens and other non-mammalian vertebrates.

Efficient Identification of Tembusu Virus CTL Epitopes in Inbred HBW/B4 Ducks Using a Novel MHC Class I-Restricted Epitope Screening Scheme

The identification of MHC class I-restricted CTL epitopes in certain species, particularly nonmammals, remains a challenge. In this study, we developed a four-step identification scheme and confirmed its efficiency by identifying the Anpl-UAA*76-restricted CTL epitopes of Tembusu virus (TMUV) in inbred haplotype ducks HBW/B4. First, the peptide binding motif of Anpl-UAA*76 was determined by random peptide library in de novo liquid chromatography-tandem mass spectrometry, a novel nonbiased, data-independent acquisition method that we previously established. Second, a total of 38 TMUV peptides matching the motif were screened from the viral proteome, among which 11 peptides were conserved across the different TMUV strains. Third, the conserved TMUV peptides were refolded in vitro with Anpl-UAA*76 and Anpl-β₂-microglobulin to verify the results from the previous two steps. To clarify the structural basis of the obtained motif, we resolved the crystal structure of Anpl-UAA*76 with the TMUV NS3 peptide LRKRQLTVL and found that Asp³⁴ is critical for the preferential binding of the B pocket to bind the second residue to arginine as an anchor residue. Fourth, the immunogenicity of the conserved TMUV peptides was tested in vivo using specific pathogen-free HBW/B4 ducks immunized with the attenuated TMUV vaccine. All 11 conserved TMUV epitopes could bind stably to Anpl-UAA*76 in vitro and stimulate the secretion of IFN-γ and lymphocyte proliferation, and three conserved and one nonconserved peptides were selected to evaluate the CTL responses in vivo by flow cytometry and their tetramers. We believe that this new scheme could improve the identification of MHC class I-restricted CTL epitopes, and our data provide a foundation for further study on duck anti-TMUV CTL immunity.

Caecal microbiota composition of experimental inbred MHC-B lines infected with IBV differs according to genetics and vaccination

Interactions between the gut microbiota and the immune system may be involved in vaccine and infection responses. In the present study, we studied the interactions between caecal microbiota composition and parameters describing the immune response in six experimental inbred chicken lines harboring different MHC haplotypes. Animals were challenge-infected with the infectious bronchitis virus (IBV), and half of them were previously vaccinated against this pathogen. We explored to what extent the gut microbiota composition and the genetic line could be related to the immune response, evaluated through flow cytometry. To do so, we characterized the caecal bacterial communities with a 16S rRNA gene amplicon sequencing approach performed one week after the IBV infectious challenge. We observed significant effects of both the vaccination and the genetic line on the microbiota after the challenge infection with IBV, with a lower bacterial richness in vaccinated chickens. We also observed dissimilar caecal community profiles among the different lines, and between the vaccinated and non-vaccinated animals. The effect of vaccination was similar in all the lines, with a reduced abundance of OTU from the Ruminococcacea UCG-014 and Faecalibacterium genera, and an increased abundance of OTU from the Eisenbergiella genus. The main association between the caecal microbiota and the immune phenotypes involved TCR_ϒδ expression on TCR_ϒδ⁺ T cells. This phenotype was negatively associated with OTU from the Escherichia-Shigella genus that were also less abundant in the lines with the highest responses to the vaccine. We proved that the caecal microbiota composition is associated with the IBV vaccine response level in inbred chicken lines, and that the TCR_ϒδ⁺ T cells (judged by TCR_ϒδ expression) may be an important component involved in this interaction, especially with bacteria from the Escherichia-Shigella genus. We hypothesized that bacteria from the Escherichia-Shigella genus increased the systemic level of bacterial lipid antigens, which subsequently mitigated poultry γδ T cells.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

The Diverse Major Histocompatibility Complex Haplotypes of a Common Commercial Chicken Line and Their Effect on Marek's Disease Virus Pathogenesis and Tumorigenesis

The major histocompatibility complex (MHC) is crucial for appropriate immune responses against invading pathogens. Chickens possess a single predominantly-expressed class I molecule with strong associations between disease resistance and MHC haplotype. For Marek's disease virus (MDV) infections of chickens, the MHC haplotype is one of the major determinants of genetic resistance and susceptibility. VALO specific pathogen free (SPF) chickens are widely used in biomedical research and vaccine production. While valuable findings originate from MDV infections of VALO SPF chickens, their MHC haplotypes and associated disease resistance remained elusive. In this study, we used several typing systems to show that VALO SPF chickens possess MHC haplotypes that include B9, B9:02, B15, B19 and B21 at various frequencies. Moreover, we associate the MHC haplotypes to MDV-induced disease and lymphoma formation and found that B15 homozygotes had the lowest tumor incidence while B21 homozygotes had the lowest number of organs with tumors. Finally, we found transmission at variable levels to all contact birds except B15/B21 heterozygotes. These data have immediate implications for the use of VALO SPF chickens and eggs in the life sciences and add another piece to the puzzle of the chicken MHC complex and its role in infections with this oncogenic herpesvirus.

Dissociation of β2m from MHC class I Triggers formation of Noncovalent, transient heavy chain dimers

At the plasma membrane of mammalian cells, major histocompatibility complex class I molecules (MHC-I) present antigenic peptides to cytotoxic T cells. Following the loss of the peptide and the light chain beta-2 microglobulin (β2m), the resulting free heavy chains (FHCs) can associate into homotypic complexes in the plasma membrane. Here, we investigate the stoichiometry and dynamics of MHC-I FHCs assemblies by combining a micropattern assay with fluorescence recovery after photobleaching (FRAP) and with single molecule co-tracking. We identify non-covalent MHC-I FHC dimers mediated by the α3 domain as the prevalent species at the plasma membrane, leading a moderate decrease in the diffusion coefficient. MHC-I FHC dimers show increased tendency to cluster into higher order oligomers as concluded from an increased immobile fraction with higher single molecule co-localization. In vitro studies with isolated proteins in conjunction with molecular docking and dynamics simulations suggest that in the complexes, the α3 domain of one FHC binds to another FHC in a manner similar to the β2m light chain.

Evolutionary pressures rendered by animal husbandry practices for avian influenza viruses to adapt to humans

Commercial poultry operations produce and crowd billions of birds every year, which is a source of inexpensive animal protein. Commercial poultry is intensely bred for desirable production traits, and currently presents very low variability at the major histocompatibility complex. This situation dampens the advantages conferred by the MHC’s high genetic variability, and crowding generates immunosuppressive stress. We address the proteins of influenza A viruses directly and indirectly involved in host specificities. We discuss how mutants with increased virulence and/or altered host specificity may arise if few class I alleles are the sole selective pressure on avian viruses circulating in immunocompromised poultry. This hypothesis is testable with peptidomics of MHC ligands. Breeding strategies for commercial poultry can easily and inexpensively include high variability of MHC as a trait of interest, to help save billions of dollars as a disease burden caused by influenza and decrease the risk of selecting highly virulent strains.

Genome-wide association study of disease resilience traits from a natural polymicrobial disease challenge model in pigs identifies the importance of the major histocompatibility complex region

Infectious diseases cause tremendous financial losses in the pork industry, emphasizing the importance of disease resilience, which is the ability of an animal to maintain performance under disease. Previously, a natural polymicrobial disease challenge model was established, in which pigs were challenged in the late nursery phase by multiple pathogens to maximize expression of genetic differences in disease resilience. Genetic analysis found that performance traits in this model, including growth rate, feed and water intake, and carcass traits, as well as clinical disease phenotypes, were heritable and could be selected for to increase disease resilience of pigs. The objectives of the current study were to identify genomic regions that are associated with disease resilience in this model, using genome-wide association studies and fine-mapping methods, and to use gene set enrichment analyses to determine whether genomic regions associated with disease resilience are enriched for previously published quantitative trait loci, functional pathways, and differentially expressed genes subject to physiological states. Multiple quantitative trait loci were detected for all recorded performance and clinical disease traits. The major histocompatibility complex region was found to explain substantial genetic variance for multiple traits, including for growth rate in the late nursery (12.8%) and finisher (2.7%), for several clinical disease traits (up to 2.7%), and for several feeding and drinking traits (up to 4%). Further fine mapping identified 4 quantitative trait loci in the major histocompatibility complex region for growth rate in the late nursery that spanned the subregions for class I, II, and III, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects, explaining 0.8–27.1% of genetic variance for growth rate and for multiple clinical disease traits. This single-nucleotide polymorphism was located in the enhancer of TRIM39 gene, which is involved in innate immune response. The major histocompatibility complex region was pleiotropic for growth rate in the late nursery and finisher, and for treatment and mortality rates. Growth rate in the late nursery showed strong negative genetic correlations in the major histocompatibility complex region with treatment or mortality rates (−0.62 to −0.85) and a strong positive genetic correlation with growth rate in the finisher (0.79). Gene set enrichment analyses found genomic regions associated with resilience phenotypes to be enriched for previously identified disease susceptibility and immune capacity quantitative trait loci, for genes that were differentially expressed following bacterial or virus infection and immune response, and for gene ontology terms related to immune and inflammatory response. In conclusion, the major histocompatibility complex and other quantitative trait loci that harbor immune-related genes were identified to be associated with disease resilience traits in a large-scale natural polymicrobial disease challenge. The major histocompatibility complex region was pleiotropic for growth rate under challenge and for clinical disease traits. Four quantitative trait loci were identified across the class I, II, and III subregions of the major histocompatibility complex for nursery growth rate under challenge, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects. The major histocompatibility complex and other quantitative trait loci identified play an important role in host response to infectious diseases and can be incorporated in selection to improve disease resilience, in particular the identified single-nucleotide polymorphism in the major histocompatibility complex class I subregion.

Reconstructing Macroevolutionary Patterns in Avian MHC Architecture With Genomic Data

The Major Histocompatibility Complex (MHC) is a hyper-polymorphic genomic region, which forms a part of the vertebrate adaptive immune system and is crucial for intra- and extra-cellular pathogen recognition (MHC-I and MHC-IIA/B, respectively). Although recent advancements in high-throughput sequencing methods sparked research on the MHC in non-model species, the evolutionary history of MHC gene structure is still poorly understood in birds. Here, to explore macroevolutionary patterns in the avian MHC architecture, we retrieved contigs with antigen-presenting MHC and MHC-related genes from available genomes based on third-generation sequencing. We identified: 1) an ancestral avian MHC architecture with compact size and tight linkage between MHC-I, MHC-IIA/IIB and MHC-related genes; 2) three major patterns of MHC-IIA/IIB unit organization in different avian lineages; and 3) lineage-specific gene translocation events (e.g., separation of the antigen-processing TAP genes from the MHC-I region in passerines), and 4) the presence of a single MHC-IIA gene copy in most taxa, showing evidence of strong purifying selection (low dN/dS ratio and low number of positively selected sites). Our study reveals long-term macroevolutionary patterns in the avian MHC architecture and provides the first evidence of important transitions in the genomic arrangement of the MHC region over the last 100 million years of bird evolution.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

Immunological Diversity Is a Cornerstone of Organismal Defense and Allorecognition across Metazoa

The ongoing arms race between hosts and microbes has fueled the evolution of novel strategies for diversifying the molecules involved in immune responses. Characterization of immune systems from an ever-broadening phylogenetic range of organisms reveals that there are many mechanisms by which this diversity can be generated and maintained. Diversification strategies operate at the level of populations, genomes, genes, and even individual transcripts. Lineage-specific innovations have been cataloged within the immune systems of both invertebrates and vertebrates. Furthermore, somatic diversification of immune receptor genes has now been described in jawless vertebrates and some invertebrate species. In addition to pathogen detection, immunological diversity plays important roles in several distinct allorecognition systems. In this Brief Review, we highlight some of the evolutionary innovations employed by a variety of metazoan species to generate the molecular diversity required to detect a vast array of molecules in the context of both immune response and self/nonself-recognition.

Approaching Genetics Through the MHC Lens: Tools and Methods for HLA Research

The current SARS-CoV-2 pandemic era launched an immediate and broad response of the research community with studies both about the virus and host genetics. Research in genetics investigated HLA association with COVID-19 based on in silico, population, and individual data. However, they were conducted with variable scale and success; convincing results were mostly obtained with broader whole-genome association studies. Here, we propose a technical review of HLA analysis, including basic HLA knowledge as well as available tools and advice. We notably describe recent algorithms to infer and call HLA genotypes from GWAS SNPs and NGS data, respectively, which opens the possibility to investigate HLA from large datasets without a specific initial focus on this region. We thus hope this overview will empower geneticists who were unfamiliar with HLA to run MHC-focused analyses following the footsteps of the Covid-19|HLA & Immunogenetics Consortium.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

Evidence of MHC class I and II influencing viral and helminth infection via the microbiome in a non-human primate

Until recently, the study of major histocompability complex (MHC) mediated immunity has focused on the direct link between MHC diversity and susceptibility to parasite infection. However, MHC genes can also influence host health indirectly through the sculpting of the bacterial community that in turn shape immune responses. We investigated the links between MHC class I and II gene diversity gut microbiome diversity and micro- (adenovirus, AdV) and macro- (helminth) parasite infection probabilities in a wild population of non-human primates, mouse lemurs of Madagascar. This setup encompasses a plethora of underlying interactions between parasites, microbes and adaptive immunity in natural populations. Both MHC classes explained shifts in microbiome composition and the effect was driven by a few select microbial taxa. Among them were three taxa (Odoribacter, Campylobacter and Prevotellaceae-UCG-001) which were in turn linked to AdV and helminth infection status, correlative evidence of the indirect effect of the MHC via the microbiome. Our study provides support for the coupled role of MHC diversity and microbial flora as contributing factors of parasite infection.

The selective pressure of the major histocompatibility complex (MHC) on microbial communities, and the potential role of this interaction in driving parasite resistance has been largely neglected. Using a natural population of the primate Microcebus griseorufus, we provide correlative evidence of two outstanding findings: that MHCI and MHCII diversity shapes the composition of the gut microbiota; and that select taxa associated with MHC diversity predicted adenovirus and helminth infection status. Our study highlights the importance of incorporating the microbiome when investigating parasite-mediated MHC selection.

Innate immune genes of the chicken MHC and related regions

Compared to the major histocompatibility complex (MHC) of typical mammals, the chicken BF/BL region is small and simple, with most of the genes playing central roles in the adaptive immune response. However, some genes of the chicken MHC are almost certainly involved in innate immunity, such as the complement component C4 and the lectin-like receptor/ligand gene pair BNK and Blec. The poorly expressed classical class I molecule BF1 is known to be recognised by natural killer (NK) cells and, analogous to mammalian immune responses, the classical class I molecules BF1 and BF2, the CD1 homologs and the butyrophilin homologs called BG may be recognised by adaptive immune lymphocytes with semi-invariant receptors in a so-called adaptate manner. Moreover, the TRIM and BG regions next to the chicken MHC, along with the genetically unlinked Y and olfactory/scavenger receptor regions on the same chromosome, have multigene families almost certainly involved in innate and adaptate responses. On this chicken microchromosome, the simplicity of the adaptive immune gene systems contrasts with the complexity of the gene systems potentially involved in innate immunity.

Comparative Analysis of Susceptibility and Severity of COVID-19 in Countries from the Eastern and the Western World Till March '21

Majority of the world’s human population today is affected by Covid-19. The disease has not only exhibited differences in susceptibility among people of different countries, but also the mortality rate. In general, Western world has been reporting a greater number of infected cases than eastern countries. Even the mortality rates are quite high there. The aim of this study was to analyse the data available on the infectivity and mortality rates of Covid-19 in different countries till March’21 and then reviewed the literature to find reasons for the differences in susceptibility and severity in eastern and western countries. The reasons for the observed differences may be: (i) Eastern countries followed stricter modalities and got grace period to create better healthcare facilities to tackle COVID-19. This probably also slowed the transmission of virus and its evolution, (ii) Vaccination policies in the east may have provided some immunity due to cross reactivity, (iii) Frequent exposure to infections at young age in eastern countries might be helping in better immunity, (iv) Mutations in viral genome may be geography based and (v) Genetic differences in the immune system of the hosts with respect to ACE receptors and MHC may be playing an important role. In this article, an attempt has been made to put forth and discuss these plausible reasons along with suitable evidences. These findings may help in future research on the diagnosis, treatment and prevention of Covid-19.

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.

Genetic Diversity of MHC B-F/B-L Region in 21 Chicken Populations

The chicken major histocompatibility complex (MHC) on chromosome 16 is the most polymorphic region across the whole genome, and also an ideal model for genetic diversity investigation. The MHC B-F/B-L region is 92 kb in length with high GC content consisting of 18 genes and one pseudogene (Blec4), which plays important roles in immune response. To evaluate polymorphism of the Chinese indigenous chickens as well as to analyze the effect of selection to genetic diversity, we used WaferGen platform to identify sequence variants of the B-F/B-L region in 21 chicken populations, including the Red Jungle Fowl (RJF), Cornish (CS), White Leghorns (WLs), 16 Chinese domestic breeds, and two well-known inbred lines 63 and 72. A total of 3,319 single nucleotide polymorphism (SNPs) and 181 INDELs in the B-F/B-L region were identified among 21 populations, of which 2,057 SNPs (62%) and 159 INDELs (88%) were novel. Most of the variants were within the intron and the flanking regions. The average variation density was 36 SNPs and 2 INDELs per kb, indicating dramatical high diversity of this region. Furthermore, BF2 was identified as the hypervariable genes with 67 SNPs per kb. Chinese domestic populations showed higher diversity than the WLs and CS. The indigenous breeds, Nandan Yao (NY), Xishuangbanna Game (XG), Gushi (GS), and Xiayan (XY) chickens, were the top four with the highest density of SNPs and INDELs. The highly inbred lines 63 and 72 have the lowest diversity, which might be resulted from a long-term intense selection for decades. Collectively, we refined the genetic map of chicken MHC B-F/B-L region, and illustrated genetic diversity of 21 chicken populations. Abundant genetic variants were identified, which not only strikingly expanded the current Ensembl SNP database, but also provided comprehensive data for researchers to further investigate association between variants in MHC and immune traits.

Coevolution between MHC class I and Antigen Processing Genes in salamanders

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Differential Frequencies of HLA-DRB1, DQA1, and DQB1 Alleles and Haplotypes Are Observed in the Arbovirus-Related Neurological Syndromes

BACKGROUND

We took advantage of the 2015-2016 Brazilian arbovirus outbreak (Zika [ZIKV]/dengue/chikungunya viruses) associated with neurological complications to type HLA-DRB1/DQA1/DQB1 variants in patients exhibiting neurological complications and in bone marrow donors from the same endemic geographical region.

METHODS

DRB1/DQA1/DQB1 loci were typed using sequence-specific oligonucleotides. In silico studies were performed using X-ray resolved dimer constructions.

RESULTS

The DQA1*01, DQA1*05, DQB1*02, or DQB1*06 genotypes/haplotypes and DQA1/DQB1 haplotypes that encode the putative DQA1/DQB1 dimers were overrepresented in the whole group of patients and in patients exhibiting peripheral neurological spectrum disorders (PSD) or encephalitis spectrum disorders (ESD). The DRB1*04, DRB1*13, and DQA1*03 allele groups protected against arbovirus neurological manifestation, being underrepresented in whole group of patients and ESD and PSD groups. Genetic and in silico studies revealed that DQA1/DQB1 dimers (1) were primarily associated with susceptibility to arbovirus infections; (2) can bind to a broad range of ZIKV peptides (235 of 1878 peptides, primarily prM and NS2A); and (3) exhibited hydrophilic and highly positively charged grooves when compared to the DRA1/DRB1 cleft. The protective dimer (DRA1/DRB1*04) bound a limited number of ZIKV peptides (40 of 1878 peptides, primarily prM).

CONCLUSION

Protective haplotypes may recognize arbovirus peptides more specifically than susceptible haplotypes.

Immune Equilibrium Depends on the Interaction Between Recognition and Presentation Landscapes

In this review, we described the structure and organization of antigen-recognizing repertoires of B and T cells from the standpoint of modern immunology. We summarized the latest advances in bioinformatics analysis of sequencing data from T and B cell repertoires and also presented contemporary ideas about the mechanisms of clonal diversity formation at different stages of organism development. At the same time, we focused on the importance of the allelic variants of the HLA genes and spectra of presented antigens for the formation of T-cell receptors (TCR) landscapes. The main idea of this review is that immune equilibrium and proper functioning of immunity are highly dependent on the interaction between the recognition and the presentation landscapes of antigens. Certain changes in these landscapes can occur during life, which can affect the protective function of adaptive immunity. We described some mechanisms associated with these changes, for example, the conversion of effector cells into regulatory cells and vice versa due to the trans-differentiation or bystander effect, changes in the clonal organization of the general TCR repertoire due to homeostatic proliferation or aging, and the background for the altered presentation of some antigens due to SNP mutations of MHC, or the alteration of the presenting antigens due to post-translational modifications. The authors suggest that such alterations can lead to an increase in the risk of the development of oncological and autoimmune diseases and influence the sensitivity of the organism to different infectious agents.

Know thy immune self and non-self: Proteomics informs on the expanse of self and non-self, and how and where they arise

T cells play an important role in the adaptive immune response to a variety of infections and cancers. Initiation of a T cell mediated immune response requires antigen recognition in a process termed MHC (major histocompatibility complex) restri ction. A T cell antigen is a composite structure made up of a peptide fragment bound within the antigen‐binding groove of an MHC‐encoded class I or class II molecule. Insight into the precise composition and biology of self and non‐self immunopeptidomes is essential to harness T cell mediated immunity to prevent, treat, or cure infectious diseases and cancers. T cell antigen discovery is an arduous task! The pioneering work in the early 1990s has made large‐scale T cell antigen discovery possible. Thus, advancements in mass spectrometry coupled with proteomics and genomics technologies make possible T cell antigen discovery with ease, accuracy, and sensitivity. Yet we have only begun to understand the breadth and the depth of self and non‐self immunopeptidomes because the molecular biology of the cell continues to surprise us with new secrets directly related to the source, and the processing and presentation of MHC ligands. Focused on MHC class I molecules, this review, therefore, provides a brief historic account of T cell antigen discovery and, against a backdrop of key advances in molecular cell biologic processes, elaborates on how proteogenomics approaches have revolutionised the field.

Allele-Specific Thresholds of Eluted Ligands for T-Cell Epitope Prediction

A common strategy for predicting candidate human leukocyte antigen class I T-cell epitopes is to use an affinity-based threshold of 500 nM. Although a 500 nM threshold across alleles effectively identifies most epitopes across alleles, findings showing that major histocompatibility complex repertoire sizes vary by allele indicate that using thresholds specific to individual alleles may improve epitope identification. In this work, we compare different strategies utilizing common and allele-specific thresholds to identify an optimal approach for T-cell epitope prediction. First, we confirmed previous observations that different human leukocyte antigen class I alleles correspond with varying repertoire sizes. Here, we define general and allele-specific thresholds that capture 80% of eluted ligands and a different set of thresholds associated with capturing 9-mer T-cell epitopes at 80% sensitivity. Our analysis revealed that allele-specific threshold performance was roughly equivalent to that of a common threshold when considering a relatively large number of alleles. However, when predicting epitopes for only a few alleles, the use of allele-specific thresholds would be preferable. Finally, we present here for public use a set of allele-specific thresholds that may be used to identify T-cell epitopes at 80% sensitivity.

•

Confirmed findings that different HLA class I alleles have varying repertoire sizes.

•

Defined common and allele-specific thresholds that capture 80% of eluted ligands.

•

Defined common and allele-specific thresholds that capture 80% of T-cell epitopes.

•

Allele-specific thresholds perform more consistently when analyzing a few alleles.

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.

MHC-B variation in maternal and paternal synthetic lines of the Argentinian Campero INTA chicken

The Campero-INTA chicken of Argentina was developed to provide a robust bird that can survive under Argentinian pasture conditions with no significant additional nutrition, producing a source of animal protein for small producers or low-income families. In previous work, we described the AH paternal line of Campero and its Major Histocompatibility Complex B region (MHC-B) variation. In this work we analyzed the three remaining synthetic lines used to produce the Campero-INTA production bird: lines AS, A, and E. Because of the association between variation within the MHC of chickens and disease resistance, MHC variation within this breed is of particular interest. MHC variability within the lines used to produce the Campero-INTA chicken was examined using a 90 SNP panel encompassing the chicken MHC-B region plus the VNTR, LEI0258, located within the chicken MHC. Across all 4lines 12 haplotypes were found, with 7 of these being previously reported in North America/European breeds, reflecting the original breed sources for these birds. Three Campero unique haplotypes were found, 2 of which likely originated from MHC recombination events. MHC-B variation for all lines involved with production of the final Campero-INTA bird has now been determined.

Chickens as a simple system for scientific discovery: The example of the MHC

Chickens have played many roles in human societies over thousands of years, most recently as an important model species for scientific discovery, particularly for embryology, virology and immunology. In the last few decades, biomedical models like mice have become the most important model organism for understanding the mechanisms of disease, but for the study of outbred populations, they have many limitations. Research on humans directly addresses many questions about disease, but frank experiments into mechanisms are limited by practicality and ethics. For research into all levels of disease simultaneously, chickens combine many of the advantages of humans and of mice, and could provide an independent, integrated and overarching system to validate and/or challenge the dogmas that have arisen from current biomedical research. Moreover, some important systems are simpler in chickens than in typical mammals. An example is the major histocompatibility complex (MHC) that encodes the classical MHC molecules, which play crucial roles in the innate and adaptive immune systems. Compared to the large and complex MHCs of typical mammals, the chicken MHC is compact and simple, with single dominantly-expressed MHC molecules that can determine the response to infectious pathogens. As a result, some fundamental principles have been easier to discover in chickens, with the importance of generalist and specialist MHC alleles being the latest example.

Research Note: Rous sarcoma growth differs among congenic lines containing major histocompatibility (B) complex recombinants

New arrangements of chicken major histocompatibility complex (MHC) class I BF and class IV BG genes are created through recombination. Characterizing the immune responses of such recombinants reveals genes or gene regions that contribute to immunity. Inbred Line UCD 003 (B17B17) served as the genetic background for congenic lines, each containing a unique MHC recombinant. After an initial cross to introduce a specific recombinant, 10 backcrosses to the inbred line produced lines with 99.9% genetic uniformity. The current study compared Rous sarcoma virus (RSV) tumor growth in 5 congenic lines homozygous for MHC recombinants (003.R1 = BF24-BG23, 003.R2 = BF2-BG23, 003.R4 = BF2-BG23, 003.R5 = BF21-BG19, and 003.R13 = BF17-BG23). Two experiments used a total of 70 birds from the 5 congenic lines inoculated with 20 pock forming units of RSV subgroup C at 6 wk of age. Tumor size was scored 6 times over 10 wk postinoculation followed by assignment of a tumor profile index (TPI) based on the tumor size scores. Tumor growth over time and rank transformed TPI values were analyzed by least squares ANOVA. Tumor size increased over the experimental period in all genotypes through 4 wk postinoculation. After this time, tumor size increased in Lines 003.R1, plateaued in Lines 003.R2, 003.R4, and 003.R13, and declined in 003.R5. Tumor growth over time was significantly lower in Line 003.R5 compared with all other genotypes. In addition, Line 003.R5 chickens had significantly lower TPI values compared with Lines 003.R2, 003.R4, and 003.R13. The TPI of Line 003.R1 did not differ significantly from any of the other genotypes. The BF21 in Line 003.R5 produced a greater response against subgroup C RSV tumors than did BF24, found in 003.R1; BF2 found in 003.R2 and R4 as well as BF17 found in 003.R13.