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

Haplotype Analysis Reveals Pleiotropic Disease Associations in the HLA Region

The human leukocyte antigen (HLA) region plays an important role in human health through involvement in immune cell recognition and maturation. While genetic variation in the HLA region is associated with many diseases, the pleiotropic patterns of these associations have not been systematically investigated. Here, we developed a haplotype approach to investigate disease associations phenome-wide for 412,181 Finnish individuals and 2,459 traits. Across the 1,035 diseases with a GWAS association, we found a 17-fold average per-SNP enrichment of hits in the HLA region. Altogether, we identified 7,649 HLA associations across 647 traits, including 1,750 associations uncovered by haplotype analysis. We find some haplotypes show trade-offs between diseases, while others consistently increase risk across traits, indicating a complex pleiotropic landscape involving a range of diseases. This study highlights the extensive impact of HLA variation on disease risk, and underscores the importance of classical and non-classical genes, as well as non-coding variation.

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

A journey to your self: The vague definition of immune self and its practical implications

The identification of immunogenic peptides has become essential in an increasing number of fields in immunology, ranging from tumor immunotherapy to vaccine development. The nature of the adaptive immune response is shaped by the similarity between foreign and self-protein sequences, a concept extensively applied in numerous studies. Can we precisely define the degree of similarity to self? Furthermore, do we accurately define immune self? In the current work, we aim to unravel the conceptual and mechanistic vagueness hindering the assessment of self-similarity. Accordingly, we demonstrate the remarkably low consistency among commonly employed measures and highlight potential avenues for future research.

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

Immune-driven clonal cell selection at the intersection among cancer, infections, autoimmunity and senescence

Immune surveillance mechanisms play a crucial role in maintaining lifelong immune homeostasis in response to pathologic stimuli and aberrant cell states. However, their persistence, especially in the context of chronic antigenic exposure, can create a fertile ground for immune evasion. These escaping cell phenotypes, harboring a variety of genomic and transcriptomic aberrances, chiefly in human leukocyte antigen (HLA) and antigen presentation machinery genes, may survive and proliferate, featuring a scenario of clonal cell expansion with immune failure characteristics. While well characterized in solid and, to some extent, hematological malignancies, little is known about their occurrence and significance in other disease contexts. Historical literature highlights the role for escaping HLA-mediated recognition as a strategy adopted by virus to evade from the immune system, hinting at the potential for immune aberrant cell expansion in the context of chronic infections. Additionally, unmasked in idiopathic aplastic anemia as a mechanism able to rescue failing hematopoiesis, HLA clonal escape may operate in autoimmune disorders, particularly in tissues targeted by aberrant immune responses. Furthermore, senescent cell status emerging as immunogenic phenotypes stimulating T cell responses, may act as a bottleneck for the selection of such immune escaping clones, blurring the boundaries between neoplastic transformation, aging and inflammation. Here we provide a fresh overview and perspective on this immune-driven clonal cell expansion, linking pathophysiological features of neoplastic, autoimmune, infectious and senescence processes exposed to immune surveillance.

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.

Host Transcriptional Regulatory Genes and Microbiome Networks Crosstalk through Immune Receptors Establishing Normal and Tumor Multiomics Metafirm of the Oral-Gut-Lung Axis

The microbiome has shown a correlation with the diet and lifestyle of each population in health and disease, the ability to communicate at the cellular level with the host through innate and adaptative immune receptors, and therefore an important role in modulating inflammatory process related to the establishment and progression of cancer. The oral cavity is one of the most important interaction windows between the human body and the environment, allowing the entry of an important number of microorganisms and their passage across the gastrointestinal tract and lungs. In this review, the contribution of the microbiome network to the establishment of systemic diseases like cancer is analyzed through their synergistic interactions and bidirectional crosstalk in the oral-gut-lung axis as well as its communication with the host cells. Moreover, the impact of the characteristic microbiota of each population in the formation of the multiomics molecular metafirm of the oral-gut-lung axis is also analyzed through state-of-the-art sequencing techniques, which allow a global study of the molecular processes involved of the flow of the microbiota environmental signals through cancer-related cells and its relationship with the establishment of the transcription factor network responsible for the control of regulatory processes involved with tumorigenesis.

