Multistage genomewide association study identifies a locus at 1q41 associated with rate of HIV-1 disease progression to clinical AIDS

Host variation in type I interferon signaling genes (MX1), C-C chemokine receptor type 5 gene, and major histocompatibility complex class I alleles in treated HIV+ noncontrollers predict viral reservoir size

OBJECTIVE

Prior genomewide association studies have identified variation in major histocompatibility complex (MHC) class I alleles and C-C chemokine receptor type 5 gene (CCR5Δ32) as genetic predictors of viral control, especially in 'elite' controllers, individuals who remain virally suppressed in the absence of therapy.

DESIGN

Cross-sectional genomewide association study.

METHODS

We analyzed custom whole exome sequencing and direct human leukocyte antigen (HLA) typing from 202 antiretroviral therapy (ART)-suppressed HIV+ noncontrollers in relation to four measures of the peripheral CD4+ T-cell reservoir: HIV intact DNA, total (t)DNA, unspliced (us)RNA, and RNA/DNA. Linear mixed models were adjusted for potential covariates including age, sex, nadir CD4+ T-cell count, pre-ART HIV RNA, timing of ART initiation, and duration of ART suppression.

RESULTS

Previously reported 'protective' host genetic mutations related to viral setpoint (e.g. among elite controllers) were found to predict smaller HIV reservoir size. The HLA 'protective' B∗57:01 was associated with significantly lower HIV usRNA (q = 3.3 × 10-3), and among the largest subgroup, European ancestry individuals, the CCR5Δ32 deletion was associated with smaller HIV tDNA (P = 4.3 × 10-3) and usRNA (P = 8.7 × 10-3). In addition, genomewide analysis identified several single nucleotide polymorphisms in MX1 (an interferon stimulated gene) that were significantly associated with HIV tDNA (q = 0.02), and the direction of these associations paralleled MX1 gene eQTL expression.

CONCLUSIONS

We observed a significant association between previously reported 'protective' MHC class I alleles and CCR5Δ32 with the HIV reservoir size in noncontrollers. We also found a novel association between MX1 and HIV total DNA (in addition to other interferon signaling relevant genes, PPP1CB, DDX3X). These findings warrant further investigation in future validation studies.

A Novel Multitasking Ant Colony Optimization Method for Detecting Multiorder SNP Interactions

MOTIVATION

Linear or nonlinear interactions of multiple single-nucleotide polymorphisms (SNPs) play an important role in understanding the genetic basis of complex human diseases. However, combinatorial analytics in high-dimensional space makes it extremely challenging to detect multiorder SNP interactions. Most classic approaches can only perform one task (for detecting k-order SNP interactions) in each run. Since prior knowledge of a complex disease is usually not available, it is difficult to determine the value of k for detecting k-order SNP interactions.

METHODS

A novel multitasking ant colony optimization algorithm (named MTACO-DMSI) is proposed to detect multiorder SNP interactions, and it is divided into two stages: searching and testing. In the searching stage, multiple multiorder SNP interaction detection tasks (from 2nd-order to kth-order) are executed in parallel, and two subpopulations that separately adopt the Bayesian network-based K2-score and Jensen-Shannon divergence (JS-score) as evaluation criteria are generated for each task to improve the global search capability and the discrimination ability for various disease models. In the testing stage, the G test statistical test is adopted to further verify the authenticity of candidate solutions to reduce the error rate.

RESULT

Three multiorder simulated disease models with different interaction effects and three real age-related macular degeneration (AMD), rheumatoid arthritis (RA) and type 1 diabetes (T1D) datasets were used to investigate the performance of the proposed MTACO-DMSI. The experimental results show that the MTACO-DMSI has a faster search speed and higher discriminatory power for diverse simulation disease models than traditional single-task algorithms. The results on real AMD data and RA and T1D datasets indicate that MTACO-DMSI has the ability to detect multiorder SNP interactions at a genome-wide scale. Availability and implementation: https://github.com/shouhengtuo/MTACO-DMSI/.

Genome-wide association study reveals genetic variants associated with HIV-1C infection in a Botswana study population

Although there have been many studies of gene variant association with different stages of HIV/AIDS progression in United States and European cohorts, few gene-association studies have assessed genic determinants in sub-Saharan African populations, which have the highest density of HIV infections worldwide. We carried out genome-wide association studies on 766 study participants at risk for HIV-1 subtype C (HIV-1C) infection in Botswana. Three gene associations (AP3B1, PTPRA, and NEO1) were shown to have significant association with HIV-1C acquisition. Each gene association was replicated within Botswana or in the United States-African American or United States-European American AIDS cohorts or in both. Each associated gene has a prior reported influence on HIV/AIDS pathogenesis. Thirteen previously discovered AIDS restriction genes were further replicated in the Botswana cohorts, extending our confidence in these prior AIDS restriction gene reports. This work presents an early step toward the identification of genetic variants associated with and affecting HIV acquisition or AIDS progression in the understudied HIV-1C afflicted Botswana population.

Host Genomics of the HIV-1 Reservoir Size and Its Decay Rate During Suppressive Antiretroviral Treatment

BACKGROUND

The primary hurdle for the eradication of HIV-1 is the establishment of a latent viral reservoir early after primary infection. Here, we investigated the potential influence of human genetic variation on the HIV-1 reservoir size and its decay rate during suppressive antiretroviral treatment.

SETTING

Genome-wide association study and exome sequencing study to look for host genetic determinants of HIV-1 reservoir measurements in patients enrolled in the Swiss HIV Cohort Study, a nation-wide prospective observational study.

METHODS

We measured total HIV-1 DNA in peripheral blood mononuclear cells from study participants, as a proxy for the reservoir size at 3 time points over a median of 5.4 years, and searched for associations between human genetic variation and 2 phenotypic readouts: the reservoir size at the first time point and its decay rate over the study period. We assessed the contribution of common genetic variants using genome-wide genotyping data from 797 patients with European ancestry enrolled in the Swiss HIV Cohort Study and searched for a potential impact of rare variants and exonic copy number variants using exome sequencing data generated in a subset of 194 study participants.

RESULTS

Genome-wide and exome-wide analyses did not reveal any significant association with the size of the HIV-1 reservoir or its decay rate on suppressive antiretroviral treatment.

CONCLUSIONS

Our results point to a limited influence of human genetics on the size of the HIV-1 reservoir and its long-term dynamics in successfully treated individuals.

Minding the gap in HIV host genetics: opportunities and challenges

Genome-wide association studies (GWAS) have been successful in identifying and confirming novel genetic variants that are associated with diverse HIV phenotypes. However, these studies have predominantly focused on European cohorts. HLA molecules have been consistently associated with HIV outcomes, some of which have been found to be population specific, underscoring the need for diversity in GWAS. Recently, there has been a concerted effort to address this gap that leads to health care (disease prevention, diagnosis, treatment) disparities with marginal improvement. As precision medicine becomes more utilized, non-European individuals will be more and more disadvantaged, as the genetic variants identified in genomic research based on European populations may not accurately reflect that of non-European individuals. Leveraging pre-existing, large, multiethnic cohorts, such as the UK Biobank, 23andMe, and the National Institute of Health's All of Us Research Program, can contribute in raising genomic research in non-European populations and ultimately lead to better health outcomes.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

Large-Scale "OMICS" Studies to Explore the Physiopatholgy of HIV-1 Infection

