Common genetic variation and the control of HIV-1 in humans

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

HIV Diagnosis and Treatment through Advanced Technologies

Human immunodeficiency virus (HIV) is the chief contributor to global burden of disease. In 2010, HIV was the fifth leading cause of disability-adjusted life years in people of all ages and leading cause for people aged 30-44 years. It is classified as a member of the family Retroviridae and genus Lentivirus based on the biological, morphological, and genetic properties. It infects different cells of the immune system, such as CD4+ T cells (T-helper cells), dendritic cells, and macrophages. HIV has two subtypes: HIV-1 and HIV-2. Among these strains, HIV-1 is the most virulent and pathogenic. Advanced diagnostic methods are exploring new ways of treatment and contributing in the reduction of HIV cases. The diagnostic techniques like PCR, rapid test, EIA, p24 antigen, and western blot have markedly upgraded the diagnosis of HIV. Antiretroviral therapy and vaccines are promising candidates in providing therapeutic and preventive regimes, respectively. Invention of CRISPR/Cas9 is a breakthrough in the field of HIV disease management.

In Vivo and in Vitro Proteome Analysis of Human Immunodeficiency Virus (HIV)-1-infected, Human CD4+ T Cells

Host-directed therapies against HIV-1 are thought to be critical for long term containment of the HIV-1 pandemic but remain elusive. Because HIV-1 infects and manipulates important effectors of both the innate and adaptive immune system, identifying modulations of the host cell systems in humans during HIV-1 infection may be crucial for the development of immune based therapies. Here, we quantified the changes of the proteome in human CD4⁺ T cells upon HIV-1 infection, both in vitro and in vivo A SWATH-MS approach was used to measure the proteome of human primary CD4⁺ T cells infected with HIV-1 in vitro as well as CD4⁺ T cells from HIV-1-infected patients with paired samples on and off antiretroviral treatment. In the in vitro experiment, the proteome of CD4⁺ T cells was quantified over a time course following HIV-1 infection. 1,725 host cell proteins and 4 HIV-1 proteins were quantified, with 145 proteins changing significantly during the time course. Changes in the proteome peaked 24 h after infection, concomitantly with significant HIV-1 protein production. In the in vivo branch of the study, CD4⁺ T cells from viremic patients and those with no detectable viral load after treatment were sorted, and the proteomes were quantified. We consistently detected 895 proteins, 172 of which were considered to be significantly different between the viremic patients and patients undergoing successful treatment. The proteome of the in vitro-infected CD4⁺ T cells was modulated on multiple functional levels, including TLR-4 signaling and the type 1 interferon signaling pathway. Perturbations in the type 1 interferon signaling pathway were recapitulated in CD4⁺ T cells from patients. The study shows that proteome maps generated by SWATH-MS indicate a range of functionally significant changes in the proteome of HIV-infected human CD4⁺ T cells. Exploring these perturbations in more detail may help identify new targets for immune based interventions.

Activating Killer Immunoglobulin Receptors and HLA-C: a successful combination providing HIV-1 control

Several studies demonstrated a relevant role of polymorphisms located within the HLA-B and -C loci and the Killer Immunoglobulin Receptors (KIRs) 3DL1 and 3DS1 in controlling HIV-1 replication. KIRs are regulatory receptors expressed at the surface of NK and CD8+ T-cells that specifically bind HLA-A and -B alleles belonging to the Bw4 supratype and all the -C alleles expressing the C1 or C2 supratype. We here disclose a novel signature associated with the Elite Controller but not with the long-term nonprogressor status concerning 2DS activating KIRs and HLA-C2 alleles insensitive to miRNA148a regulation. Overall, our findings support a crucial role of NK cells in the control of HIV-1 viremia.

Estimating the Respective Contributions of Human and Viral Genetic Variation to HIV Control

We evaluated the fraction of variation in HIV-1 set point viral load attributable to viral or human genetic factors by using joint host/pathogen genetic data from 541 HIV infected individuals. We show that viral genetic diversity explains 29% of the variation in viral load while host factors explain 8.4%. Using a joint model including both host and viral effects, we estimate a total of 30% heritability, indicating that most of the host effects are reflected in viral sequence variation.

Viral loads of Human Immunodeficiency Virus infections are correlated between the donor and the recipient of the transmission pair. Similarly, human genetic factors may modulate viral load. In this study we estimate the extents to which viral load is heritable either via the viral genotype (from donor to recipient) or via the host’s Human Leukocyte Antigen (HLA) genotype. We find that a major fraction of inter individual variability is explained by the similarity of the viral genotypes, and that human genetic variation in the HLA region provide little additional explanatory power.

RICH2 is implicated in viraemic control of HIV-1 in black South African individuals

An intronic single nucleotide polymorphism (SNP) in RICH2 (rs2072255; 255ⁱ), in complete linkage disequilibrium (LD) with an exonic SNP (rs2072254; 254^e), has been identified in a genome wide association study to be associated with progression to AIDS in Caucasian individuals. RICH2 links tetherin to the cortical actin network and the RICH2/tetherin interaction has been shown to be important for the downstream activation of NF-κβ and the consequential promotion of proinflammatory responses. We investigated the role of these two SNPs in natural control of HIV-1 in black South Africans including healthy controls (HCs; N=102) and antiretroviral-naive HIV-1-infected controllers (HICs; N=52) and progressors (N=74). HICs were stratified as elite controllers (ECs; N=11), viraemic controllers (VCs; N=30), high viral load (VL) long term non-progressors (HVL LTNPs; N=11) and also according to VL<400RNA copies/ml (HICs VL<400; N=20) and VL>400RNA copies/ml (HICs VL>400; N=32). Results showed that in contrast to Caucasians who had very strong LD between these SNPs (r²=0.97), black populations exhibited low LD (r²=0.11-0.27), however a 254^e minor allele was always present with a 255ⁱ minor allele but not vice versa. The SNPs did not show significant over- or underrepresentation in any particular group, however the combination of 254^e major allele homozygosity and 255ⁱ heterozygosity (254^eAA/255ⁱGA) was underrepresented in HICs (OR=3.26; P=0.04) and VCs (OR=7.77; P=0.02) compared to HCs, and in HICs VL>400 compared to both HCs (P=0.002) and progressors (P=0.02). A lower CD4⁺ T-cell count was associated with 254^eAA/255ⁱGA and 255ⁱ (GA+AA) in the total HIV-1-infected group (P=0.043) and progressors (P=0.017), respectively. In silico analysis predicted loss of an exonic splice enhancer site with the 254^e-G allele. We postulate that altered splicing of RICH2 will affect levels of RICH2 expression and consequently NF-κβ activation. These findings point to a role for RICH2 and tetherin in viraemic natural control of HIV-1.

Increased T cell breadth and antibody response elicited in prime-boost regimen by viral vector encoded homologous SIV Gag/Env in outbred CD1 mice

BACKGROUND

A major obstacle for the development of HIV vaccines is the virus' worldwide sequence diversity. Nevertheless, the presence of T cell epitopes within conserved regions of the virus' structural Gag protein and conserved structures in the envelope (env) sequence raises the possibility that cross-reactive responses may be induced by vaccination. In this study, the aim was to investigate the importance of antigenic match on immunodominance and breadth of obtainable T cell responses.

METHODS

