Polymorphisms of large effect explain the majority of the host genetic contribution to variation of HIV-1 virus load

Significance A proportion of the variation in HIV-1 viral load in the infected population is influenced by host genetics. Using a large sample of infected individuals (n = 6,315) with genome-wide genotype data, we sought to map genomic regions that influence HIV viral load and quantify their impact. We identified amino acid positions located in the binding groove of class I HLA proteins (HLA-A and -B) and SNPs in the chemokine (C-C motif) receptor 5 gene region that together explain 14.5% of the observed variation in HIV viral load. Controlling for these signals, we estimated that an additional 5.5% can be explained by common, additive genetic variation. Thus, we demonstrate that common variants of large effect explain the majority of the host genetic component of HIV viral load. Previous genome-wide association studies (GWAS) of HIV-1–infected populations have been underpowered to detect common variants with moderate impact on disease outcome and have not assessed the phenotypic variance explained by genome-wide additive effects. By combining the majority of available genome-wide genotyping data in HIV-infected populations, we tested for association between ∼8 million variants and viral load (HIV RNA copies per milliliter of plasma) in 6,315 individuals of European ancestry. The strongest signal of association was observed in the HLA class I region that was fully explained by independent effects mapping to five variable amino acid positions in the peptide binding grooves of the HLA-B and HLA-A proteins. We observed a second genome-wide significant association signal in the chemokine (C-C motif) receptor (CCR) gene cluster on chromosome 3. Conditional analysis showed that this signal could not be fully attributed to the known protective CCR5Δ32 allele and the risk P1 haplotype, suggesting further causal variants in this region. Heritability analysis demonstrated that common human genetic variation—mostly in the HLA and CCR5 regions—explains 25% of the variability in viral load. This study suggests that analyses in non-European populations and of variant classes not assessed by GWAS should be priorities for the field going forward.