In this review, we present the main large-scale experimental studies that have been performed in the HIV/AIDS field. These “omics” studies are based on several technologies including genotyping, RNA interference, and transcriptome or epigenome analysis. Due to the direct connection with disease evolution, there has been a large focus on genotyping cohorts of well-characterized patients through genome-wide association studies (GWASs), but there have also been several invitro studies such as small interfering RNA (siRNA) interference or transcriptome analyses of HIV-1–infected cells. After describing the major results obtained with these omics technologies—including some with a high relevance for HIV-1 treatment—we discuss the next steps that the community needs to embrace in order to derive new actionable therapeutic or diagnostic targets. Only integrative approaches that combine all big data results and consider their complex interactions will allow us to capture the global picture of HIV molecular pathogenesis. This novel challenge will require large collaborative efforts and represents a huge open field for innovative bioinformatics approaches.

Association of a single nucleotide polymorphism in the ubxn6 gene with long-term non-progression phenotype in HIV-positive individuals

OBJECTIVES

The long-term non-progressors (LTNPs) are a heterogeneous group of HIV-positive individuals characterized by their ability to maintain high CD4⁺ T-cell counts and partially control viral replication for years in the absence of antiretroviral therapy. The present study aims to identify host single nucleotide polymorphisms (SNPs) associated with non-progression in a cohort of 352 individuals.

METHODS

DNA microarrays and exome sequencing were used for genotyping about 240 000 functional polymorphisms throughout more than 20 000 human genes. The allele frequencies of 85 LTNPs were compared with a control population. SNPs associated with LTNPs were confirmed in a population of typical progressors. Functional analyses in the affected gene were carried out through knockdown experiments in HeLa-P4, macrophages and dendritic cells.

RESULTS

Several SNPs located within the major histocompatibility complex region previously related to LTNPs were confirmed in this new cohort. The SNP rs1127888 (UBXN6) surpassed the statistical significance of these markers after Bonferroni correction (q = 2.11 × 10^-6). An uncommon allelic frequency of rs1127888 among LTNPs was confirmed by comparison with typical progressors and other publicly available populations. UBXN6 knockdown experiments caused an increase in CAV1 expression and its accumulation in the plasma membrane. In vitro infection of different cell types with HIV-1 replication-competent recombinant viruses caused a reduction of the viral replication capacity compared with their corresponding wild-type cells expressing UBXN6.

CONCLUSIONS

A higher prevalence of Ala31Thr in UBXN6 was found among LTNPs within its N-terminal region, which is crucial for UBXN6/VCP protein complex formation. UBXN6 knockdown affected CAV1 turnover and HIV-1 replication capacity.

Polymorphisms in the 3'-UTR of SCD5 gene are associated with hepatocellular carcinoma in Korean population

The purpose of the study was to assess the relationship between polymorphisms of the SCD5 and MMP1 gene and hepatocellular carcinoma (HCC). The gene polymorphisms with a minor allele frequency (MAF) > 0.05 were selected eight SNPs (rs6840, rs1065403, rs3821974, and rs3733230 in 3'-UTR; rs4693472, rs3733227, rs1848067, and rs6535374 in intron region) of SCD5 gene and two SNPs (rs1799750 and rs1144393 in promoter region) of MMP1 gene. The genotype of SCD5 and MMP1 gene SNPs were determined by direct sequencing and pyrosequencing, respectively. One hundred sixty-two patients with HCC and two hundred twenty-five control subjects were recruited in Korean male population. In terms of genotype frequencies, SCD5 genotype TC, GA, AG, and CG of rs6840, rs1065403, rs3821974, and rs3733230, respectively were higher in control group than HCC. In addition, these genotype decreased the risk (rs6840; OR 0.55, 95% CI 0.31-0.99; rs1065403; OR 0.46, 95% CI 0.26-0.83; rs3821974; OR 0.56, 95% CI 0.31-0.99; rs3733230; OR 0.62, 95% CI 0.34-1.12) of HCC, which may work as a prevention of HCC. At least one minor allele carrier of SCD5 gene polymorphisms were related to decreased risk of HCC for AFP cut-point levels > 200 or > 400 ng/ml, respectively. Our results indicate that polymorphisms in the 3'-UTR of the SCD5 gene may associated with HCC in the Korean male population.

Molecular epidemiology of infectious diseases

Molecular epidemiology (ME) is a branch of epidemiology developed by merging molecular biology into epidemiological studies. In this paper, the authors try to discuss the ways that molecular epidemiology studies identify infectious diseases' causation and pathogenesis, and unravel infectious agents' sources, reservoirs, circulation pattern, transmission pattern, transmission probability, and transmission order. They bring real-world examples of research works in each area to make each study design more understandable. They also address some research areas and study design aspects that need further attention in future. They close with some thoughts about future directions in this field and emphasize on the need for training competent molecular epidemiology specialists that are capable of dealing with rapid advances in the field.

Host genetic variation and HIV disease: from mapping to mechanism

This review aims to provide a summary of current knowledge of host genetic effects on human immunodeficiency virus (HIV) disease. Mapping of simple single nucleotide polymorphisms (SNP) has been largely successful in HIV, but more complex genetic associations involving haplotypic or epigenetic variation, for example, remain elusive. Mechanistic insights explaining SNP associations are incomplete, but continue to be forthcoming. The number of robust immunogenetic correlates of HIV is modest and their discovery mostly predates the genome-wide era. Nevertheless, genome-wide evaluations have nicely validated the impact of HLA and CCR5 variants on HIV disease, and importantly, made clear the many false positive associations that were previously suggested by studies using the candidate gene approach. We describe how multiple HIV outcome measures such as acquisition, viral control, and immune decline have been studied in adults and in children, but that collectively these identify only the two replicable loci responsible for modifying HIV disease, CCR5, and HLA. Recent heritability estimates in this disease corroborate the modest impact of genetic determinants and their oligogenic nature. While the mechanism of protection afforded by genetic variants that diminish CCR5 expression is clear, new aspects of HLA class I-mediated protection continue to be uncovered. We describe how these genetic findings have enhanced insights into immunobiology, been clinically translated into CCR5 antagonists, allowed prioritization of antigens for vaccination efforts, and identified targets for genome-editing interventions. Finally, we describe how studies of genetically complex parts of the genome using new tools may begin revealing additional correlates.

Psoriasis risk SNPs and their association with HIV-1 control

Human evolution has resulted in selection for genetic polymorphisms beneficial in the defense against pathogens. However, such polymorphisms may have the potential to heighten the risk of autoimmune disease. Here, we investigated whether psoriasis-associated single nucleotide polymorphisms influence host control of HIV-1 infection. We studied psoriasis and viral immune response variants in three HIV-positive cohorts: (1) HIV-1 controllers and non-controllers in the Study of the Consequences of the Protease Inhibitor Era (SCOPE) cohort (n=366), (2) Individuals with primary HIV infection in the Options cohort (n=675), and (3) HIV-positive injection drug users from the Urban Health Study (UHS) (n=987). We found a strong association of two psoriasis MHC variants, rs9264942 and rs3021366, with both HIV-1 controller status and viral load, and identified another Class III MHC variant rs9368699 to be strongly associated with viral load. A number of genetic variants outside the MHC (SOX5, TLR9, SDC4, PROX1, IL12B, TLR4, MBL-2, TYK2, IFIH1) demonstrated nominal significance. Overall, our data suggest that several psoriasis variants within the MHC have a robust impact on HIV-1 control, while variants outside the MHC require further investigation.

Jointly determining significance levels of primary and replication studies by controlling the false discovery rate in two-stage genome-wide association studies

In genome-wide association studies, we normally discover associations between genetic variants and diseases/traits in primary studies, and validate the findings in replication studies. We consider the associations identified in both primary and replication studies as true findings. An important question under this two-stage setting is how to determine significance levels in both studies. In traditional methods, significance levels of the primary and replication studies are determined separately. We argue that the separate determination strategy reduces the power in the overall two-stage study. Therefore, we propose a novel method to determine significance levels jointly. Our method is a reanalysis method that needs summary statistics from both studies. We find the most powerful significance levels when controlling the false discovery rate in the two-stage study. To enjoy the power improvement from the joint determination method, we need to select single nucleotide polymorphisms for replication at a less stringent significance level. This is a common practice in studies designed for discovery purpose. We suggest this practice is also suitable in studies with validation purpose in order to identify more true findings. Simulation experiments show that our method can provide more power than traditional methods and that the false discovery rate is well-controlled. Empirical experiments on datasets of five diseases/traits demonstrate that our method can help identify more associations. The R-package is available at: http://bioinformatics.ust.hk/RFdr.html .

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.

Impact of polymorphisms in the HCP5 and HLA-C, and ZNRD1 genes on HIV viral load

AIMS

Single nucleotide polymorphisms (SNPs) in the human leucocyte antigen (HLA) complex P5 (HCP5), HLA-C, and near the zinc ribbon domain containing 1 (ZNRD1) have been shown to influence viral load (VL) set point in HIV-infected individuals with a known seroconversion onset. We aimed to determine the influence of HCP5 rs2395029, HLA-C rs9264942, and ZNRD1 rs3869068 on VL in antiretroviral-naïve individuals and on time to the first VL<51 copies/ml and on CD4(+) T-cell recovery after initiation of combination antiretroviral therapy (cART).

MATERIAL AND METHODS

We genotyped the rs2395029 (A>C), rs9264942 (T>C), and rs3869068 (C>T) SNPs in 1897 Caucasians from The Danish HIV Cohort Study - a prospective, nationwide, population-based study of HIV-infected individuals in Denmark. General linear models evaluated the effect of SNPs on VL in antiretroviral-naïve individuals 0-18months after diagnosis and on CD4(+) T-cell recovery during cART. Cox proportional hazard regression analysis assessed the association with time to first VL<51 copies/ml. All models were assuming additive genetic effects.

RESULTS

The rs2395029, rs9264942, and rs3869068 minor alleles were associated with lower VL in antiretroviral-naïve individuals (rs2395029: [mean VL (copies/ml)], A/A: 70,795 [61,660-79,433], A/C: 33,884 [19,498-58,884], P=0.002; rs9264942: TT: 81,283 [67,608-97,724], T/C: 63,096 [54,954-75,858], CC: 38,905 [25,119-58,884], P<0.0001; rs3869068, CC: 72,444 [63,096-83,176], C/T: 45,709 [33,113-64,565], TT: 58,884 [20,417-169,824], P=0.01). Moreover, the C-alleles of rs2395029 and rs9264942 were associated with shorter time to VL<51 copies/ml: (HR [95% confidence interval], 1.67 [1.09-1.72], P=0.008; 1.16 [1.06-1.28], P=0.002; 1.30 [1.08-1.53], P=0.005, respectively, adjusted for last VL before cART). None of the SNPs predicted CD4(+) T-cell recovery during cART.

CONCLUSIONS

The minor alleles of rs2395029, rs9264942, and rs3689068 associate with lower VL among antiretroviral-naïve individuals and with shorter time to first VL<51copies/ml during cART even after adjustment for VL before cART.

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.

GWATCH: a web platform for automated gene association discovery analysis

Background

As genome-wide sequence analyses for complex human disease determinants are expanding, it is increasingly necessary to develop strategies to promote discovery and validation of potential disease-gene associations.

Findings

Here we present a dynamic web-based platform – GWATCH – that automates and facilitates four steps in genetic epidemiological discovery: 1) Rapid gene association search and discovery analysis of large genome-wide datasets; 2) Expanded visual display of gene associations for genome-wide variants (SNPs, indels, CNVs), including Manhattan plots, 2D and 3D snapshots of any gene region, and a dynamic genome browser illustrating gene association chromosomal regions; 3) Real-time validation/replication of candidate or putative genes suggested from other sources, limiting Bonferroni genome-wide association study (GWAS) penalties; 4) Open data release and sharing by eliminating privacy constraints (The National Human Genome Research Institute (NHGRI) Institutional Review Board (IRB), informed consent, The Health Insurance Portability and Accountability Act (HIPAA) of 1996 etc.) on unabridged results, which allows for open access comparative and meta-analysis.