Pathogens are an important driving force for the rapid spread of symbionts in an insect host

One of the biggest challenges for pathogens invading hosts is microbial symbionts but the role of pathogens in symbionts in nature is unknown. By tracking the dynamics of the entomopathogenic fungal Cordyceps javanica and symbionts in natural populations of the whitefly Bemisia tabaci from 2016 to 2021 across China, we reveal that Rickettsia, a newly invaded symbiont, is positively correlated with the pathogen in both frequency and density. We confirm that applying pathogen pressure can selectively drive Rickettsia to sudden fixation in whiteflies both in the laboratory and in the field. Furthermore, the driving force is elucidated by the Rickettsia-conferred suppression of pathogen infection quantity, proliferation and sporulation, acting as a potential barrier of onward transmission of the pathogen. These results show that pathogens are an important driving force for rapid shifts in host symbionts in the natural niche.

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.

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.

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.

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.

Large registry-based analysis of genetic predisposition to tuberculosis identifies genetic risk factors at HLA

Tuberculosis is a significant public health concern resulting in the death of over 1 million individuals each year worldwide. While treatment options and vaccines exist, a substantial number of infections still remain untreated or are caused by treatment resistant strains. Therefore, it is important to identify mechanisms that contribute to risk and prognosis of tuberculosis as this may provide tools to understand disease mechanisms and provide novel treatment options for those with severe infection. Our goal was to identify genetic risk factors that contribute to the risk of tuberculosis and to understand biological mechanisms and causality behind the risk of tuberculosis. A total of 1895 individuals in the FinnGen study had International Classification of Diseases-based tuberculosis diagnosis. Genome-wide association study analysis identified genetic variants with statistically significant association with tuberculosis at the human leukocyte antigen (HLA) region (P < 5e-8). Fine mapping of the HLA association provided evidence for one protective haplotype tagged by HLA DQB1*05:01 (P = 1.82E-06, OR = 0.81 [CI 95% 0.74-0.88]), and predisposing alleles tagged by HLA DRB1*13:02 (P = 0.00011, OR = 1.35 [CI 95% 1.16-1.57]). Furthermore, genetic correlation analysis showed association with earlier reported risk factors including smoking (P < 0.05). Mendelian randomization supported smoking as a risk factor for tuberculosis (inverse-variance weighted P < 0.05, OR = 1.83 [CI 95% 1.15-2.93]) with no significant evidence of pleiotropy. Our findings indicate that specific HLA alleles associate with the risk of tuberculosis. In addition, lifestyle risk factors such as smoking contribute to the risk of developing tuberculosis.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

Incidence of respiratory virus illness and hospitalizations in a Panama and El Salvador birth cohort, 2014-2018

Background

Respiratory viruses remain a key cause of early childhood illness, hospitalization, and death globally.The recent pandemic has rekindled interest in the control of respiratory viruses among paediatric populations. We estimate the burden of such viruses among children <2 years.

Methods

Enrolled neonates were followed until two years of age. Weekly active symptom monitoring for the development of acute respiratory illnesses (ARI) defined as cough, rhinorrhoea, difficulty breathing, asthenia, anorexia, irritability, or vomiting was conducted. When the child had ARI and fever, nasopharyngeal swabbing was performed, and samples were tested through singleplex RT-PCR. Incidence of respiratory viruses was calculated by dividing the number of laboratory-confirmed detections by the person-time accrued during weeks when that virus was detectable through national surveillance then corrected for under-ascertainment among untested children.

Findings

During December 2014-November 2017, 1567 enrolled neonates contributed 2,186.9 person-years (py). Six in ten (64·4%) children developed ARI (total 2493 episodes). Among children <2 years, incidence of respiratory syncytial virus (RSV)-associated ARI episodes (21·0, 95%CI 19·3-22·8, per 100py) and rhinovirus-associated (20·5, 95%CI 20·4-20·7) were similar and higher than parainfluenza 1-3-associated (14·2, 95%CI 12·2-16·1), human metapneumovirus-associated (9·2, 95%CI 7·7-10·8), influenza-associated (5·9, 95%CI 4·4-7·5), and adenovirus-associated ARI episodes (5·1, 95%CI 5·0-5·2). Children aged <3 months had the highest rates of RSV ARI (49·1, 95%CI 44·0-54·1 per 100py) followed by children aged 3-5 (25·1, 95%CI 20·1-30·0), 6-11 (17·6, 95%CI 13·2-21·9), and 12-23 months (11·9, 95%CI 10·8-12·9). One in ten children with RSV was referred to the hospital (2·5, 95%CI 2·1-2·8, per 100py).

