HLA and HIV-1: heterozygote advantage and B*35-Cw*04 disadvantage

Predicted Epitope Abundance Supports Vaccine-Induced Cytotoxic Protection Against SARS-CoV-2 Variants of Concern

The effect of emerging SARS-CoV-2 variants on vaccine efficacy is of critical importance. In this study, the potential impact of mutations that facilitate escape from the cytotoxic cellular immune response in these new virus variants for the 551 most abundant HLA class I alleles was analyzed. Computational prediction showed that most of these alleles, that cover >90% of the population, contain enough epitopes without escape mutations in the principal SARS-CoV-2 variants. These data suggest that the cytotoxic cellular immune protection elicited by vaccination is not greatly affected by emerging SARS-CoV-2 variants.

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.

Collaboration of a Detrimental HLA-B*35:01 Allele with HLA-A*24:02 in Coevolution of HIV-1 with T Cells Leading to Poorer Clinical Outcomes

Although mutant-specific T cells are elicited in some individuals infected with HIV-1 mutant viruses, the detailed characteristics of these T cells remain unknown. A recent study showed that the accumulation of strains expressing Nef135F, which were selected by HLA-A*24:02-restricted T cells, was associated with poor outcomes in individuals with the detrimental HLA-B*35:01 allele and that HLA-B*35:01-restricted NefYF9 (Nef135-143)-specific T cells failed to recognize target cells infected with Nef135F mutant viruses. Here, we investigated HLA-B*35:01-restricted T cells specific for the NefFF9 epitope incorporating the Nef135F mutation. Longitudinal T-cell receptor (TCR) clonotype analysis demonstrated that 3 types of HLA-B*35:01-restricted T cells (wild-type [WT] specific, mutant specific, and cross-reactive) with different T cell repertoires were elicited during the clinical course. HLA-B*35:01⁺ individuals possessing wild-type-specific T cells had a significantly lower plasma viral load (pVL) than those with mutant-specific and/or cross-reactive T cells, even though the latter T cells effectively recognized the mutant virus-infected cells. These results suggest that mutant-specific and cross-reactive T cells could only partially suppress HIV-1 replication in vivo. An ex vivo analysis of the T cells showed higher expression of PD-1 on cross-reactive T cells and lower expression of CD160/2B4 on the mutant-specific T cells than other T cells, implying that these inhibitory and stimulatory molecules are key to the reduced function of these T cells. In the present study, we demonstrate that mutant-specific and cross-reactive T cells do not contribute to the suppression of HIV-1 replication in HIV-1-infected individuals, even though they have the capacity to recognize mutant virus-infected cells. Thus, the collaboration of HLA-A*24:02 with the detrimental allele HLA-B*35:01 resulted in the coevolution of HIV-1 alongside virus-specific T cells, leading to poorer clinical outcomes.

IMPORTANCE HIV-1 escape mutations are selected under pressure from HIV-1-specific CD8⁺ T cells. Accumulation of these mutations in circulating viruses impairs the control of HIV-1 by HIV-1-specific T cells. Although it is known that HIV-1-specific T cells recognizing mutant virus were elicited in some individuals infected with a mutant virus, the role of these T cells remains unclear. Accumulation of phenylalanine at HIV-1 Nef135 (Nef135F), which is selected by HLA-A*24:02-restricted T cells, led to poor clinical outcome in individuals carrying the detrimental HLA-B*35:01 allele. In the present study, we found that HLA-B*35:01-restricted mutant-specific and cross-reactive T cells were elicited in HLA-B*35:01⁺ individuals infected with the Nef135F mutant virus. These T cells could not effectively suppress HIV-1 replication in vivo even though they could recognize mutant virus-infected cells in vitro. Mutant-specific and cross-reactive T cells expressed lower levels of stimulatory molecules and higher levels of inhibitory molecules, respectively, suggesting a potential mechanism whereby these T cells fail to suppress HIV-1 replication in HIV-1-infected individuals.

Immune Profile of SARS-CoV-2 Variants of Concern

The spread of the current Sars-Cov-2 pandemics leads to the development of mutations that are constantly monitored because they could affect the efficacy of vaccines. Three recently identified mutated strains, known as variants of concern, are rapidly spreading worldwide. Here, we study possible effects of these mutations on the immune response to Sars-Cov-2 infection using NetTepi a computational method based on artificial neural networks that considers binding and stability of peptides obtained by proteasome degradation for widely represented HLA class I alleles present in human populations as well as the T-cell propensity of viral peptides that measures their immune response. Our results show variations in the number of potential highly ranked peptides ranging between 0 and 20% depending on the specific HLA allele. The results can be useful to design more specific vaccines.

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.

Heterozygosity of the major histocompatibility complex predicts later self-reported pubertal maturation in men

Individual variation in the age of pubertal onset is linked to physical and mental health, yet the factors underlying this variation are poorly understood. Life history theory predicts that individuals at higher risk of mortality due to extrinsic causes such as infectious disease should sexually mature and reproduce earlier, whereas those at lower risk can delay puberty and continue to invest resources in somatic growth. We examined relationships between a genetic predictor of infectious disease resistance, heterozygosity of the major histocompatibility complex (MHC), referred to as the human leukocyte antigen (HLA) gene in humans, and self-reported pubertal timing. In a combined sample of men from Canada (n = 137) and the United States (n = 43), MHC heterozygosity predicted later self-reported pubertal development. These findings suggest a genetic trade-off between immunocompetence and sexual maturation in human males.

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.

Exome Sequencing Reveals a Putative Role for HLA-C*03:02 in Control of HIV-1 in African Pediatric Populations

Human leucocyte antigen (HLA) class I molecules present endogenously processed antigens to T-cells and have been linked to differences in HIV-1 disease progression. HLA allelotypes show considerable geographical and inter-individual variation, as does the rate of progression of HIV-1 disease, with long-term non-progression (LTNP) of disease having most evidence of an underlying genetic contribution. However, most genetic analyses of LTNP have occurred in adults of European ancestry, limiting the potential transferability of observed associations to diverse populations who carry the burden of disease. This is particularly true of HIV-1 infected children. Here, using exome sequencing (ES) to infer HLA allelotypes, we determine associations with HIV-1 LTNP in two diverse African pediatric populations. We performed a case-control association study of 394 LTNPs and 420 rapid progressors retrospectively identified from electronic medical records of pediatric HIV-1 populations in Uganda and Botswana. We utilized high-depth ES to perform high-resolution HLA allelotyping and assessed evidence of association between HLA class I alleles and LTNP. Sixteen HLA alleles and haplotypes had significantly different frequencies between Uganda and Botswana, with allelic differences being more prominent in HLA-A compared to HLA-B and C allelotypes. Three HLA allelotypes showed association with LTNP, including a novel association in HLA-C (HLA-B∗57:03, aOR 3.21, Pc = 0.0259; B∗58:01, aOR 1.89, Pc = 0.033; C∗03:02, aOR 4.74, Pc = 0.033). Together, these alleles convey an estimated population attributable risk (PAR) of non-progression of 16.5%. We also observed novel haplotype associations with HLA-B∗57:03-C∗07:01 (aOR 5.40, Pc = 0.025) and HLA-B∗58:01-C∗03:02 (aOR 4.88, Pc = 0.011) with a PAR of 9.8%, as well as a previously unreported independent additive effect and heterozygote advantage of HLA-C∗03:02 with B∗58:01 (aOR 4.15, Pc = 0.005) that appears to limit disease progression, despite weak LD (r 2 = 0.18) between these alleles. These associations remained irrespective of gender or country. In one of the largest studies of HIV in Africa, we find evidence of a protective effect of canonical HLA-B alleles and a novel HLA-C association that appears to augment existing HIV-1 control alleles in pediatric populations. Our findings outline the value of using multi-ethnic populations in genetic studies and offer a novel HIV-1 association of relevance to ongoing vaccine studies.