Conclusions

GWATCH is suitable for both GWAS and whole genome sequence association datasets. We illustrate the utility of GWATCH with three large genome-wide association studies for HIV-AIDS resistance genes screened in large multicenter cohorts; however, association datasets from any study can be uploaded and analyzed by GWATCH.

Global, regional, and national incidence and mortality for HIV, tuberculosis, and malaria during 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013

BACKGROUND

The Millennium Declaration in 2000 brought special global attention to HIV, tuberculosis, and malaria through the formulation of Millennium Development Goal (MDG) 6. The Global Burden of Disease 2013 study provides a consistent and comprehensive approach to disease estimation for between 1990 and 2013, and an opportunity to assess whether accelerated progress has occured since the Millennium Declaration.

METHODS

To estimate incidence and mortality for HIV, we used the UNAIDS Spectrum model appropriately modified based on a systematic review of available studies of mortality with and without antiretroviral therapy (ART). For concentrated epidemics, we calibrated Spectrum models to fit vital registration data corrected for misclassification of HIV deaths. In generalised epidemics, we minimised a loss function to select epidemic curves most consistent with prevalence data and demographic data for all-cause mortality. We analysed counterfactual scenarios for HIV to assess years of life saved through prevention of mother-to-child transmission (PMTCT) and ART. For tuberculosis, we analysed vital registration and verbal autopsy data to estimate mortality using cause of death ensemble modelling. We analysed data for corrected case-notifications, expert opinions on the case-detection rate, prevalence surveys, and estimated cause-specific mortality using Bayesian meta-regression to generate consistent trends in all parameters. We analysed malaria mortality and incidence using an updated cause of death database, a systematic analysis of verbal autopsy validation studies for malaria, and recent studies (2010-13) of incidence, drug resistance, and coverage of insecticide-treated bednets.

FINDINGS

Globally in 2013, there were 1·8 million new HIV infections (95% uncertainty interval 1·7 million to 2·1 million), 29·2 million prevalent HIV cases (28·1 to 31·7), and 1·3 million HIV deaths (1·3 to 1·5). At the peak of the epidemic in 2005, HIV caused 1·7 million deaths (1·6 million to 1·9 million). Concentrated epidemics in Latin America and eastern Europe are substantially smaller than previously estimated. Through interventions including PMTCT and ART, 19·1 million life-years (16·6 million to 21·5 million) have been saved, 70·3% (65·4 to 76·1) in developing countries. From 2000 to 2011, the ratio of development assistance for health for HIV to years of life saved through intervention was US$4498 in developing countries. Including in HIV-positive individuals, all-form tuberculosis incidence was 7·5 million (7·4 million to 7·7 million), prevalence was 11·9 million (11·6 million to 12·2 million), and number of deaths was 1·4 million (1·3 million to 1·5 million) in 2013. In the same year and in only individuals who were HIV-negative, all-form tuberculosis incidence was 7·1 million (6·9 million to 7·3 million), prevalence was 11·2 million (10·8 million to 11·6 million), and number of deaths was 1·3 million (1·2 million to 1·4 million). Annualised rates of change (ARC) for incidence, prevalence, and death became negative after 2000. Tuberculosis in HIV-negative individuals disproportionately occurs in men and boys (versus women and girls); 64·0% of cases (63·6 to 64·3) and 64·7% of deaths (60·8 to 70·3). Globally, malaria cases and deaths grew rapidly from 1990 reaching a peak of 232 million cases (143 million to 387 million) in 2003 and 1·2 million deaths (1·1 million to 1·4 million) in 2004. Since 2004, child deaths from malaria in sub-Saharan Africa have decreased by 31·5% (15·7 to 44·1). Outside of Africa, malaria mortality has been steadily decreasing since 1990.

