The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control.

Introduction: We have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers. Methods: We performed a genome-wide Association Study comparing 543 HIV Elite Controllers with 3,272 uninfected controls of European descent. Using the latest database for imputation, we analyzed 35,552 Single Nucleotide Polymorphisms (SNPs) within the Major Histocompatibility Complex (MHC) region. Results: Our analysis identified 2,626 SNPs significantly associated (p<5. 10-8) with elite control of HIV-1 infection, including well-established MHC signals such as the rs2395029-G allele which tags HLA-B*57:01. A thorough investigation of SNPs in linkage disequilibrium with rs2395029 revealed an extensive haploblock spanning 1.9 megabases in the MHC region tagging HLA-B*57:01, comprising 379 SNP alleles impacting 72 genes. This haploblock contains damaging variations in proteins like NOTCH4 and DXO and is also associated with a strong differential pattern of expression of multiple MHC genes such as HLA-B, MICB, and ZBTB12. The study was expanded to include two cohorts of seropositive African-American individuals, where a haploblock tagging the HLA-B*57:03 allele was similarly associated with control of viral load. The mRNA expression profile of this haploblock in African Americans closely mirrored that in the European cohort. Discussion: These findings suggest that additional molecular mechanisms beyond the conventional antigen-presenting role of class I HLA molecules may contribute to the observed influence of HLA-B*57:01/B*57:03 alleles on HIV-1 elite control. Overall, this study has uncovered a large haploblock associated with HLA-B*57 alleles, providing novel insights into their massive effect on HIV-1 elite control.

Identification of New Biological Pathways Involved in Skin Aging From the Analysis of French Women Genome-Wide Data.

Skin aging is an ineluctable process leading to the progressive loss of tissue integrity and is characterized by various outcomes such as wrinkling and sagging. Researchers have identified impacting environmental factors (sun exposure, smoking, etc.) and several molecular mechanisms leading to skin aging. We have previously performed genome-wide association studies (GWAS) in 502 very-well characterized French women, looking for associations with four major outcomes of skin aging, namely, photoaging, solar lentigines, wrinkling, and sagging, and this has led to new insights into the molecular mechanisms of skin aging. Since individual SNP associations in GWAS explain only a small fraction of the genetic impact in complex polygenic phenotypes, we have made the integration of these genotypes into the reference Kegg biological pathways and looked for associations by the gene set enrichment analysis (GSEA) approach. 106 pathways were tested for association with the four outcomes of skin aging. This biological pathway analysis revealed new relevant pathways and genes, some likely specific of skin aging such as the WNT7B and PRKCA genes in the "melanogenesis" pathway and some likely involved in global aging such as the DDB1 gene in the "nucleotide excision repair" pathway, not picked up in the previously published GWAS. Overall, our results suggest that the four outcomes of skin aging possess specific molecular mechanisms such as the "proteasome" and "mTOR signaling pathway" but may also share common molecular mechanisms such as "nucleotide excision repair."

Comparative evaluation of shape retrieval methods on macromolecular surfaces: an application of computer vision methods in structural bioinformatics.

MOTIVATION: The investigation of the structure of biological systems at the molecular level gives insights about their functions and dynamics. Shape and surface of biomolecules are fundamental to molecular recognition events. Characterizing their geometry can lead to more adequate predictions of their interactions. In the present work, we assess the performance of reference shape retrieval methods from the computer vision community on protein shapes. RESULTS: Shape retrieval methods are efficient in identifying orthologous proteins and tracking large conformational changes. This work illustrates the interest for the protein surface shape as a higher-level representation of the protein structure that (i) abstracts the underlying protein sequence, structure or fold, (ii) allows the use of shape retrieval methods to screen large databases of protein structures to identify surficial homologs and possible interacting partners and (iii) opens an extension of the protein structure-function paradigm toward a protein structure-surface(s)-function paradigm. AVAILABILITYAND IMPLEMENTATION: All data are available online at http://datasetmachat.drugdesign.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discriminating Agonist from Antagonist Ligands of the Nuclear Receptors Using Different Chemoinformatics Approaches.

Nuclear receptors (NRs) constitute an important class of therapeutic targets. During the last 4 years, we tackled the pharmacological profile assessment of NR ligands for which we constructed the NRLiSt BDB. We evaluated and compared the performance of different virtual screening approaches: mean of molecular descriptor distribution values, molecular docking and 3D pharmacophore models. The simple comparison of the distribution profiles of 4885 molecular descriptors between the agonist and antagonist datasets didn't provide satisfying results. We obtained an overall good performance with the docking method we used, Surflex-Dock which was able to discriminate agonist from antagonist ligands. But the availability of PDB structures in the "pharmacological-profile-to-predict-bound-state" (agonist-bound or antagonist-bound) and the availability of enough ligands of both pharmacological profiles constituted limits to generalize this protocol for all NRs. Finally, the 3D pharmacophore modeling approach, allowed us to generate selective agonist pharmacophores and selective antagonist pharmacophores that covered more than 99 % of the whole NRLiSt BDB. This study allowed a better understanding of the pharmacological modulation of NRs with small molecules and could be extended to other therapeutic classes.

