HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification

The distinct spatiotemporal evolutionary landscape of HBV and HDV largely determines the unique epidemic features of HDV globally

Chronic infection of hepatitis B virus (HBV) and hepatitis D virus (HDV) causes the most severe form of viral hepatitis. Due to the dependence on HBV, HDV was deemed to co-evolve and co-migrate with HBV. However, we previously found that the naturally occurred HDV/HBV combinations do not always reflect the most efficient virological adaptation (Wang et al., 2021). Moreover, regions with heavy HBV burden do not always correlate with high HDV prevalence (e.g., East Asia), and vice versa (e.g., Central Asia). Herein, we systematically elucidated the spatiotemporal evolutionary landscape of HDV to understand the unique epidemic features of HDV. We found that the MRCA of HDV was from South America around the late 13th century, was globally dispersed mainly via Central Asia, and evolved into eight genotypes from the 19th to 20th century. In contrast, the MRCA of HBV was from Europe ∼23.7 thousand years ago (Kya), globally dispersed mainly via Africa and East Asia, and evolved into eight genotypes ∼1100 years ago. When HDV stepped in, all present-day HBV genotypes had already formed and its global genotypic distribution had stayed stable geographically. Nevertheless, regionalized HDV adapted to local HBV genotypes and human lineages, contributing to the global geographical separation of HDV genotypes. Additionally, a sharp increase in HDV infections was observed after the 20th century. In conclusion, HDV exhibited a distinct spatiotemporal distribution path compared with HBV. This unique evolutionary relationship largely fostered the unique epidemic features we observe nowadays. Moreover, HDV infections may continue to ramp up globally, thus more efforts are urgently needed to combat this disease.

A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

Disentangling archaic introgression and genomic signatures of selection at human immunity genes

Pathogens and infectious diseases have imposed exceptionally strong selective pressure on ancient and modern human genomes and contributed to the current variation in many genes. There is evidence that modern humans acquired immune variants through interbreeding with ancient hominins, but the impact of such variants on human traits is not fully understood. The main objectives of this research were to infer the genetic signatures of positive selection that may be involved in adaptation to infectious diseases and to investigate the function of Neanderthal alleles identified within a set of 50 Lithuanian genomes. Introgressed regions were identified using the machine learning tool ArchIE. Recent positive selection signatures were analysed using iHS. We detected high-scoring signals of positive selection at innate immunity genes (EMB, PARP8, HLAC, and CDSN) and evaluated their interactions with the structural proteins of pathogens. Interactions with human immunodeficiency virus (HIV) 1 and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were identified. Overall, genomic regions introgressed from Neanderthals were shown to be enriched in genes related to immunity, keratinocyte differentiation, and sensory perception.

Genomic Markers Associated with Cytomegalovirus DNAemia in Kidney Transplant Recipients

Human cytomegalovirus (CMV) is a major pathogen after solid organ transplantation, leading to high morbidity and mortality. Transplantation from a CMV-seropositive donor to a CMV-seronegative recipient (D+/R-) is associated with high risk of CMV disease. However, that risk is not uniform, suggesting a role for host factors in immune control of CMV. To identify host genetic factors that control CMV DNAemia post transplantation, we performed a whole-exome association study in two cohorts of D+/R- kidney transplant recipients. Quantitative CMV DNA was measured for at least one year following transplantation. Several CMV-protective single-nucleotide polymorphisms (SNPs) were identified in the first cohort (72 patients) but were not reproducible in the second cohort (126 patients). A meta-analysis of both cohorts revealed several SNPs that were significantly associated with protection from CMV DNAemia. The copy number variation of several genes was significantly different between recipients with and without CMV DNAemia. Amongst patients with CMV DNAemia in the second cohort, several variants of interest (p < 5 × 10-5), the most common of which was NLRC5, were associated with peak viral load. We provide new predictive genetic markers for protection of CMV DNAemia. These markers should be validated in larger cohorts.

Persistence of SARS-CoV-2 IgG Antibody Response among South African Adults: A Prospective Cohort Study

This study assesses the durability of severe acute respiratory coronavirus-2 (SARS-CoV-2) anti-nucleocapsid (anti-N) immunoglobulin G (IgG) after infection and examines its association with established risk factors among South African healthcare workers (HCWs). Blood samples were obtained from 390 HCWs with diagnosis of coronavirus disease 2019 (COVID-19) for assay of the SARS-CoV-2 anti-N IgG at two time points (Phase 1 and 2) between November 2020 and February 2021. Out of 390 HCWs with a COVID-19 diagnosis, 267 (68.5%) had detectable SARS-CoV-2 anti-N IgG antibodies at the end of Phase I. These antibodies persisted for 4-5 and 6-7 months in 76.4% and 16.1%, respectively. In the multivariate logistic regression model analysis, Black participants were more likely to sustain SARS-CoV-2 anti-N IgG for 4-5 months. However, participants who were HIV positive were less likely to sustain SARS-CoV-2 anti-N IgG antibodies for 4-5 months. In addition, individuals who were <45 years of age were more likely to sustain SARS-CoV-2 anti-N IgG for 6-7 months. Of the 202 HCWs selected for Phase 2, 116 participants (57.4%) had persistent SARS-CoV-2 anti-N IgG for an extended mean period of 223 days (7.5 months). Findings support the longevity of vaccine responses against SARS-CoV-2 in Black Africans.

HLA-B*13, B*35 and B*39 Alleles Are Closely Associated With the Lack of Response to ART in HIV Infection: A Cohort Study in a Population of Northern Brazil

Introduction