An immunogenetic view of COVID-19

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

SNPs in Mammary Gland Epithelial Cells Unraveling Potential Difference in Milk Production Between Jersey and Kashmiri Cattle Using RNA Sequencing

Deep RNA sequencing experiment was employed to detect putative single nucleotide polymorphisms (SNP) in mammary epithelial cells between two diverse cattle breeds (Jersey and Kashmiri) to understand the variations in the coding regions that reflect differences in milk production traits. The low milk-producing Kashmiri cattle are being replaced by crossbreeding practices with Jersey cattle with the aim of improving milk production. However, crossbred animals are prone to infections and various other diseases resulting in unsustainable milk production. In this study, we tend to identify high-impact SNPs from Jersey and Kashmiri cows (utilizing RNA-Seq data) to delineate key pathways mediating milk production traits in both breeds. A total of 607 (442 SNPs and 169 INDELs) and 684 (464 SNPs and 220 INDELs) high-impact variants were found specific to Jersey and Kashmir cattle, respectively. Based on our results, we conclude that in Jersey cattle, genes with high-impact SNPs were enriched in nucleotide excision repair pathway, ABC transporter, and metabolic pathways like glycerolipid metabolism, pyrimidine metabolism, and amino acid synthesis (glycine, serine, and threonine). Whereas, in Kashmiri cattle, the most enriched pathways include endocytosis pathway, innate immunity pathway, antigen processing pathway, insulin resistance pathway, and signaling pathways like TGF beta and AMPK which could be a possible defense mechanism against mammary gland infections. A varied set of SNPs in both breeds, suggests a clear differentiation at the genomic level; further analysis of high-impact SNPs are required to delineate their effect on these pathways.

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.

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.

Characteristics of the MACS/WIHS Combined Cohort Study: Opportunities for Research on Aging With HIV in the Longest US Observational Study of HIV

In 2019, the National Institutes of Health combined the Multicenter AIDS Cohort Study (MACS) and the Women's Interagency HIV Study (WIHS) into the MACS/WIHS Combined Cohort Study (MWCCS). In this paper, participants who made a study visit during October 2018-September 2019 (targeted for MWCCS enrollment) are described by human immunodeficiency virus (HIV) serostatus and compared with people living with HIV (PLWH) in the United States. Participants include 2,115 women and 1,901 men with a median age of 56 years (interquartile range, 48-63); 62% are PLWH. Study sites encompass the South (18%), the Mid-Atlantic/Northeast (45%), the West Coast (22%), and the Midwest (15%). Participant race/ethnicity approximates that of PLWH throughout the United States. Longitudinal data and specimens collected for 35 years (men) and 25 years (women) were combined. Differences in data collection and coding were reviewed, and key risk factor and comorbidity data were harmonized. For example, recent use of alcohol (62%) and tobacco (28%) are common, as are dyslipidemia (64%), hypertension (56%), obesity (42%), mildly or severely impaired daily activities (31%), depressive symptoms (28%), and diabetes (22%). The MWCCS repository includes serum, plasma, peripheral blood mononuclear cells, cell pellets, urine, cervicovaginal lavage samples, oral samples, B-cell lines, stool, and semen specimens. Demographic differences between the MACS and WIHS can confound analyses by sex. The merged MWCCS is both an ongoing observational cohort study and a valuable resource for harmonized longitudinal data and specimens for HIV-related research.

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.

Individual HLA-A, -B, -C, and -DRB1 Genotypes Are No Major Factors Which Determine COVID-19 Severity

HLA molecules are key restrictive elements to present intracellular antigens at the crossroads of an effective T-cell response against SARS-CoV-2. To determine the impact of the HLA genotype on the severity of SARS-CoV-2 courses, we investigated data from 6,919 infected individuals. HLA-A, -B, and -DRB1 allotypes grouped into HLA supertypes by functional or predicted structural similarities of the peptide-binding grooves did not predict COVID-19 severity. Further, we did not observe a heterozygote advantage or a benefit from HLA diplotypes with more divergent physicochemical peptide-binding properties. Finally, numbers of in silico predicted viral T-cell epitopes did not correlate with the severity of SARS-CoV-2 infections. These findings suggest that the HLA genotype is no major factor determining COVID-19 severity. Moreover, our data suggest that the spike glycoprotein alone may allow for abundant T-cell epitopes to mount robust T-cell responses not limited by the HLA genotype.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele–pathogen interactions. Specific MHC allele–pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

MiDAS-Meaningful Immunogenetic Data at Scale

Human immunogenetic variation in the form of HLA and KIR types has been shown to be strongly associated with a multitude of immune-related phenotypes. However, association studies involving immunogenetic loci most commonly involve simple analyses of classical HLA allelic diversity, resulting in limitations regarding the interpretability and reproducibility of results. We here present MiDAS, a comprehensive R package for immunogenetic data transformation and statistical analysis. MiDAS recodes input data in the form of HLA alleles and KIR types into biologically meaningful variables, allowing HLA amino acid fine mapping, analyses of HLA evolutionary divergence as well as experimentally validated HLA-KIR interactions. Further, MiDAS enables comprehensive statistical association analysis workflows with phenotypes of diverse measurement scales. MiDAS thus closes the gap between the inference of immunogenetic variation and its efficient utilization to make relevant discoveries related to immune and disease biology. It is freely available under a MIT license.

Using Linear Discriminant Analysis to Characterize Novel Single Nucleotide Polymorphisms and Expression Profile Changes in Genes of Three Breeds of Rabbit (Oryctolagus cuniculus)

The objectives of this study were to investigate polymorphisms and changes in expression patterns of the genes FGF5, PGAM2, TLR2 and IL10 in V-line, Baladi Black and Baladi Red rabbits. Blood samples were collected from 180 healthy rabbits (n = 60 for each breed) for DNA extraction and DNA sequencing. At 3 mo of age, 20 randomly selected females from each breed were euthanized for gene expression quantification in muscle and spleen samples. PCR-DNA sequencing revealed single nucleotide polymorphisms (SNPs) among the 3 breeds that provided a monomorphic pattern for 3 of the 4 genes analyzed. Linear discriminant analysis (LDA) was used to classify the SNPs of these genes in the 3 breeds. The overall percentage of correctly classified cases for the model was 75%, with percentages of 100% for FGF5, 63% for IL10, and 100% for TLR2. Breed was a significant predictor for gene classification with estimation (1.00). Expression profiles of the genes were higher in V-line as compared with Baladi Black or Baladi Red. The LDA discriminated the 3 breeds using results of the gene expression profile as predictors for classification. Overall, 73% of the cases were correctly classified by gene expression. The identified SNPs, along with changes in mRNA levels of FGF5, PGAM2, TLR2, and IL10, could provide a biomarker for efficient characterization of rabbit breeds and could thus help develop marker assisted selection for growth and immune traits in rabbits.