Outbred CD1 mice were immunized with either heterologous (SIVmac239 and HIV-1 clade B consensus) or homologous (SIVmac239) gag sequences using adenovirus (Ad5) and MVA vectors. Env (SIVmac239) was co-encoded in the vectors to study the induction of antibodies, which is a primary target of current HIV vaccine designs. All three vaccines were designed as virus-encoded virus-like particle vaccines. Antibody responses were analysed by ELISA, avidity ELISA, and neutralization assay. T cell responses were determined by intracellular cytokine staining of splenocytes.

RESULTS

The homologous Env/Gag prime-boost regimen induced higher Env binding antibodies, and induced stronger and broader Gag specific CD8+ T cell responses than the homologous Env/heterologous Gag prime-boost regimen. Homologous Env/heterologous Gag immunization resulted in selective boosting of Env specific CD8+ T cell responses and consequently a paradoxical decreased recognition of variant sequences including conserved elements of p24 Gag.

CONCLUSIONS

These results contrast with related studies using Env or Gag as the sole antigen and suggest that prime-boost immunizations based on homologous SIVmac239 Gag inserts is an efficient component of genetic VLP vaccines-both for induction of potent antibody responses and cross-reactive CD8+ T cell responses.

Clinical Applications of Genome Editing to HIV Cure

Despite significant advances in HIV drug treatment regimens, which grant near-normal life expectancies to infected individuals who have good virological control, HIV infection itself remains incurable. In recent years, novel gene- and cell-based therapies have gained increasing attention due to their potential to provide a functional or even sterilizing cure for HIV infection with a one-shot treatment. A functional cure would keep the infection in check and prevent progression to AIDS, while a sterilizing cure would eradicate all HIV viruses from the patient. Genome editing is the most precise form of gene therapy, able to achieve permanent genetic disruption, modification, or insertion at a predesignated genetic locus. The most well-studied candidate for anti-HIV genome editing is CCR5, an essential coreceptor for the majority of HIV strains, and the lack of which confers HIV resistance in naturally occurring homozygous individuals. Genetic disruption of CCR5 to treat HIV has undergone clinical testing, with seven completed or ongoing trials in T cells and hematopoietic stem and progenitor cells, and has shown promising safety and potential efficacy profiles. Here we summarize clinical findings of CCR5 editing for HIV therapy, as well as other genome editing-based approaches under pre-clinical development. The anticipated development of more sophisticated genome editing technologies should continue to benefit HIV cure efforts.

Wide variation in susceptibility of transmitted/founder HIV-1 subtype C Isolates to protease inhibitors and association with in vitro replication efficiency

The gag gene is highly polymorphic across HIV-1 subtypes and contributes to susceptibility to protease inhibitors (PI), a critical class of antiretrovirals that will be used in up to 2 million individuals as second-line therapy in sub Saharan Africa by 2020. Given subtype C represents around half of all HIV-1 infections globally, we examined PI susceptibility in subtype C viruses from treatment-naïve individuals. PI susceptibility was measured in a single round infection assay of full-length, replication competent MJ4/gag chimeric viruses, encoding the gag gene and 142 nucleotides of pro derived from viruses in 20 patients in the Zambia-Emory HIV Research Project acute infection cohort. Ten-fold variation in susceptibility to PIs atazanavir and lopinavir was observed across 20 viruses, with EC50s ranging 0.71-6.95 nM for atazanvir and 0.64-8.54 nM for lopinavir. Ten amino acid residues in Gag correlated with lopinavir EC50 (p < 0.01), of which 380 K and 389I showed modest impacts on in vitro drug susceptibility. Finally a significant relationship between drug susceptibility and replication capacity was observed for atazanavir and lopinavir but not darunavir. Our findings demonstrate large variation in susceptibility of PI-naïve subtype C viruses that appears to correlate with replication efficiency and could impact clinical outcomes.

CD4/CD8 ratio and CD8 counts predict CD4 response in HIV-1-infected drug naive and in patients on cART

Plasma HIV viral load is related to declining CD4 lymphocytes. The extent to which CD8 cells, in addition to RNA viral load, predict the depletion of CD4 cells is not well characterized so far. We examine if CD8 cell count is a prognostic factor for CD4 cell counts during an HIV infection.A longitudinal analysis is conducted using data from the Swiss HIV cohort study collected between January 2000 and October 2014. Linear mixed regression models were applied to observations from HIV-1-infected treatment naive patients (NAIVE) and cART-treated patients to predict the short-term evolution of CD4 cell counts. For each subgroup, it was quantified to which extent CD8 cell counts or CD4/CD8 ratios are prognostic factors for disease progression.In both subgroups, 2500 NAIVE and 8902 cART patients, past CD4 cells are positively (P < 0.0001) and past viral load is negatively (P < 0.0001) associated with the outcome. Including additionally past CD8 cell counts improves the fit significantly (P < 0.0001) and increases the marginal explained variation 31.7% to 40.7% for the NAIVE and from 44.1% to 50.7% for the cART group. The past CD4/CD8 ratio (instead of the past CD8 level) is positively associated with the outcome, increasing the explained variation further to 41.8% for NAIVE and 51.9% for cART.

A genome-wide association meta-analysis of diarrhoeal disease in young children identifies FUT2 locus and provides plausible biological pathways

More than a million childhood diarrhoeal episodes occur worldwide each year, and in developed countries a considerable part of them are caused by viral infections. In this study, we aimed to search for genetic variants associated with diarrhoeal disease in young children by meta-analyzing genome-wide association studies, and to elucidate plausible biological mechanisms. The study was conducted in the context of the Early Genetics and Lifecourse Epidemiology (EAGLE) consortium. Data about diarrhoeal disease in two time windows (around 1 year of age and around 2 years of age) was obtained via parental questionnaires, doctor interviews or medical records. Standard quality control and statistical tests were applied to the 1000 Genomes imputed genotypic data. The meta-analysis (N = 5758) followed by replication (N = 3784) identified a genome-wide significant association between rs8111874 and diarrhoea at age 1 year. Conditional analysis suggested that the causal variant could be rs601338 (W154X) in the FUT2 gene. Children with the A allele, which results in a truncated FUT2 protein, had lower risk of diarrhoea. FUT2 participates in the production of histo-blood group antigens and has previously been implicated in the susceptibility to infections, including Rotavirus and Norovirus Gene-set enrichment analysis suggested pathways related to the histo-blood group antigen production, and the regulation of ion transport and blood pressure. Among others, the gastrointestinal tract, and the immune and neuro-secretory systems were detected as relevant organs. In summary, this genome-wide association meta-analysis suggests the implication of the FUT2 gene in diarrhoeal disease in young children from the general population.

CD8+ T Cell Breadth and Ex Vivo Virus Inhibition Capacity Distinguish between Viremic Controllers with and without Protective HLA Class I Alleles

The mechanisms of viral control and loss of viral control in chronically infected individuals with or without protective HLA class I alleles are not fully understood. We therefore characterized longitudinally the immunological and virological features that may explain divergence in disease outcome in 70 HIV-1 C-clade-infected antiretroviral therapy (ART)-naive South African adults, 35 of whom possessed protective HLA class I alleles. We demonstrate that, over 5 years of longitudinal study, 35% of individuals with protective HLA class I alleles lost viral control compared to none of the individuals without protective HLA class I alleles (P = 0.06). Sustained HIV-1 control in patients with protective HLA class I alleles was characteristically related to the breadth of HIV-1 CD8(+) T cell responses against Gag and enhanced ability of CD8(+) T cells to suppress viral replication ex vivo In some cases, loss of virological control was associated with reduction in the total breadth of CD8(+) T cell responses in the absence of differences in HIV-1-specific CD8(+) T cell polyfunctionality or proliferation. In contrast, viremic controllers without protective HLA class I alleles possessed reduced breadth of HIV-1-specific CD8(+) T cell responses characterized by reduced ability to suppress viral replication ex vivo These data suggest that the control of HIV-1 in individuals with protective HLA class I alleles may be driven by broad CD8(+) T cell responses with potent viral inhibitory capacity while control among individuals without protective HLA class I alleles may be more durable and mediated by CD8(+) T cell-independent mechanisms.

