MHC heterozygosity confers a selective advantage against multiple-strain infections

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

A comparative evaluation of the role of olfaction in attachment

Olfactory information plays an important role in the attachment and bonding processes for both humans and non-human animals. Odor cues obtained via individual body odor facilitate attachment and bonding processes across species with regard to both mate selection and mother-infant bonding. The purpose of the present paper is to summarize the role of odor as it pertains to bond formation and maintenance in the mother-infant bond for human infants and non-human animal infants, and for mate selection among human adults and non-human animals. We then synthesize this summary with literature on attachment and existing evidence for the relationships between olfaction and attachment processes. Finally, we suggest avenues for areas of future research.

MHC heterozygosity limits T cell receptor variability in CD4 T cells

αβ T cell receptor (TCR) V(D)J genes code for billions of TCR combinations. However, only some appear on peripheral T cells in any individual because, to mature, thymocytes must react with low affinity but not high affinity with thymus expressed major histocompatibility (MHC)/peptides. MHC proteins are very polymorphic. Different alleles bind different peptides. Therefore, any individual might express many different MHC alleles to ensure that some peptides from an invader are bound to MHC and activate T cells. However, most individuals express limited numbers of MHC alleles. To explore this, we compared the TCR repertoires of naïve CD4 T cells in mice expressing one or two MHC alleles. Unexpectedly, the TCRs in heterozygotes were less diverse that those in the sum of their MHC homozygous relatives. Our results suggest that thymus negative selection cancels out the advantages of increased thymic positive selection in the MHC heterozygotes.

The DRB3 gene of the bovine major histocompatibility complex: discovery, diversity and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine Major Histocompatibility Complex (MHC), also known as the Bovine Leucocyte Antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

Organisation and evolution of the major histocompatibility complex class I genes in cetaceans

A quarter of marine mammals are at risk of extinction, with disease and poor habitat quality contributing to population decline. Investigation of the Major Histocompatibility Complex (MHC) provides insight into species' capacity to respond to immune and environmental challenges. The eighteen available cetacean chromosome level genomes were used to annotate MHC Class I loci, and to reconstruct the phylogenetic relationship of the described loci. The highest number of loci was observed in the striped dolphin (Stenella coeruleoalba), while the least was observed in the pygmy sperm whale (Kogia breviceps) and rough toothed dolphin (Steno bredanensis). Of the species studied, Mysticetes had the most pseudogenes. Evolutionarily, MHC Class I diverged before the speciation of cetaceans. Yet, locus one was genomically and phylogenetically similar in many species, persisting over evolutionary time. This characterisation of MHC Class I in cetaceans lays the groundwork for future population genetics and MHC expression studies.

An immunogenetic basis for lung cancer risk

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response

BACKGROUND

The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer.

METHODS

RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles.

RESULTS

The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls.

CONCLUSIONS

The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.

Genomic heterozygosity is associated with a lower risk of osteoarthritis

BACKGROUND

Genomic heterozygosity has been shown to confer a health advantage in humans and play a protective role in complex diseases. Given osteoarthritis (OA) is a highly polygenic disease, we set out to determine if an association exists between OA and genomic heterozygosity.

RESULTS

End-stage knee and hip OA patients and healthy controls were recruited from the Newfoundland and Labrador (NL) population. The Arthritis Research UK Osteoarthritis Genetics (arcOGEN) consortium database was utilized as a replication cohort. DNA was extracted from blood samples and genotyped. Individual rates of observed heterozygosity (HetRate) and heterozygosity excess (HetExcess) relative to the expected were mathematically derived, and standardized to a z-score. Logistic regression modeling was used to examine the association between OA and HetRate or HetExcess. A total of 559 knee and hip OA patients (mean age 66.5 years, body mass index (BMI) 33.7 kg/m2, and 55% females) and 118 healthy controls (mean age 56.4 years, BMI 29.5 kg/m2, and 59% female) were included in the NL cohort analysis. We found that OA had an inverse relationship with HetRate and HetExcess with odds ratios of 0.64 (95% CI: 0.45-0.91) and 0.65 (95% CI: 0.45-0.93) per standard deviation (SD), respectively. The arcOGEN data included 2,019 end-stage knee and hip OA patients and 2,029 healthy controls, validating our findings with HetRate and HetExcess odds ratios of 0.60 (95% CI: 0.56-0.64) and 0.44 (95% CI: 0.40-0.47) per SD, respectively.

CONCLUSIONS

Our results are the first to clearly show evidence, from two separate cohorts, that reduced genomic heterozygosity confers a risk for the future development of OA.

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.

Population-wide genetic analysis of Ovar-DQA1 and DQA2 loci across sheep breeds in India revealed their evolutionary importance and fitness of sheep in a tropical climate

Genetic variability at the major histocompatibility complex (MHC) is important in any species due to significant role played by MHC for antigen presentation. DQA locus has not been studied for its genetic variability across sheep population in India. In the present study, MHC of sheep at DQA1 and DQA2 loci were evaluated across 17 Indian sheep breeds. Results revealed high degree of heterozygosity (10.34% to 100% for DQA1 and 37.39 to 100% for DQA2). 18 DQA1 alleles and 22 DQA2 alleles were isolated in different breeds. Nucleotide content for DQA region revealed richness of AT content (54.85% for DQA1 and 53.89% for DQA2). DQA1 and DQA2 sequences clustered independently. We could see evidence of divergence of DQA as DQA1 and DQA2 across sheep breeds. Wu-Kabat variability index revealed vast genetic variation across DQA1 and DQA2, specifically at peptide binding sites (PBS) that consisted 21 residues for DQA1 and 17 residues for DQA2. Evolutionary analysis revealed the presence of positive and balancing selection for DQA1 locus, however DQA2 was under purifying selection across sheep breeds. Higher heterozygosity and large diversity at both loci especially at PBS indicated the fitness of the sheep population for evading pathogens and adapt to the harsh tropical climate.