Optimal Maturation of the SIV-Specific CD8+ T Cell Response after Primary Infection Is Associated with Natural Control of SIV: ANRS SIC Study

Highly efficient CD8⁺ T cells are associated with natural HIV control, but it has remained unclear how these cells are generated and maintained. We have used a macaque model of spontaneous SIVmac251 control to monitor the development of efficient CD8⁺ T cell responses. Our results show that SIV-specific CD8⁺ T cells emerge during primary infection in all animals. The ability of CD8⁺ T cells to suppress SIV is suboptimal in the acute phase but increases progressively in controller macaques before the establishment of sustained low-level viremia. Controller macaques develop optimal memory-like SIV-specific CD8⁺ T cells early after infection. In contrast, a persistently skewed differentiation phenotype characterizes memory SIV-specific CD8⁺ T cells in non-controller macaques. Accordingly, the phenotype of SIV-specific CD8⁺ T cells defined early after infection appears to favor the development of protective immunity in controllers, whereas SIV-specific CD8⁺ T cells in non-controllers fail to gain antiviral potency, feasibly as a consequence of early defects imprinted in the memory pool.

Systematic comparative study of computational methods for HLA typing from next-generation sequencing

The human leukocyte antigen (HLA) system plays an important role in hematopoietic stem cell transplantation (HSCT) and organ transplantations, immune disorders as well as oncological immunotherapy. However, HLA typing remains a challenging task due to the high level of polymorphism and homology among HLA genes. Based on the high-throughput next-generation sequencing data, new HLA typing algorithms and software tools were developed. But there is still a deficit of systematic comparative studies to assist in the selection of the optimal analytical approaches under different conditions. Here, we present a detailed comparison of eight software tools for HLA typing on different real datasets (whole-genome sequencing, whole-exome sequencing and transcriptomic sequencing data) and in-silico samples with different sequencing lengths, depths, and error rates. We figure out the algorithms with the best efficiency in different scenarios, and demonstrate the effect of different raw reads on analytical performances. Our results provide a comprehensive picture of specifications and performances of the eight existing HLA genotyping algorithms, which could assist researchers in selecting the most appropriate tool for specific raw datasets.

The association of human leukocyte antigen alleles with clinical disease progression in HIV-positive cohorts with varied treatment strategies

OBJECTIVES

The Strategic Timing of AntiRetroviral Treatment (START) and Strategies for Management of Antiretroviral Therapy (SMART) trials demonstrated that ART can partly reverse clinically defined immune dysfunction induced by HIV replication. As control of HIV replication is influenced by the HLA region, we explored whether HLA alleles independently influence the risk of clinical events in HIV+ individuals.

DESIGN

Cohort study.

METHODS

In START and SMART participants, associations between imputed HLA alleles and AIDS, infection-related cancer, herpes virus-related AIDS events, chronic inflammation-related conditions, and bacterial pneumonia were assessed. Cox regression was used to estimate hazard ratios for the risk of events among allele carriers versus noncarriers. Models were adjusted for sex, age, geography, race, time-updated CD4+ T-cell counts and HIV viral load and stratified by treatment group within trials. HLA class I and II alleles were analyzed separately. The Benjamini--Hochberg procedure was used to limit the false discovery rate to less than 5% (i.e. q value <0.05).

RESULTS

Among 4829 participants, there were 132 AIDS events, 136 chronic inflammation-related conditions, 167 bacterial pneumonias, 45 infection-related cancers, and 49 herpes virus-related AIDS events. Several associations with q value less than 0.05 were found: HLA-DQB1∗06:04 and HLA-DRB1∗13:02 with AIDS (adjusted HR [95% CI] 2.63 [1.5-4.6] and 2.25 [1.4-3.7], respectively), HLA-B∗15:17 and HLA-DPB1∗15:01 with bacterial pneumonia (4.93 [2.3-10.7] and 4.33 [2.0-9.3], respectively), and HLA-A∗69:01 with infection-related cancer (15.26 [3.5-66.7]). The carriage frequencies of these alleles were 10% or less.

CONCLUSION

This hypothesis-generating study suggests that certain HLA alleles may influence the risk of immune dysfunction-related events irrespective of viral load and CD4+ T-cell count.

A Novel Sample Selection Approach to Aid the Identification of Factors That Correlate With the Control of HIV-1 Infection

Individuals infected with HIV display varying rates of viral control and disease progression, with a small percentage of individuals being able to spontaneously control infection in the absence of treatment. In attempting to define the correlates associated with natural protection against HIV, extreme heterogeneity in the datasets generated from systems methodologies can be further complicated by the inherent variability encountered at the population, individual, cellular and molecular levels. Furthermore, such studies have been limited by the paucity of well-characterised samples and linked epidemiological data, including duration of infection and clinical outcomes. To address this, we selected 10 volunteers who rapidly and persistently controlled HIV, and 10 volunteers each, from two control groups who failed to control (based on set point viral loads) from an acute and early HIV prospective cohort from East and Southern Africa. A propensity score matching approach was applied to control for the influence of five factors (age, risk group, virus subtype, gender, and country) known to influence disease progression on causal observations. Fifty-two plasma proteins were assessed at two timepoints in the 1st year of infection. We independently confirmed factors known to influence disease progression such as the B^*57 HLA Class I allele, and infecting virus Subtype. We demonstrated associations between circulating levels of MIP-1α and IL-17C, and the ability to control infection. IL-17C has not been described previously within the context of HIV control, making it an interesting target for future studies to understand HIV infection and transmission. An in-depth systems analysis is now underway to fully characterise host, viral and immunological factors contributing to control.

A Comparison of Cell Activation, Exhaustion, and Expression of HIV Coreceptors and Restriction Factors in HIV-1- and HIV-2-Infected Nonprogressors

Human immunodeficiency viruses induce rare attenuated diseases due either to HIV-1 in the exceptional long-term nonprogressors (LTNPs) or to HIV-2 in West Africa. To better understand characteristics of these two disease types we performed a multiplex comparative analysis of cell activation, exhaustion, and expression of coreceptors and restriction factors in CD4 T cells susceptible to harbor those viruses. We analyzed by flow cytometry the expression of HLA-DR, PD1, CCR5, CXCR6, SAMHD1, Blimp-1, and TRIM5α on CD4 T cell subsets from 10 HIV-1⁺ LTNPs and 14 HIV-2⁺ (12 nonprogressors and 2 progressors) of the ANRS CO-15 and CO-5 cohorts, respectively, and 12 HIV^- healthy donors (HD). The V3 loop of the HIV-1 envelope from 6 HIV-1+ LTNPs was sequenced to determine the CXCR6-binding capacity. Proportions of HLA-DR⁺ and PD1+ cells were higher in memory CD4 T subsets from HIV-1 LTNPs compared with HIV-2 and HD. Similar findings were observed for CCR5+ cells although limited to central-memory CD4 T cell (TCM) and follicular helper T cell subsets, whereas all major subsets from HIV-1 LTNPs contained less CXCR6+ cells compared with HIV-2. All six V3 loop sequences from HIV-1 LTNPs contained a proline at position 326. Proportions of SAMHD1+ cells were higher in all resting CD4 T subsets from HIV-1 LTNPs compared with the other groups, whereas Blimp-1+ and Trim5α+ cells did not differ. The CD4 T cell subsets from HIV-1 LTNPs differ from those of HIV-2-infected subjects by higher levels of activation, exhaustion, and SAMHD1 expression that can reflect the distinct patterns of host/virus relationships.