INTERPRETATION

Our estimates of the number of people living with HIV are 18·7% smaller than UNAIDS's estimates in 2012. The number of people living with malaria is larger than estimated by WHO. The number of people living with HIV, tuberculosis, or malaria have all decreased since 2000. At the global level, upward trends for malaria and HIV deaths have been reversed and declines in tuberculosis deaths have accelerated. 101 countries (74 of which are developing) still have increasing HIV incidence. Substantial progress since the Millennium Declaration is an encouraging sign of the effect of global action.

FUNDING

Bill & Melinda Gates Foundation.

Host genetics and immune control of HIV-1 infection: fine mapping for the extended human MHC region in an African cohort

Multiple MHC loci encoding human leukocyte antigens (HLA) have allelic variants unequivocally associated with differential immune control of HIV-1 infection. Fine mapping based on single nucleotide polymorphisms (SNPs) in the extended MHC (xMHC) region is expected to reveal causal or novel factors and to justify a search for functional mechanisms. We have tested the utility of a custom fine-mapping platform (the ImmunoChip) for 172 HIV-1 seroconverters (SCs) and 449 seroprevalent individuals (SPs) from Lusaka, Zambia, with a focus on more than 6,400 informative xMHC SNPs. When conditioned on HLA and non-genetic factors previously associated with HIV-1 viral load (VL) in the study cohort, penalized approaches (HyperLasso models) identified an intergenic SNP (rs3094626 between RPP21 and HLA-E) and an intronic SNP (rs3134931 in NOTCH4) as novel correlates of early set-point VL in SCs. The minor allele of rs2857114 (downstream from HLA-DOB) was an unfavorable factor in SPs. Joint models based on demographic features, HLA alleles and the newly identified SNP variants could explain 29% and 15% of VL variance in SCs and SPs, respectively. These findings and bioinformatics strongly suggest that both classic and non-classic MHC genes deserve further investigation, especially in Africans with relatively short haplotype blocks.

Evidence for an epistatic effect between TP53 R72P and MDM2 T309G SNPs in HIV infection: a cross-sectional study in women from South Brazil

OBJECTIVE

To investigate the associations of TP53 R72P and MDM2 T309G SNPs with HPV infection status, HPV oncogenic risk and HIV infection status.

DESIGN

Cross-sectional study combining two groups (150 HIV-negative and 100 HIV-positive) of women.

METHODS

Data was collected using a closed questionnaire. DNA was extracted from cervical samples. HPV infection status was determined by nested-PCR, and HPV oncogenic risk group by Sanger sequencing. Both SNPS were genotyped by PCR-RFLP. Crude and adjusted associations involving each exposure (R72P and T309G SNPs, as well as 13 models of epistasis) and each outcome (HPV status, HPV oncogenic risk group and HIV infection) were assessed using logistic regression.

RESULTS

R72P SNP was protectively associated with HPV status (overdominant model), as well as T309G SNP with HPV oncogenic risk (strongest in the overdominant model). No epistatic model was associated with HPV status, but a dominant (R72P over T309G) protective epistatic effect was observed for HPV oncogenic risk. HIV status was strongly associated (risk factor) with different epistatic models, especially in models based on a visual inspection of the results. Moreover, HIV status was evidenced to be an effect mediator of the associations involving HPV oncogenic risk.

CONCLUSIONS

We found evidence for a role of R72P and T309G SNPs in HPV status and HPV oncogenic risk (respectively), and strong associations were found for an epistatic effect in HIV status. Prospective studies in larger samples are warranted to validate our findings, which point to a novel role of these SNPs in HIV infection.

Genetic predictors of cervical dysplasia in African American HIV-infected women: ACTG DACS 268

OBJECTIVE

To examine genome-wide associations in HIV-infected women with a history of cervical dysplasia compared with HIV-infected women with no history of abnormal Papanicolaou (Pap) tests.

DESIGN

Case-control study using data from women analyzed for the HIV Controllers Study and enrolled in HIV treatment-naïve studies in the AIDS Clinical Trials Group (ACTG).

METHODS

Genotyping utilized Illumina HumanHap 650 Y or 1MDuo platforms. After quality control and principal component analysis, ~610,000 significant single nucleotide polymorphisms (SNPs) were tested for association. Threshold for significance was P < 5 × 10(-8) for genome-wide associations.

RESULTS

No significant genomic association was observed between women with low-grade dysplasia and controls. The genome-wide association study (GWAS) analysis between women with high-grade dysplasia or invasive cervical cancer and normal controls identified significant SNPs. In the analyses limited to African American women, 11 SNPs were significantly associated with the development of high-grade dysplasia or cancer after correcting for multiple comparisons. The model using significant SNPs alone had improved accuracy in predicting high-grade dysplasia in African American women compared to the use of clinical data (area under the receiver operating characteristic curve for genetic and clinical model = 0.9 and 0.747, respectively).

CONCLUSIONS

These preliminary data serve as proof of concept that there may be a genetic predisposition to developing high-grade cervical dysplasia in African American HIV-infected women. Given the small sample size, the results need to be validated in a separate cohort.

Genetic variation and HIV-associated neurologic disease

HIV-associated neurologic disease continues to be a significant complication in the era of highly active antiretroviral therapy. A substantial subset of the HIV-infected population shows impaired neuropsychological performance as a result of HIV-mediated neuroinflammation and eventual central nervous system (CNS) injury. CNS compartmentalization of HIV, coupled with the evolution of genetically isolated populations in the CNS, is responsible for poor prognosis in patients with AIDS, warranting further investigation and possible additions to the current therapeutic strategy. This chapter reviews key advances in the field of neuropathogenesis and studies that have highlighted how molecular diversity within the HIV genome may impact HIV-associated neurologic disease. We also discuss the possible functional implications of genetic variation within the viral promoter and possibly other regions of the viral genome, especially in the cells of monocyte-macrophage lineage, which are arguably key cellular players in HIV-associated CNS disease.

Insights from host genomics into HIV infection and disease: Identification of host targets for drug development

HIV susceptibility and disease progression show a substantial degree of individual heterogeneity, ranging from fast progressors to long-term non progressors or elite controllers, that is, subjects that control infection in the absence of therapy. Recent years have seen a significant increase in understanding of the host genetic determinants of susceptibility to HIV infection and disease progression, driven in large part by candidate gene studies, genome-wide association studies, genome-wide transcriptome analyses, and large-scale functional screens. These studies have identified common variants in host loci that clearly influence disease progression, characterized the scale and dynamics of gene and protein expression changes in response to infection, and provided the first comprehensive catalogue of genes and pathways involved in viral replication. This review highlights the potential of host genomic influences in antiviral therapy by pointing to promising novel drug targets but also providing the basis of the identification and validation of host mechanisms that might be susceptible targets for novel antiviral therapies.

Association study of common genetic variants and HIV-1 acquisition in 6,300 infected cases and 7,200 controls