A new 3p25 locus is associated with liver fibrosis progression in human immunodeficiency virus/hepatitis C virus-coinfected patients.

There is growing evidence that human genetic variants contribute to liver fibrosis in subjects with hepatitis C virus (HCV) monoinfection, but this aspect has been little investigated in patients coinfected with HCV and human immunodeficiency virus (HIV). We performed the first genome-wide association study of liver fibrosis progression in patients coinfected with HCV and HIV, using the well-characterized French National Agency for Research on AIDS and Viral Hepatitis CO13 HEPAVIH cohort. Liver fibrosis was assessed by elastography (FibroScan), providing a quantitative fibrosis score. After quality control, a genome-wide association study was conducted on 289 Caucasian patients, for a total of 8,426,597 genotyped (Illumina Omni2.5 BeadChip) or reliably imputed single-nucleotide polymorphisms. Single-nucleotide polymorphisms with P values <10-6 were investigated in two independent replication cohorts of European patients infected with HCV alone. Two signals of genome-wide significance (P < 5 × 10-8 ) were obtained. The first, on chromosome 3p25 and corresponding to rs61183828 (P = 3.8 × 10-9 ), was replicated in the two independent cohorts of patients with HCV monoinfection. The cluster of single-nucleotide polymorphisms in linkage disequilibrium with rs61183828 was located close to two genes involved in mechanisms affecting both cell signaling and cell structure (CAV3) or HCV replication (RAD18). The second signal, obtained with rs11790131 (P = 9.3 × 10-9 ) on chromosome region 9p22, was not replicated. CONCLUSION: This genome-wide association study identified a new locus associated with liver fibrosis severity in patients with HIV/HCV coinfection, on chromosome 3p25, a finding that was replicated in patients with HCV monoinfection; these results provide new relevant hypotheses for the pathogenesis of liver fibrosis in patients with HIV/HCV coinfection that may help define new targets for drug development or new prognostic tests, to improve patient care. (Hepatology 2016;64:1462-1472).

A genome-wide association study in Caucasian women suggests the involvement of HLA genes in the severity of facial solar lentigines.

Solar lentigines are a common feature of sun-induced skin ageing. Little is known, however, about the genetic factors contributing to their development. In this genome-wide association study, we aimed to identify genetic loci associated with solar lentigines on the face in 502 middle-aged French women. Nine SNPs, gathered in two independent blocks on chromosome 6, exhibited a false discovery rate below 25% when looking for associations with the facial lentigine score. The first block, in the 6p22 region, corresponded to intergenic SNPs and also exhibited a significant association with forehead lentigines (P = 1.37 × 10(-8) ). The second block, within the 6p21 HLA region, was associated with decreased HLA-C expression according to several eQTL databases. Interestingly, these SNPs were also in high linkage disequilibrium with the HLA-C*0701 allele (r(2)  = 0.95). We replicated an association recently found by GWAS in the IRF4 gene. Finally, a complementary study on 44 selected candidate SNPs revealed novel associations in the MITF gene. Overall, our results point to several mechanisms involved in the severity of facial lentigines, including HLA/immunity and the melanogenesis pathway.

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

BACKGROUND: Many genome-wide association studies have been performed on progression towards the acquired immune deficiency syndrome (AIDS) and they mainly identified associations within the HLA loci. In this study, we demonstrate that the integration of biological information, namely gene expression data, can enhance the sensitivity of genetic studies to unravel new genetic associations relevant to AIDS. METHODS: We collated the biological information compiled from three databases of expression quantitative trait loci (eQTLs) involved in cells of the immune system. We derived a list of single nucleotide polymorphisms (SNPs) that are functional in that they correlate with differential expression of genes in at least two of the databases. We tested the association of those SNPs with AIDS progression in two cohorts, GRIV and ACS. Tests on permuted phenotypes of the GRIV and ACS cohorts or on randomised sets of equivalent SNPs allowed us to assess the statistical robustness of this method and to estimate the true positive rate. RESULTS: Eight genes were identified with high confidence (p = 0.001, rate of true positives 75%). Some of those genes had previously been linked with HIV infection. Notably, ENTPD4 belongs to the same family as CD39, whose expression has already been associated with AIDS progression; while DNAJB12 is part of the HSP90 pathway, which is involved in the control of HIV latency. Our study also drew our attention to lesser-known functions such as mitochondrial ribosomal proteins and a zinc finger protein, ZFP57, which could be central to the effectiveness of HIV infection. Interestingly, for six out of those eight genes, down-regulation is associated with non-progression, which makes them appealing targets to develop drugs against HIV.

Evidence after imputation for a role of MICA variants in nonprogression and elite control of HIV type 1 infection.