IMPORTANCE

Host mechanisms of natural HIV-1 control are not fully understood. In a longitudinal study of antiretroviral therapy (ART)-naive individuals, we show that those with protective HLA class I alleles subsequently experienced virologic failure compared to those without protective alleles. Among individuals with protective HLA class I alleles, viremic control was associated with broad CD8(+) T cells that targeted the Gag protein, and CD8(+) T cells from these individuals exhibited superior virus inhibition capacity. In individuals without protective HLA class I alleles, HIV-1-specific CD8(+) T cell responses were narrow and poorly inhibited virus replication. These results suggest that broad, highly functional cytotoxic T cells (cytotoxic T lymphocytes [CTLs]) against the HIV-1 Gag protein are associated with control among those with protective HLA class I alleles and that loss of these responses eventually leads to viremia. A subset of individuals appears to have alternative, non-CTL mechanisms of viral control. These controllers may hold the key to an effective HIV vaccine.

Models of Viral Population Dynamics

Models of viral population dynamics have contributed enormously to our understanding of the pathogenesis and transmission of several infectious diseases, the coevolutionary dynamics of viruses and their hosts, the mechanisms of action of drugs, and the effectiveness of interventions. In this chapter, we review major advances in the modeling of the population dynamics of the human immunodeficiency virus (HIV) and briefly discuss adaptations to other viruses.

Association between KIR genotypes and HLA-B alleles on viral load in Southern Brazilian individuals infected by HIV-1 subtypes B and C

There is a great variety of HIV-1 subtypes circulating in Brazil, including subtype C, whose prevalence is on the rise, particularly in the southern region. Many host and viral genetic factors may be involved in this trend. We evaluated the influence of human leukocyte antigen (HLA) class I alleles and killer-cell immunoglobulin-like receptor (KIR) genotypes on viral set point and T-CD4(+) parameters in 84 treatment-naïve HIV-1-positive individuals. Frequency data in the infected group were compared to data of 548 healthy control subjects. Individuals with the KIR AA genotype had a higher viral load (VL) than individuals with the KIR Bx genotype. The HIV-1 group was subdivided into three subgroups according to HLA-B allele presence: those with protection to disease alleles (HLA-B(+)), accelerated disease progression alleles (HLA-B(-)), or neither (HLA-B(o)) were grouped. We observed a significant effect of the HLA-B allele presence on VL. The HLA-B(+) group had significantly lower VL than the HLA-B(-) group and trended toward a lower VL than the HLA-B(o) group. There were significant differences between groups expressing extreme VL values: KIR-AA+HLA-B(-) vs. KIR Bx+HLA-B(+) and KIR-AA+HLA-B(o)vs. KIR Bx+HLA-B(+). The relationship of KIR/HLA host genetics with slow HIV disease progression in southern Brazil may be useful for vaccine developers, epidemiologists, and clinicians.

Role of APOBEC3F Gene Variation in HIV-1 Disease Progression and Pneumocystis Pneumonia

Human APOBEC3 cytidine deaminases are intrinsic resistance factors to HIV-1. However, HIV-1 encodes a viral infectivity factor (Vif) that degrades APOBEC3 proteins. In vitro APOBEC3F (A3F) anti-HIV-1 activity is weaker than A3G but is partially resistant to Vif degradation unlike A3G. It is unknown whether A3F protein affects HIV-1 disease in vivo. To assess the effect of A3F gene on host susceptibility to HIV- acquisition and disease progression, we performed a genetic association study in six well-characterized HIV-1 natural cohorts. A common six-Single Nucleotide Polymorphism (SNP) haplotype of A3F tagged by a codon-changing variant (p. I231V, with allele (V) frequency of 48% in European Americans) was associated with significantly lower set-point viral load and slower rate of progression to AIDS (Relative Hazards (RH) = 0.71, 95% CI: 0.56, 0.91) and delayed development of pneumocystis pneumonia (PCP) (RH = 0.53, 95% CI: 0.37-0.76). A validation study in the International Collaboration for the Genomics of HIV (ICGH) showed a consistent association with lower set-point viral load. An in vitro assay revealed that the A3F I231V variant may influence Vif mediated A3F degradation. Our results provide genetic epidemiological evidence that A3F modulates HIV-1/AIDS disease progression.

New challenges and promises in solid organ transplantation pharmacogenetics: the genetic variability of proteins involved in the pharmacodynamics of immunosuppressive drugs

Interindividual variability in immunosuppressive drug responses might be partly explained by genetic variants in proteins involved in the immune response or associated with IS pharmacodynamics. On a general basis, the pharmacogenetics of drug target proteins is less known and understood than that of proteins involved in drug disposition pathways. The aim of this review is to facilitate research related to the pharmacodynamics of the main immunosuppressive drugs used in solid organ transplantation. We elaborated a quality of evidence grading system based on a literature review and identified 'highly recommended', 'recommended' or 'potential' candidates for further research. It is likely that a number of additional rare variants might further explain drug response phenotypes in transplantation, and particularly the most severe ones. The advent of next-generation sequencing will help to identify those variants.

Role of HLA Adaptation in HIV Evolution

Killing of HIV-infected cells by CD8⁺ T-cells imposes strong selection pressure on the virus toward escape. The HLA class I molecules that are successful in mediating some degree of control over the virus are those that tend to present epitopes in conserved regions of the proteome, such as in p24 Gag, in which escape also comes at a significant cost to viral replicative capacity (VRC). In some instances, compensatory mutations can fully correct for the fitness cost of such an escape variant; in others, correction is only partial. The consequences of these events within the HIV-infected host, and at the population level following transmission of escape variants, are discussed. The accumulation of escape mutants in populations over the course of the epidemic already shows instances of protective HLA molecules losing their impact, and in certain cases, a modest decline in HIV virulence in association with population-level increase in mutants that reduce VRC.

The Characteristics of Heterozygous Protein Truncating Variants in the Human Genome

Sequencing projects have identified large numbers of rare stop-gain and frameshift variants in the human genome. As most of these are observed in the heterozygous state, they test a gene’s tolerance to haploinsufficiency and dominant loss of function. We analyzed the distribution of truncating variants across 16,260 autosomal protein coding genes in 11,546 individuals. We observed 39,893 truncating variants affecting 12,062 genes, which significantly differed from an expectation of 12,916 genes under a model of neutral de novo mutation (p<10⁻⁴). Extrapolating this to increasing numbers of sequenced individuals, we estimate that 10.8% of human genes do not tolerate heterozygous truncating variants. An additional 10 to 15% of truncated genes may be rescued by incomplete penetrance or compensatory mutations, or because the truncating variants are of limited functional impact. The study of protein truncating variants delineates the essential genome and, more generally, identifies rare heterozygous variants as an unexplored source of diversity of phenotypic traits and diseases.