Interpretation

Children frequently developed viral ARI and a substantive proportion required hospital care. Such findings suggest the importance of exploring the value of new interventions and increasing uptake of existing prevention measures to mitigate burden of epidemic-prone respiratory viruses.

Funding

The study was supported by the Centers for Disease Control and Prevention.

The immunogenetic impact of European colonization in the Americas

The introduction of pathogens originating from Eurasia into the Americas during early European contact has been associated with high mortality rates among Indigenous peoples, likely contributing to their historical and precipitous population decline. However, the biological impacts of imported infectious diseases and resulting epidemics, especially in terms of pathogenic effects on the Indigenous immunity, remain poorly understood and highly contentious to this day. Here, we examine multidisciplinary evidence underpinning colonization-related immune genetic change, providing contextualization from anthropological studies, paleomicrobiological evidence of contrasting host-pathogen coevolutionary histories, and the timings of disease emergence. We further summarize current studies examining genetic signals reflecting post-contact Indigenous population bottlenecks, admixture with European and other populations, and the putative effects of natural selection, with a focus on ancient DNA studies and immunity-related findings. Considering current genetic evidence, together with a population genetics theoretical approach, we show that post-contact Indigenous immune adaptation, possibly influenced by selection exerted by introduced pathogens, is highly complex and likely to be affected by multifactorial causes. Disentangling putative adaptive signals from those of genetic drift thus remains a significant challenge, highlighting the need for the implementation of population genetic approaches that model the short time spans and complex demographic histories under consideration. This review adds to current understandings of post-contact immunity evolution in Indigenous peoples of America, with important implications for bettering our understanding of human adaptation in the face of emerging infectious diseases.

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.

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.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Allelic and haplotypic data of MHC class II Alu insertions in Ngazidja (Comoros archipelago) and insight on its historical biology

This study aimed to analyse, for the first time, five MHC class II polymorphic Alu insertions in a population with a strong Sub-Saharan African genetic background: the Ngazidja islanders and compare its allelic and haplotypic data with Worldwide populations. The genotyping was performed in 80 individuals, using simple PCR and agarose gel electrophoresis methods. Allele and haplotype frequencies, genetic diversity, Hardy-Weinberg equilibrium deviations, normalized deviate of homozygotes and pairwise linkage disequilibrium were estimated. The phylogenetic analyses included the available population data. In Ngazidja, the MHC class II Alu insertion frequencies ranged from 0.119 for AluORF10 to 0.588 for AluDPB2. Concerning haplotypes, the most predominant were the ones with only the AluDPB2 insertion allele (AluDPB2*2-AluDQA2*1-AluDQA1*1-AluDRB1*1-AluORF10*1), followed by the theoretical ancestral haplotype with no Alu insertions (AluDPB2*1-AluDQA2*1-AluDQA1*1-AluDRB1*1-AluORF10*1) and finally the haplotype with the AluDPB2 and AluDQA1 Alu insertions (AluDPB2*2-AluDQA2*1-AluDQA1*2-AluDRB1*1-AluORF10*1) with frequencies of 19.2%, 15% and 12.9%, respectively. In the phylogenetic analyses, our results indicate that the Ngazidja people are genetically differentiated from the other populations of the analysis; we found also a new haplotype that can be probably characteristic of Sub-Saharans and finally confirm the usefulness of these markers as genetic and evolutionary tools for studying genetic variations among populations of different origins.

Comprehensive Analysis of CD4+ T Cell Response Cross-Reactive to SARS-CoV-2 Antigens at the Single Allele Level of HLA Class II