Integration of the immune memory into the pathogen-driven MHC polymorphism hypothesis

Major histocompatibility complex (MHC) genes (referred to as human leukocyte antigen or HLA in humans) are a key component of vertebrate immune systems, coding for proteins which present antigens to T-cells. These genes are outstanding in their degree of polymorphism, with important consequences for human and animal health. The polymorphism is thought to arise from selection pressures imposed by pathogens on MHC allomorphs, which differ in their antigen-binding capacity. However, the existing theory has not considered MHC selection in relation to the formation of immune memory. In this paper, we argue that this omission limits our understanding of the evolution of MHC polymorphism and its role in disease. We review recent evidence that has emerged from the massive research effort related to the SARS-CoV-2 pandemics, and which provides new evidence for the role of MHC in shaping immune memory. We then discuss why the inclusion of immune memory within the existing theory may have non-trivial consequence for our understanding of the evolution of MHC polymorphism. Finally, we will argue that neglecting immune memory hinders our interpretation of empirical findings, and postulate that future studies focusing on pathogen-driven MHC selection would benefit from stratifying the available data according to the history of infection (and vaccination, if relevant).

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

Genomic inbreeding estimation, runs of homozygosity, and heterozygosity-enriched regions uncover signals of selection in the Quarter Horse racing line

With the advent of genomics, significant progress has been made in the genetic improvement of livestock species, particularly through increased accuracy in predicting breeding values for selecting superior animals and the possibility of performing a high-resolution genetic scan throughout the genome of an individual. The main objectives of this study were to estimate the individual genomic inbreeding coefficient based on runs of homozygosity (FROH ), to identify and characterize runs of homozygosity and heterozygosity (ROH and ROHet, respectively; length and distribution) throughout the genome, and to map selection signatures in relevant chromosomal regions in the Quarter Horse racing line. A total of 336 animals registered with the Brazilian Association of Quarter Horse Breeders (ABQM) were genotyped. One hundred and twelve animals were genotyped using the Equine SNP50 BeadChip (Illumina, USA), with 54,602 single nucleotide polymorphisms (SNPs; 54K). The remaining 224 samples were genotyped using the Equine SNP70 BeadChip (Illumina, USA) with 65,157 SNPs (65K). To ensure data quality, we excluded animals with a call rate below 0.9. We also excluded SNPs located on non-autosomal chromosomes, as well as those with a call rate below 0.9 or a p-value below 1 × 10-5 for Hardy-Weinberg equilibrium. The results indicate moderate to high genomic inbreeding, with 46,594 ROH and 16,101 ROHet detected. In total, 30 and 14 candidate genes overlap with ROH and ROHet regions, respectively. The ROH islands showed genes linked to crucial biological processes, such as cell differentiation (CTBP1, WNT5B, and TMEM120B), regulation of glucose metabolic process (MAEA and NKX1-1), heme transport (PGRMC2), and negative regulation of calcium ion import (VDAC1). In ROHet, the islands showed genes related to respiratory capacity (OR7D19, OR7D4G, OR7D4E, and OR7D4J) and muscle repair (EGFR and BCL9). These findings could aid in selecting animals with greater regenerative capacity and developing treatments for muscle disorders in the QH breed. This study serves as a foundation for future research on equine breeds. It can contribute to developing reproductive strategies in animal breeding programs to improve and preserve the Quarter Horse breed.

Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

Heterozygote advantage at HLA class I and II loci and reduced risk of colorectal cancer

Objective

Reduced diversity at Human Leukocyte Antigen (HLA) loci may adversely affect the host's ability to recognize tumor neoantigens and subsequently increase disease burden. We hypothesized that increased heterozygosity at HLA loci is associated with a reduced risk of developing colorectal cancer (CRC).

Methods

We imputed HLA class I and II four-digit alleles using genotype data from a population-based study of 5,406 cases and 4,635 controls from the Molecular Epidemiology of Colorectal Cancer Study (MECC). Heterozygosity at each HLA locus and the number of heterozygous genotypes at HLA class -I (A, B, and C) and HLA class -II loci (DQB1, DRB1, and DPB1) were quantified. Logistic regression analysis was used to estimate the risk of CRC associated with HLA heterozygosity. Individuals with homozygous genotypes for all loci served as the reference category, and the analyses were adjusted for sex, age, genotyping platform, and ancestry. Further, we investigated associations between HLA diversity and tumor-associated T cell repertoire features, as measured by tumor infiltrating lymphocytes (TILs; N=2,839) and immunosequencing (N=2,357).

Results

Individuals with all heterozygous genotypes at all three class I genes had a reduced odds of CRC (OR: 0.74; 95% CI: 0.56-0.97, p= 0.031). A similar association was observed for class II loci, with an OR of 0.75 (95% CI: 0.60-0.95, p= 0.016). For class-I and class-II combined, individuals with all heterozygous genotypes had significantly lower odds of developing CRC (OR: 0.66, 95% CI: 0.49-0.87, p= 0.004) than those with 0 or one heterozygous genotype. HLA class I and/or II diversity was associated with higher T cell receptor (TCR) abundance and lower TCR clonality, but results were not statistically significant.