Association of HLA Class I Genotypes With Severity of Coronavirus Disease-19

Human leukocyte antigen (HLA) class I molecules play a crucial role in the development of a specific immune response to viral infections by presenting viral peptides at the cell surface where they will be further recognized by T cells. In the present manuscript, we explored whether HLA class I genotypes can be associated with the critical course of Coronavirus Disease-19 by searching possible connections between genotypes of deceased patients and their age at death. HLA-A, HLA-B, and HLA-C genotypes of n = 111 deceased patients with COVID-19 (Moscow, Russia) and n = 428 volunteers were identified with next-generation sequencing. Deceased patients were split into two groups according to age at the time of death: n = 26 adult patients aged below 60 and n = 85 elderly patients over 60. With the use of HLA class I genotypes, we developed a risk score (RS) which was associated with the ability to present severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) peptides by the HLA class I molecule set of an individual. The resulting RS was significantly higher in the group of deceased adults compared to elderly adults [p = 0.00348, area under the receiver operating characteristic curve (AUC ROC = 0.68)]. In particular, presence of HLA-A*01:01 allele was associated with high risk, while HLA-A*02:01 and HLA-A*03:01 mainly contributed to low risk. The analysis of patients with homozygosity strongly highlighted these results: homozygosity by HLA-A*01:01 accompanied early deaths, while only one HLA-A*02:01 homozygote died before 60 years of age. Application of the constructed RS model to an independent Spanish patients cohort (n = 45) revealed that the score was also associated with the severity of the disease. The obtained results suggest the important role of HLA class I peptide presentation in the development of a specific immune response to COVID-19.

State-of-the-art genome inference in the human MHC

The Major Histocompatibility Complex (MHC) on the short arm of chromosome 6 is associated with more diseases than any other region of the genome; it encodes the antigen-presenting Human Leukocyte Antigen (HLA) proteins and is one of the key immunogenetic regions of the genome. Accurate genome inference and interpretation of MHC association signals have traditionally been hampered by the region's uniquely complex features, such as high levels of polymorphism; inter-gene sequence homologies; structural variation; and long-range haplotype structures. Recent algorithmic and technological advances have, however, significantly increased the accessibility of genetic variation in the MHC; these developments include (i) accurate SNP-based HLA type imputation; (ii) genome graph approaches for variation-aware genome inference from next-generation sequencing data; (iii) long-read-based diploid de novo assembly of the MHC; (iv) cost-effective targeted MHC sequencing methods. Applied to hundreds of thousands of samples over the last years, these technologies have already enabled significant biological discoveries, for example in the field of autoimmune disease genetics. Remaining challenges concern the development of integrated methods that leverage haplotype-resolved de novo assembly of the MHC for the development of improved MHC genotyping methods for short reads and the construction of improved reference panels for SNP-based imputation. Improved genome inference in the MHC can crucially contribute to an improved genetic and functional understanding of many immune-related phenotypes and diseases.

Levels of Human Immunodeficiency Virus DNA Are Determined Before ART Initiation and Linked to CD8 T-Cell Activation and Memory Expansion

Initiation of antiretroviral therapy (ART) in early compared with chronic human immunodeficiency virus (HIV) infection is associated with a smaller HIV reservoir. This longitudinal analysis of 60 individuals who began ART during primary HIV infection (PHI) investigates which pre- and posttherapy factors best predict HIV DNA levels (a correlate of reservoir size) after treatment initiation during PHI. The best predictor of HIV DNA at 1 year was pre-ART HIV DNA, which was in turn significantly associated with CD8 memory T-cell differentiation (effector memory, naive, and T-bet⁻Eomes⁻ subsets), CD8 T-cell activation (CD38 expression) and T-cell immunoglobulin and mucin-domain containing-3 (Tim-3) expression on memory T cells. No associations were found for any immunological variables after 1 year of ART. Levels of HIV DNA are determined around the time of ART initiation in individuals treated during PHI. CD8 T-cell activation and memory expansion are linked to HIV DNA levels, suggesting the importance of the initial host-viral interplay in eventual reservoir size.

HIV-1 and human genetic variation

Over the past four decades, research on the natural history of HIV infection has described how HIV wreaks havoc on human immunity and causes AIDS. HIV host genomic research, which aims to understand how human genetic variation affects our response to HIV infection, has progressed from early candidate gene studies to recent multi-omic efforts, benefiting from spectacular advances in sequencing technology and data science. In addition to invading cells and co-opting the host machinery for replication, HIV also stably integrates into our own genome. The study of the complex interactions between the human and retroviral genomes has improved our understanding of pathogenic mechanisms and suggested novel preventive and therapeutic approaches against HIV infection.

Critical effect of Pol escape mutations associated with detrimental allele HLA-C*15: 05 on clinical outcome in HIV-1 subtype A/E infection

OBJECTIVE

The mechanism explaining the role of detrimental HLA alleles in HIV-1 infections has been investigated in very few studies. HLA-A*29:01-B*07:05-C*15:05 is a detrimental haplotype in HIV-1 subtype A/E-infected Vietnamese individuals. The accumulation of mutations at Pol 653/657 is associated with a poor clinical outcome in these individuals. However, the detrimental HLA allele and the mechanism responsible for its detrimental effect remains unknown. Therefore, in this current study we identified the detrimental HLA allele and investigated the mechanism responsible for the detrimental effect.

DESIGN AND METHODS

A T-cell epitope including Pol 653/657 and its HLA restriction were identified by using overlapping HIV-1 peptides and cell lines expressing a single HLA. The effect of the mutations on the T-cell recognition of HIV-1-infected cells was investigated by using target cells infected with the mutant viruses. The effect of these mutations on the clinical outcome was analyzed in 74 HLA-C*15:05 Vietnamese infected with the subtype A/E virus.

RESULTS

We identified HLA-C*15:05-restricted SL9 epitope including Pol 653/657. PolS653A/T/L mutations within this epitope critically impaired the T-cell recognition of HIV-1-infected cells, indicating that these mutations had escaped from the T cells. T-cell responders infected with these mutants showed significantly lower CD4 T-cell counts than those with the wild-type virus or Pol S653K/Q mutants, which are not associated with HLA-C*15:05.

CONCLUSION

The accumulation of Pol S653A/T/L escape mutants critically affected the control of HIV-1 by SL9-specific T cells and led to a poor clinical outcome in the subtype A/E-infected individuals having the detrimental HLA-C*15:05 allele.

[In vivo protective mechanisms of neutralizing antibodies against simian immunodeficiency virus replicatio]

