The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control.

Introduction: We have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers. Methods: We performed a genome-wide Association Study comparing 543 HIV Elite Controllers with 3,272 uninfected controls of European descent. Using the latest database for imputation, we analyzed 35,552 Single Nucleotide Polymorphisms (SNPs) within the Major Histocompatibility Complex (MHC) region. Results: Our analysis identified 2,626 SNPs significantly associated (p<5. 10-8) with elite control of HIV-1 infection, including well-established MHC signals such as the rs2395029-G allele which tags HLA-B*57:01. A thorough investigation of SNPs in linkage disequilibrium with rs2395029 revealed an extensive haploblock spanning 1.9 megabases in the MHC region tagging HLA-B*57:01, comprising 379 SNP alleles impacting 72 genes. This haploblock contains damaging variations in proteins like NOTCH4 and DXO and is also associated with a strong differential pattern of expression of multiple MHC genes such as HLA-B, MICB, and ZBTB12. The study was expanded to include two cohorts of seropositive African-American individuals, where a haploblock tagging the HLA-B*57:03 allele was similarly associated with control of viral load. The mRNA expression profile of this haploblock in African Americans closely mirrored that in the European cohort. Discussion: These findings suggest that additional molecular mechanisms beyond the conventional antigen-presenting role of class I HLA molecules may contribute to the observed influence of HLA-B*57:01/B*57:03 alleles on HIV-1 elite control. Overall, this study has uncovered a large haploblock associated with HLA-B*57 alleles, providing novel insights into their massive effect on HIV-1 elite control.

Transcriptomics and chromatin accessibility in multiple African population samples.

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

Community Insights in Phylogenetic HIV Research: The CIPHR Project Protocol.

Inferring HIV transmission networks from HIV sequences is gaining popularity in the field of HIV molecular epidemiology. However, HIV sequences are often analyzed at distance from those affected by HIV epidemics, namely without the involvement of communities most affected by HIV. These remote analyses often mean that knowledge is generated in absence of lived experiences and socio-economic realities that could inform the ethical application of network-derived information in 'real world' programmes. Procedures to engage communities are noticeably absent from the HIV molecular epidemiology literature. Here we present our team's protocol for engaging community activists living in Nairobi, Kenya in a knowledge exchange process - The CIPHR Project (Community Insights in Phylogenetic HIV Research). Drawing upon a community-based participatory approach, our team will (1) explore the possibilities and limitations of HIV molecular epidemiology for key population programmes, (2) pilot a community-based HIV molecular study, and (3) co-develop policy guidelines on conducting ethically safe HIV molecular epidemiology. Critical dialogue with activist communities will offer insight into the potential uses and abuses of using such information to sharpen HIV prevention programmes. The outcome of this process holds importance to the development of policy frameworks that will guide the next generation of the global response.

Host genetic variation at a locus near CHD1L impacts HIV sequence diversity in a South African population.

IMPORTANCE: It has been previously shown that genetic variants near CHD1L on chromosome 1 are associated with reduced HIV VL in African populations. However, the impact of these variants on viral diversity and how they restrict viral replication are unknown. We report on a regional association analysis in a South African population and show evidence of selective pressure by variants near CHD1L on HIV RT and gag. Our findings provide further insight into how genetic variability at this locus contributes to host control of HIV in a South African population.

Africa-specific human genetic variation near CHD1L associates with HIV-1 load.

HIV-1 remains a global health crisis1, highlighting the need to identify new targets for therapies. Here, given the disproportionate HIV-1 burden and marked human genome diversity in Africa2, we assessed the genetic determinants of control of set-point viral load in 3,879 people of African ancestries living with HIV-1 participating in the international collaboration for the genomics of HIV3. We identify a previously undescribed association signal on chromosome 1 where the peak variant associates with an approximately 0.3 log10-transformed copies per ml lower set-point viral load per minor allele copy and is specific to populations of African descent. The top associated variant is intergenic and lies between a long intergenic non-coding RNA (LINC00624) and the coding gene CHD1L, which encodes a helicase that is involved in DNA repair4. Infection assays in iPS cell-derived macrophages and other immortalized cell lines showed increased HIV-1 replication in CHD1L-knockdown and CHD1L-knockout cells. We provide evidence from population genetic studies that Africa-specific genetic variation near CHD1L associates with HIV replication in vivo. Although experimental studies suggest that CHD1L is able to limit HIV infection in some cell types in vitro, further investigation is required to understand the mechanisms underlying our observations, including any potential indirect effects of CHD1L on HIV spread in vivo that our cell-based assays cannot recapitulate.

Cytokine/chemokine profiles in people with recent infection by Mycobacterium tuberculosis.