Conclusion

Our findings support a heterozygote advantage for the HLA class-I and -II loci, indicating an important role for HLA genetic variability in the etiology of CRC.

CEGA: a method for inferring natural selection by comparative population genomic analysis across species

We developed maximum likelihood method for detecting positive selection or balancing selection using multilocus or genomic polymorphism and divergence data from two species. The method is especially useful for investigating natural selection in noncoding regions. Simulations demonstrate that the method outperforms existing methods in detecting both positive and balancing selection. We apply the method to population genomic data from human and chimpanzee. The list of genes identified under selection in the noncoding regions is prominently enriched in pathways related to the brain and nervous system. Therefore, our method will serve as a useful tool for comparative population genomic analysis.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

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

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

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

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

Facial attractiveness is only weakly linked to genome-wide heterozygosity

Introduction

It has been frequently suggested that overall genomic heterozygosity and, particularly, heterozygosity of loci on the so-called major histocompatibility complex (MHC), which is responsible for the recognition of foreign substances/ pathogens and the recognition of self and non-self, is associated with better health and better resistance to infections and parasites. It has further been speculated that such a potentially beneficial heterozygosity can be detected through body odor and facial attractiveness.

Methods

On the basis of genome wide SNP data (713,014 SNPs) of participants from the Wisconsin Longitudinal Study, we therefore investigated whether homozygosity either on the MHC (measured as inbreeding coefficient) or genome-wide (measured as runs of homozygosity and inbreeding coefficient) is associated with rated facial attractiveness.

Results

Although we found that the genome-wide average length of homozygous segments and the genome-wide inbreeding coefficient are significantly negatively associated with some measures of facial attractiveness, if corrected for multiple testing, any significant association was no longer formally significant after correction. In addition, the variance in facial attractiveness explained by the genome wide homozygosity is very low (<0.15%). We did not find any significant association between the inbreeding coefficient on the MHC and facial attractiveness.

Discussion

We only find a weak association of genome- wide heterozygosity and facial attractiveness.

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

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

Heterogeneity in disease resistance and the impact of antibiotics in the US

We hypothesize that the impact of antibiotics is moderated by a population's inherent (genetic) resistance to infectious disease. Using the introduction of sulfa drugs in 1937, we show that US states that are more genetically susceptible to infectious disease saw larger declines in their bacterial mortality rates following the introduction of sulfa drugs in 1937. This suggests area-level genetic endowments of disease resistance and the discovery of medical technologies have acted as substitutes in determining levels of health across the US. We also document immediate effects of sulfa drug exposure to the age of the workforce and cumulative effects on educational attainment for cohorts exposed to sulfa drugs in early life.

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

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

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

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

A single genomic region involving a putative chromosome rearrangement in flat oyster (Ostrea edulis) is associated with differential host resilience to the parasite Bonamia ostreae