Multiple genome-wide association studies (GWAS) have been performed in HIV-1 infected individuals, identifying common genetic influences on viral control and disease course. Similarly, common genetic correlates of acquisition of HIV-1 after exposure have been interrogated using GWAS, although in generally small samples. Under the auspices of the International Collaboration for the Genomics of HIV, we have combined the genome-wide single nucleotide polymorphism (SNP) data collected by 25 cohorts, studies, or institutions on HIV-1 infected individuals and compared them to carefully matched population-level data sets (a list of all collaborators appears in Note S1 in Text S1). After imputation using the 1,000 Genomes Project reference panel, we tested approximately 8 million common DNA variants (SNPs and indels) for association with HIV-1 acquisition in 6,334 infected patients and 7,247 population samples of European ancestry. Initial association testing identified the SNP rs4418214, the C allele of which is known to tag the HLA-B*57:01 and B*27:05 alleles, as genome-wide significant (p = 3.6×10⁻¹¹). However, restricting analysis to individuals with a known date of seroconversion suggested that this association was due to the frailty bias in studies of lethal diseases. Further analyses including testing recessive genetic models, testing for bulk effects of non-genome-wide significant variants, stratifying by sexual or parenteral transmission risk and testing previously reported associations showed no evidence for genetic influence on HIV-1 acquisition (with the exception of CCR5Δ32 homozygosity). Thus, these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Comparing the frequency differences between common DNA variants in disease-affected cases and in unaffected controls has been successful in uncovering the genetic component of multiple diseases. This approach is most effective when large samples of cases and controls are available. Here we combine information from multiple studies of HIV infected patients, including more than 6,300 HIV+ individuals, with data from 7,200 general population samples of European ancestry to test nearly 8 million common DNA variants for an impact on HIV acquisition. With this large sample we did not observe any single common genetic variant that significantly associated with HIV acquisition. We further tested 22 variants previously identified by smaller studies as influencing HIV acquisition. With the exception of a deletion polymorphism in the CCR5 gene (CCR5Δ32) we found no convincing evidence to support these previous associations. Taken together these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Unravelling the mechanisms of durable control of HIV-1

Untreated HIV-1 infection typically progresses to AIDS within 10 years, but less than 1% of infected individuals remain healthy and have normal CD4(+) T cell counts and undetectable viral loads; some individuals have remained this way for 35 years and counting. Through a combination of large population studies of cohorts of these 'HIV-1 controllers' and detailed studies of individual patients, a heterogeneous picture has emerged regarding the basis for this remarkable resistance to AIDS progression. In this Review, we highlight the host genetic factors, the viral genetic factors and the immunological factors that are associated with the controller phenotype, we discuss emerging methodological approaches that could facilitate a better understanding of spontaneous HIV-1 immune control in the future, and we delineate implications for a 'functional cure' of HIV-1 infection.

MicroRNA regulation and its effects on cellular transcriptome in human immunodeficiency virus-1 (HIV-1) infected individuals with distinct viral load and CD4 cell counts

Background

Disease progression in the absence of therapy varies significantly in HIV-1 infected individuals. Both viral and host cellular molecules are implicated; however, the exact role of these factors and/or the mechanism involved remains elusive. To understand how microRNAs (miRNAs), which are regulators of transcription and translation, influence host cellular gene expression (mRNA) during HIV-1 infection, we performed a comparative miRNA and mRNA microarray analysis using PBMCs obtained from infected individuals with distinct viral load and CD4 counts.

Methods

RNA isolated from PBMCs obtained from HIV-1 seronegative and HIV-1 positive individuals with distinct viral load and CD4 counts were assessed for miRNA and mRNA profile. Selected miRNA and mRNA transcripts were validated using in vivo and in vitro infection model.

Results

Our results indicate that HIV-1 positive individuals with high viral load (HVL) showed a dysregulation of 191 miRNAs and 309 mRNA transcripts compared to the uninfected age and sex matched controls. The miRNAs miR-19b, 146a, 615-3p, 382, 34a, 144 and 155, that are known to target innate and inflammatory factors, were significantly upregulated in PBMCs with high viral load, as were the inflammatory molecules CXCL5, CCL2, IL6 and IL8, whereas defensin, CD4, ALDH1, and Neurogranin (NRGN) were significantly downregulated. Using the transcriptome profile and predicted target genes, we constructed the regulatory networks of miRNA-mRNA pairs that were differentially expressed between control, LVL and HVL subjects. The regulatory network revealed an inverse correlation of several miRNA-mRNA pair expression patterns, suggesting HIV-1 mediated transcriptional regulation is in part likely through miRNA regulation.

Conclusions

Results from our studies indicate that gene expression is significantly altered in PBMCs in response to virus replication. It is interesting to note that the infected individuals with low or undetectable viral load exhibit a gene expression profile very similar to control or uninfected subjects. Importantly, we identified several new mRNA targets (Defensin, Neurogranin, AIF) as well as the miRNAs that could be involved in regulating their expression through the miRNA-mRNA interaction.

Immunogenetics: Genome-Wide Association of Non-Progressive HIV and Viral Load Control: HLA Genes and Beyond

Very early after the identification of the human immunodeficiency virus (HIV), host genetics factors were anticipated to play a role in viral control and disease progression. As early as the mid-1990s, candidate gene studies demonstrated a central role for the chemokine co-receptor/ligand (e.g., CCR5) and human leukocyte antigen (HLA) systems. In the last decade, the advent of genome-wide arrays opened a new era for unbiased genetic exploration of the genome and brought big expectations for the identification of new unexpected genes and pathways involved in HIV/AIDS. More than 15 genome-wide association studies targeting various HIV-linked phenotypes have been published since 2007. Surprisingly, only the two HIV-chemokine co-receptors and HLA loci have exhibited consistent and reproducible statistically significant genetic associations. In this chapter, we will review the findings from the genome-wide studies focusing especially on non-progressive and HIV control phenotypes, and discuss the current perspectives.

Host genomic influences on HIV/AIDS

The AIDS era has seen multiple advances in the power of genetics research; scores of host genetic protective factors have been nominated and several have translated to the bedside. We discuss how genomics may inform HIV/AIDS prevention, treatment and eradication.

Immunogenetics of human echinococcosis

Susceptibility and resistance to human Echinococcus infection and disease, although poorly understood, appear to reflect a complex interaction of parasite and host immunological and genetic factors. Disease stage, progression, and prognosis following treatment appear to be strongly influenced by cytokine and antibody profiles, and more recent evidence has suggested an important role of dendritic cells (DCs) and T regulatory cells (Tregs) in immunomodulation. Microarrays have supported these findings, highlighting both known and novel pathways involved in chronic murine disease. Genetic studies to date have been few and with limited success. Advanced genomic approaches, such as genome-wide association studies (GWAS), may provide further insight to identify the relevant pathways involved, thereby facilitating a new approach for the development of new clinical therapies.

Genome-wide association studies on HIV susceptibility, pathogenesis and pharmacogenomics

Susceptibility to HIV-1 and the clinical course after infection show a substantial heterogeneity between individuals. Part of this variability can be attributed to host genetic variation. Initial candidate gene studies have revealed interesting host factors that influence HIV infection, replication and pathogenesis. Recently, genome-wide association studies (GWAS) were utilized for unbiased searches at a genome-wide level to discover novel genetic factors and pathways involved in HIV-1 infection. This review gives an overview of findings from the GWAS performed on HIV infection, within different cohorts, with variable patient and phenotype selection. Furthermore, novel techniques and strategies in research that might contribute to the complete understanding of virus-host interactions and its role on the pathogenesis of HIV infection are discussed.

Host genes important to HIV replication and evolution

Recent years have seen a significant increase in understanding of the host genetic and genomic determinants of susceptibility to HIV-1 infection and disease progression, driven in large part by candidate gene studies, genome-wide association studies, genome-wide transcriptome analyses, and large-scale in vitro genome screens. These studies have identified common variants in some host loci that clearly influence disease progression, characterized the scale and dynamics of gene and protein expression changes in response to infection, and provided the first comprehensive catalogs of genes and pathways involved in viral replication. Experimental models of AIDS and studies in natural hosts of primate lentiviruses have complemented and in some cases extended these findings. As the relevant technology continues to progress, the expectation is that such studies will increase in depth (e.g., to include host whole exome and whole genome sequencing) and in breadth (in particular, by integrating multiple data types).