Identifying protective adaptive immune responses against human immunodeficiency virus type 1 (HIV-1), mainly comprising CD8+ cytotoxic T lymphocyte (CTL) and neutralizing antibody (NAb) responses, is crucial for understanding in vivo mechanisms of viral persistence and developing prophylactic/intervention strategies. In HIV-1 and pathogenic simian immunodeficiency virus (SIV) infections, CTL responses play the canonical role in primary viral replication control, whereas NAb responses are impaired. This NAb impairment in early infection conversely highlights the necessity of elucidating anti-HIV/SIV antibody defense/induction mechanisms, and one approach to analyze the impact of NAbs on HIV/SIV infection is passive immunization. We have analyzed a simian AIDS model of highly pathogenic SIVmac239-infected rhesus macaques, and characterized that a single acute-phase passive infusion of SIV-specific polyclonal NAbs drives a synergistic qualitative boosting of virus-specific T-cell responses, resulting in sustained SIV replication control. This in vivo functional augmentation of virus-specific T cells by NAbs in the SIV model provides insights into the design of protective immunity against HIV-1 infection.

MHC Class II Heterozygosity Associated With Attractiveness of Men and Women

The genes of the Major Histocompatibility Complex (MHC), which plays a fundamental role in the immune system, are some of the most diverse genes in vertebrates and have been connected to mate choice in several species, including humans. While studies suggest a positive relationship between MHC diversity and male facial attractiveness, the connection of MHC diversity to other visual traits and female attractiveness is still unclear. The purpose of this study was to investigate further whether MHC heterozygosity, indicating genetic quality, is associated with visual traits affecting mate preferences in humans. In total 74 Latvian men and 49 women were genotyped for several MHC loci and rated for facial and, in men, also body attractiveness. The results indicate a preference for MHC heterozygous female and male faces. However, the initially positive relationship between MHC heterozygosity and facial attractiveness becomes non-significant in females, when controlling for multiple testing, and in males, when age and fat content is taken into account, referring to the importance of adiposity in immune function and thus also attractiveness. Thus overall the effect of MHC heterozygosity on attractiveness seems weak. When considering separate loci, we show that the main gene related to facial attractiveness is the MHC class II DQB1; a gene important also in viral infections and autoimmune diseases. Indeed, in our study, heterozygous individuals are rated significantly more attractive than their homozygous counterparts, only in relation to gene DQB1. This study is the first to indicate a link between DQB1 and attractiveness in humans.

Depletion of MHC supertype during domestication can compromise immunocompetence

The major histocompatibility complex (MHC) plays an important role in infectious disease resistance. The presence of certain MHC alleles and functionally similar groups of MHC alleles (i.e., supertypes) has been associated with resistance to particular parasite species. Farmed and domesticated fish stocks are often depleted in their MHC alleles and supertype diversity, possibly as a consequence of artificial selection for desirable traits, inbreeding (loss of heterozygosity), genetic drift (loss of allelic diversity) and/or reduced parasite biodiversity. Here we quantify the effects of depletion of MHC class II genotype and supertype variation on resistance to the parasite Gyrodactylus turnbulli in guppies (Poecilia reticulata). Compared to the descendants of wild‐caught guppies, ornamental fish had a significantly reduced MHC variation (i.e., the numbers of MHC alleles and supertypes per individual, and per population). In addition, ornamental fish were significantly more susceptible to G. turnbulli infections, accumulating peak intensity 10 times higher than that of their wildtype counterparts. Four out of 13 supertypes were associated with a significantly reduced parasite load, and the presence of some supertypes had a dramatic effect on the intensity of infection. Remarkably, the ornamental and wildtype fish differed in the supertypes that were associated with parasite resistance. Analysis with a genetic algorithm showed that resistance‐conferring supertypes of the wildtype and ornamental fish shared two unique amino acids in the peptide‐binding region of the MHC that were not found in any other alleles. These data show that the supertype demarcation captures some, but not all, of the variation in the immune function of the alleles. This study highlights the importance of managing functional MHC diversity in livestock, and suggests there might be some immunological redundancy among MHC supertypes.

Host genetics and infectious disease: new tools, insights and translational opportunities

Understanding how human genetics influence infectious disease susceptibility offers the opportunity for new insights into pathogenesis, potential drug targets, risk stratification, response to therapy and vaccination. As new infectious diseases continue to emerge, together with growing levels of antimicrobial resistance and an increasing awareness of substantial differences between populations in genetic associations, the need for such work is expanding. In this Review, we illustrate how our understanding of the host–pathogen relationship is advancing through holistic approaches, describing current strategies to investigate the role of host genetic variation in established and emerging infections, including COVID-19, the need for wider application to diverse global populations mirroring the burden of disease, the impact of pathogen and vector genetic diversity and a broad array of immune and inflammation phenotypes that can be mapped as traits in health and disease. Insights from study of inborn errors of immunity and multi-omics profiling together with developments in analytical methods are further advancing our knowledge of this important area.

Infectious diseases are an ever-present global threat. In this Review, Kwok, Mentzer and Knight discuss our latest understanding of how human genetics influence susceptibility to disease. Furthermore, they discuss emerging progress in the interplay between host and pathogen genetics, molecular responses to infection and vaccination, and opportunities to bring these aspects together for rapid responses to emerging diseases such as COVID-19.

The HLA diversity of the Anthony Nolan register

While the success of allogeneic stem cell transplantation depends on a high degree of HLA compatibility between donor and patient, finding a suitable donor remains challenging due to the hyperpolymorphic nature of HLA genes. We calculated high-resolution allele, haplotype and phenotype frequencies for HLA-A, -C, -B, -DRB1 and -DQB1 for 10 subpopulations of the Anthony Nolan (AN) register using an in-house expectation-maximisation (EM) algorithm run on mixed resolution HLA data, covering 676 155 individuals. Sample sizes range from 599 410 for British/Irish North West European (BINWE) individuals, the largest subpopulation in the United Kingdom to 1105 for the British Bangladeshi population. Calculation of genetic distance between the subpopulations based on haplotype frequencies shows three broad clusters, each following a major continental group: European, African and Asian. We further analysed the HLA haplotype and phenotype diversity of each subpopulation, and found that 35.52% of BINWE individuals ranging to 98.34% of Middle Eastern individuals on the register had a unique phenotype within their subpopulation. These analyses and the allele, haplotype and phenotype frequency data of the subpopulation on the AN register are a valuable resource in understanding the HLA diversity in the United Kingdom and can be used to improve the accuracy of match likelihoods and to inform future donor recruitment strategies.

HLA tapasin independence: broader peptide repertoire and HIV control

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

Haplobanking induced pluripotent stem cells for clinical use

The development of induced pluripotent stem cells (iPSCs) by Shinya Yamanaka and colleagues in 2006 has led to a potential new paradigm in cellular therapeutics, including the possibility of producing patient-specific, disease-specific and immune matched allogeneic cell therapies. One can envisage two routes to immunologically compatible iPSC therapies: using genetic modification to generate a 'universal donor' with reduced expression of Human Leukocyte Antigens (HLA) and other immunological targets or developing a haplobank containing iPSC lines specifically selected to provide HLA matched products to large portions of the population. HLA matched lines can be stored in a designated physical or virtual global bank termed a 'haplobank'. The process of 'iPSC haplobanking' refers to the banking of iPSC cell lines, selected to be homozygous for different HLA haplotypes, from which therapeutic products can be derived and matched immunologically to patient populations. By matching iPSC and derived products to a patient's HLA class I and II molecules, one would hope to significantly reduce the risk of immune rejection and the use of immunosuppressive medication. Immunosuppressive drugs are used in several conditions (including autoimmune disease and in transplantation procedures) to reduce rejection of infused cells, or transplanted tissue and organs, due to major and minor histocompatibility differences between donor and recipient. Such regimens can lead to immune compromise and pathological consequences such as opportunistic infections or malignancies due to decreased cancer immune surveillance. In this article, we will discuss what is practically involved if one is developing and executing an iPSC haplobanking strategy.