Introduction: The risk of progression to tuberculosis disease is highest within the first year after M. tuberculosis infection (TBI). We hypothesize that people with newly acquired TBI have a unique cytokine/chemokine profile that could be used as a potential biomarker. Methods: We evaluated socio-demographic variables and 18 cytokines/chemokines in plasma samples from a cohort of people deprived of liberty (PDL) in two Colombian prisons: 47 people diagnosed with pulmonary TB, 24 with new TBI, and 47 non-infected individuals. We performed a multinomial regression to identify the immune parameters that differentiate the groups. Results: The concentration of immune parameters changed over time and was affected by the time of incarceration. The concentration of sCD14, IL-18 and IP-10 differed between individuals with new TBI and short and long times of incarceration. Among people with short incarceration, high concentrations of MIP-3α were associated with a higher risk of a new TBI, and higher concentrations of Eotaxin were associated with a lower risk of a new TBI. Higher concentrations of sCD14 and TNF-α were associated with a higher risk of TB disease, and higher concentrations of IL-18 and MCP-1 were associated with a lower risk of TB disease. Conclusions: There were cytokines/chemokines associated with new TBI and TB disease. However, the concentration of immune mediators varies by the time of incarceration among people with new TBI. Further studies should evaluate the changes of these and other cytokines/chemokines over time to understand the immune mechanisms across the spectrum of TB.

Multiancestry sex-stratified genomic associations with HIV viral load and controller status from the ICGH.

Biological sex and host genetics influence HIV pathogenesis. Females have a higher likelihood of spontaneous viral control and lower set point viral load (spVL). No prior studies have assessed sex-specific genetics of HIV. To address this, we performed a sex-stratified genome-wide association study using data from the ICGH. Although it is the largest collection of genomic data in HIV, this multiethnic sample of 9,705 people is 81.3% male. We sought to identify sex-specific genetic variants and genes associated with HIV spVL and control. We confirmed associations in the HLA and CCR5 regions in males and HLA in females. Gene-based analyses detected associations between HIV spVL and PET100, PCP2, XAB2, and STXBP2 only in males. We detected variants with a significant sex-differential effect on spVL in SDC3 and PUM1 (rs10914268) and PSORS1C2 (rs1265159) and on HIV control in SUB1 (rs687659), AL158151.3, PTPA, and IER5L (rs4387067). Those variants have epigenetic and genetic interactions with relevant genes with both cis and trans effects. In summary, we identified sex-shared associations at the single-variant level, sex-specific associations at the gene-based level, and genetic variants with significant differential effects between the sexes.

Gene expression profiling identifies candidate biomarkers for new latent tuberculosis infections. A cohort study.

OBJECTIVE: To determine the gene expression profile in individuals with new latent tuberculosis infection (LTBI), and to compare them with people with active tuberculosis (TB) and those exposed to TB but not infected. DESIGN: A prospective cohort study. Recruitment and follow-up were conducted between September 2016 to December 2018. Gene expression and data processing and analysis from April 2019 to April 2021. SETTING: Two male Colombian prisons. PARTICIPANTS: 15 new tuberculin skin test (TST) converters (negative TST at baseline that became positive during follow-up), 11 people that continued with a negative TST after two years of follow-up, and 10 people with pulmonary ATB. MAIN OUTCOME MEASURES: Gene expression profile using RNA sequencing from PBMC samples. The differential expression was assessed using the DESeq2 package in Bioconductor. Genes with |logFC| >1.0 and an adjusted p-value < 0.1 were differentially expressed. We analyzed the differences in the enrichment of KEGG pathways in each group using InterMiner. RESULTS: The gene expression was affected by the time of incarceration. We identified group-specific differentially expressed genes between the groups: 289 genes in people with a new LTBI and short incarceration (less than three months of incarceration), 117 in those with LTBI and long incarceration (one or more years of incarceration), 26 in ATB, and 276 in the exposed but non-infected individuals. Four pathways encompassed the largest number of down and up-regulated genes among individuals with LTBI and short incarceration: cytokine signaling, signal transduction, neutrophil degranulation, and innate immune system. In individuals with LTBI and long incarceration, the only enriched pathway within up-regulated genes was Emi1 phosphorylation. CONCLUSIONS: Recent infection with MTB is associated with an identifiable RNA pattern related to innate immune system pathways that can be used to prioritize LTBI treatment for those at greatest risk for developing active TB.

Expansion of cytotoxic tissue-resident CD8+ T cells and CCR6+CD161+ CD4+ T cells in the nasal mucosa following mRNA COVID-19 vaccination.