Immune reconstitution failure after HIV treatment is a multifactorial phenomenon that may also be associated with a single polymorphism of human leukocyte antigen (HLA); however, few reports include patients from the Brazilian Amazon. Our objective was to evaluate the association of the immunogenic profile of the “classical” HLA-I and HLA-II loci with treatment nonresponse in a regional cohort monitored over 24 months since HIV diagnosis.

Materials and Methods

Treatment-free participants from reference centers in the state of Pará, Brazil, were enrolled. Infection screening was performed using enzyme immunoassays (Murex AG/AB Combination DiaSorin, UK) and confirmed by immunoblots (Bio-Manguinhos, FIOCRUZ). Plasma viral load was quantified by real-time PCR (ABBOTT, Chicago, Illinois, USA). CD4⁺/CD8⁺ T lymphocyte quantification was performed by immunophenotyping and flow cytometry (BD Biosciences, San Jose, CA, USA). Infection was monitored via test and logistics platforms (SISCEL and SICLOM). Therapeutic response failure was inferred based on CD4⁺ T lymphocyte quantification after 1 year of therapy. Loci A, B and DRB1 were genotyped using PCR-SSO (One Lambda Inc., Canoga Park, CA, USA). Statistical tests were applied using GENEPOP, GraphPad Prism 8.4.3 and BioEstat 5.3.

Results

Of the 270 patients monitored, 134 responded to treatment (CD4⁺ ≥ 500 cells/µL), and 136 did not respond to treatment (CD4⁺ < 500 cells/µL). The allele frequencies of the loci were similar to heterogeneous populations. The allelic profile of locus B was statistically associated with treatment nonresponse, and the B*13, B*35 and B*39 alleles had the greatest probabilistic influence. The B*13 allele had the highest risk of treatment nonresponse, and carriers of the allele had a detectable viral load and a CD4+ T lymphocyte count less than 400 cells/µL with up to 2 years of therapy. The B*13 allele was associated with a switch in treatment regimens, preferably to efavirenz (EFZ)-based regimens, and among those who switched regimens, half had a history of coinfection with tuberculosis.

Conclusions

The allelic variants of the B locus are more associated with non-response to therapy in people living with HIV (PLHIV) from a heterogeneous population in the Brazilian Amazon.

Epitope length variants balance protective immune responses and viral escape in HIV-1 infection

Cytotoxic T lymphocyte (CTL) and natural killer (NK) cell responses to a single optimal 10-mer epitope (KK10) in the human immunodeficiency virus type-1 (HIV-1) protein p24Gag are associated with enhanced immune control in patients expressing human leukocyte antigen (HLA)-B∗27:05. We find that proteasomal activity generates multiple length variants of KK10 (4-14 amino acids), which bind TAP and HLA-B∗27:05. However, only epitope forms ≥8 amino acids evoke peptide length-specific and cross-reactive CTL responses. Structural analyses reveal that all epitope forms bind HLA-B∗27:05 via a conserved N-terminal motif, and competition experiments show that the truncated epitope forms outcompete immunogenic epitope forms for binding to HLA-B∗27:05. Common viral escape mutations abolish (L136M) or impair (R132K) production of KK10 and longer epitope forms. Peptide length influences how well the inhibitory NK cell receptor KIR3DL1 binds HLA-B∗27:05 peptide complexes and how intraepitope mutations affect this interaction. These results identify a viral escape mechanism from CTL and NK responses based on differential antigen processing and peptide competition.

Neutralization diversity of HIV-1 Indian subtype C envelopes obtained from cross sectional and followed up individuals against broadly neutralizing monoclonal antibodies having distinct gp120 specificities

BACKGROUND

The potential use of the broadly neutralizing monoclonal antibodies (bnAbs) towards prophylaxis and treatment to HIV-1 is currently being explored. While a number of promising bnAbs have been discovered and a few of them have progressed towards clinical development, their extent of neutralization coverage with respect to global HIV-1 variants given the existence of genetically distinct subtypes and recombinants circulating globally is not clearly known. In the present study, we examined the variation in the neutralization susceptibility of pseudoviruses expressing 71 full length primary HIV-1 subtype C envs obtained from limited cross-sectional individuals over different time points against four bnAbs that target gp120 with distinct specificities: VRC01, CAP256-VRC26.25, PGDM1400 and PGT121.

RESULTS

We found significant variations in the susceptibility of Indian clade C to these four bnAbs. These variations were found to be distinct to that observed in African subtype C based on the existing datasets and concordant with their sequence diversity. Trend analysis indicated an increasing neutralization resistance observed over time with CAP25-VRC26.25, PGDM1400 and PGT121 when tested on pseudoviruses expressing envs obtained from 1999 to 2016. However, inconsistent trend in neutralization susceptibility was observed, when pseudoviruses expressing envs obtained from three followed up individuals were examined. Finally, through predictive analysis of the 98 Indian subtype C including those assessed in the present study by employing additive model implemented in CombiNAber ( http://www.hiv.lanl.gov ), we observed two possibilities where combinations of three bnAbs (VRC01/CAP56-VRC26.25/PGT121 and PGDM1400/CAP256-VRC26.25/PGT121) could achieve near 100% neutralization coverage.

CONCLUSIONS

Our findings not only indicate disparate intra-clade C genetic vis-à-vis neutralization diversities but also warrant the need for more comprehensive study using additional isolates towards comparing inter and intra-clade neutralization diversities which will be necessary for selecting the bnAb combinations suitable for optimal coverage of the region-specific HIV-1 circulating subtypes. Expanding these efforts is imperative for designing efficacious bnAb based intervention strategies for India as well as subtype C in general.