Genome sequencing provides evidence for large numbers of putative protein truncating variants in humans. Most truncating variants are only observed in few individuals but are collectively prevalent and widely distributed across the coding genome. Most of the truncating variants are so rare that they are only observed in heterozygosis. The current study identifies 10% of genes where heterozygous truncations are not observed and describes their biological characteristics. In addition, for genes where rare truncations are observed, we argue that these are an unexplored source of diversity of phenotypic traits and diseases.

Polymorphisms of large effect explain the majority of the host genetic contribution to variation of HIV-1 virus 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.

Identification of Genes Whose Expression Profile Is Associated with Non-Progression towards AIDS Using eQTLs

BACKGROUND

Many genome-wide association studies have been performed on progression towards the acquired immune deficiency syndrome (AIDS) and they mainly identified associations within the HLA loci. In this study, we demonstrate that the integration of biological information, namely gene expression data, can enhance the sensitivity of genetic studies to unravel new genetic associations relevant to AIDS.

METHODS

We collated the biological information compiled from three databases of expression quantitative trait loci (eQTLs) involved in cells of the immune system. We derived a list of single nucleotide polymorphisms (SNPs) that are functional in that they correlate with differential expression of genes in at least two of the databases. We tested the association of those SNPs with AIDS progression in two cohorts, GRIV and ACS. Tests on permuted phenotypes of the GRIV and ACS cohorts or on randomised sets of equivalent SNPs allowed us to assess the statistical robustness of this method and to estimate the true positive rate.

RESULTS

Eight genes were identified with high confidence (p = 0.001, rate of true positives 75%). Some of those genes had previously been linked with HIV infection. Notably, ENTPD4 belongs to the same family as CD39, whose expression has already been associated with AIDS progression; while DNAJB12 is part of the HSP90 pathway, which is involved in the control of HIV latency. Our study also drew our attention to lesser-known functions such as mitochondrial ribosomal proteins and a zinc finger protein, ZFP57, which could be central to the effectiveness of HIV infection. Interestingly, for six out of those eight genes, down-regulation is associated with non-progression, which makes them appealing targets to develop drugs against HIV.

Predicting functional regulatory SNPs in the human antimicrobial peptide genes DEFB1 and CAMP in tuberculosis and HIV/AIDS

Single nucleotide polymorphisms (SNPs) in transcription factor binding sites (TFBSs) within gene promoter region or enhancers can modify the transcription rate of genes related to complex diseases. These SNPs can be called regulatory SNPs (rSNPs). Data compiled from recent projects, such as the 1000 Genomes Project and ENCODE, has revealed essential information used to perform in silico prediction of the molecular and biological repercussions of SNPs within TFBS. However, most of these studies are very limited, as they only analyze SNPs in coding regions or when applied to promoters, and do not integrate essential biological data like TFBSs, expression profiles, pathway analysis, homotypic redundancy (number of TFBSs for the same TF in a region), chromatin accessibility and others, which could lead to a more accurate prediction. Our aim was to integrate different data in a biologically coherent method to analyze the proximal promoter regions of two antimicrobial peptide genes, DEFB1 and CAMP, that are associated with tuberculosis (TB) and HIV/AIDS. We predicted SNPs within the promoter regions that are more likely to interact with transcription factors (TFs). We also assessed the impact of homotypic redundancy using a novel approach called the homotypic redundancy weight factor (HWF). Our results identified 10 SNPs, which putatively modify the binding affinity of 24 TFs previously identified as related to TB and HIV/AIDS expression profiles (e.g. KLF5, CEBPA and NFKB1 for TB; FOXP2, BRCA1, CEBPB, CREB1, EBF1 and ZNF354C for HIV/AIDS; and RUNX2, HIF1A, JUN/AP-1, NR4A2, EGR1 for both diseases). Validating with the OregAnno database and cell-specific functional/non functional SNPs from additional 13 genes, our algorithm performed 53% sensitivity and 84.6% specificity to detect functional rSNPs using the DNAseI-HUP database. We are proposing our algorithm as a novel in silico method to detect true functional rSNPs in antimicrobial peptide genes. With further improvement, this novel method could be applied to other promoters in order to design probes and to discover new drug targets for complex diseases.

Low-cost simultaneous detection of CCR5-delta32 and HLA-B*5701 alleles in human immunodeficiency virus type 1 infected patients by selective multiplex endpoint PCR

Host genetic traits impact susceptibility to human immunodeficiency virus type 1 (HIV-1) infection, disease progression as well as antiretroviral drug pharmacokinetics and toxicity. Remarkable examples include a 32-bp deletion in the CCR5 coreceptor molecule (CCR5-delta32) impairing attachment of monocytotropic HIV-1 to the host cell membrane and the HLA-B*5701 allele, strongly associated with a potentially fatal hypersensitivity reaction triggered by abacavir, a nucleoside inhibitor of HIV reverse transcriptase. We developed a simple selective multiplex endpoint PCR method for simultaneous analysis of both genetic traits. Two primers were designed for amplification of a region surrounding the CCR5 32-bp deletion site. One common forward primer and two reverse primers with different 3' termini targeting the HLA-B*570101 and HLA-B*570102 alleles were designed for HLA-B*5701 analysis. A panel of 110 reference DNA samples typed in the HLA-B locus was used for development and blind validation of the assay. All the 45 HLA-B*5701 positive and the 55 HLA-B*5701 negative samples were correctly identified. The CCR5-delta32 allele was readily detected in 7 samples and did not interfere with detection of HLA-B*5701 while providing an internal amplification control. Multiplex PCR products were easily identified in agarose gels with no background noise. This simple and low-cost end-point selective multiplex PCR can conveniently screen HIV patients for the protective CCR5-delta32 allele and the risk of developing abacavir hypersensitivity reaction.

CD4:CD8 lymphocyte ratio as a quantitative measure of immunologic health in HIV-1 infection: findings from an African cohort with prospective data

In individuals with human immunodeficiency virus type 1 (HIV-1) infection, CD4:CD8 lymphocyte ratio is often recognized as a quantitative outcome that reflects the critical role of both CD4(+) and CD8(+) T-cells in HIV-1 pathogenesis or disease progression. Our work aimed to first establish the dynamics and clinical relevance of CD4:CD8 ratio in a cohort of native Africans and then to examine its association with viral and host factors, including: (i) length of infection, (ii) demographics, (iii) HIV-1 viral load (VL), (iv) change in CD4(+) T-lymphocyte count (CD4 slope), (v) HIV-1 subtype, and (vi) host genetics, especially human leukocyte antigen (HLA) variants. Data from 499 HIV-1 seroconverters with frequent (monthly to quarterly) follow-up revealed that CD4:CD8 ratio was stable in the first 3 years of infection, with a modest correlation with VL and CD4 slope. A relatively normal CD4:CD8 ratio (>1.0) in early infection was associated with a substantial delay in disease progression to severe immunodeficiency (<350 CD4 cells/μl), regardless of other correlates of HIV-1 pathogenesis (adjusted hazards ratio (HR) = 0.43, 95% confidence interval (CI) = 0.29-0.63, P < 0.0001). Low VL (<10,000 copies/ml) and HLA-A*74:01 were the main predictors of CD4:CD8 ratio >1.0, but HLA variants (e.g., HLA-B*57 and HLA-B*81) previously associated with VL and/or CD4 trajectories in eastern and southern Africans had no obvious impact on CD4:CD8 ratio. Collectively, these findings suggest that CD4:CD8 ratio is a robust measure of immunologic health with both clinical and epidemiological implications.