The effect of human leukocyte antigen A1 and B35-Cw4 on sustained BK polyomavirus DNAemia after renal transplantation

Human leukocyte antigen (HLA) class I presentation pathway plays a central role in natural killer (NK) cell and cytotoxic T-cell activities against BK polyomavirus (BKPyV) DNAemia. We determined the risk of sustained BKPyV DNAemia in 175 consecutive renal transplant recipients considering the simultaneous effect of donor/recipient HLA class I antigens and pre- or post-transplant variables. Median (IQR) age was 53 (44-64) years, and 37% of patients were female. 40 patients (22.9%) developed sustained BKPyV DNAemia [median (IQR) viral load: 9740 (4350-17 125) copies/ml]. In the Cox proportional hazard analysis, HLA-A1 (HR: 3.06, 95% CI: 1.51-6.17) and HLA-B35-Cw4 (HR: 4.63, 95% CI: 2.12-10.14) significantly increased the risk of sustained BKPyV DNAemia, while 2 HLA-C mismatches provided a marginally protective effect (HR: 0.32, 95% CI: 0.10-0.98). HLA-Cw4 is a ligand for NK cell inhibitory receptor, and HLA-B35 is in strong linkage disequilibrium with the HLA-Cw4 allele. The association between HLA-B35-Cw4 expression and sustained BKPyV DNAemia supports the important role of cytotoxic T cells and NK cells that would normally control BKPyV activation through engagement with immunoglobulin-like killer receptors (KIRs). Further studies are required to investigate the effect of HLA-C alleles along with NK cell activity against BKPyV DNAemia.

Genomics-based immuno-oncology: bridging the gap between immunology and tumor biology

The first hypotheses about how the immune system affects cancers were proposed in the early 20th century. These early concepts about cancer immunosurveillance were further developed in the decades that followed, but a detailed understanding of cancer immunity remained elusive. It was only recently, through the advent of high-throughput technologies, that scientists gained the ability to profile tumors with a resolution that allowed for granular assessment of both tumor cells and the tumor microenvironment. The advent of immune checkpoint inhibitors (ICIs), which have proven to be effective cancer therapies in many malignancies, has spawned great interest in developing biomarkers for efficacy, an endeavor that highlighted the value of dissecting tumor immunity using large-scale methods. Response to ICI therapy has been shown to be a highly complex process, where the dynamics of tumor and immune cells is key to success. The need to understand the biologic mechanisms at the tumor-immune interface has given rise to the field of cancer immunogenomics, a discipline that aims to bridge the gap between cancer genomics and classical immunology. We provide a broad overview of this emerging branch of translational science, summarizing common platforms used and recent discoveries in the field, which are having direct clinical implications. Our discussion will be centered around the genetic foundations governing tumor immunity and molecular determinants associated with clinical benefit from ICI therapy. We emphasize the importance of molecular diversity as a driver of anti-tumor immunity and discuss how these factors can be probed using genomic approaches.

Specific Class I HLA Supertypes but Not HLA Zygosity or Expression Are Associated with Outcomes following HLA-Matched Allogeneic Hematopoietic Cell Transplant: HLA Supertypes Impact Allogeneic HCT Outcomes

Maximizing the probability of antigen presentation to T cells through diversity in HLAs can enhance immune responsiveness and translate into improved clinical outcomes, as evidenced by the association of heterozygosity and supertypes at HLA class I loci with improved survival in patients with advanced solid tumors treated with immune checkpoint inhibitors. We investigated the impact of HLA heterozygosity, supertypes, and surface expression on outcomes in adult and pediatric patients with acute myeloid leukemia (AML), myelodysplastic syndrome, acute lymphoblastic leukemia, and non-Hodgkin lymphoma who underwent 8/8 HLA-matched, T cell replete, unrelated, allogeneic hematopoietic cell transplant (HCT) from 2000 to 2015 using patient data reported to the Center for International Blood and Marrow Transplant Research. HLA class I heterozygosity and HLA expression were not associated with overall survival, relapse, transplant-related mortality (TRM), disease-free survival (DFS), and acute graft-versus-host disease following HCT. The HLA-B62 supertype was associated with decreased TRM in the entire patient cohort (hazard ratio [HR], 0.79; 95% CI, 0.69 to 0.90; P = .00053). The HLA-B27 supertype was associated with worse DFS in patients with AML (HR = 1.21; 95% CI, 1.10 to 1.32; P = .00005). These findings suggest that the survival benefit of HLA heterozygosity seen in solid tumor patients receiving immune checkpoint inhibitors does not extend to patients undergoing allogeneic HCT. Certain HLA supertypes, however, are associated with TRM and DFS, suggesting that similarities in peptide presentation between supertype members play a role in these outcomes. Beyond implications for prognosis following HCT, these findings support the further investigation of these HLA supertypes and the specific immune peptides important for transplant outcomes.

Role of Escape Mutant-Specific T Cells in Suppression of HIV-1 Replication and Coevolution with HIV-1

Escape mutant-specific CD8⁺ T cells were elicited in some individuals infected with escape mutants, but it is still unknown whether these CD8⁺ T cells can suppress HIV-1 replication. We clarified that Gag280V mutation were selected by HLA-B*52:01-restricted CD8⁺ T cells specific for the GagRI8 protective epitope, whereas the Gag280V virus could frequently elicit GagRI8-6V mutant-specific CD8⁺ T cells. GagRI8-6V mutant-specific T cells had a strong ability to suppress the replication of the Gag280V mutant virus both in vitro and in vivo. In addition, these T cells contributed to the selection of wild-type virus in HLA-B*52:01⁺ Japanese individuals. We for the first time demonstrated that escape mutant-specific CD8⁺ T cells can suppress HIV-1 replication and play an important role in the coevolution with HIV-1. Thus, the present study highlighted an important role of escape mutant-specific T cells in the control of HIV-1 and coevolution with HIV-1.

The accumulation of HIV-1 escape mutations affects HIV-1 control by HIV-1-specific T cells. Some of these mutations can elicit escape mutant-specific T cells, but it still remains unclear whether they can suppress the replication of HIV-1 mutants. It is known that HLA-B*52:01-restricted RI8 (Gag 275 to 282; RMYSPTSI) is a protective T cell epitope in HIV-1 subtype B-infected Japanese individuals, though 3 Gag280A/S/V mutations are found in 26% of them. Gag280S and Gag280A were HLA-B*52:01-associated mutations, whereas Gag280V was not, implying a different mechanism for the accumulation of Gag280 mutations. In this study, we investigated the coevolution of HIV-1 with RI8-specific T cells and suppression of HIV-1 replication by its escape mutant-specific T cells both in vitro and in vivo. HLA-B*52:01⁺ individuals infected with Gag280A/S mutant viruses failed to elicit these mutant epitope-specific T cells, whereas those with the Gag280V mutant one effectively elicited RI8-6V mutant-specific T cells. These RI8-6V-specific T cells suppressed the replication of Gag280V virus and selected wild-type virus, suggesting a mechanism affording no accumulation of the Gag280V mutation in the HLA-B*52:01⁺ individuals. The responders to wild-type (RI8-6T) and RI8-6V mutant peptides had significantly higher CD4 counts than nonresponders, indicating that the existence of not only RI8-6T-specific T cells but also RI8-6V-specific ones was associated with a good clinical outcome. The present study clarified the role of escape mutant-specific T cells in HIV-1 evolution and in the control of HIV-1.