Common human coronaviruses have been circulating undiagnosed worldwide. These common human coronaviruses share partial sequence homology with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2); therefore, T cells specific to human coronaviruses are also cross-reactive with SARS-CoV-2 antigens. Herein, we defined CD4⁺ T cell responses that were cross-reactive with SARS-CoV-2 antigens in blood collected in 2016–2018 from healthy donors at the single allele level using artificial antigen-presenting cells (aAPC) expressing a single HLA class II allotype. We assessed the allotype-restricted responses in the 42 individuals using the aAPCs matched 22 HLA-DR alleles, 19 HLA-DQ alleles, and 13 HLA-DP alleles. The response restricted by the HLA-DR locus showed the highest magnitude, and that by HLA-DP locus was higher than that by HLA-DQ locus. Since two alleles of HLA-DR, -DQ, and -DP loci are expressed co-dominantly in an individual, six different HLA class II allotypes can be used to the cross-reactive T cell response. Of the 16 individuals who showed a dominant T cell response, five, one, and ten showed a dominant response by a single allotype of HLA-DR, -DQ, and -DP, respectively. The single allotype-restricted T cells responded to only one antigen in the five individuals and all the spike, membrane, and nucleocapsid proteins in the six individuals. In individuals heterozygous for the HLA-DPA and HLA-DPB loci, four combinations of HLA-DP can be expressed, but only one combination showed a dominant response. These findings demonstrate that cross-reactive T cells to SARS-CoV-2 respond with single-allotype dominance.

Allelic and haplotypic HLA diversity in indigenous Malaysian populations explored using Next Generation Sequencing

The heterogenous population of Malaysia includes more than 50 indigenous groups, and characterizing their HLA diversity would not only provide insights to their ancestry, but also on the effects of natural selection on their genome. We utilized hybridization-based sequence capture and short-read sequencing on the HLA region of 172 individuals representing seven indigenous groups in Malaysia (Jehai, Kintaq, Temiar, Mah Meri, Seletar, Temuan, Bidayuh). Allele and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DQA1, -DQB1, -DPA1, and -DPB1 revealed several ancestry-informative markers. Using SNP-based heterozygosity and pairwise Fst, we observed signals of natural selection, particularly in HLA-A, -C and -DPB1 genes. Consequently, we showed the impact of natural selection on phylogenetic inference using HLA and non-HLA SNPs. We demonstrate the utility of Next Generation Sequencing for generating unambiguous, high-throughput, high-resolution HLA data that adds to our knowledge of HLA diversity and natural selection in indigenous minority groups.

Self-mediated positive selection of T cells sets an obstacle to the recognition of nonself

Adaptive immune recognition is mediated by the binding of peptide-human leukocyte antigen complexes by T cells. Positive selection of T cells in the thymus is a fundamental step in the generation of a responding T cell repertoire: only those T cells survive that recognize human peptides presented on the surface of cortical thymic epithelial cells. We propose that while this step is essential for optimal immune function, the process results in a defective T cell repertoire because it is mediated by self-peptides. To test our hypothesis, we focused on amino acid motifs of peptides in contact with T cell receptors. We found that motifs rarely or not found in the human proteome are unlikely to be recognized by the immune system just like the ones that are not expressed in cortical thymic epithelial cells or not presented on their surface. Peptides carrying such motifs were especially dissimilar to human proteins. Importantly, we present our main findings on two independent T cell activation datasets and directly demonstrate the absence of naïve T cells in the repertoire of healthy individuals. We also show that T cell cross-reactivity is unable to compensate for the absence of positively selected T cells. Additionally, we show that the proposed mechanism could influence the risk for different infectious diseases. In sum, our results suggest a side effect of T cell positive selection, which could explain the nonresponsiveness to many nonself peptides and could improve the understanding of adaptive immune recognition.

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.

Nanovaccines against Animal Pathogens: The Latest Findings

Nowadays, safe and efficacious vaccines represent powerful and cost-effective tools for global health and economic growth. In the veterinary field, these are undoubtedly key tools for improving productivity and fighting zoonoses. However, cases of persistent infections, rapidly evolving pathogens having high variability or emerging/re-emerging pathogens for which no effective vaccines have been developed point out the continuing need for new vaccine alternatives to control outbreaks. Most licensed vaccines have been successfully used for many years now; however, they have intrinsic limitations, such as variable efficacy, adverse effects, and some shortcomings. More effective adjuvants and novel delivery systems may foster real vaccine effectiveness and timely implementation. Emerging vaccine technologies involving nanoparticles such as self-assembling proteins, virus-like particles, liposomes, virosomes, and polymeric nanoparticles offer novel, safe, and high-potential approaches to address many vaccine development-related challenges. Nanotechnology is accelerating the evolution of vaccines because nanomaterials having encapsulation ability and very advantageous properties due to their size and surface area serve as effective vehicles for antigen delivery and immunostimulatory agents. This review discusses the requirements for an effective, broad-coverage-elicited immune response, the main nanoplatforms for producing it, and the latest nanovaccine applications for fighting animal pathogens.

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.

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.

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.

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

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

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.

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

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