Disease progression despite protective HLA expression in an HIV-infected transmission pair

Background

The precise immune responses mediated by HLA class I molecules such as HLA-B*27:05 and HLA-B*57:01 that protect against HIV disease progression remain unclear. We studied a CRF01_AE clade HIV infected donor-recipient transmission pair in which the recipient expressed both HLA-B*27:05 and HLA-B*57:01.

Results

Within 4.5 years of diagnosis, the recipient had progressed to meet criteria for antiretroviral therapy initiation. We employed ultra-deep sequencing of the full-length virus genome in both donor and recipient as an unbiased approach by which to identify specific viral mutations selected in association with progression. Using a heat map method to highlight differences in the viral sequences between donor and recipient, we demonstrated that the majority of the recipient’s mutations outside of Env were within epitopes restricted by HLA-B*27:05 and HLA-B*57:01, including the well-studied Gag epitopes. The donor, who also expressed HLA alleles associated with disease protection, HLA-A*32:01/B*13:02/B*14:01, showed selection of mutations in parallel with disease progression within epitopes restricted by these protective alleles.

Conclusions

These studies of full-length viral sequences in a transmission pair, both of whom expressed protective HLA alleles but nevertheless failed to control viremia, are consistent with previous reports pointing to the critical role of Gag-specific CD8+ T cell responses restricted by protective HLA molecules in maintaining immune control of HIV infection. The transmission of subtype CRF01_AE clade infection may have contributed to accelerated disease progression in this pair as a result of clade-specific sequence differences in immunodominant epitopes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0179-z) contains supplementary material, which is available to authorized users.

[Prevalence study of the genetic markers associated with slow progression of human inmunodefiency virus type 1 in the Galician population (Northwest of Spain)]

INTRODUCTION

The deletion in the CCR5 gene (CCR5Δ32), the HLA-B*27:05, and polymorphisms rs2395029 and rs9264942 have been associated with slower progression of HIV-1.

METHODS

An analysis was performed on 408 patients on follow-up. The analysis of viral load, CD4+ Tlymphocytes and other clinical variables since the diagnosis of the infection were collected.

RESULTS

The prevalence of the genetic markers rs9264942, CCR5wt/Δ32, rs2395029, HLA-B*27:05 was 17.9%, 11.5%, 7.6%, and 6.4%, respectively. Of all the patients, 354 were classified as progressors and 46 as long-term non-progressors (LTNPs). Except for the HLA-B*27:05 allele, other genetic markers were associated with slower progression: CCR5wt/Δ32 (P=.011) and SNPs rs2395029 and rs9264942 (P<.0001), as well as their association (P<.0001).

CONCLUSION

The prevalence of the HLA-B*57:01 allele was higher than described nationally. No association could be found between the HLA-B*27:05 allele and the presence of slower disease progression.

Protective Effect of HLA-B*5701 and HLA-C -35 Genetic Variants in HIV-Positive Caucasians from Northern Poland

Aim of the Study

Association of two HLA class I variants with HIV-1 pretreatment viremia, CD4+ T cell count at the care-entry and CD4+ T cell nadir.

Methods

414 HIV-positive Caucasians (30% women) aged 19-73 years were genotyped for HLA-C -35 (rs9264942) and HLA-B*5701 variants. HIV-1 viral load, as well as CD4+ T cell count at care-entry and nadir, were compared across alleles, genotypes and haplotypes.

Results

HLA-C -35 C/C genotype was found in 17.6% patients, C/T genotype in 48.1%, and T/T genotype in 34.3% patients. HLA-B*5701 variant was present in 5.8% of studied population. HIV plasma viremia in the group with C allele was significantly lower (p=0.0002) compared to T/T group [mean:4.66 log (SD:1.03) vs. 5.07 (SD:0.85) log HIV-RNA copies/ml, respectively], while CD4+ T cell count at baseline was notably higher among C allele carriers compared to T/T homozygotes [median: 318 (IQR:127-537) cells/μl vs. median: 203 (IQR:55-410) cells/μl, respectively] (p=0.0007). Moreover, CD4+ T cell nadir among patients with C allele [median: 205 (IQR:83.5-390) cells/μl] was significantly higher compared to T/T group [median: 133 (IQR:46-328) cells/μl] (p=0.006). Among cases with HLA-B*5701 allele, significantly lower pretreatment viremia and higher baseline CD4+ T cell count were found (mean: 4.08 [SD: 1.2] vs. mean: 4.84 [SD:0.97] log HIV-RNA copies/ml, p=0.003 and 431 vs. 270 cells/μl, p=0.04, respectively) compared to HLA-B*5701 negative individuals. The lowest viremia (mean: 3.85 log [SD:1.3]) HIV-RNA copies/ml and the highest baseline and nadir CD4+ T cell [median: 476 (IQR:304-682) vs. median: 361 (IQR: 205-574) cells/μl, respectively) were found in individuals with HLA-B*5701(+)/HLA-C –35 C/C haplotype.

Conclusions

HLA-C -35 C and HLA-B*5701 allele exert a favorable effect on the immunological (higher baseline and nadir CD4+ T cell count) and virologic (lower pretreatment HIV viral load) variables. This protective effect is additive for the compound HLA-B*5701(+)/HLA-C -35 C/C haplotype.

Genome-Wide Association Studies of HIV-1 Host Control in Ethnically Diverse Chinese Populations

Genome-wide association studies (GWASs) have revealed several genetic loci associated with HIV-1 outcome following infection (e.g., HLA-C at 6p21.33) in multi-ethnic populations with genetic heterogeneity and racial/ethnic differences among Caucasians, African-Americans, and Hispanics. To systematically investigate the inherited predisposition to modulate HIV-1 infection in Chinese populations, we performed GWASs in three ethnically diverse HIV-infected patients groups (i.e., HAN, YUN, and XIN, N = 538). The reported loci at 6p21.33 was validated in HAN (e.g., rs9264942, P = 0.0018). An independent association signal (rs2442719, P = 7.85 × 10⁻⁷, HAN group) in the same region was observed. Imputation results suggest that haplotype HLA-B*13:02/C*06:02, which can partially account for the GWAS signal, is associated with lower viral load in Han Chinese. Moreover, several novel loci were identified using GWAS approach including the top association signals at 6q13 (KCNQ5, rs947612, P = 2.15 × 10⁻⁶), 6p24.1 (PHACTR1, rs202072, P = 3.8 × 10⁻⁶), and 11q12.3 (SCGB1D4, rs11231017, P = 7.39 × 10⁻⁷) in HAN, YUN, and XIN groups, respectively. Our findings imply shared or specific mechanisms for host control of HIV-1 in ethnically diverse Chinese populations, which may shed new light on individualized HIV/AIDS therapy in China.

Discordant Impact of HLA on Viral Replicative Capacity and Disease Progression in Pediatric and Adult HIV Infection