Multicohort genomewide association study reveals a new signal of protection against HIV-1 acquisition

BACKGROUND

To date, only mutations in CCR5 have been shown to confer resistance to human immunodeficiency virus type 1 (HIV-1) infection, and these explain only a small fraction of the observed variability in HIV susceptibility.

METHODS

We performed a meta-analysis between 2 independent European genomewide association studies, each comparing HIV-1 seropositive cases with normal population controls known to be HIV uninfected, to identify single-nucleotide polymorphisms (SNPs) associated with the HIV-1 acquisition phenotype. SNPs exhibiting P < 10(-5) in this first stage underwent second-stage analysis in 2 independent US cohorts of European descent.

RESULTS

After the first stage, a single highly significant association was revealed for the chromosome 8 rs6996198 with HIV-1 acquisition and was replicated in both second-stage cohorts. Across the 4 groups, the rs6996198-T allele was consistently associated with a significant reduced risk of HIV-1 infection, and the global meta-analysis reached genomewide significance: P(combined) = 7.76 × 10(-8).

CONCLUSIONS

We provide strong evidence of association for a common variant with HIV-1 acquisition in populations of European ancestry. This protective signal against HIV-1 infection is the first identified outside the CCR5 nexus. First clues point to a potential functional role for a nearby candidate gene, CYP7B1, but this locus warrants further investigation.

Incorporation of innate immune effector mechanisms in the formulation of a vaccine against HIV-1

The realization of a major role for events that occur during acute viremia that dictate the course of disease both in HIV-1 infected humans and susceptible SIV infected non-human primates has prompted an intense interest in studies of the contribution of innate immune effector mechanisms. It is reasoned that findings from such studies may be important and need to be incorporated into the design and formulation of potential candidate vaccines against HIV-1. This review serves to outline the various non-human primate models that can best serve to address this issue, a summary of our knowledge on the various subsets of NK cells (one of the major innate immune cell lineage) that have an impact on the course of disease, the potential pathways that regulate their function and the potential role of the KIRs on SIV-induced disease course. Finally, the major points from this report and the data presented on similar subjects by other investigators is utilized to provide a summary of the potential future directions that we need to take in efforts to move this field forward.

Genome-wide association studies and susceptibility to infectious diseases

Progress in genomics and the associated technological, statistical and bioinformatics advances have facilitated the successful implementation of genome-wide association studies (GWAS) towards understanding the genetic basis of common diseases. Infectious diseases contribute significantly to the global burden of disease and there is robust epidemiological evidence that host genetic factors are important determinants of the outcome of interactions between host and pathogen. Indeed, infectious diseases have exerted profound selective pressure on human evolution. However, the application of GWAS to infectious diseases has been relatively limited compared with non-communicable diseases. Here we review GWAS findings for important infectious diseases, including malaria, tuberculosis and HIV. We highlight some of the pitfalls recognized more generally for GWAS, as well as issues specific to infection, including the role of the pathogen which also has a genome. We also discuss the challenges encountered when studying African populations which are genetically more ancient and more diverse that other populations and disproportionately bear the main global burden of serious infectious diseases.

Role of exonic variation in chemokine receptor genes on AIDS: CCRL2 F167Y association with pneumocystis pneumonia

Chromosome 3p21-22 harbors two clusters of chemokine receptor genes, several of which serve as major or minor coreceptors of HIV-1. Although the genetic association of CCR5 and CCR2 variants with HIV-1 pathogenesis is well known, the role of variation in other nearby chemokine receptor genes remain unresolved. We genotyped exonic single nucleotide polymorphisms (SNPs) in chemokine receptor genes: CCR3, CCRL2, and CXCR6 (at 3p21) and CCR8 and CX3CR1 (at 3p22), the majority of which were non-synonymous. The individual SNPs were tested for their effects on disease progression and outcomes in five treatment-naïve HIV-1/AIDS natural history cohorts. In addition to the known CCR5 and CCR2 associations, significant associations were identified for CCR3, CCR8, and CCRL2 on progression to AIDS. A multivariate survival analysis pointed to a previously undetected association of a non-conservative amino acid change F167Y in CCRL2 with AIDS progression: 167F is associated with accelerated progression to AIDS (RH = 1.90, P = 0.002, corrected). Further analysis indicated that CCRL2-167F was specifically associated with more rapid development of pneumocystis pneumonia (PCP) (RH = 2.84, 95% CI 1.28-6.31) among four major AIDS-defining conditions. Considering the newly defined role of CCRL2 in lung dendritic cell trafficking, this atypical chemokine receptor may affect PCP through immune regulation and inducing inflammation.

Polymorphism in HIV-1 dependency factor PDE8A affects mRNA level and HIV-1 replication in primary macrophages

Four genome-wide RNAi screens have recently identified hundreds of HIV-1 dependency factors (HDFs). Previously, we reported a large variation in the ability of HIV-1 to replicate in monocyte-derived macrophages (MDM) derived from >400 healthy seronegative blood donors. Here we determined whether SNPs in genes encoding newly identified HDFs were associated with this variation in HIV-1 replication. We found a significant association between the minor allele of SNP rs2304418 in phosphodiesterase 8A (PDE8A) and lower HIV-1 replication (p=2.4×10(-6)). The minor allele of SNP rs2304418 was also significantly associated with lower PDE8A mRNA levels in MDM (p=8.3×10(-5)). In accordance with this, overexpression of PDE8A in HEK293T cells resulted in increased HIV-1 replication, while subsequent knock-down of PDE8A decreased replication. This study links host genetic variation in a newly identified HDF to variation in HIV-1 replication in a relevant primary target cell for HIV-1 and may provide new leads for treatment of this infection.

Comparative expression profile of miRNA and mRNA in primary peripheral blood mononuclear cells infected with human immunodeficiency virus (HIV-1)

Host cells respond to exogenous infectious agents such as viruses, including HIV-1. Studies have evaluated the changes associated with virus infection at the transcriptional and translational levels of the cellular genes involved in specific pathways. While this approach is useful, in our view it provides only a partial view of genome-wide changes. Recently, technological advances in the expression profiling at the microRNA (miRNA) and mRNA levels have made it possible to evaluate the changes in the components of multiple pathways. To understand the role of miRNA and its interplay with host cellular gene expression (mRNA) during HIV-1 infection, we performed a comparative global miRNA and mRNA microarray using human PBMCs infected with HIV-1. The PBMCs were derived from multiple donors and were infected with virus generated from the molecular clone pNL4-3. The results showed that HIV-1 infection led to altered regulation of 21 miRNAs and 444 mRNA more than 2-fold, with a statistical significance of p<0.05. Furthermore, the differentially regulated miRNA and mRNA were shown to be associated with host cellular pathways involved in cell cycle/proliferation, apoptosis, T-cell signaling, and immune activation. We also observed a number of inverse correlations of miRNA and mRNA expression in infected PBMCs, further confirming the interrelationship between miRNA and mRNA regulation during HIV-1 infection. These results for the first time provide evidence that the miRNA profile could be an early indicator of host cellular dysfunction induced by HIV-1.

Genome-wide association scan in HIV-1-infected individuals identifying variants influencing disease course

Background

AIDS develops typically after 7–11 years of untreated HIV-1 infection, with extremes of very rapid disease progression (<2 years) and long-term non-progression (>15 years). To reveal additional host genetic factors that may impact on the clinical course of HIV-1 infection, we designed a genome-wide association study (GWAS) in 404 participants of the Amsterdam Cohort Studies on HIV-1 infection and AIDS.