IMPORTANCE Escape mutant-specific CD8⁺ T cells were elicited in some individuals infected with escape mutants, but it is still unknown whether these CD8⁺ T cells can suppress HIV-1 replication. We clarified that Gag280V mutation were selected by HLA-B*52:01-restricted CD8⁺ T cells specific for the GagRI8 protective epitope, whereas the Gag280V virus could frequently elicit GagRI8-6V mutant-specific CD8⁺ T cells. GagRI8-6V mutant-specific T cells had a strong ability to suppress the replication of the Gag280V mutant virus both in vitro and in vivo. In addition, these T cells contributed to the selection of wild-type virus in HLA-B*52:01⁺ Japanese individuals. We for the first time demonstrated that escape mutant-specific CD8⁺ T cells can suppress HIV-1 replication and play an important role in the coevolution with HIV-1. Thus, the present study highlighted an important role of escape mutant-specific T cells in the control of HIV-1 and coevolution with HIV-1.

Estimating the Binding of Sars-CoV-2 Peptides to HLA Class I in Human Subpopulations Using Artificial Neural Networks

Epidemiological studies show that SARS-CoV-2 infection leads to severe symptoms only in a fraction of patients, but the determinants of individual susceptibility to the virus are still unknown. The major histocompatibility complex (MHC) class I exposes viral peptides in all nucleated cells and is involved in the susceptibility to many human diseases. Here, we use artificial neural networks to analyze the binding of SARS-CoV-2 peptides with polymorphic human MHC class I molecules. In this way, we identify two sets of haplotypes present in specific human populations: the first displays weak binding with SARS-CoV-2 peptides, while the second shows strong binding and T cell propensity. Our work offers a useful support to identify the individual susceptibility to COVID-19 and illustrates a mechanism underlying variations in the immune response to SARS-CoV-2. A record of this paper’s transparent peer review process is included in the Supplemental Information.

•

Binding of SARS-CoV-2 peptides to HLA molecules is computed

•

Weakly or strongly binding haplotypes are identified in human populations

•

Results explain variations in the individual immune response to SARS-CoV-2

Genetic and environmental determinants of human TCR repertoire diversity

T cell discrimination of self and non-self is the foundation of the adaptive immune response, and is orchestrated by the interaction between T cell receptors (TCRs) and their cognate ligands presented by major histocompatibility (MHC) molecules. However, the impact of host immunogenetic variation on the diversity of the TCR repertoire remains unclear. Here, we analyzed a cohort of 666 individuals with TCR repertoire sequencing. We show that TCR repertoire diversity is positively associated with polymorphism at the human leukocyte antigen class I (HLA-I) loci, and diminishes with age and cytomegalovirus (CMV) infection. Moreover, our analysis revealed that HLA-I polymorphism and age independently shape the repertoire in healthy individuals. Our data elucidate key determinants of human TCR repertoire diversity, and suggest a mechanism underlying the evolutionary fitness advantage of HLA-I heterozygosity.

Innate, non-cytolytic CD8+ T cell-mediated suppression of HIV replication by MHC-independent inhibition of virus transcription

MHC-I-restricted, virus-specific cytotoxic CD8+ T cells (CTLs) may control human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) replication via the recognition and killing of productively infected CD4+ T cells. Several studies in SIV-infected macaques suggest that CD8+ T cells may also decrease virus production by suppressing viral transcription. Here, we show that non-HIV-specific, TCR-activated non-cytolytic CD8+ T cells suppress HIV transcription via a virus- and MHC-independent immunoregulatory mechanism that modulates CD4+ T cell proliferation and activation. We also demonstrate that this CD8+ T cell-mediated effect promotes the survival of infected CD4+ T cells harboring integrated, inducible virus. Finally, we used RNA sequencing and secretome analyses to identify candidate cellular pathways that are involved in the virus-silencing mediated by these CD8+ T cells. This study characterizes a previously undescribed mechanism of immune-mediated HIV silencing that may be involved in the establishment and maintenance of the reservoir under antiretroviral therapy and therefore represent a major obstacle to HIV eradication.

arcasHLA: high-resolution HLA typing from RNAseq

MOTIVATION

The human leukocyte antigen (HLA) locus plays a critical role in tissue compatibility and regulates the host response to many diseases, including cancers and autoimmune di3orders. Recent improvements in the quality and accessibility of next-generation sequencing have made HLA typing from standard short-read data practical. However, this task remains challenging given the high level of polymorphism and homology between HLA genes. HLA typing from RNA sequencing is further complicated by post-transcriptional modifications and bias due to amplification.

RESULTS

Here, we present arcasHLA: a fast and accurate in silico tool that infers HLA genotypes from RNA-sequencing data. Our tool outperforms established tools on the gold-standard benchmark dataset for HLA typing in terms of both accuracy and speed, with an accuracy rate of 100% at two-field resolution for Class I genes, and over 99.7% for Class II. Furthermore, we evaluate the performance of our tool on a new biological dataset of 447 single-end total RNA samples from nasopharyngeal swabs, and establish the applicability of arcasHLA in metatranscriptome studies.

AVAILABILITY AND IMPLEMENTATION

arcasHLA is available at https://github.com/RabadanLab/arcasHLA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Selective Upregulation of CTLA-4 on CD8+ T Cells Restricted by HLA-B*35Px Renders them to an Exhausted Phenotype in HIV-1 infection

HLA-B*35Px is associated with HIV-1 disease rapid progression to AIDS. However, the mechanism(s) underlying this deleterious effect of this HLA allele on HIV-1 infection outcome has not fully understood. CD8⁺ T cells play a crucial role to control the viral replication but impaired CD8⁺ T cells represent a major hallmark of HIV-1 infection. Here, we examined the effector functions of CD8⁺ T cells restricted by HLA-B*35Px (HLA-B*35:03 and HLA-B*35:02), HLA-B*27/B57 and non-HLA-B*27/B57 (e.g. HLA-A*01, A*02, A*03, A*11, A*24, A*26, B*40, B*08, B*38, B*44). CD8⁺ T cells restricted by HLA-B*35Px exhibited an impaired phenotype compared with those restricted by HLA-B*27/B57 and even non-HLA-B*27/B57. CD8⁺ T cells restricted by non-HLA-B*27/B57 when encountered their cognate epitopes upregulated TIM-3 and thus became suppressed by regulatory T cells (Tregs) via TIM-3: Galectin-9 (Gal-9). Strikingly, CD8⁺ T cells restricted by HLA-B*35Px expressed fewer TIM-3 and therefore did not get suppressed by Tregs, which was similar to CD8⁺ T cells restricted by HLA-B*27/B57. Instead, CD8⁺ T cells restricted by HLA-B*35Px upon recognition of their cognate epitopes upregulated CTLA-4. The transcriptional and impaired phenotype (e.g. poor effector functions) of HIV-specific CD8⁺ T cells restricted by HLA-B*35 was related to persistent CTLA-4, elevated Eomes and blimp-1 but poor T-bet expression. As such, anti-CTLA-4 antibody, Ipilimumab, reversed the impaired proliferative capacity of antigen-specific CD8⁺ T cells restricted by HLA-B*35Px but not others. This study supports the concept that CD8⁺ T resistance to Tregs-mediated suppression is related to allele restriction rather than the epitope specificity. Our results aid to explain a novel mechanism for the inability of HIV-specific CD8⁺ T cells restricted by HLA-B*35Px to control viral replication.