HLA class I polymorphism has a major influence on adult HIV disease progression. An important mechanism mediating this effect is the impact on viral replicative capacity (VRC) of the escape mutations selected in response to HLA-restricted CD8+ T-cell responses. Factors that contribute to slow progression in pediatric HIV infection are less well understood. We here investigate the relationship between VRC and disease progression in pediatric infection, and the effect of HLA on VRC and on disease outcome in adult and pediatric infection. Studying a South African cohort of >350 ART-naïve, HIV-infected children and their mothers, we first observed that pediatric disease progression is significantly correlated with VRC. As expected, VRCs in mother-child pairs were strongly correlated (p = 0.004). The impact of the protective HLA alleles, HLA-B*57, HLA-B*58:01 and HLA-B*81:01, resulted in significantly lower VRCs in adults (p<0.0001), but not in children. Similarly, in adults, but not in children, VRCs were significantly higher in subjects expressing the disease-susceptible alleles HLA-B*18:01/45:01/58:02 (p = 0.007). Irrespective of the subject, VRCs were strongly correlated with the number of Gag CD8+ T-cell escape mutants driven by HLA-B*57/58:01/81:01 present in each virus (p = 0.0002). In contrast to the impact of VRC common to progression in adults and children, the HLA effects on disease outcome, that are substantial in adults, are small and statistically insignificant in infected children. These data further highlight the important role that VRC plays both in adult and pediatric progression, and demonstrate that HLA-independent factors, yet to be fully defined, are predominantly responsible for pediatric non-progression.

The impact of host genetic variation on infection with HIV-1

The outcome after infection with the human immunodeficiency virus type 1 (HIV-1) is a complex phenotype determined by interactions among the pathogen, the human host and the surrounding environment. An impact of host genetic variation on HIV-1 susceptibility was identified early in the pandemic, with a major role attributed to the genes encoding class I human leukocyte antigens (HLA) and the chemokine receptor CCR5. Studies using genome-wide data sets have underscored the strength of these associations relative to variants located throughout the rest of the genome. However, the extent to which additional polymorphisms influence HIV-1 disease progression, and how much of the variability in outcome can be attributed to host genetics, remain largely unclear. Here we discuss findings concerning the functional impact of associated variants, outline methods for quantifying the host genetic component and examine how available genome-wide data sets may be leveraged to discover gene variants that affect the outcome of HIV-1 infection.

The genetic architecture of the human immune system: a bioresource for autoimmunity and disease pathogenesis

Despite recent discoveries of genetic variants associated with autoimmunity and infection, genetic control of the human immune system during homeostasis is poorly understood. We undertook a comprehensive immunophenotyping approach, analyzing 78,000 immune traits in 669 female twins. From the top 151 heritable traits (up to 96% heritable), we used replicated GWAS to obtain 297 SNP associations at 11 genetic loci, explaining up to 36% of the variation of 19 traits. We found multiple associations with canonical traits of all major immune cell subsets and uncovered insights into genetic control for regulatory T cells. This data set also revealed traits associated with loci known to confer autoimmune susceptibility, providing mechanistic hypotheses linking immune traits with the etiology of disease. Our data establish a bioresource that links genetic control elements associated with normal immune traits to common autoimmune and infectious diseases, providing a shortcut to identifying potential mechanisms of immune-related diseases.

The Interplay Between Host Genetic Variation, Viral Replication, and Microbial Translocation in Untreated HIV-Infected Individuals

Systemic immune activation, a major determinant of human immunodeficiency virus (HIV) disease progression, is the result of a complex interplay between viral replication, dysregulation of the immune system, and microbial translocation due to gut mucosal damage. Although human genetic variants influencing HIV load have been identified, it is unknown how much the host genetic background contributes to interindividual differences in other determinants of HIV pathogenesis such as gut damage and microbial translocation. Using samples and data from 717 untreated participants in the Swiss HIV Cohort Study and a genome-wide association study design, we searched for human genetic determinants of plasma levels of intestinal fatty acid-binding protein (I-FABP/FABP2), a marker of gut damage, and of soluble CD14 (sCD14), a marker of lipopolysaccharide bioactivity and microbial translocation. We also assessed the correlations between HIV load, sCD14, and I-FABP. Although we found no genome-wide significant determinant of the tested plasma markers, we observed strong associations between sCD14 and both HIV load and I-FABP, shedding new light on the relationships between processes that drive progression of untreated HIV infection.

High heritability is compatible with the broad distribution of set point viral load in HIV carriers

Set point viral load in HIV patients ranges over several orders of magnitude and is a key determinant of disease progression in HIV. A number of recent studies have reported high heritability of set point viral load implying that viral genetic factors contribute substantially to the overall variation in viral load. The high heritability is surprising given the diversity of host factors associated with controlling viral infection. Here we develop an analytical model that describes the temporal changes of the distribution of set point viral load as a function of heritability. This model shows that high heritability is the most parsimonious explanation for the observed variance of set point viral load. Our results thus not only reinforce the credibility of previous estimates of heritability but also shed new light onto mechanisms of viral pathogenesis.

Following an initial peak in viremia, the viral load in HIV infected patients settles down to a set point which remains more or less stable during chronic HIV infection. This set point viral load is one of the key factors determining the rate of disease progression. The extent to which it is determined by the virus versus host genetics is thus central to developing a better understanding of disease progression. Here we develop an analytical model that describes the changes of the distribution of set point viral load in the HIV carrier population over a full cycle of transmission. Applying this model to patient data we find that the most parsimonious explanation for the observed large variation of set point viral load across HIV patients is that set point viral load is highly heritable from donors to recipients. This implies that set point viral load is to a considerable extent under the genetic control of the virus.

Transmitted virus fitness and host T cell responses collectively define divergent infection outcomes in two HIV-1 recipients

Control of virus replication in HIV-1 infection is critical to delaying disease progression. While cellular immune responses are a key determinant of control, relatively little is known about the contribution of the infecting virus to this process. To gain insight into this interplay between virus and host in viral control, we conducted a detailed analysis of two heterosexual HIV-1 subtype A transmission pairs in which female recipients sharing three HLA class I alleles exhibited contrasting clinical outcomes: R880F controlled virus replication while R463F experienced high viral loads and rapid disease progression. Near full-length single genome amplification defined the infecting transmitted/founder (T/F) virus proteome and subsequent sequence evolution over the first year of infection for both acutely infected recipients. T/F virus replicative capacities were compared in vitro, while the development of the earliest cellular immune response was defined using autologous virus sequence-based peptides. The R880F T/F virus replicated significantly slower in vitro than that transmitted to R463F. While neutralizing antibody responses were similar in both subjects, during acute infection R880F mounted a broad T cell response, the most dominant components of which targeted epitopes from which escape was limited. In contrast, the primary HIV-specific T cell response in R463F was focused on just two epitopes, one of which rapidly escaped. This comprehensive study highlights both the importance of the contribution of the lower replication capacity of the transmitted/founder virus and an associated induction of a broad primary HIV-specific T cell response, which was not undermined by rapid epitope escape, to long-term viral control in HIV-1 infection. It underscores the importance of the earliest CD8 T cell response targeting regions of the virus proteome that cannot mutate without a high fitness cost, further emphasizing the need for vaccines that elicit a breadth of T cell responses to conserved viral epitopes.

Immuno-epidemiological modeling of HIV-1 predicts high heritability of the set-point virus load, while selection for CTL escape dominates virulence evolution

It has been suggested that HIV-1 has evolved its set-point virus load to be optimized for transmission. Previous epidemiological models and studies into the heritability of set-point virus load confirm that this mode of adaptation within the human population is feasible. However, during the many cycles of replication between infection of a host and transmission to the next host, HIV-1 is under selection for escape from immune responses, and not transmission. Here we investigate with computational and mathematical models how these two levels of selection, within-host and between-host, are intertwined. We find that when the rate of immune escape is comparable to what has been observed in patients, immune selection within hosts is dominant over selection for transmission. Surprisingly, we do find high values for set-point virus load heritability, and argue that high heritability estimates can be caused by the 'footprints' left by differing hosts' immune systems on the virus.