European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long-term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F ST > 0.25) unlike neutral loci (F ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE-C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE-C8. A range of relevant immune functions were specifically enriched among genes annotated on OE-C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker-assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.

HLA homozygosity is associated with Non-Hodgkin lymphoma

The "heterozygote advantage" hypothesis has been postulated regarding the role of human leukocyte antigen (HLA) in non-Hodgkin lymphoma (NHL), where homozygous loci are associated with an increased risk of disease. In this retrospective study, we analyzed the HLA homozygosity of 3789 patients with aplastic anemia (AA), acute lymphocytic leukemia (ALL), acute myeloblastic leukemia (AML), chronic lymphocytic leukemia (CLL), chronic myeloid leukemia (CML), myelodysplastic syndrome (MDS), multiple myeloma (MM), and non-Hodgkin lymphoma (NHL) at HLA-A, B, C, DRB1 and DQB1 loci compared to 169,964 normal controls. HLA homozygosity at one or more loci was only associated with an increased risk in NHL patients (OR = 1.28, 95% CI [1.09, 1.50], p = 0.002). This association was not seen in any of the other hematologic diseases. Homozygosity at HLA-A alone, HLA-B + C only, and HLA-DRB1 + DQB1 only was also significantly associated with NHL. Finally, we observed a 17% increased risk of NHL with each additional homozygous locus (OR per locus = 1.17, 95% CI [1.08, 1.25], p trend = 2.4 × 10^-5). These results suggest that reduction of HLA diversity could predispose individuals to an increased risk of developing NHL.

Susceptibility and Severity of COVID-19 Are Both Associated With Lower Overall Viral-Peptide Binding Repertoire of HLA Class I Molecules, Especially in Younger People

An important number of studies have been conducted on the potential association between human leukocyte antigen (HLA) genes and COVID-19 susceptibility and severity since the beginning of the pandemic. However, case-control and peptide-binding prediction methods tended to provide inconsistent conclusions on risk and protective HLA alleles, whereas some researchers suggested the importance of considering the overall capacity of an individual's HLA Class I molecules to present SARS-CoV-2-derived peptides. To close the gap between these approaches, we explored the distributions of HLA-A, -B, -C, and -DRB1 1st-field alleles in 142 Iranian patients with COVID-19 and 143 ethnically matched healthy controls, and applied in silico predictions of bound viral peptides for each individual's HLA molecules. Frequency comparison revealed the possible predisposing roles of HLA-A*03, B*35, and DRB1*16 alleles and the protective effect of HLA-A*32, B*58, B*55, and DRB1*14 alleles in the viral infection. None of these results remained significant after multiple testing corrections, except HLA-A*03, and no allele was associated with severity, either. Compared to peptide repertoires of individual HLA molecules that are more likely population-specific, the overall coverage of virus-derived peptides by one's HLA Class I molecules seemed to be a more prominent factor associated with both COVID-19 susceptibility and severity, which was independent of affinity index and threshold chosen, especially for people under 60 years old. Our results highlight the effect of the binding capacity of different HLA Class I molecules as a whole, and the more essential role of HLA-A compared to HLA-B and -C genes in immune responses against SARS-CoV-2 infection.

HLA Diversity in Saudi Population: High Frequency of Homozygous HLA Alleles and Haplotypes

Human leukocyte antigens (HLA) diversity has a tremendous impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions. Here, we conducted a retrospective study of HLA class I and class II homozygosity at allelic and haplotype levels in unrelated individuals genotyped from 2012 to 2016 in a tertiary hospital in the capital of Saudi Arabia. Among 5,000 individuals, 2,773 individuals meet inclusion criteria and were retrospectively analyzed for HLA-A, -B, -C-DRB1, and -DQB1 homozygosity at allelic and haplotype levels. HLA molecular typing was performed using a commercial reverse sequence-specific oligonucleotide (rSSO) kit. We were able to identify 15 HLA-A, 20 HLA-B, 11 HLA-C, 13 HLA-DRB1, and five HLA-DQB1 homozygous alleles demonstrating a very low genetic diversity in the Saudi population. The highest homozygosity in HLA class I was found in locus C followed by A and B (20.3% > 16.1% > 15.5%; p < 0.001) where the most homozygote alleles were A*02 (9.2%), B*51 and B*50 (5.7% and 3.7%), and C*07, C*06, and C*15 (7.2%, 5.48%, and 3.3%) and in HLA class II, the highest homozygosity was found in locus DQB1 compared to DRB1 (31.71% > 19.2%; p < 0.001), with the most common homozygote alleles being DRB1*07 and DRB1*04 (5.33% and 4.2%) and DQB1*02, DQB1*06, and DQB1*03 (13.55%, 7.92%, and 7.64%). The frequency of finding an individual with one homozygote allele was (24.6%), two homozygote alleles (13.5%), three homozygote alleles (4.7%), four homozygote alleles (3.4%), and five alleles were (4.8%). The most frequent homozygote haplotypes are A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02. This study shows low diversity of both class I and II alleles and haplotypes in the Saudi population, which would have a significant impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

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.

HLA Zygosity Increases Risk of Hepatitis B Virus-Associated Hepatocellular Carcinoma

BACKGROUND

Diversity in the HLA genes might be associated with disease outcomes-the heterozygote advantage hypothesis. We tested this hypothesis in relation to hepatitis B virus (HBV)-associated hepatocellular carcinoma (HCC).

METHODS

We utilized DNA from > 10 000 Taiwanese individuals with current or past HBV infection to examine the association between HLA diversity and critical natural history steps in the progression from HBV infection to HCC. Individuals were classified as homozygotes at a given locus when imputed to carry the same 4-digit allele for the 2 HLA alleles at that locus.

RESULTS

Increase in number of homozygous HLA class II loci was associated with an increased risk of chronic HBV infection (Ptrend = 1.18 × 10-7). Among chronic HBV carriers, increase in number of homozygous HLA class II loci was also associated with an increased risk of HBV-associated HCC (Ptrend = .031). For individual HLA loci, HLA-DQB1 homozygosity was significantly associated with HCC risk (adjusted hazard ratio = 1.40; 95% confidence interval, 1.06-1.84). We also found that zygosity affects risk of HCC through its ability to affect viral control.

CONCLUSIONS

Homozygosity at HLA class II loci, particularly HLA-DQB1, is associated with a higher risk of HBV-associated HCC.

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.

Associations between MHC class II variation and phenotypic traits in a free-living sheep population

Pathogen-mediated selection (PMS) is thought to maintain the high level of allelic diversity observed in the major histocompatibility complex (MHC) class II genes. A comprehensive way to demonstrate contemporary selection is to examine associations between MHC variation and individual fitness. As individual fitness is hard to measure, many studies examine associations between MHC variation and phenotypic traits, including direct or indirect measures of adaptive immunity thought to contribute to fitness. Here, we tested associations between MHC class II variation and five phenotypic traits measured in free-living sheep captured in August: weight, strongyle faecal egg count, and plasma IgA, IgE and IgG immunoglobulin titres against the gastrointestinal nematode parasite Teladorsagia circumcincta. We found no association between MHC class II variation and weight or strongyle faecal egg count. We did, however, find associations between MHC class II variation and immunoglobulin levels which varied with isotype, age and sex. Our results suggest associations between MHC and phenotypic traits are more likely to be found for traits more closely associated with pathogen defence than integrative traits such as bodyweight and highlight the association between MHC variation and antibodies in wild populations.

Germline HLA-B evolutionary divergence influences the efficacy of immune checkpoint blockade therapy in gastrointestinal cancer

BACKGROUND

The human leukocyte antigen class I (HLA-I) genotype has been linked with differential immune responses to infectious disease and cancer. However, the clinical relevance of germline HLA-mediated immunity in gastrointestinal (GI) cancer remains elusive.

METHODS

This study retrospectively analyzed the genomic profiling data from 84 metastatic GI cancer patients treated with immune checkpoint blockade (ICB) recruited from Peking University Cancer Hospital (PUCH). A publicly available dataset from the Memorial Sloan Kettering (MSK) Cancer Center (MSK GI cohort) was employed as the validation cohort. For the PUCH cohort, we performed HLA genotyping by whole exome sequencing (WES) analysis on the peripheral blood samples from all patients. Tumor tissues from 76 patients were subjected to WES analysis and immune oncology-related RNA profiling. We studied the associations of two parameters of germline HLA as heterozygosity and evolutionary divergence (HED, a quantifiable measure of HLA-I evolution) with the clinical outcomes of patients in both cohorts.

RESULTS

Our data showed that neither HLA heterozygosity nor HED at the HLA-A/HLA-C locus correlated with the overall survival (OS) in the PUCH cohort. Interestingly, in both the PUCH and MSK GI cohorts, patients with high HLA-B HED showed a better OS compared with low HLA-B HED subgroup. Of note, a combinatorial biomarker of HLA-B HED and tumor mutational burden (TMB) may better stratify potential responders. Furthermore, patients with high HLA-B HED were characterized with a decreased prevalence of multiple driver gene mutations and an immune-inflamed phenotype.

CONCLUSIONS

Our results unveil how HLA-B evolutionary divergence influences the ICB response in patients with GI cancers, supporting its potential utility as a combinatorial biomarker together with TMB for patient stratification in the future.

Integrated analysis reveals prognostic value of HLA-I LOH in triple-negative breast cancer

Background

Triple-negative breast cancers (TNBCs), especially those non-immune-inflamed tumors, have a poor prognosis and limited therapies. Human leukocyte antigen (HLA)-I not only contributes to antitumor immune response and the phenotype of the tumor microenvironment, but also is a negative predictor of outcomes after immunotherapy. However, the importance of HLA functional status in TNBCs remains poorly understood.

Methods

Using the largest original multiomics datasets on TNBCs, we systematically characterized the HLA-Ⅰ status of TNBCs from the perspective of HLA-Ⅰ homogeneity and loss of heterozygosity (LOH). The prognostic significance of HLA-I status was measured. To explain the potential mechanism of prognostic value in HLA-Ⅰ status, the mutational signature, copy number alteration, neoantigen and intratumoral heterogeneity were measured. Furthermore, the correlation between HLA-Ⅰ functional status and the tumor immune microenvironment was analyzed.

Results

LOH and homogeneity in HLA-I accounted for 18% and 21% of TNBCs, respectively. HLA-I LOH instead of HLA-I homogeneity was an independent prognostic biomarker in TNBCs. In particular, for patients with non-immune-inflamed tumors, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH. Furthermore, integrated genomic and transcriptomic analysis showed that HLA-I LOH was accompanied by upregulated scores of mutational signature 3 and homologous recombination deficiency scores, which implied the failure of DNA double-strand break repair. Moreover, HLA-I LOH had higher mutation and neoantigen loads and more subclones than HLA-I non-LOH. These results indicated that although HLA-I LOH tumors with failure of DNA double-strand break repair were prone to produce neoantigens, their limited capacity for antigen presentation finally contributed to poor immune selection pressure.

Conclusion

Our study illustrates the genomic landscape of HLA-I functional status and stresses the prognostic significance of HLA-I LOH in TNBCs. For “cold” tumors in TNBCs, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH.

Preen gland microbiota covary with major histocompatibility complex genotype in a songbird

Pathogen-mediated selection at the major histocompatibility complex (MHC) is thought to promote MHC-based mate choice in vertebrates. Mounting evidence implicates odour in conveying MHC genotype, but the underlying mechanisms remain uncertain. MHC effects on odour may be mediated by odour-producing symbiotic microbes whose community structure is shaped by MHC genotype. In birds, preen oil is a primary source of body odour and similarity at MHC predicts similarity in preen oil composition. Hypothesizing that this relationship is mediated by symbiotic microbes, we characterized MHC genotype, preen gland microbial communities and preen oil chemistry of song sparrows (Melospiza melodia). Consistent with the microbial mediation hypothesis, pairwise similarity at MHC predicted similarity in preen gland microbiota. Counter to this hypothesis, overall microbial similarity did not predict chemical similarity of preen oil. However, permutation testing identified a maximally predictive set of microbial taxa that best reflect MHC genotype, and another set of taxa that best predict preen oil chemical composition. The relative strengths of relationships between MHC and microbes, microbes and preen oil, and MHC and preen oil suggest that MHC may affect host odour both directly and indirectly. Thus, birds may assess MHC genotypes based on both host-associated and microbially mediated odours.

Donor HLA Class 1 Evolutionary Divergence Is a Major Predictor of Liver Allograft Rejection : A Retrospective Cohort Study

BACKGROUND

The HLA evolutionary divergence (HED), a continuous metric quantifying the peptidic differences between 2 homologous HLA alleles, reflects the breadth of the immunopeptidome presented to T lymphocytes.

OBJECTIVE

To assess the potential effect of donor or recipient HED on liver transplant rejection.

DESIGN

Retrospective cohort study.

SETTING

Liver transplant units.

PATIENTS

1154 adults and 113 children who had a liver transplant between 2004 and 2018.

MEASUREMENTS

Liver biopsies were done 1, 2, 5, and 10 years after the transplant and in case of liver dysfunction. Donor-specific anti-HLA antibodies (DSAs) were measured in children at the time of biopsy. The HED was calculated using the physicochemical Grantham distance for class I (HLA-A or HLA-B) and class II (HLA-DRB1 or HLA-DQB1) alleles. The influence of HED on the incidence of liver lesions was analyzed through the inverse probability weighting approach based on covariate balancing, generalized propensity scores.

RESULTS

In adults, class I HED of the donor was associated with acute rejection (hazard ratio [HR], 1.09 [95% CI, 1.03 to 1.16]), chronic rejection (HR, 1.20 [CI, 1.10 to 1.31]), and ductopenia of 50% or more (HR, 1.33 [CI, 1.09 to 1.62]) but not with other histologic lesions. In children, class I HED of the donor was also associated with acute rejection (HR, 1.16 [CI, 1.03 to 1.30]) independent of the presence of DSAs. There was no effect of either donor class II HED or recipient class I or class II HED on the incidence of liver lesions in adults and children.

LIMITATION

The DSAs were measured only in children.

CONCLUSION

Class I HED of the donor predicts acute or chronic rejection of liver transplant. This novel and accessible prognostic marker could orientate donor selection and guide immunosuppression.

PRIMARY FUNDING SOURCE

Institut National de la Santé et de la Recherche Médicale.

Comparison of major histocompatibility complex-B variability in Sri Lankan indigenous chickens with five global chicken populations using MHC-B SNP panel

In the present study, we investigated the major histocompatibility complex (MHC)-B haplotypes diversity of Sri Lankan indigenous chickens from three different geographical sites consisting of highly mixed populations using 90 SNPs in the MHC-B region. A total of 48 haplotypes were identified. Those included 37 novel haplotypes and 11 previously identified 'standard' haplotypes. The MHC-linked marker, LEI0258, had 23 alleles showing less diversity than defined by MHC-B SNP haplotypes. Among those identified haplotypes, five standard haplotypes-BSNP-O02, BSNP-M01, BSNP-A04, BSNP-K03, BSNP-T04-were most commonly observed, suggesting past introgression of imported breeds. Comparison of the MHC-B haplotypes of Sri Lankan and four other global populations with previously defined haplotypes indicated the sharing of 23 standard haplotypes with common origins. Novel haplotypes are population-specific and not shared among the geographical boundaries. Backyard indigenous chickens are unselected, highly crossbred, and generally thrive under dynamic environmental conditions. Hence free-range production systems may be responsible for maintaining high diversity in the MHC-B region with novel haplotypes.

Haematological Indicators of Response to Erythropoietin Therapy in Chronic Renal Failure Patients on Haemodialysis: Impact of Angiotensin-Converting Enzyme rs4343 Gene Polymorphism

PURPOSE

This is the first cross-sectional study studying the changes in haematological indicators of the response to recombinant human erythropoietin (rHuEPO) therapy in chronic renal failure (CRF) patients on haemodialysis (HD) stratified according to ACE G2350A (rs4343) gene polymorphism.

DESIGN

An observational cross-sectional study.

SETTING

Nephrology department and Biochemistry and molecular biology department, faculty of medicine, Cairo University.

PATIENTS

A total of 256 CRF patients on HD for at least six months (162 male and 103 female) and 160 healthy subjects (122 male and 38 female) were recruited in the current study after signing a consent form. ACE G2350A (rs4343) Insertion/Deletion (I/D) was tested, the association between ACE G2350A (RS4343) gene polymorphisms and patients response to rHuEpo was evaluated.

RESULTS

ACE G2350A (rs4343) I/D was the most prevalent genotype, while I/I genotype was the lowest prevalent among patient or control subjects included in the study. D allele is the most prevalent allele, either among patients or the control group. Hemoglobin (Hb) level in patients with I/I and Deletion/Deletion (D/D) genotype was significantly higher compared to those with I/D genotype (P = 0.012 and P = 0.005, respectively). Serum iron in the I/D genotype was significantly higher than those with either I/I or D/D genotype (P = 0.045 and P = 0.018, respectively). Angiotensin-converting enzyme (ACE) content, total leukocytic count (TLC), and soluble erythropoietin receptor (sEpoR) were independent predictors of Hb level. The ACE gene, TLC, and serum iron were the independent factors that may affect the Haematocrit (Hct) level. ACE G2350A (rs4343) gene polymorphisms may affect the HD patient's responses to rHuEPOs.

CONCLUSION

In HD patients, screening for ACE G2350A (rs4343) gene polymorphisms before rHuEpo administration may help predict patient response.

Functional immunogenetic variation, rather than local adaptation, predicts ectoparasite infection intensity in a model fish species

Natural host populations differ in their susceptibility to infection by parasites, and these intrapopulation differences are still an incompletely understood component of host‐parasite dynamics. In this study, we used controlled infection experiments with wild‐caught guppies (Poecilia reticulata) and their ectoparasite Gyrodactylus turnbulli to investigate the roles of local adaptation and host genetic composition (immunogenetic and neutral) in explaining differences in susceptibility to infection. We found differences between our four study host populations that were consistent between two parasite source populations, with no indication of local adaptation by either host or parasite at two tested spatial scales. Greater values of host population genetic variability metrics broadly aligned with lower population mean infection intensity, with the best alignments associated with major histocompatibility complex (MHC) “supertypes”. Controlling for intrapopulation differences and potential inbreeding variance, we found a significant negative relationship between individual‐level functional MHC variability and infection: fish carrying more MHC supertypes experienced infections of lower severity, with limited evidence for supertype‐specific effects. We conclude that population‐level differences in host infection susceptibility probably reflect variation in parasite selective pressure and/or host evolutionary potential, underpinned by functional immunogenetic variation.

Major histocompatibility complex class II DR and DQ evolution and variation in wild capuchin monkey species (Cebinae)

The major histocompatibility complex (MHC) is an important gene complex contributing to adaptive immunity. Studies of platyrrhine MHC have focused on identifying experimental models of immune system function in the equivalent Human Leukocyte Antigen (HLA). These genes have thus been explored primarily in captive platyrrhine individuals from research colonies. However, investigations of standing MHC variation and evolution in wild populations are essential to understanding its role in immunity, sociality and ecology. Capuchins are a promising model group exhibiting the greatest habitat diversity, widest diet breadth and arguably the most social complexity among platyrrhines, together likely resulting in varied immunological challenges. We use high-throughput sequencing to characterize polymorphism in four Class II DR and DQ exons for the first time in seven capuchin species. We find evidence for at least three copies for DQ genes and at least five for DRB, with possible additional unrecovered diversity. Our data also reveal common genotypes that are inherited across our most widely sampled population, Cebus imitator in Sector Santa Rosa, Costa Rica. Notably, phylogenetic analyses reveal that platyrrhine DQA sequences form a monophyletic group to the exclusion of all Catarrhini sequences examined. This result is inconsistent with the trans-species hypothesis for MHC evolution across infraorders in Primates and provides further evidence for the independent origin of current MHC genetic diversity in Platyrrhini. Identical allele sharing across cebid species, and more rarely genera, however, does underscore the complexity of MHC gene evolution and the need for more comprehensive assessments of allelic diversity and genome structure.

Association Between Human Leukocyte Antigen Class I and II Diversity and Non-virus-associated Solid Tumors

Homozygosity at human leukocyte antigen (HLA) loci might lead to reduced immunosurveillance and increased disease risk, including cancers caused by infection or of hematopoietic origin. To investigate the association between HLA zygosity and risk of non-virus-associated solid tumors, we leveraged genome-wide association study (GWAS) data from over 28,000 individuals of European ancestry who participated in studies of 12 cancer sites (bladder, brain, breast, colon, endometrial, kidney, lung, ovary, pancreas, prostate, skin, and testis). Information on HLA zygosity was obtained by imputation; individuals were classified as homozygotes at a given locus when imputed to carry the same four-digit allele at that locus. We observed no evidence for an association between zygosity at six HLA loci and all cancers combined. Increase in number of homozygous at HLA class I loci, class II loci, or class I and II loci was also not associated with cancer overall (P trend = 0.28), with adjusted odds ratios (ORs) for risk-per-locus of 1.00 [95% confidence intervals (CIs) = 0.97, 1.03], 1.02 (0.99, 1.04), and 1.01 (0.99, 1.02), respectively. This study does not support a strong role for HLA zygosity on risk of non-virus-associated solid tumors.

HLA-A homozygosis is associated with susceptibility to COVID-19

The purpose of this single center retrospective study was to investigate the relationship between HLA and ABO polymorphisms and COVID-19 susceptibility and severity in kidney transplant recipients. It included 720 recipients who had COVID-19 and 1680 controls composed by recipients in follow-up who did not contact the transplantation center for COVID-19 symptoms, up to the moment of their inclusion in the study. HLA-A, -B, and -DRB1 allele groups and ABO frequencies were compared between recipients with COVID-19 (all cases, or separately mild/moderate and severe disease) and controls. The HLA association study was conducted in two case-control series and only associations that showed a p-value <0.05 in both series were considered. No HLA association regarding COVID-19 occurrence or severity met this criterion. Homozygosity at HLA-A locus was associated with COVID-19 susceptibility (odds ratio 1.4) but not severity. Blood groups A and O were associated with susceptibility and resistance to COVID-19, respectively. COVID-19 severity was associated only with older age and cardiac disease, in a multivariate analysis. We conclude that an influence of HLA on COVID-19 susceptibility is supported by the association with homozygosity at HLA-A locus but that there is no evidence for a role of any particular HLA-A, -B, or -DRB1 polymorphism. Thus, we suggest that what matters is the overall capability of an individual's HLA molecules to present SARS-CoV-2 peptides to T cells, a factor that might have a great influence on the breadth of the immune response.

Ancient DNA Study in Medieval Europeans Shows an Association Between HLA-DRB1*03 and Paratyphoid Fever

Outbreaks of infectious diseases repeatedly affected medieval Europe, leaving behind a large number of dead often inhumed in mass graves. Human remains interred in two burial pits from 14^th century CE Germany exhibited molecular evidence of Salmonella enterica Paratyphi C (S. Paratyphi C) infection. The pathogen is responsible for paratyphoid fever, which was likely the cause of death for the buried individuals. This finding presented the unique opportunity to conduct a paratyphoid fever association study in a European population. We focused on HLA-DRB1*03:01 that is a known risk allele for enteric fever in present-day South Asians. We generated HLA profiles for 29 medieval S. Paratyphi C cases and 24 contemporaneous controls and compared these to a modern German population. The frequency of the risk allele was higher in the medieval cases (29.6%) compared to the contemporaneous controls (13%; p = 0.189), albeit not significantly so, possibly because of small sample sizes. Indeed, in comparison with the modern controls (n = 39,689; 10.2%; p = 0.005) the frequency difference became statistically significant. This comparison also suggested a slight decrease in the allele’s prevalence between the medieval and modern controls. Up to now, this is the first study on the genetic predisposition to Salmonella infection in Europeans and the only association analysis on paratyphoid fever C. Functional investigation using computational binding prediction between HLA variants and S. Paratyphi and S. Typhi peptides supported a reduced recognition capacity of bacterial proteins by DRB1*03:01 relative to other common DRB1 variants. This pattern could potentially explain the disease association. Our results suggest a slightly reduced predisposition to paratyphoid fever in modern Europeans. The causative allele, however, is still common today, which can be explained by a trade-off, as DRB1*03:01 is protective against infectious respiratory diseases such as severe respiratory syndrome (SARS). It is thus possible that the allele also provided resistance to corona-like viruses in the past.

The scent of attraction and the smell of success: crossmodal influences on person perception

In recent decades, there has been an explosion of research into the crossmodal influence of olfactory cues on multisensory person perception. Numerous peer-reviewed studies have documented that a variety of olfactory stimuli, from ambient malodours through to fine fragrances, and even a range of chemosensory body odours can influence everything from a perceiver's judgments of another person's attractiveness, age, affect, health/disease status, and even elements of their personality. The crossmodal and multisensory contributions to such effects are reviewed and the limitations/peculiarities of the research that have been published to date are highlighted. At the same time, however, it is important to note that the presence of scent (and/or the absence of malodour) can also influence people's (i.e., a perceiver's) self-confidence which may, in turn, affect how attractive they appear to others. Several potential cognitive mechanisms have been put forward to try and explain such crossmodal/multisensory influences, and some of the neural substrates underpinning these effects have now been characterized. At the end of this narrative review, a number of the potential (and actual) applications for, and implications of, such crossmodal/multisensory phenomena involving olfaction are outlined briefly.