Vaccines against SARS-CoV-2 have shown high efficacy in clinical trials, yet a full immunologic characterization of these vaccines, particularly within the human upper respiratory tract, is less well known. Here, we enumerate and phenotype T cells in nasal mucosa and blood using flow cytometry before and after vaccination with the Pfizer-BioNTech COVID-19 vaccine (n = 21). Tissue-resident memory (Trm) CD8+ T cells expressing CD69+CD103+ increase in number ~12 days following the first and second doses, by 0.31 and 0.43 log10 cells per swab respectively (p = 0.058 and p = 0.009 in adjusted linear mixed models). CD69+CD103+CD8+ T cells in the blood decrease post-vaccination. Similar increases in nasal CD8+CD69+CD103- T cells are observed, particularly following the second dose. CD4+ cells co-expressing CCR6 and CD161 are also increased in abundance following both doses. Stimulation of nasal CD8+ T cells with SARS-CoV-2 spike peptides elevates expression of CD107a at 2- and 6-months (p = 0.0096) post second vaccine dose, with a subset of donors also expressing increased cytokines. These data suggest that nasal T cells may be induced and contribute to the protective immunity afforded by this vaccine.

Increased prevalence of clonal hematopoiesis of indeterminate potential amongst people living with HIV.

People living with human immunodeficiency virus (PLWH) have significantly increased risk for cardiovascular disease in part due to inflammation and immune dysregulation. Clonal hematopoiesis of indeterminate potential (CHIP), the age-related acquisition and expansion of hematopoietic stem cells due to leukemogenic driver mutations, increases risk for both hematologic malignancy and coronary artery disease (CAD). Since increased inflammation is hypothesized to be both a cause and consequence of CHIP, we hypothesized that PLWH have a greater prevalence of CHIP. We searched for CHIP in multi-ethnic cases from the Swiss HIV Cohort Study (SHCS, n = 600) and controls from the Atherosclerosis Risk in the Communities study (ARIC, n = 8111) from blood DNA-derived exome sequences. We observed that HIV is associated with a twofold increase in CHIP prevalence, both in the whole study population and in a subset of 230 cases and 1002 matched controls selected by propensity matching to control for demographic imbalances (SHCS 7%, ARIC 3%, p = 0.005). We also observed that ASXL1 is the most commonly mutated CHIP-associated gene in PLWH. Our results suggest that CHIP may contribute to the excess cardiovascular risk observed in PLWH.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response.

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

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.

The influence of human genetic variation on Epstein-Barr virus sequence diversity.

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Rectal microbiota diversity in Kenyan MSM is inversely associated with frequency of receptive anal sex, independent of HIV status.

OBJECTIVE: Both HIV infection and identifying as MSM have been linked to altered rectal microbiota composition, but few studies have studied sexual behavioural associations with rectal microbiota within MSM. In addition, most rectal microbiota studies in MSM have been limited geographically to Europe and North America, and replication of findings in lower and middle-income countries is lacking. DESIGN: A cross-sectional study. METHODS: We enrolled MSM from Nairobi, Kenya, and determined their HIV/sexually transmitted infection status. Rectal specimens were obtained for 16s rRNA sequencing of the rectal microbiota, and sexual behaviour was characterized using a standardized questionnaire. Microbiome differences were modelled using nonparametric statistics, Bray-Curtis ecological distance metrics and analyses of differential taxa abundance. Multivariable linear regression was used to model HIV status and recent sexual activity as predictors of alpha diversity, controlling for a range of covariates. RESULTS: Alpha diversity was consistently lower in Kenyan HIV-infected MSM (n = 80), including those on antiretroviral therapy (ART) compared with HIV-uninfected MSM. A statistical trend was observed for clustering of HIV status by Prevotella or Bacteroides dominance (P = 0.13). Several taxa were enriched in HIV-positive men, including Roseburia, Lachnospira, Streptococcus and Granulicatella. Receptive anal sex with several types of sexual partners (paying, regular, casual) was associated with lower Chao1 and Simpson diversity, independent of HIV status, while HIV infection was associated lower Chao1 (P = 0.030) but not Simpson diversity (P = 0.49). CONCLUSION: Both HIV infection and sexual behaviour were associated with rectal microflora alpha diversity, in particular richness, but not Prevotella spp. dominance, in Kenyan MSM. Associations were more robust for sexual behaviour.

Dissemination of Mycobacterium tuberculosis is associated to a SIGLEC1 null variant that limits antigen exchange via trafficking extracellular vesicles.

The identification of individuals with null alleles enables studying how the loss of gene function affects infection. We previously described a non-functional variant in SIGLEC1, which encodes the myeloid-cell receptor Siglec-1/CD169 implicated in HIV-1 cell-to-cell transmission. Here we report a significant association between the SIGLEC1 null variant and extrapulmonary dissemination of Mycobacterium tuberculosis (Mtb) in two clinical cohorts comprising 6,256 individuals. Local spread of bacteria within the lung is apparent in Mtb-infected Siglec-1 knockout mice which, despite having similar bacterial load, developed more extensive lesions compared to wild type mice. We find that Siglec-1 is necessary to induce antigen presentation through extracellular vesicle uptake. We postulate that lack of Siglec-1 delays the onset of protective immunity against Mtb by limiting antigen exchange via extracellular vesicles, allowing for an early local spread of mycobacteria that increases the risk for extrapulmonary dissemination.