Genome-wide association studies in Africans and African Americans: expanding the framework of the genomics of human traits and disease

Genomic research is one of the tools for elucidating the pathogenesis of diseases of global health relevance and paving the research dimension to clinical and public health translation. Recent advances in genomic research and technologies have increased our understanding of human diseases, genes associated with these disorders, and the relevant mechanisms. Genome-wide association studies (GWAS) have proliferated since the first studies were published several years ago and have become an important tool in helping researchers comprehend human variation and the role genetic variants play in disease. However, the need to expand the diversity of populations in GWAS has become increasingly apparent as new knowledge is gained about genetic variation. Inclusion of diverse populations in genomic studies is critical to a more complete understanding of human variation and elucidation of the underpinnings of complex diseases. In this review, we summarize the available data on GWAS in recent African ancestry populations within the western hemisphere (i.e. African Americans and peoples of the Caribbean) and continental African populations. Furthermore, we highlight ways in which genomic studies in populations of recent African ancestry have led to advances in the areas of malaria, HIV, prostate cancer, and other diseases. Finally, we discuss the advantages of conducting GWAS in recent African ancestry populations in the context of addressing existing and emerging global health conditions.

Host genetic factors associated with symptomatic primary HIV infection and disease progression among Argentinean seroconverters

BACKGROUND

Variants in HIV-coreceptor C-C chemokine receptor type 5 (CCR5) and Human leukocyte antigen (HLA) genes are the most important host genetic factors associated with HIV infection and disease progression. Our aim was to analyze the association of these genetic factors in the presence of clinical symptoms during Primary HIV Infection (PHI) and disease progression within the first year.

METHODS

Seventy subjects diagnosed during PHI were studied (55 symptomatic and 15 asymptomatic). Viral load (VL) and CD4 T-cell count were evaluated. HIV progression was defined by presence of B or C events and/or CD4 T-cell counts <350 cell/mm3. CCR5 haplotypes were characterized by polymerase chain reaction and SDM-PCR-RFLP. HLA-I characterization was performed by Sequencing.

RESULTS

