A cohort-based study of host gene expression: tumor suppressor and innate immune/inflammatory pathways associated with the HIV reservoir size
0
by Ashok K. Dwivedi, Germán G. Gornalusse, David A. Siegel, Alton Barbehenn, Cassandra Thanh, Rebecca Hoh, Kristen S. Hobbs, Tony Pan, Erica A. Gibson, Jeffrey Martin, Frederick Hecht, Christopher Pilcher, Jeffrey Milush, Michael P. Busch, Mars Stone, Meei-Li Huang, Julieta Reppetti, Phuong M. Vo, Claire N. Levy, Pavitra Roychoudhury, Keith R. Jerome, Florian Hladik, Timothy J. Henrich, Steven G. Deeks, Sulggi A. Lee
The major barrier to an HIV cure is the HIV reservoir: latently-infected cells that persist despite effective antiretroviral therapy (ART). There have been few cohort-based studies evaluating host genomic or transcriptomic predictors of the HIV reservoir. We performed host RNA sequencing and HIV reservoir quantification (total DNA [tDNA], unspliced RNA [usRNA], intact DNA) from peripheral CD4+ T cells from 191 ART-suppressed people with HIV (PWH). After adjusting for nadir CD4+ count, timing of ART initiation, and genetic ancestry, we identified two host genes for which higher expression was significantly associated with smaller total DNA viral reservoir size, P3H3 and NBL1, both known tumor suppressor genes. We then identified 17 host genes for which lower expression was associated with higher residual transcription (HIV usRNA). These included novel associations with membrane channel (KCNJ2, GJB2), inflammasome (IL1A, CSF3, TNFAIP5, TNFAIP6, TNFAIP9, CXCL3, CXCL10), and innate immunity (TLR7) genes (FDR-adjusted q<0.05). Gene set enrichment analyses further identified significant associations of HIV usRNA with TLR4/microbial translocation (q = 0.006), IL-1/NRLP3 inflammasome (q = 0.008), and IL-10 (q = 0.037) signaling. Protein validation assays using ELISA and multiplex cytokine assays supported these observed inverse host gene correlations, with P3H3, IL-10, and TNF-α protein associations achieving statistical significance (p<0.05). Plasma IL-10 was also significantly inversely associated with HIV DNA (p = 0.016). HIV intact DNA was not associated with differential host gene expression, although this may have been due to a large number of undetectable values in our study. To our knowledge, this is the largest host transcriptomic study of the HIV reservoir. Our findings suggest that host gene expression may vary in response to the transcriptionally active reservoir and that changes in cellular proliferation genes may influence the size of the HIV reservoir. These findings add important data to the limited host genetic HIV reservoir studies to date.