Computation of haplotypes on SNPs subsets: advantage of the "global method"

BACKGROUND: Genetic association studies aim at finding correlations between a disease state and genetic variations such as SNPs or combinations of SNPs, termed haplotypes. Some haplotypes have a particular biological meaning such as the ones derived from SNPs located in the promoters, or the ones derived from non synonymous SNPs. All these haplotypes are "subhaplotypes" because they refer only to a part of the SNPs found in the gene. Until now, subhaplotypes were directly computed from the very SNPs chosen to constitute them, without taking into account the rest of the information corresponding to the other SNPs located in the gene. In the present work, we describe an alternative approach, called the "global method", which takes into account all the SNPs known in the region and compare the efficacy of the two "direct" and "global" methods. RESULTS: We used empirical haplotypes data sets from the GH1 promoter and the APOE gene, and 10 simulated datasets, and randomly introduced in them missing information (from 0% up to 20%) to compare the 2 methods. For each method, we used the PHASE haplotyping software since it was described to be the best. We showed that the use of the "global method" for subhaplotyping leads always to a better error rate than the classical direct haplotyping. The advantage provided by this alternative method increases with the percentage of missing genotyping data (diminution of the average error rate from 25% to less than 10%). We applied the global method software on the GRIV cohort for AIDS genetic associations and some associations previously identified through direct subhaplotyping were found to be erroneous. CONCLUSION: The global method for subhaplotyping can reduce, sometimes dramatically, the error rate on patient resolutions and haplotypes frequencies. One should thus use this method in order to minimise the risk of a false interpretation in genetic studies involving subhaplotypes. In practice the global method is always more efficient than the direct method, but a combination method taking into account the level of missing information in each subject appears to be even more interesting when the level of missing information becomes larger (>10%)

ISHAPE: new rapid and accurate software for haplotyping

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licens

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 infectio

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

IntroductionWe have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers.MethodsWe 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.ResultsOur 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.DiscussionThese 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

Filtering Redundancies For Sequence Similarity Search Programs.

Database scanning programs such as BLAST and FASTA are used nowadays by most biologists for the post-genomic processing of DNA or protein sequence information (in particular to retrieve the structure/function of uncharacterized proteins). Unfortunately, their results can be polluted by identical alignments (called redundancies) coming from the same protein or DNA sequences present in different entries of the database. This makes the efficient use of the listed alignments difficult. Pretreatment of databases has been proposed to suppress strictly identical entries. However, there still remain many identical alignments since redundancies may occur locally for entries corresponding to various fragments of the same sequence or for entries corresponding to very homologous sequences but differing at the level of a few residues such as ortholog proteins. In the present work, we show that redundant alignments can be indeed numerous even when working with a pretreated non-redundant data bank, going as high as 60% of the output results according to the query and the bank. Therefore the accuracy and the efficiency of the post-genomic work will be greatly increased if these redundancies are removed. To solve this up to now unaddressed problem, we have developed an algorithm that allows for the efficient and safe suppression of all the redundancies with no loss of information. This algorithm is based on various filtering steps that we describe here in the context of the Automat similarity search program, and such an algorithm should also be added to the other similarity search programs (BLAST, FASTA, etc...)

Measure of the impact of the number of bootstrap samples on the size and the relevance of the candidate haplotypes space for the APOE gene

<p><b>Copyright information:</b></p><p>Taken from "ISHAPE: new rapid and accurate software for haplotyping"</p><p>http://www.biomedcentral.com/1471-2105/8/205</p><p>BMC Bioinformatics 2007;8():205-205.</p><p>Published online 15 Jun 2007</p><p>PMCID:PMC1919397.</p><p></p> Black line and left scale are for ICR (capture rate of true haplotypes configurations). Dashed line and right scale are for ANCR (average number of candidates per genotype)

Exploration of associations between phospholipase A2 gene family polymorphisms and AIDS progression using the SNPlex method.

Members of the secreted phospholipase A2 (PLA2) protein family can inhibit HIV-1 virus replication in vitro. To evaluate the impact of PLA2 gene polymorphisms on AIDS disease development, we studied 12 family members using SNPlextrade mark technology that permitted simultaneous typing of 70 tagging Single Nucleotide Polymorphisms (tagSNPs). The study utilized HIV-1 seropositive donors with slow progressor (n=168) or rapid progressor (n=54) status, plus 355 control subjects. All donors were Caucasian (total 577 individuals). Genetic associations yielded mainly 0.0

Improvement of collagen-induced arthritis by active immunization against murine IL-1 beta peptides designed by molecular modelling

Interleukin-1 beta (IL-1 beta) is a crucial cytokine in inflammation processes and has been implicated in the pathogenesis of several chronic inflammatory diseases. Strategies designed to blocking IL- 10 by passive administration of inhibitors (mAbs, IL-1 receptor antagonist) have previously demonstrated efficacy in rheumatoid arthritis (RA). Using molecular modelling, we have defined three murine IL-1 beta peptide regions characterized by their close proximity to the receptor. Synthetic peptides corresponding to these regions, in cyclic and linear form, were delivered as immunogens in Swiss mice, resulting insignificant levels of autoantibodies directed against the native murine IL-1 beta cytokine as determined by ELISA and by an assay for neutralization of IL-1 beta biological activity. More importantly, one of the cyclic peptides showed a protective effect against inflammation and articular destruction in DBA/1 mouse collagen-induced arthritis, a model of RA. The high rate of success observed foractive immunization against cytokine peptides in vivo suggests that the in silico approach to autoantigen design may be a promising avenue for the development of anti-cytokine immunotherapeutics. (c) 2005 Elsevier Ltd. All rights reserved

Gene expression profiling of peripheral blood mononuclear cells from women with cervical lesions reveals new markers of cancer

International audienceCervical cancer (CC) is a multifactorial disease of which human papillomavirus (HPV) is the main etiological agent. Despite cervical Pap smear screening and anti‑HPV vaccination, CC remains a major public health issue. Identification of specific gene expression signatures in the blood could allow better insight into the immune response of CC and could provide valuable information for the development of novel biomarkers. The present study performed a transcriptomic analysis of peripheral blood mononuclear cells (PBMCs) from Senegalese patients with CC (n=31), low‑grade cervical intraepithelial neoplasia (CIN1; n=27) and from healthy control (CTR) subjects (n=29). Individuals in the CIN1 and CTR groups exhibited similar patterns in gene expression. A total of 182 genes were revealed to be differentially expressed in patients with CC compared with individuals in the CIN1 and CTR groups. The IL1R2, IL18R1, MMP9 and FKBP5 genes were the most upregulated, whereas the T‑cell receptor α gene TRA was the most downregulated in the CC group compared with in the CIN1 and CTR groups. The pathway enrichment analysis of the differentially expressed genes revealed pathways directly and indirectly linked to inflammation. To the best of our knowledge, the present study is the first large transcriptomic study on CC performed using PBMCs from African women; the results revealed the involvement of genes and pathways related to inflammation, most notably the IL‑1 pathway, and the involvement of downregulation of the T‑cell receptor α, a key component of the immune response. Several of the stated genes have already been reported in other cancer studies as putative blood biomarkers, thus reinforcing the requirement for deeper investigation. These findings may aid in the development of innovative clinical biomarkers for CC prevention and should be further replicated in other populations