A rare group of HIV-infected individuals with HLA-B*35Px rapidly progress to AIDS but those with HLA-B*27 and HLA-B*57 spare disease progression. Previous studies have suggested that viral mutation may prevent a robust immune response against the virus in these with HLA-B*35Px. However, the functionality of HIV-specific CD8⁺ T cells restricted by HLA-B*35Px remains unclear. In this study, we demonstrate that HIV-specific CD8⁺ T cells restricted by HLA-B*35Px (HLA-B*35:03 and HLA-B*35:02) exhibit an impaired phenotype (e.g. low proliferative capacity, poor cytotoxic molecules expression and, poor cytokine production ability). Interestingly, CD8⁺ T cells restricted by HLA-B*27/B*57 evade regulatory T cells (Tregs) suppression but not those restricted by non-HLA-B*27/B*57. CD8⁺ T cells restricted by non-HLA-B*27/B*57 when encountering their epitopes upregulate TIM-3 but not those restricted by HLA-B*27/B*57 and HLA-B*35Px. As a result, CD8⁺ T cells restricted by non-HLA-B*27/B*57 become suppressed by Tregs via TIM-3: Galectin-9 interactions. Strikingly, CD8⁺ T cells restricted by HLA-B*35Px upregulate CTLA-4 when encountering their epitopes, which render them to an exhausted phenotype. This differential response is linked to the up-regulation of Eomes, Blimp-1 but low T-bet expression in CD8⁺ T cells restricted by HLA-B*35Px. These results implicate that reinvigoration of these cells might be feasible using an anti-CTLA-4 antibody.

MHC-associated Baylisascaris schroederi load informs the giant panda reintroduction program

Reintroducing captive giant pandas (Ailuropoda melanoleuca) to the wild is the ultimate goal of their ex situ conservation. Choosing higher fitness candidates to train prior to release is the first step in the giant panda reintroduction program. Disease resistance is one important index of individual fitness and presumed to be related to variation at major histocompatibility complex genes (MHC). Here, we used seven polymorphic functional MHC genes (Aime-C, Aime-I, Aime-L, Aime-DQA1, Aime-DQA2, Aime-DQB1 and Aime-DRB3) and estimate their relationship with Baylisascaris schroederi (Ascarididae) infection in giant panda. We found that DQA1 heterozygous pandas were less frequently infected than homozygotes. The presence of one MHC genotype and one MHC allele were also associated with B. schroederi infection: Aime-C*0203 and Aime-L*08 were both associated with B. schroederi resistance. Our results indicate that both heterozygosity and certain MHC variants are important for panda disease resistance, and should therefore be considered in future reintroduction programs for this species alongside conventional selection criteria (such as physical condition and pedigree-based information).

•

MHC heterozygous pandas were less frequently infected by Baylisascaris schroederi than homozygotes.

•

Presence of Aime-C*0203 and Aime-L*08 are associated with Baylisascaris schroederi resistance.

•

MHC types are important for panda parasite resistance.

Association Between Single-Nucleotide Polymorphisms in HLA Alleles and Human Immunodeficiency Virus Type 1 Viral Load in Demographically Diverse, Antiretroviral Therapy-Naive Participants From the Strategic Timing of AntiRetroviral Treatment Trial

The impact of variation in host genetics on replication of human immunodeficiency virus type 1 (HIV-1) in demographically diverse populations remains uncertain. In the current study, we performed a genome-wide screen for associations of single-nucleotide polymorphisms (SNPs) to viral load (VL) in antiretroviral therapy-naive participants (n = 2440) with varying demographics from the Strategic Timing of AntiRetroviral Treatment (START) trial. Associations were assessed using genotypic data generated by a customized SNP array, imputed HLA alleles, and multiple linear regression. Genome-wide significant associations between SNPs and VL were observed in the major histocompatibility complex class I region (MHC I), with effect sizes ranging between 0.14 and 0.39 log10 VL (copies/mL). Supporting the SNP findings, we identified several HLA alleles significantly associated with VL, extending prior observations that the (MHC I) is a major host determinant of HIV-1 control with shared genetic variants across diverse populations and underscoring the limitations of genome-wide association studies as being merely a screening tool.

Hepatitis B infection outcome is associated with novel human leukocyte antigen variants in Ghanaian cohort

IMPACT STATEMENT

Genetic association studies can determine the effect size of gene loci on disease outcomes. In the arena of HBV infections, HLA alleles that associate with HBV outcomes can be used in clinical management decisions. This potential translational utility can shape the future management of HBV infections by identifying at-risk individuals and tailoring medical interventions accordingly. This precision medicine motif is currently only a nascent idea. However, it has stakes that may well override the current "wait and see" approach of clinical management of HBV infections. Here, we have identified HLA alleles associated with HBV outcome in a Ghanaian cohort. Our findings support the motif that HLA alleles associate with HBV outcome along geo-ethnic lines. This buttresses the need for further population pivoted studies. In the long term, our findings add to efforts towards the development of an HLA molecular-based algorithm for predicting HBV infection outcomes.

HIV Transmission Chains Exhibit Greater HLA-B Homogeneity Than Randomly Expected

BACKGROUND

HIV's capacity to escape immune recognition by human leukocyte antigen (HLA) is a core component of HIV pathogenesis. A better understanding of the distribution of HLA class I in HIV-infected patients would improve our knowledge of pathogenesis in relation to the host HLA type and could better improve therapeutic strategies against HIV.

MATERIALS AND METHODS

Three hundred one to 325 transmission pairs and 469-496 clusters were identified for analysis among Swiss HIV Cohort Study (SHCS) participants using HIV pol sequences from the drug resistance database. HLA class I data were compiled at 3 specificity levels: 4-digit, 2-digit alleles, and HLA-B supertype. The analysis tabulated HLA-I homogeneity as 2 measures: the proportion of transmission pairs, which are HLA concordant, and the average percentage of allele matches within all clusters. These measures were compared with the mean value across randomizations with randomly assorted individuals.

RESULTS

We repeated the analysis for different HLA classification levels and separately for HLA-A, -B, and -C. Subanalyses by the risk group were performed for HLA-B. HLA-B showed significantly greater homogeneity in the transmission chains (2-digit clusters: 0.291 vs. 0.251, P value = 0.009; supertype clusters: 0.659 vs. 0.611, P value = 0.002; supertype pairs: 0.655 vs. 0.608, P value = 0.014). Risk group restriction caused the effect to disappear for men-who-have-sex-with-men but not for other risk groups. We also examined if protective HLA alleles B27 and B57 were under- or overrepresented in the transmission chains, although this yielded no significant pattern.

CONCLUSIONS

The HLA-B alleles of patients within HIV-1 transmission chains segregate in homogenous clusters/pairs, potentially indicating preferential transmission among HLA-B concordant individuals.