Cédric Coulonges

Genomewide association study of an AIDS-nonprogression cohort emphasizes the role played by HLA genes (ANRS genomewide association study 02)

To elucidate the genetic factors predisposing to AIDS progression, we analyzed a unique cohort of 275 human immunodeficiency virus (HIV) type 1-seropositive nonprogressor patients in relation to a control group of 1352 seronegative individuals in a genomewide association study (GWAS). The strongest association was obtained for HCP5 rs2395029 (P=6.79x10(-10); odds ratio, 3.47) and was possibly linked to an effect of sex. Interestingly, this single-nucleotide polymorphism (SNP) was in high linkage disequilibrium with HLA-B, MICB, TNF, and several other HLA locus SNPs and haplotypes. A meta-analysis of our genomic data combined with data from the previously conducted Euro-CHAVI (Center for HIV/AIDS Vaccine Immunology) GWAS confirmed the HCP5 signal (P=3.02x10(-19)) and identified several new associations, all of them involving HLA genes: MICB, TNF, RDBP, BAT1-5, PSORS1C1, and HLA-C. Finally, stratification by HCP5 rs2395029 genotypes emphasized an independent role for ZNRD1, also in the HLA locus, and this finding was confirmed by experimental data. The present study, the first GWAS of HIV-1 nonprogressors, underscores the potential for some HLA genes to control disease progression soon after infection.

Genomewide Association Study of a Rapid Progression Cohort Identifies New Susceptibility Alleles for AIDS (ANRS Genomewide Association Study 03)

BACKGROUND Previous genomewide association studies (GWASs) of AIDS have targeted end points based on the control of viral load and disease nonprogression. The discovery of genetic factors that predispose individuals to rapid progression to AIDS should also reveal new insights into the molecular etiology of the pathology. METHODS We undertook a case-control GWAS of a unique cohort of 85 human immunodeficiency virus type 1 (HIV-1)-infected patients who experienced rapid disease progression, using Illumina HumanHap300 BeadChips. The case group was compared with a control group of 1352 individuals for the 291,119 autosomal single-nucleotide polymorphisms (SNPs) passing the quality control tests, using the false-discovery rate (FDR) statistical method for multitest correction. RESULTS Novel associations with rapid progression (FDR, < or = 25%) were identified for PRMT6 (P = 6.1 x 10(-7); odds ratio [OR], 0.24), SOX5 (P = 1.8 x 10(-6); OR, 0.45), RXRG (P = 3.9 x 10(-6); OR, 3.29), and TGFBRAP1 (P = 7 x 10(-6); OR, 0.34). The haplotype analysis identified exonic and promoter SNPs potentially important for PRMT6 and TGFBRAP1 function. CONCLUSIONS The statistical and biological relevance of these associations and their high ORs underscore the power of extreme phenotypes for GWASs, even with a modest sample size. These genetic results emphasize the role of the transforming growth factor beta pathway in the pathogenesis of HIV-1 disease. Finally, the wealth of information provided by this study should help unravel new diagnostic and therapeutic targets.

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−11). 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.

Exhaustive Genotyping of the CEM15 (APOBEC3G) Gene and Absence of Association with AIDS Progression in a French Cohort

CEM15 (or APOBEC3G) has recently been identified as an inhibitor of human immunodeficiency virus type 1 (HIV-1) replication in vitro. To evaluate the impact of its genetic variations on the progression of acquired immunodeficiency syndrome (AIDS), we have performed an extensive genetic analysis of CEM15. We have sequenced CEM15 in a cohort of 327 HIV-1-seropositive patients with extreme disease progression phenotypes--either slow progression or rapid progression--and in 446 healthy control subjects, all of white descent. We have identified 29 polymorphisms with allele frequencies >1%, 14 of which were newly characterized. There were no significant associations between the polymorphisms or haplotypes of CEM15 and a disease progression phenotype in our cohort.

Shape-IT: new rapid and accurate algorithm for haplotype inference

BackgroundWe have developed a new computational algorithm, Shape-IT, to infer haplotypes under the genetic model of coalescence with recombination developed by Stephens et al in Phase v2.1. It runs much faster than Phase v2.1 while exhibiting the same accuracy. The major algorithmic improvements rely on the use of binary trees to represent the sets of candidate haplotypes for each individual. These binary tree representations: (1) speed up the computations of posterior probabilities of the haplotypes by avoiding the redundant operations made in Phase v2.1, and (2) overcome the exponential aspect of the haplotypes inference problem by the smart exploration of the most plausible pathways (ie. haplotypes) in the binary trees.ResultsOur results show that Shape-IT is several orders of magnitude faster than Phase v2.1 while being as accurate. For instance, Shape-IT runs 50 times faster than Phase v2.1 to compute the haplotypes of 200 subjects on 6,000 segments of 50 SNPs extracted from a standard Illumina 300 K chip (13 days instead of 630 days). We also compared Shape-IT with other widely used software, Gerbil, PL-EM, Fastphase, 2SNP, and Ishape in various tests: Shape-IT and Phase v2.1 were the most accurate in all cases, followed by Ishape and Fastphase. As a matter of speed, Shape-IT was faster than Ishape and Fastphase for datasets smaller than 100 SNPs, but Fastphase became faster -but still less accurate- to infer haplotypes on larger SNP datasets.ConclusionShape-IT deserves to be extensively used for regular haplotype inference but also in the context of the new high-throughput genotyping chips since it permits to fit the genetic model of Phase v2.1 on large datasets. This new algorithm based on tree representations could be used in other HMM-based haplotype inference software and may apply more largely to other fields using HMM.

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

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

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.

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−6), 4 SNPs showed statistical evidence of association with the slow progressor phenotype. The best result was for HCP5 rs2395029 [P = 8.54 × 10−15, 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−10, OR = 2.9) and NOTCH4 (9.08 × 10−07, OR = 2.32). The fourth association corresponded to rs2072255 located in RICH2 (P = 3.30 × 10−06, 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.

A Genome-Wide Association Study in Caucasian Women Points Out a Putative Role of the STXBP5L Gene in Facial Photoaging

A genome-wide association study (GWAS) was conducted on 502 French middle-aged Caucasian women to identify genetic factors that may affect skin aging severity. A high-throughput Illumina Human Omni1-Quad beadchip was used. After single-nucleotide polymorphism (SNP) quality controls, 795,063 SNPs remained for analysis purposes. Possible stratification was first examined using the Eigenstrat method, and then the relationships between genotypes and four skin aging indicators (global photoaging, lentigines, wrinkles, and sagging) were investigated separately by linear regressions adjusted on age, smoking habits, lifetime sun exposure, hormonal status, and the two main Eigen vectors. One signal passed the Bonferroni threshold (P=1.53 × 10(-8)) and was significantly associated with global photoaging. It was also correlated with the wrinkling score and the sagging score. According to HapMap, this SNP, rs322458, was in linkage disequilibrium (LD) with intronic SNPs of the STXBP5L gene, which is expressed in the skin. In addition, it was also in LD with another SNP that increases the expression of the FBXO40 gene in the skin. These two genes, which were not previously described in the context of aging, may constitute good candidates for the investigation of molecular mechanisms of skin photoaging.

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.