Methods

The association of SNP genotypes with the clinical course of HIV-1 infection was tested in Cox regression survival analyses using AIDS-diagnosis and AIDS-related death as endpoints.

Results

Multiple, not previously identified SNPs, were identified to be strongly associated with disease progression after HIV-1 infection, albeit not genome-wide significant. However, three independent SNPs in the top ten associations between SNP genotypes and time between seroconversion and AIDS-diagnosis, and one from the top ten associations between SNP genotypes and time between seroconversion and AIDS-related death, had P-values smaller than 0.05 in the French Genomics of Resistance to Immunodeficiency Virus cohort on disease progression.

Conclusions

Our study emphasizes that the use of different phenotypes in GWAS may be useful to unravel the full spectrum of host genetic factors that may be associated with the clinical course of HIV-1 infection.

CD39/adenosine pathway is involved in AIDS progression

HIV-1 infection is characterized by a chronic activation of the immune system and suppressed function of T lymphocytes. Regulatory CD4+ CD25(high) FoxP3+CD127(low) T cells (Treg) play a key role in both conditions. Here, we show that HIV-1 positive patients have a significant increase of Treg-associated expression of CD39/ENTPD1, an ectoenzyme which in concert with CD73 generates adenosine. We show in vitro that the CD39/adenosine axis is involved in Treg suppression in HIV infection. Treg inhibitory effects are relieved by CD39 down modulation and are reproduced by an adenosine-agonist in accordance with a higher expression of the adenosine A2A receptor on patients' T cells. Notably, the expansion of the Treg CD39+ correlates with the level of immune activation and lower CD4+ counts in HIV-1 infected patients. Finally, in a genetic association study performed in three different cohorts, we identified a CD39 gene polymorphism that was associated with down-modulated CD39 expression and a slower progression to AIDS.

Replication-competent simian immunodeficiency virus (SIV) Gag escape mutations archived in latent reservoirs during antiretroviral treatment of SIV-infected macaques

In response to pressure exerted by major histocompatibility complex (MHC) class I-mediated CD8(+) T cell control, human immunodeficiency virus (HIV) escape mutations often arise in immunodominant epitopes recognized by MHC class I alleles. While the current standard of care for HIV-infected patients is treatment with highly active antiretroviral therapy (HAART), suppression of viral replication in these patients is not absolute and latently infected cells persist as lifelong reservoirs. To determine whether HIV escape from MHC class I-restricted CD8(+) T cell control develops during HAART treatment and then enters latent reservoirs in the periphery and central nervous system (CNS), with the potential to emerge as replication-competent virus, we tracked the longitudinal development of the simian immunodeficiency virus (SIV) Gag escape mutation K165R in HAART-treated SIV-infected pigtailed macaques. Key findings of these studies included: (i) SIV Gag K165R escape mutations emerged in both plasma and cerebrospinal fluid (CSF) during the decaying phase of viremia after HAART initiation before suppression of viral replication, (ii) SIV K165R Gag escape mutations were archived in latent proviral DNA reservoirs, including the brain in animals receiving HAART that suppressed viral replication, and (iii) replication-competent SIV Gag K165R escape mutations were present in the resting CD4(+) T cell reservoir in HAART-treated SIV-infected macaques. Despite early administration of aggressive antiretroviral treatment, HIV immune escape from CD8(+) T cell control can still develop during the decaying phases of viremia and then persist in latent reservoirs, including the brain, with the potential to emerge if HAART therapy is interrupted.

Genome-wide association study implicates PARD3B-based AIDS restriction

BACKGROUND

Host genetic variation influences human immunodeficiency virus (HIV) infection and progression to AIDS. Here we used clinically well-characterized subjects from 5 pretreatment HIV/AIDS cohorts for a genome-wide association study to identify gene associations with rate of AIDS progression.

METHODS

European American HIV seroconverters (n = 755) were interrogated for single-nucleotide polymorphisms (SNPs) (n = 700,022) associated with progression to AIDS 1987 (Cox proportional hazards regression analysis, co-dominant model).

RESULTS

Association with slower progression was observed for SNPs in the gene PARD3B. One of these, rs11884476, reached genome-wide significance (relative hazard = 0.3; P =3. 370 × 10(-9)) after statistical correction for 700,022 SNPs and contributes 4.52% of the overall variance in AIDS progression in this study. Nine of the top-ranked SNPs define a PARD3B haplotype that also displays significant association with progression to AIDS (hazard ratio, 0.3; P = 3.220 × 10(-8)). One of these SNPs, rs10185378, is a predicted exonic splicing enhancer; significant alteration in the expression profile of PARD3B splicing transcripts was observed in B cell lines with alternate rs10185378 genotypes. This SNP was typed in European cohorts of rapid progressors and was found to be protective for AIDS 1993 definition (odds ratio, 0.43, P = .025).

CONCLUSIONS

These observations suggest a potential unsuspected pathway of host genetic influence on the dynamics of AIDS progression.

Progression to AIDS in South Africa is associated with both reverting and compensatory viral mutations

We lack the understanding of why HIV-infected individuals in South Africa progress to AIDS. We hypothesised that in end-stage disease there is a shifting dynamic between T cell imposed immunity and viral immune escape, which, through both compensatory and reverting viral mutations, results in increased viral fitness, elevated plasma viral loads and disease progression. We explored how T cell responses, viral adaptation and viral fitness inter-relate in South African cohorts recruited from Bloemfontein, the Free State (n = 278) and Durban, KwaZulu-Natal (n = 775). Immune responses were measured by γ-interferon ELISPOT assays. HLA-associated viral polymorphisms were determined using phylogenetically corrected techniques, and viral replication capacity (VRC) was measured by comparing the growth rate of gag-protease recombinant viruses against recombinant NL4-3 viruses. We report that in advanced disease (CD4 counts <100 cells/µl), T cell responses narrow, with a relative decline in Gag-directed responses (p<0.0001). This is associated with preserved selection pressure at specific viral amino acids (e.g., the T242N polymorphism within the HLA-B*57/5801 restricted TW10 epitope), but with reversion at other sites (e.g., the T186S polymorphism within the HLA-B*8101 restricted TL9 epitope), most notably in Gag and suggestive of "immune relaxation". The median VRC from patients with CD4 counts <100 cells/µl was higher than from patients with CD4 counts ≥ 500 cells/µl (91.15% versus 85.19%, p = 0.0004), potentially explaining the rise in viral load associated with disease progression. Mutations at HIV Gag T186S and T242N reduced VRC, however, in advanced disease only the T242N mutants demonstrated increasing VRC, and were associated with compensatory mutations (p = 0.013). These data provide novel insights into the mechanisms of HIV disease progression in South Africa. Restoration of fitness correlates with loss of viral control in late disease, with evidence for both preserved and relaxed selection pressure across the HIV genome. Interventions that maintain viral fitness costs could potentially slow progression.

Strenuous resistance to natural HIV-1 disease progression: viral controllers and long-term nonprogressors

HIV-1 infection leads to AIDS and death within 8–10 years for most individuals in the absence of antiretroviral therapy (ART). However, a minority of infected individuals show the unusual capacity to spontaneously control disease progression after infection in the absence of any ART. So-called ‘long-term nonprogressors’ are defined by maintenance of peripheral CD4+ T-cell counts >500 cells/µl and good health without ART for >7 years since infection. More recently, ART-naive individuals who spontaneously control their viremia levels at either <50 or <2000 copies of RNA/ml for at least 12 months in the absence of ART have been named ‘elite controllers’ and ‘HIV controllers’, respectively. The overlap between long-term nonprogressors and elite controllers/HIV controllers is partial, and both groups collectively account for <5% of all infected individuals. Unraveling the nature of their relative resistance to HIV-1 disease progression would be of great value for HIV-prevention strategies.