Functional biology in its natural context: A search for emergent simplicity

The immeasurable complexity at every level of biological organization creates a daunting task for understanding biological function. Here, we highlight the risks of stripping it away at the outset and discuss a possible path toward arriving at emergent simplicity of understanding while still embracing the ever-changing complexity of biotic interactions that we see in nature.

An Immunoinformatics Approach for SARS-CoV-2 in Latam Populations and Multi-Epitope Vaccine Candidate Directed towards the World's Population

The coronavirus pandemic is a major public health crisis affecting global health systems with dire socioeconomic consequences, especially in vulnerable regions such as Latin America (LATAM). There is an urgent need for a vaccine to help control contagion, reduce mortality and alleviate social costs. In this study, we propose a rational multi-epitope candidate vaccine against SARS-CoV-2. Using bioinformatics, we constructed a library of potential vaccine peptides, based on the affinity of the most common major human histocompatibility complex (HLA) I and II molecules in the LATAM population to predict immunological complexes among antigenic, non-toxic and non-allergenic peptides extracted from the conserved regions of 92 proteomes. Although HLA-C, had the greatest antigenic peptide capacity from SARS-CoV-2, HLA-B and HLA-A, could be more relevant based on COVID-19 risk of infection in LATAM countries. We also used three-dimensional structures of SARS-CoV-2 proteins to identify potential regions for antibody production. The best HLA-I and II predictions (with increased coverage in common alleles and regions evoking B lymphocyte responses) were grouped into an optimized final multi-epitope construct containing the adjuvants Beta defensin-3, TpD, and PADRE, which are recognized for invoking a safe and specific immune response. Finally, we used Molecular Dynamics to identify the multi-epitope construct which may be a stable target for TLR-4/MD-2. This would prove to be safe and provide the physicochemical requirements for conducting experimental tests around the world.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

Comprehensive Analysis of CD4+ T Cell Responses to CMV pp65 Antigen Restricted by Single HLA-DR, -DQ, and -DP Allotype Within an Individual

Within an individual, six different HLA class II heterodimers are expressed co-dominantly by two alleles of HLA-DR, -DQ, and -DP loci. However, it remained unclear which HLA allotypes were used in T cell responses to a given antigen. For the measurement of the CD4⁺ T cell responses restricted by a single HLA allotype, we established a panel of artificial antigen-presenting cells (aAPCs) expressing each single HLA allele of 20 HLA-DRB1, 16 HLA-DQ, and 13 HLA-DP alleles. CD4⁺ T cell responses to cytomegalovirus (CMV) pp65 restricted by single HLA class II allotype defined in 45 healthy donors. The average magnitude of CD4⁺ T cell responses by HLA-DR allotypes was higher than HLA-DQ and HLA-DP allotypes. CD4⁺ T cell responses by DRA*01:01/DRB1*04:06, DQA1*01:02/DQB1*06:02, DPA1*02:02/DPB1*05:01 were higher among the other alleles in each HLA-DR, -DQ, and -DP locus. Interestingly, the frequencies of HLA-DR alleles and the positivity of specific allotypes showed an inverse correlation. One allotype within individuals is dominantly used in CD4⁺ T cell response in 49% of donors, and two allotypes showed that in 7% of donors, and any positive response was detected in 44% of donors. Even if one individual had several dominant alleles, CD4⁺ T cell responses tended to be restricted by only one of them. Furthermore, CD8⁺ and CD4⁺ T cell responses by HLA class I and class II were correlated. Our results demonstrate that the CD4⁺ T cell preferentially use a few dominant HLA class II allotypes within individuals, similar to CD8⁺ T cell response to CMV pp65.