Symptoms during PHI were significantly associated with lower frequency of CCR5-CF1 (1.8% vs. 26.7%, p = 0.006). Rapid progression was significantly associated with higher frequency of CCR5-CF2 (16.7% vs. 0%, p = 0.024) and HLA-A*11 (16.7% vs. 1.2%, p = 0.003) and lower frequency of HLA-C*3 (2.8% vs. 17.5%, p = 0.035). Higher baseline VL was significantly associated with presence of HLA-A*11, HLA-A*24, and absence of HLA-A*31 and HLA-B*57. Higher 6-month VL was significantly associated with presence of CCR5-HHE, HLA-A*24, HLA-B*53, and absence of HLA-A*31 and CCR5-CF1. Lower baseline CD4 T-cell count was significantly associated with presence of HLA-A*24/*33, HLA-B*53, CCR5-CF2 and absence of HLA-A*01/*23 and CCR5-HHA. Lower 6-month CD4 T-cell count was associated with presence of HLA-A*24 and HLA-B*53, and absence of HLA-A*01 and HLA-B*07/*39. Moreover, lower 12-month CD4 T-cell count was significantly associated with presence of HLA-A*33, HLA-B*14, HLA-C*08, CCR5-CF2, and absence of HLA-B*07 and HLA-C*07.

CONCLUSION

Several host factors were significantly associated with disease progression in PHI subjects. Most results agree with previous studies performed in other groups. However, some genetic factor associations are being described for the first time, highlighting the importance of genetic studies at a local level.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22 981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10⁻⁸). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Testing for risk and protective trends in genetic analyses of HIV acquisition

Host genetics studies of HIV-1 acquisition are critically important for the identification of new targets for drug and vaccine development. Analyses of such studies typically focus on pairwise comparisons of three different groups: HIV-1 positive individuals, HIV-1 high-risk seronegative individuals, and population controls. Because there is a clear expectation of how gene frequencies of risk or protective alleles would be ordered in the three groups, we are able to construct a statistical framework that offers a consistent increase in power over a wide-range of the magnitude of risk/protective effects. In this paper, we develop tests that constrain the alternative hypothesis to appropriately reflect risk or protective trends jointly across the three groups and show that they lead to a substantial increase in power over the naive pairwise approach. We develop both likelihood-ratio and score statistics that test for genetic effects across the three groups while constraining the alternative hypothesis to reflect biologically motivated trends of risk or protection. The asymptotic distribution of both statistics (likelihood ratio and score) is derived. We investigate the performance of our approach via extensive simulation studies using a biologically motivated model of HIV-1 acquisition, and find that our proposed approach leads to an increase in power of roughly 10-28%. We illustrate our approach with an analysis of the effect of the CCR5Δ32 mutation on HIV acquisition.

Differential microRNA expression profile between stimulated PBMCs from HIV-1 infected elite controllers and viremic progressors

Background

The emerging relationship between microRNAs (miRNA) and viral-control is a topic of interest in the field of HIV. Host-genome might play an important role in the control of viremia. The aim of this study was to assess the specific miRNA profile that could contribute to the control of HIV replication in Elite Controllers

Results

After adequate normalization, expression profile of 286 human miRNAs (hsa-miR) was evaluated in phytohaemagglutinin-stimulated PBMCs from 29 individuals classified in 4 groups: 8 elite controllers (EC; viral load <50 cp/ml without treatment), 8 viremic progressors (VP; VL>5000 cp/ml without treatment), 8 patients under antiretroviral treatment (ART; VL<200 cp/ml) and 5 uninfected individuals (HIV-) through TaqMan Array Human microRNA Cards v3.0. A differential expression pattern consisting of 23 miRNAs became significantly different when comparing EC and VP. Profiling analysis segregated the population in two different blocks: while EC and HIV- clustered together in the same block (EC/HIV-_block 1), VP and ART individuals clustered together in a second block (VP/ART_block 2). Two inversely expressed miRNA patterns were determined within those two blocks: a set of 4 miRNAs (hsa-miR-221, -27a, -27b and -29b) was up-expressed in EC/HIV-_block and down-expressed in VP/ART_block while 19 miRNAs were down-expressed in block 1 and up-expressed in block 2. Differential miRNAs were successfully validated through individual RT-qPCR assays.

Conclusions

Profile in EC resembled HIV- and differentially clusters with VP and ART. Therefore, differential clustering does not rely on undetectable viremia.

Disentangling human tolerance and resistance against HIV

Title: Human tolerance against HIV An evolutionary ecology perspective on clinical data reveals that human traits can affect how well an individual tolerates HIV infection, and identifies host immunity factors associated with disease tolerance.

In ecology, “disease tolerance” is defined as an evolutionary strategy of hosts against pathogens, characterized by reduced or absent pathogenesis despite high pathogen load. To our knowledge, tolerance has to date not been quantified and disentangled from host resistance to disease in any clinically relevant human infection. Using data from the Swiss HIV Cohort Study, we investigated if there is variation in tolerance to HIV in humans and if this variation is associated with polymorphisms in the human genome. In particular, we tested for associations between tolerance and alleles of the Human Leukocyte Antigen (HLA) genes, the CC chemokine receptor 5 (CCR5), the age at which individuals were infected, and their sex. We found that HLA-B alleles associated with better HIV control do not confer tolerance. The slower disease progression associated with these alleles can be fully attributed to the extent of viral load reduction in carriers. However, we observed that tolerance significantly varies across HLA-B genotypes with a relative standard deviation of 34%. Furthermore, we found that HLA-B homozygotes are less tolerant than heterozygotes. Lastly, tolerance was observed to decrease with age, resulting in a 1.7-fold difference in disease progression between 20 and 60-y-old individuals with the same viral load. Thus, disease tolerance is a feature of infection with HIV, and the identification of the mechanisms involved may pave the way to a better understanding of pathogenesis.

When confronted with pathogens, hosts can either evolve to fight them or learn to live with them. The first of these two strategies is called “resistance” and the second “tolerance”. In the context of HIV, many genes conferring resistance have been identified, but no tolerance genes are known. Using statistical techniques originating from plant ecology, we analyzed data from an HIV cohort to look for differences in tolerance between HIV-infected individuals and tested whether they go hand in hand with genetic differences. We found that younger people are more tolerant to HIV infection. We also observed that individuals who carry two different alleles of HLA-B, an important immunity gene, are more tolerant. These findings add to our understanding of how hosts tolerate infections and could open new avenues for treating infections.

Mitochondrial DNA variation and virologic and immunological HIV outcomes in African Americans

OBJECTIVE

To evaluate the impact of mitochondrial DNA (mtDNA) haplogroups on virologic and immunological outcomes of HIV infection.

DESIGN

HAART-naive African American adolescent participants to the Reaching for Excellence in Adolescent Care and Health study.

METHODS

The mtDNA haplogroups were inferred from sequenced mtDNA hypervariable regions HV1 and HV2 and their predictive value on HIV outcomes were evaluated in linear mixed models, controlled for human leukocyte antigen (HLA)-B27, HLA-B57 and HLA-B35-Px alleles and other covariates.

RESULTS

We report data showing that the mtDNA L2 lineage, a group composed of L2a, L2b and L2e mtDNA haplogroups in the studied population, is significantly associated (beta  = -0.08; Bonferroni-adjusted P = 0.004) with decline of CD4 T cells (median loss of 8 ± 1 cells per month) in HAART-naive HIV-infected individuals of African American descent (n = 133). No significant association (P < 0.05) with set-point viral load was observed with any of the tested mtDNA haplogroups. The present data concur with previous findings in the AIDS Clinical Trials Group study 384, implicating the L2 lineage with slower CD4 T-cell recovery after antiretroviral therapy in African Americans.

CONCLUSIONS

Whereas the L2 lineage showed an association with unfavorable immunological outcomes of HIV infection, its phylogenetic divergence from J and U5a, two lineages associated with accelerated HIV progression in European Americans, raises the possibility that interactions with common nucleus-encoded variants drive HIV progression. Disentangling the effects of mitochondrial and nuclear gene variants on the outcomes of HIV infection is an important step to be taken toward a better understanding of HIV/AIDS pathogenesis and pharmacogenomics.

Ring finger protein 39 genetic variants associate with HIV-1 plasma viral loads and its replication in cell culture

BACKGROUND

The human immunodeficiency virus (HIV-1) exploits host proteins to complete its life cycle. Genome-wide siRNA approaches suggested that host proteins affect HIV-1 replication. However, the results barely overlapped. RING finger protein 39 (RNF39) has been identified from genome-wide association studies. However, its function during HIV-1 replication remains unclear.

METHODS AND RESULTS

We investigated the relationship between common RNF39 genetic variants and HIV-1 viral loads. The effect of RNF39 protein knockdown or overexpression on HIV-1 replication was then investigated in different cell lines. Two genetic variants were associated with HIV-1 viral loads. Patients with the ht1-GG/GG haplotype presented lower RNF39 expression levels and lower HIV-1 viral load. RNF39 knockdown inhibited HIV-1 expression.

CONCLUSIONS

RNF39 protein may be involved in HIV-1 replication as observed in genetic studies on patients with HIV-1 and in in vitro cell cultures.

HLA-C -35kb expression SNP is associated with differential control of β-HPV infection in squamous cell carcinoma cases and controls

A single nucleotide polymorphism (SNP) 35 kb upstream of the HLA-C gene is associated with HLA-C expression, and the high expressing genotype (CC) has been associated with HIV-I control. HLA-C is unique among the classical MHC class I molecules for its role in the control of viral infections and recognition of abnormal or missing self. This immunosurveillance is central to the pathogenesis of non-melanoma skin cancer (NMSC), and of squamous cell carcinoma (SCC) in particular. While sun exposure is a major risk factor for these cancers, cutaneous infections with genus β-HPV have been implicated in the development of SCC. We hypothesized that the high expression HLA-C genotype is associated with β-HPV infections. Therefore, we investigated the association between β-HPV serology and the -35 kb SNP (rs9264942) in a population-based case-control study of 510 SCC cases and 608 controls. Among controls, the high expression -35 kb SNP genotype (CC) reduced the likelihood of positive serology for multiple (≥2) β-HPV infections (OR = 0.49, 95% CI: 0.25-0.97), and β-HPV species 2 infection (OR = 0.43, 95% CI: 0.23-0.79). However, no association with β-HPV status was observed among SCC cases. Our findings suggest that underlying immunogenotype plays an important role in differential control of β-HPV in SCC cases and controls.

The dissection of complex susceptibility to infectious disease: bacterial, viral and parasitic infections

Infectious diseases are the result of the exposure of susceptible hosts to pathogenic microbes. Genetic factors are important determinants of host susceptibility and efforts are being made to establish the molecular identity of such genetic susceptibility variants by genome-wide association studies. Results obtained to date partly confirm already known genetic vulnerabilities, but also point to new and unexpected mechanisms of susceptibility that extend from classical innate and acquired immunity to weaknesses in constitutional resistance. These studies also revealed an overlap in genetic control between infectious disease and other common immune and inflammatory disorders.

Large scale analysis of phenotype-pathway relationships based on GWAS results

The widely used pathway-based approach for interpreting Genome Wide Association Studies (GWAS), assumes that since function is executed through the interactions of multiple genes, different perturbations of the same pathway would result in a similar phenotype. This assumption, however, was not systemically assessed on a large scale. To determine whether SNPs associated with a given complex phenotype affect the same pathways more than expected by chance, we analyzed 368 phenotypes that were studied in >5000 GWAS. We found 216 significant phenotype-pathway associations between 70 of the phenotypes we analyzed and known pathways. We also report 391 strong phenotype-phenotype associations between phenotypes that are affected by the same pathways. While some of these associations confirm previously reported connections, others are new and could shed light on the molecular basis of these diseases. Our findings confirm that phenotype-associated SNPs cluster into pathways much more than expected by chance. However, this is true for <20% (70/368) of the phenotypes. Different types of phenotypes show markedly different tendencies: Virtually all autoimmune phenotypes show strong clustering of SNPs into pathways, while most cancers and metabolic conditions, and all electrophysiological phenotypes, could not be significantly associated with any pathway despite being significantly associated with a large number of SNPs. While this may be due to missing data, it may also suggest that these phenotypes could result only from perturbations of specific genes and not from other perturbations of the same pathway. Further analysis of pathway-associated versus gene-associated phenotypes is, therefore, needed in order to understand disease etiology and in order to promote better drug target selection.