A strong case for viral genetic factors in HIV virulence

HIV infections show great variation in the rate of progression to disease, and the role of viral genetic factors in this variation had remained poorly characterized until recently. Now a series of four studies [1-4] published within a year has filled this important gap and has demonstrated a robust effect of the viral genotype on HIV virulence.

Screening low-frequency SNPS from genome-wide association study reveals a new risk allele for progression to AIDS

BACKGROUND

Seven genome-wide association studies (GWAS) have been published in AIDS, and only associations in the HLA region on chromosome 6 and CXCR6 have passed genome-wide significance.

METHODS

We reanalyzed the data from 3 previously published GWAS, targeting specifically low-frequency SNPs (minor allele frequency <5%). Two groups composed of 365 slow progressors and 147 rapid progressors from Europe and the United States were compared with a control group of 1394 seronegative individuals using Eigenstrat corrections.

RESULTS

Of the 8584 SNPs with minor allele frequency <5% in cases and controls (Bonferroni threshold = 5.8 × 10⁻⁶), 4 SNPs showed statistical evidence of association with the slow progressor phenotype. The best result was for HCP5 rs2395029 [P = 8.54 × 10⁻¹⁵, odds ratio (OR) = 3.41] in the HLA locus, in partial linkage disequilibrium with 2 additional chromosome 6 associations in C6orf48 (P = 3.03 × 10⁻¹⁰, OR = 2.9) and NOTCH4 (9.08 × 10⁻⁰⁷, OR = 2.32). The fourth association corresponded to rs2072255 located in RICH2 (P = 3.30 × 10⁻⁰⁶, OR = 0.43) in chromosome 17. Using HCP5 rs2395029 as a covariate, the C6orf48 and NOTCH4 signals disappeared, but the RICH2 signal still remained significant.

CONCLUSIONS

Besides the already known chromosome 6 associations, the analysis of low-frequency SNPs brought up a new association in the RICH2 gene. Interestingly, RICH2 interacts with BST-2 known to be a major restriction factor for HIV-1 infection. Our study has thus identified a new candidate gene for AIDS molecular etiology and confirms the interest of singling out low-frequency SNPs to exploit GWAS data.

The search for host genetic factors of HIV/AIDS pathogenesis in the post-genome era: progress to date and new avenues for discovery

Though pursuit of host genetic factors that influence the pathogenesis of HIV began over two decades ago, progress has been slow. Initial genome-level searches for variations associated with HIV-related traits have yielded interesting candidates, but less in the way of novel pathways to be exploited for therapeutic targets. More recent genome-wide association studies (GWAS) that include different phenotypes, novel designs, and that have examined different population characteristics suggest novel targets and affirm the utility of additional searches. Recent findings from these GWAS are reviewed, new directions for research are identified, and the promise of systems biology to yield novel insights is discussed.

Genome-wide association study identifies single nucleotide polymorphism in DYRK1A associated with replication of HIV-1 in monocyte-derived macrophages

BACKGROUND

HIV-1 infected macrophages play an important role in rendering resting T cells permissive for infection, in spreading HIV-1 to T cells, and in the pathogenesis of AIDS dementia. During highly active anti-retroviral treatment (HAART), macrophages keep producing virus because tissue penetration of antiretrovirals is suboptimal and the efficacy of some is reduced. Thus, to cure HIV-1 infection with antiretrovirals we will also need to efficiently inhibit viral replication in macrophages. The majority of the current drugs block the action of viral enzymes, whereas there is an abundance of yet unidentified host factors that could be targeted. We here present results from a genome-wide association study identifying novel genetic polymorphisms that affect in vitro HIV-1 replication in macrophages.

METHODOLOGY/PRINCIPAL FINDINGS

Monocyte-derived macrophages from 393 blood donors were infected with HIV-1 and viral replication was determined using Gag p24 antigen levels. Genomic DNA from individuals with macrophages that had relatively low (n = 96) or high (n = 96) p24 production was used for SNP genotyping with the Illumina 610 Quad beadchip. A total of 494,656 SNPs that passed quality control were tested for association with HIV-1 replication in macrophages, using linear regression. We found a strong association between in vitro HIV-1 replication in monocyte-derived macrophages and SNP rs12483205 in DYRK1A (p = 2.16 × 10(-5)). While the association was not genome-wide significant (p<1 × 10(-7)), we could replicate this association using monocyte-derived macrophages from an independent group of 31 individuals (p = 0.0034). Combined analysis of the initial and replication cohort increased the strength of the association (p = 4.84 × 10(-6)). In addition, we found this SNP to be associated with HIV-1 disease progression in vivo in two independent cohort studies (p = 0.035 and p = 0.0048).

CONCLUSIONS/SIGNIFICANCE

These findings suggest that the kinase DYRK1A is involved in the replication of HIV-1, in vitro in macrophages as well as in vivo.

Leprosy and the human genome

Despite the availability of effective treatment for several decades, leprosy remains an important medical problem in many regions of the world. Infection with Mycobacterium leprae can produce paucibacillary disease, characterized by well-formed granulomas and a Th1 T-cell response, or multibacillary disease, characterized by poorly organized cellular infiltrates and Th2 cytokines. These diametric immune responses confer states of relative resistance or susceptibility to leprosy, respectively, and have well-defined clinical manifestations. As a result, leprosy provides a unique opportunity to dissect the genetic basis of human in vivo immunity. A series of studies over the past 40 years suggests that host genes influence the risk of leprosy acquisition and the predilection for different clinical forms of the disease. However, a comprehensive, cellular, and molecular view of the genes and variants involved is still being assembled. In this article, we review several decades of human genetic studies of leprosy, including a number of recent investigations. We emphasize genetic analyses that are validated by the replication of the same phenotype in independent studies or supported by functional experiments demonstrating biological mechanisms of action for specific polymorphisms. Identifying and functionally exploring the genetic and immunological factors that underlie human susceptibility to leprosy have yielded important insights into M. leprae pathogenesis and are likely to advance our understanding of the immune response to other pathogenic mycobacteria. This knowledge may inform new treatment or vaccine strategies for leprosy or tuberculosis.

Sample size requirements to detect gene-environment interactions in genome-wide association studies

Many complex diseases are likely to be a result of the interplay of genes and environmental exposures. The standard analysis in a genome-wide association study (GWAS) scans for main effects and ignores the potentially useful information in the available exposure data. Two recently proposed methods that exploit environmental exposure information involve a two-step analysis aimed at prioritizing the large number of SNPs tested to highlight those most likely to be involved in a GE interaction. For example, Murcray et al. ([2009] Am J Epidemiol 169:219–226) proposed screening on a test that models the G-E association induced by an interaction in the combined case-control sample. Alternatively, Kooperberg and LeBlanc ([2008] Genet Epidemiol 32:255–263) suggested screening on genetic marginal effects. In both methods, SNPs that pass the respective screening step at a pre-specified significance threshold are followed up with a formal test of interaction in the second step. We propose a hybrid method that combines these two screening approaches by allocating a proportion of the overall genomewide significance level to each test. We show that the Murcray et al. approach is often the most efficient method, but that the hybrid approach is a powerful and robust method for nearly any underlying model. As an example, for a GWAS of 1 million markers including a single true disease SNP with minor allele frequency of 0.15, and a binary exposure with prevalence 0.3, the Murcray, Kooperberg and hybrid methods are 1.90, 1.27, and 1.87 times as efficient, respectively, as the traditional case-control analysis to detect an interaction effect size of 2.0.