From Chickens to Humans: The Importance of Peptide Repertoires for MHC Class I Alleles

In humans, killer immunoglobulin-like receptors (KIRs), expressed on natural killer (NK) and thymus-derived (T) cells, and their ligands, primarily the classical class I molecules of the major histocompatibility complex (MHC) expressed on nearly all cells, are both polymorphic. The variation of this receptor-ligand interaction, based on which alleles have been inherited, is known to play crucial roles in resistance to infectious disease, autoimmunity, and reproduction in humans. However, not all the variation in response is inherited, since KIR binding can be affected by a portion of the peptide bound to the class I molecules, with the particular peptide presented affecting the NK response. The extent to which the large multigene family of chicken immunoglobulin-like receptors (ChIRs) is involved in functions similar to KIRs is suspected but not proven. However, much is understood about the two MHC-I molecules encoded in the chicken MHC. The BF2 molecule is expressed at a high level and is thought to be the predominant ligand of cytotoxic T lymphocytes (CTLs), while the BF1 molecule is expressed at a much lower level if at all and is thought to be primarily a ligand for NK cells. Recently, a hierarchy of BF2 alleles with a suite of correlated properties has been defined, from those expressed at a high level on the cell surface but with a narrow range of bound peptides to those expressed at a lower level on the cell surface but with a very wide repertoire of bound peptides. Interestingly, there is a similar hierarchy for human class I alleles, although the hierarchy is not as wide. It is a question whether KIRs and ChIRs recognize class I molecules with bound peptide in a similar way, and whether fastidious to promiscuous hierarchy of class I molecules affect both T and NK cell function. Such effects might be different from those predicted by the similarities of peptide-binding based on peptide motifs, as enshrined in the idea of supertypes. Since the size of peptide repertoire can be very different for alleles with similar peptide motifs from the same supertype, the relative importance of these two properties may be testable.

Evidence for height and immune function trade-offs among preadolescents in a high pathogen population

BACKGROUND

In an energy-limited environment, caloric investments in one characteristic should trade-off with investments in other characteristics. In high pathogen ecologies, biasing energy allocation towards immune function over growth would be predicted, given strong selective pressures against early-life mortality.

METHODOLOGY

In the present study, we use flow cytometry to examine trade-offs between adaptive immune function (T cell subsets, B cells), innate immune function (natural killer cells), adaptive to innate ratio and height-for-age z scores (HAZ) among young children (N = 344; aged 2 months-8 years) in the Bolivian Amazon, using maternal BMI and child weight-for-height z scores (WHZ) as proxies for energetic status.

RESULTS

Markers of adaptive immune function negatively associate with child HAZ, a pattern most significant in preadolescents (3+ years). In children under three, maternal BMI appears to buffer immune and HAZ associations, while child energetic status (WHZ) moderates relationships in an unexpected direction: HAZ and immune associations are greater in preadolescents with higher WHZ. Children with low WHZ maintain similar levels of adaptive immune function, but are shorter compared to high WHZ peers.

CONCLUSIONS

Reduced investment in growth in favor of immunity may be necessary for survival in high pathogen contexts, even under energetic constraints. Further, genetic and environmental factors are important considerations for understanding variation in height within this population. These findings prompt consideration of whether there may be a threshold of investment into adaptive immunity required for survival in high pathogen environments, and thus question the universal relevance of height as a marker of health.

LAY SUMMARY

Adaptive immune function is negatively associated with child height in this high pathogen environment. Further, low weight-for-height children are shorter but maintain similar immune levels. Findings question the relevance of height as a universal health marker, given that costs and benefits of height versus immunity may be calibrated to local ecology.

The complex pattern of genetic associations of leprosy with HLA class I and class II alleles can be reduced to four amino acid positions