Parasitic fauna of Eurasian beavers (Castor fiber) in Sweden (1997-1998)

BACKGROUND

The parasitic fauna of beavers (Castor fiber and C. canadensis) has been well studied in many parts of their respective areas of distribution. In Scandinavia there have, however, been limited investigations conducted on the parasites of beavers in recent times. The present study is the first quantitative survey of parasites on beavers living in Sweden and elsewhere in Scandinavia. We investigated the parasitic fauna of the Eurasian beaver (C. fiber) in a North-South gradient in Sweden. The aim of the study was to investigate parasite distribution and prevalence in particular, related to average yearly air temperature and different age groups of beavers. A total of 30 beavers were sampled at eight localities, spanning a 720 km North-South gradient during the springs of 1997 and 1998.

RESULTS

Five parasite taxa were identified. Four of these were present in all of the examined beavers, Stichorchis subtriquetrus (trematode), Travassosius rufus (nematode), Platypsyllus castoris (coleopteran), and Schizocarpus spp. (arachnid). A higher number of new infections of S. subtriquetrus, and more adults of T. rufus, were seen in beavers in southern Sweden where temperatures are higher. One-year old beavers had a higher infestation of S. subtriquetrus, but not of T. rufus, than older individuals.

CONCLUSIONS

The parasite fauna of Swedish beavers mirrored the impoverished parasite fauna of the original Norwegian population, and the high prevalence of parasites could be due to low major histocompatibility complex (MHC) polymorphism. Young beavers had a higher load of trematodes, probably depending on behavioural and ecological factors. Warmer temperatures in southern localities likely contributed to increased endoparasite loads.