Past genome-wide association studies (GWAS) involving individuals with AIDS have mainly identified associations in the HLA region. Using the latest software, we imputed 7 million single-nucleotide polymorphisms (SNPs)/indels of the 1000 Genomes Project from the GWAS-determined genotypes of individuals in the Genomics of Resistance to Immunodeficiency Virus AIDS nonprogression cohort and compared them with those of control cohorts. The strongest signals were in MICA, the gene encoding major histocompatibility class I polypeptide-related sequence A (P = 3.31 × 10(-12)), with a particular exonic deletion (P = 1.59 × 10(-8)) in full linkage disequilibrium with the reference HCP5 rs2395029 SNP. Haplotype analysis also revealed an additive effect between HLA-C, HLA-B, and MICA variants. These data suggest a role for MICA in progression and elite control of human immunodeficiency virus type 1 infection.

NRLiSt BDB, the manually curated nuclear receptors ligands and structures benchmarking database.

Nuclear receptors (NRs) constitute an important class of drug targets. We created the most exhaustive NR-focused benchmarking database to date, the NRLiSt BDB (NRs ligands and structures benchmarking database). The 9905 compounds and 339 structures of the NRLiSt BDB are ready for structure-based and ligand-based virtual screening. In the present study, we detail the protocol used to generate the NRLiSt BDB and its features. We also give some examples of the errors that we found in ChEMBL that convinced us to manually review all original papers. Since extensive and manually curated experimental data about NR ligands and structures are provided in the NRLiSt BDB, it should become a powerful tool to assess the performance of virtual screening methods on NRs, to assist the understanding of NR's function and modulation, and to support the discovery of new drugs targeting NRs. NRLiSt BDB is freely available online at http://nrlist.drugdesign.fr .

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.

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⁻¹5, 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⁻°7, OR = 2.32). The fourth association corresponded to rs2072255 located in RICH2 (P = 3.30 × 10⁻°6, 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.

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.

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.

Exhaustive genotyping of the interleukin-1 family genes and associations with AIDS progression in a French cohort.

Interleukin (IL)-1 family members are key players in inflammatory processes but have been the subject of few studies of acquired immunodeficiency syndrome (AIDS). To better evaluate the impact of the IL-1 family on AIDS development, we genotyped the IL1 alpha , IL1 beta , IL1Ra, and IL1R1 genes in 245 slow progressor (SP) and 82 rapid progressor (RP) human immunodeficiency virus type 1-seropositive patients as well as in 446 control subjects, all of whom were of white ethnicity. One hundred sixteen frequent polymorphisms were identified, of which 23 were newly characterized by our study. Many putative associations were found between single-nucleotide polymorphism (SNP) or haplotype alleles and the extreme profiles of progression. Most of them corresponded to weak associations (.01

Associations of the IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFN (gamma) R1 cytokine receptor genes with AIDS progression in a French AIDS cohort.

We have performed an extensive analysis of Th1/Th2 cytokine receptors IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1 gene polymorphisms to evaluate their impact on AIDS progression. The coding regions and promoters of these genes were sequenced in the genetics of resistance to immunodeficiency virus cohort, composed of 327 HIV-1-positive patients with extreme progression phenotypes, slow and rapid progressors, and of 446 healthy control subjects, all of them of Caucasian descent. Overall, 104 single nucleotide polymorphisms and four insertions/deletions with a minor allelic frequency higher than 1% were identified, 21 of them being newly characterized. We observed weak associations for 13 polymorphisms of IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1, and 11 haplotypes of IL2Ralpha, IL4Ralpha, and IFNgammaR1. However, we could not relate these positive signals to any relevant biological information on the gene function. To affirm these putative associations in AIDS, further confirmation on other AIDS cohorts will be needed. This complete catalog of polymorphisms in IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1 cytokine receptor genes should also be useful for investigating associations in other immune-related diseases.

Genomic approach of AIDS pathogenesis: exhaustive genotyping of the TNFR1 gene in a French AIDS cohort.

Large-scale genomic studies in cohorts have been made possible for the last few years thanks to the progress of molecular biology and bioinformatics. This systematic approach allows a better understanding of the molecular mechanisms of disease development and as a consequence can contribute to the rational design of new diagnostic and therapeutic tools. We present here the exhaustive genotyping of a candidate gene, tumor necrosis factor receptor 1 (TNFR1), in the genetic of resistance to immunodeficiency virus (GRIV) AIDS cohort. This gene was chosen because it is likely to be involved in the apoptosis pathways of CD4+ and CD8+ T-cells during human immunodeficiency virus 1 (HIV-1) infection. Seven frequent polymorphisms were characterized in 319 HIV-1 seropositive patients from the GRIV cohort with extreme disease progression phenotypes, slow progression or rapid progression, and in 427 healthy controls. The TNFR1 gene locus does not appear to be part of any haploblock and contains only a small haploblock of two successive SNPs. One promoter SNP (TNFR1_17444594, position -581) and one intronic SNP (TNFR1_27223241, position +11511) gave weak positive signals of association (resp. P=0.03 and P=0.04) as well as two haplotypes. To our knowledge, this is the first genetic association study dealing with the TNFR1 gene in AIDS and the putative associations identified will need to be validated through other AIDS cohort analyses or by further biological experimentation.