Leprosy is a chronic disease caused by Mycobacterium leprae. Worldwide, more than 200,000 new patients are affected by leprosy annually, making it the second most common mycobacterial disease after tuberculosis. The MHC/HLA region has been consistently identified as carrying major leprosy susceptibility variants in different populations at times with inconsistent results. To establish the unambiguous molecular identity of classical HLA class I and class II leprosy susceptibility factors, we applied next-generation sequencing to genotype with high-resolution 11 HLA class I and class II genes in 1,155 individuals from a Vietnamese leprosy case-control sample. HLA alleles belonging to an extended haplotype from HLA-A to HLA-DPB1 were associated with risk to leprosy. This susceptibility signal could be reduced to the HLA-DRB1*10:01~ HLA-DQA1*01:05 alleles which were in complete linkage disequilibrium (LD). In addition, haplotypes containing HLA-DRB3~ HLA-DRB1*12:02 and HLA-C*07:06~ HLA-B*44:03~ HLA-DRB1*07:01 alleles were found as two independent protective factors for leprosy. Moreover, we replicated the previously associated HLA-DRB1*15:01 as leprosy risk factor and HLA-DRB1*04:05~HLA-DQA1*03:03 as protective alleles. When we narrowed the analysis to the single amino acid level, we found that the associations of the HLA alleles were largely captured by four independent amino acids at HLA-DRβ1 positions 57 (D) and 13 (F), HLA-B position 63 (E) and HLA-A position 19 (K). Hence, analyses at the amino acid level circumvented the ambiguity caused by strong LD of leprosy susceptibility HLA alleles and identified four distinct leprosy susceptibility factors.

Despite global efforts to eliminate leprosy over the past 25 years, more than 200,000 new cases are reported annually, and leprosy still represents a major public health problem in endemic regions. Leprosy presents a strong link with the host genetic background. The most significant susceptibility factors are located in the MHC region and likely involve classical HLA genes. However, the molecular identity of the HLA class I/II-leprosy risk factor(s) has been a matter of longstanding scientific dispute. By conducting a comprehensive sequenced-based analysis of HLA class I and class II genes, we are able to provide a unifying view of the complex relationship of leprosy susceptibility and HLA alleles. In addition, we show that four amino acid polymorphisms in HLA-DRβ1, HLA-B and HLA-A are sufficient to explain the majority of leprosy-HLA associations which opens the way for select protein-HLA peptide binding studies.

Binding affinities of 438 HLA proteins to complete proteomes of seven pandemic viruses and distributions of strongest and weakest HLA peptide binders in populations worldwide

We report detailed peptide‐binding affinities between 438 HLA Class I and Class II proteins and complete proteomes of seven pandemic human viruses, including coronaviruses, influenza viruses and HIV‐1. We contrast these affinities with HLA allele frequencies across hundreds of human populations worldwide. Statistical modelling shows that peptide‐binding affinities classified into four distinct categories depend on the HLA locus but that the type of virus is only a weak predictor, except in the case of HIV‐1. Among the strong HLA binders (IC₅₀ ≤ 50), we uncovered 16 alleles (the top ones being A*02:02, B*15:03 and DRB1*01:02) binding more than 1% of peptides derived from all viruses, 9 (top ones including HLA‐A*68:01, B*15:25, C*03:02 and DRB1*07:01) binding all viruses except HIV‐1, and 15 (top ones A*02:01 and C*14:02) only binding coronaviruses. The frequencies of strongest and weakest HLA peptide binders differ significantly among populations from different geographic regions. In particular, Indigenous peoples of America show both higher frequencies of strongest and lower frequencies of weakest HLA binders. As many HLA proteins are found to be strong binders of peptides derived from distinct viral families, and are hence promiscuous (or generalist), we discuss this result in relation to possible signatures of natural selection on HLA promiscuous alleles due to past pathogenic infections. Our findings are highly relevant for both evolutionary genetics and the development of vaccine therapies. However they should not lead to forget that individual resistance and vulnerability to diseases go beyond the sole HLA allelic affinity and depend on multiple, complex and often unknown biological, environmental and other variables.

Origins of the RAG Transposome and the MHC

How innate immunity gave rise to adaptive immunity in vertebrates remains unknown. We propose an evolutionary scenario beginning with pathogen-associated molecular pattern(s) (PAMPs) being presented by molecule(s) on one cell to specific receptor(s) on other cells, much like MHC molecules and T cell receptors (TCRs). In this model, mutations in MHC-like molecule(s) that bound new PAMP(s) would not be recognized by original TCR-like molecule(s), and new MHC-like gene(s) would be lost by neutral drift. Integrating recombination activating gene (RAG) transposon(s) in a TCR-like gene would result in greater recognition diversity, with new MHC-like variants recognized and selected, along with a new RAG/TCR-like system. MHC genes would be selected to present many peptides, through multigene families, allelic polymorphism, and peptide-binding promiscuity.

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

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