Genetic Variation and Population Differentiation in the Bovine Lymphocyte Antigen DRB3.2 Locus of South African Nguni Crossbred Cattle

Simple Summary

Indigenous cattle breeds are important to their native environment as they confer significant and often unique adaptability traits. The Nguni is one such breeds that is indigenous to the Southern African region. This breed confers resistance to thermal stress and diseases, amongst other factors. The bovine major histocompatibility complex is an important region, which codes for alleles that have been associated with a plethora of diseases. In the current study, the genetic diversity within this region was assessed in Nguni crossbred cattle. This was done to detect the gene pool of the Nguni breed, and to identify genes that might be important within this breed. The populations displayed a high degree of genetic diversity, and some alleles were common throughout the populations and accounted for a significant portion of the total alleles. This high genetic diversity could account for the great adaptability of the Nguni breed to Southern Africa.

Abstract

The bovine lymphocyte antigen (BoLA-DRB3) gene is an important region that codes for glycoproteins responsible for the initiation of an immune response. BoLA-DRB3 alleles have been demonstrated to be associated with disease resistance/tolerance. Therefore, great genetic diversity is correlated with better adaptation, fitness, and robustness. The current study was conducted to assess the population genetic structure of the BoLA-DRB3 gene in Nguni crossbred cattle using polymerase chain reaction-sequence based typing (PCR-SBT). High genetic diversity was detected, with 30 alleles, 11 of which are novel to the study. Alleles DRB3*0201, DRB3*0701, DRB*0901, and DRB*1601 were present in all populations and accounted for nearly around 50% of all observed alleles. A mean genetic diversity (H_E) of 0.93 was detected. The high overall genetic diversity is possibly associated with pathogen-assisted selection and heterozygote advantage. Such high diversity might explain the hardiness of the Nguni crossbred cattle to the Southern African region. Low population genetic structure was identified (F_ST = 0.01), suggesting possible gene flow between populations and retention of similar alleles. The study was undertaken to bridge the dearth of such studies in South African breeds and it is imperative for effective sustainability of indigenous breeds and the implementation of effective breeding strategies.