Diversity of HLA Class I and Class II blocks and conserved extended haplotypes in Lacandon Mayans

Here we studied HLA blocks and haplotypes in a group of 218 Lacandon Maya Native American using a high-resolution next generation sequencing (NGS) method. We assessed the genetic diversity of HLA class I and class II in this population, and determined the most probable ancestry of Lacandon Maya HLA class I and class II haplotypes. Importantly, this Native American group showed a high degree of both HLA homozygosity and linkage disequilibrium across the HLA region and also lower class II HLA allelic diversity than most previously reported populations (including other Native American groups). Distinctive alleles present in the Lacandon population include HLA-A*24:14 and HLA-B*40:08. Furthermore, in Lacandons we observed a high frequency of haplotypes containing the allele HLA-DRB1*04:11, a relatively frequent allele in comparison with other neighboring indigenous groups. The specific demographic history of the Lacandon population including inbreeding, as well as pathogen selection, may have elevated the frequencies of a small number of HLA class II alleles and DNA blocks. To assess the possible role of different selective pressures in determining Native American HLA diversity, we evaluated the relationship between genetic diversity at HLA-A, HLA-B and HLA-DRB1 and pathogen richness for a global dataset and for Native American populations alone. In keeping with previous studies of such relationships we included distance from Africa as a covariate. After correction for multiple comparisons we did not find any significant relationship between pathogen diversity and HLA genetic diversity (as measured by polymorphism information content) in either our global dataset or the Native American subset of the dataset. We found the expected negative relationship between genetic diversity and distance from Africa in the global dataset, but no relationship between HLA genetic diversity and distance from Africa when Native American populations were considered alone.

Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds

Pathogen communities can vary substantially between geographical regions due to different environmental conditions. However, little is known about how host immune systems respond to environmental variation across macro-ecological and evolutionary scales. Here, we select 37 species of songbird that inhabit diverse environments, including African and Palaearctic residents and Afro-Palaearctic migrants, to address how climate and habitat have influenced the evolution of key immune genes, the major histocompatibility complex class I (MHC-I). Resident species living in wetter regions, especially in Africa, had higher MHC-I diversity than species living in drier regions, irrespective of the habitats they occupy. By contrast, no relationship was found between MHC-I diversity and precipitation in migrants. Our results suggest that the immune system of birds has evolved greater pathogen recognition in wetter tropical regions. Furthermore, evolving transcontinental migration appears to have enabled species to escape wet, pathogen-rich areas at key periods of the year, relaxing selection for diversity in immune genes and potentially reducing immune system costs.

Combinatorial mutagenesis of rapidly evolving residues yields super-restrictor antiviral proteins

Antagonistic interactions drive host–virus evolutionary arms races, which often manifest as recurrent amino acid changes (i.e., positive selection) at their protein–protein interaction interfaces. Here, we investigated whether combinatorial mutagenesis of positions under positive selection in a host antiviral protein could enhance its restrictive properties. We tested approximately 700 variants of human MxA, generated by combinatorial mutagenesis, for their ability to restrict Thogotovirus (THOV). We identified MxA super-restrictors with increased binding to the THOV nucleoprotein (NP) target protein and 10-fold higher anti-THOV restriction relative to wild-type human MxA, the most potent naturally occurring anti-THOV restrictor identified. Our findings reveal a means to elicit super-restrictor antiviral proteins by leveraging signatures of positive selection. Although some MxA super-restrictors of THOV were impaired in their restriction of H5N1 influenza A virus (IAV), other super-restrictor variants increased THOV restriction without impairment of IAV restriction. Thus, broadly acting antiviral proteins such as MxA mitigate breadth-versus-specificity trade-offs that could otherwise constrain their adaptive landscape.

Antagonistic interactions drive host–virus evolutionary arms-races, often manifesting as recurrent amino acid changes at their protein–protein interaction interfaces. This study shows that evolution-guided combinatorial mutagenesis of a host antiviral protein enhances its restrictive properties, revealing constraints that shape breadth-versus-specificity trade-offs in antiviral proteins.

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.