An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant.

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

Wasl is crucial to maintain microglial core activities during glioblastoma initiation stages.

Microglia actively promotes the growth of high-grade gliomas. Within the glioma microenvironment an amoeboid microglial morphology has been observed, however the underlying causes and the related impact on microglia functions and their tumor promoting activities is unclear. Using the advantages of the larval zebrafish model, we identified the underlying mechanism and show that microglial morphology and functions are already impaired during glioma initiation stages. The presence of pre-neoplastic HRasV12 expressing cells induces an amoeboid morphology of microglia, increases microglial numbers and decreases their motility and phagocytic activity. RNA sequencing analysis revealed lower expression levels of the actin nucleation promoting factor wasla in microglia. Importantly, a microglia specific rescue of wasla expression restores microglial morphology and functions. This results in increased phagocytosis of pre-neoplastic cells and slows down tumor progression. In conclusion, we identified a mechanism that de-activates core microglial functions within the emerging glioma microenvironment. Restoration of this mechanism might provide a way to impair glioma growth.

Gene expression profiling reveals a conserved microglia signature in larval zebrafish.

Microglia are the resident macrophages of the brain. Over the past decade, our understanding of the function of these cells has significantly improved. Microglia do not only play important roles in the healthy brain but are involved in almost every brain pathology. Gene expression profiling allowed to distinguish microglia from other macrophages and revealed that the full microglia signature can only be observed in vivo. Thus, animal models are irreplaceable to understand the function of these cells. One of the popular models to study microglia is the zebrafish larva. Due to their optical transparency and genetic accessibility, zebrafish larvae have been employed to understand a variety of microglia functions in the living brain. Here, we performed RNA sequencing of larval zebrafish microglia at different developmental time points: 3, 5, and 7 days post fertilization (dpf). Our analysis reveals that larval zebrafish microglia rapidly acquire the core microglia signature and many typical microglia genes are expressed from 3 dpf onwards. The majority of changes in gene expression happened between 3 and 5 dpf, suggesting that differentiation mainly takes place during these days. Furthermore, we compared the larval microglia transcriptome to published data sets of adult zebrafish microglia, mouse microglia, and human microglia. Larval microglia shared a significant number of expressed genes with their adult counterparts in zebrafish as well as with mouse and human microglia. In conclusion, our results show that larval zebrafish microglia mature rapidly and express the core microglia gene signature that seems to be conserved across species.

A genome wide association study identifies new genes potentially associated with eyelid sagging.

Sagging eyelid is considered as an outward of skin ageing and may cause medical issues. However, little is known about the factors involved in sagging eyelid. The study, which aims at determining genetic risk factors for eyelid sagging, was conducted in a cohort of 502 unrelated Caucasian women living in the Paris region. All included participants were aged between 44 and 70 years old (mean age, 57.6 years old). The severity of sagging eyelid was graded in 6 categories by a dermatologist using standardized photographs of the face. A genome wide association study adjusted on potential risk factors (including age and smoking habits) was conducted to identify genetic associations. Two single nucleotide polymorphisms in total linkage disequilibrium on chromosome 10, rs16927253 (P = 7.07 × 10-10 ) and rs4746957 (P = 1.06 × 10-8 ), were significantly associated with eyelid sagging severity. The rs16927253-T and rs4746957-A alleles showed a dominant protective effect towards eyelid sagging. These polymorphisms are located in intronic parts of the H2AFY2 gene which encodes a member of the H2A histone family and very close to the AIFM2 gene that induces apoptosis. Additionally, single nucleotide polymorphisms with a false discovery rate below 0.25 were located nearby the type XIII collagen COL13A1 gene on chromosome 10 and in the ADAMTS18 gene on chromosome 16. Several relevant genes were identified by the genome wide association study for their potential role in the sagging eyelid severity.

Argonaute proteins regulate HIV-1 multiply spliced RNA and viral production in a Dicer independent manner.

Argonaute (Ago) proteins associate with microRNAs (miRNAs) to form the core of the RNA-induced silencing complex (RISC) that mediates post-transcriptional gene silencing of target mRNAs. As key players in anti-viral defense, Ago proteins are thought to have the ability to interact with human immunodeficiency virus type 1 (HIV-1) RNA. However, the role of this interaction in regulating HIV-1 replication has been debated. Here, we used high throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to explore the interaction between Ago2 and HIV-1 RNA in infected cells. By only considering reads of 50 nucleotides length in our analysis, we identified more than 30 distinct binding sites for Ago2 along the viral RNA genome. Using reporter assays, we found four binding sites, located near splice donor sites, capable of repressing Luciferase gene expression in an Ago-dependent manner. Furthermore, inhibition of Ago1 and Ago2 levels in cells expressing HIV-1 led to an increase of viral multiply spliced transcripts and to a strong reduction in the extracellular CAp24 level. Depletion of Dicer did not affect these activities. Our results highlight a new role of Ago proteins in the control of multiply spliced HIV-1 transcript levels and viral production, independently of the miRNA pathway.

Predicting the affinity of Farnesoid X Receptor ligands through a hierarchical ranking protocol: a D3R Grand Challenge 2 case study.

The Drug Design Data Resource (D3R) Grand Challenges are blind contests organized to assess the state-of-the-art methods accuracy in predicting binding modes and relative binding free energies of experimentally validated ligands for a given target. The second stage of the D3R Grand Challenge 2 (GC2) was focused on ranking 102 compounds according to their predicted affinity for Farnesoid X Receptor. In this task, our workflow was ranked 5th out of the 77 submissions in the structure-based category. Our strategy consisted in (1) a combination of molecular docking using AutoDock 4.2 and manual edition of available structures for binding poses generation using SeeSAR, (2) the use of HYDE scoring for pose selection, and (3) a hierarchical ranking using HYDE and MM/GBSA. In this report, we detail our pose generation and ligands ranking protocols and provide guidelines to be used in a prospective computer aided drug design program.

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.

HLA-check: evaluating HLA data from SNP information.

BACKGROUND: The major histocompatibility complex (MHC) region of the human genome, and specifically the human leukocyte antigen (HLA) genes, play a major role in numerous human diseases. With the recent progress of sequencing methods (eg, Next-Generation Sequencing, NGS), the accurate genotyping of this region has become possible but remains relatively costly. In order to obtain the HLA information for the millions of samples already genotyped by chips in the past ten years, efficient bioinformatics tools, such as SNP2HLA or HIBAG, have been developed that infer HLA information from the linkage disequilibrium existing between HLA alleles and SNP markers in the MHC region. RESULTS: In this study, we first used ShapeIT and Impute2 to implement an imputation method akin to SNP2HLA and found a comparable quality of imputation on a European dataset. More importantly, we developed a new tool, HLA-check, that allows for the detection of aberrant HLA allele calling with regard to the SNP genotypes in the region. Adding this tool to the HLA imputation software increases dramatically their accuracy, especially for HLA class I genes. CONCLUSION: Overall, HLA-check was able to identify a limited number of implausible HLA typings (less than 10%) in a population, and these samples can then either be removed or be retyped by NGS for HLA association analysis.

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 general approach for haplotype phasing across the full spectrum of relatedness.

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Haplotype estimation using sequencing reads.

High-throughput sequencing technologies produce short sequence reads that can contain phase information if they span two or more heterozygote genotypes. This information is not routinely used by current methods that infer haplotypes from genotype data. We have extended the SHAPEIT2 method to use phase-informative sequencing reads to improve phasing accuracy. Our model incorporates the read information in a probabilistic model through base quality scores within each read. The method is primarily designed for high-coverage sequence data or data sets that already have genotypes called. One important application is phasing of single samples sequenced at high coverage for use in medical sequencing and studies of rare diseases. Our method can also use existing panels of reference haplotypes. We tested the method by using a mother-father-child trio sequenced at high-coverage by Illumina together with the low-coverage sequence data from the 1000 Genomes Project (1000GP). We found that use of phase-informative reads increases the mean distance between switch errors by 22% from 274.4 kb to 328.6 kb. We also used male chromosome X haplotypes from the 1000GP samples to simulate sequencing reads with varying insert size, read length, and base error rate. When using short 100 bp paired-end reads, we found that using mixtures of insert sizes produced the best results. When using longer reads with high error rates (5-20 kb read with 4%-15% error per base), phasing performance was substantially improved.

JASSA: a comprehensive tool for prediction of SUMOylation sites and SIMs.

MOTIVATION: Post-translational modification by the Small Ubiquitin-like Modifier (SUMO) proteins, a process termed SUMOylation, is involved in many fundamental cellular processes. SUMO proteins are conjugated to a protein substrate, creating an interface for the recruitment of cofactors harboring SUMO-interacting motifs (SIMs). Mapping both SUMO-conjugation sites and SIMs is required to study the functional consequence of SUMOylation. To define the best candidate sites for experimental validation we designed JASSA, a Joint Analyzer of SUMOylation site and SIMs. RESULTS: JASSA is a predictor that uses a scoring system based on a Position Frequency Matrix derived from the alignment of experimental SUMOylation sites or SIMs. Compared with existing web-tools, JASSA displays on par or better performances. Novel features were implemented towards a better evaluation of the prediction, including identification of database hits matching the query sequence and representation of candidate sites within the secondary structural elements and/or the 3D fold of the protein of interest, retrievable from deposited PDB files. AVAILABILITY AND IMPLEMENTATION: JASSA is freely accessible at http://www.jassa.fr/. Website is implemented in PHP and MySQL, with all major browsers supported. CONTACT: guillaume.beauclair@inserm.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Stress-induced neuroinflammation is mediated by GSK3-dependent TLR4 signaling that promotes susceptibility to depression-like behavior.

Most psychiatric and neurological diseases are exacerbated by stress. Because this may partially result from stress-induced inflammation, we examined factors involved in this stress response. After a paradigm of inescapable foot shock stress that causes learned helplessness depression-like behavior, eighteen cytokines and chemokines increased in mouse hippocampus, peaking 6-12h after stress. A 24h prior pre-conditioning stress accelerated the rate of stress-induced hippocampal cytokine and chemokine increases, with most reaching peak levels after 1-3h, often without altering the maximal levels. Toll-like receptor 4 (TLR4) was involved in this response because most stress-induced hippocampal cytokines and chemokines were attenuated in TLR4 knockout mice. Stress activated glycogen synthase kinase-3 (GSK3) in wild-type mouse hippocampus, but not in TLR4 knockout mice. Administration of the antidepressant fluoxetine or the GSK3 inhibitor TDZD-8 reduced the stress-induced increases of most hippocampal cytokines and chemokines. Stress increased hippocampal levels of the danger-associated molecular pattern (DAMP) protein high mobility group box 1 (HMGB1), activated the inflammatory transcription factor NF-κB, and the NLRP3 inflammasome. Knockdown of HMGB1 blocked the acceleration of cytokine and chemokine increases in the hippocampus caused by two successive stresses. Fluoxetine treatment blocked stress-induced up-regulation of HMGB1 and subsequent NF-κB activation, whereas TDZD-8 administration attenuated NF-κB activation downstream of HMGB1. To test if stress-induced cytokines and chemokines contribute to depression-like behavior, the learned helplessness model was assessed. Antagonism of TNFα modestly reduced susceptibility to learned helplessness induction, whereas TLR4 knockout mice were resistant to learned helplessness. Thus, stress-induces a broad inflammatory response in mouse hippocampus that involves TLR4, GSK3, and downstream inflammatory signaling, and these stress responses contribute to susceptibility to depression-like behavior in mice.

Screening Explorer-An Interactive Tool for the Analysis of Screening Results.

Screening Explorer is a web-based application that allows for an intuitive evaluation of the results of screening experiments using complementary metrics in the field. The usual evaluation of screening results implies the separate generation and apprehension of the ROC, predictiveness, and enrichment curves and their global metrics. Similarly, partial metrics need to be calculated repeatedly for different fractions of a data set and there exists no handy tool that allows reading partial metrics simultaneously on different charts. For a deeper understanding of the results of screening experiments, we rendered their analysis straightforward by linking these metrics interactively in an interactive usable web-based application. We also implemented simple consensus scoring methods based on scores normalization, standardization (z-scores), and compounds ranking to evaluate the enrichments that can be expected through methods combination. Two demonstration data sets allow the users to easily apprehend the functions of this tool that can be applied to the analysis of virtual and experimental screening results. Screening Explorer is freely accessible at http://stats.drugdesign.fr .

DRAM triggers lysosomal membrane permeabilization and cell death in CD4(+) T cells infected with HIV.

Productive HIV infection of CD4(+) T cells leads to a caspase-independent cell death pathway associated with lysosomal membrane permeabilization (LMP) and cathepsin release, resulting in mitochondrial outer membrane permeabilization (MOMP). Herein, we demonstrate that HIV infection induces damage-regulated autophagy modulator (DRAM) expression in a p53-dependent manner. Knocking down the expression of DRAM and p53 genes with specific siRNAs inhibited autophagy and LMP. However, inhibition of Atg5 and Beclin genes that prevents autophagy had a minor effect on LMP and cell death. The knock down of DRAM gene inhibited cytochrome C release, MOMP and cell death. However, knocking down DRAM, we increased viral infection and production. Our study shows for the first time the involvement of DRAM in host-pathogen interactions, which may represent a mechanism of defense via the elimination of infected cells.

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.

ALDsuite: Dense marker MALD using principal components of ancestral linkage disequilibrium.

BACKGROUND: Mapping by admixture linkage disequilibrium (MALD) is a whole genome gene mapping method that uses LD from extended blocks of ancestry inherited from parental populations among admixed individuals to map associations for diseases, that vary in prevalence among human populations. The extended LD queried for marker association with ancestry results in a greatly reduced number of comparisons compared to standard genome wide association studies. As ancestral population LD tends to confound the analysis of admixture LD, the earliest algorithms for MALD required marker sets sufficiently sparse to lack significant ancestral LD between markers. However current genotyping technologies routinely provide dense SNP data, which convey more information than sparse sets, if this information can be efficiently used. There are currently no software solutions that offer both local ancestry inference using dense marker data and disease association statistics. RESULTS: We present here an R package, ALDsuite, which accounts for local LD using principal components of haplotypes from surrogate ancestral population data, and includes tools for quality control of data, MALD, downstream analysis of results and visualization graphics. CONCLUSIONS: ALDsuite offers a fast, accurate estimation of global and local ancestry and comes bundled with the tools needed for MALD, from data quality control through mapping of and visualization of disease genes.

Benchmarking Data Sets for the Evaluation of Virtual Ligand Screening Methods: Review and Perspectives.

Virtual screening methods are commonly used nowadays in drug discovery processes. However, to ensure their reliability, they have to be carefully evaluated. The evaluation of these methods is often realized in a retrospective way, notably by studying the enrichment of benchmarking data sets. To this purpose, numerous benchmarking data sets were developed over the years, and the resulting improvements led to the availability of high quality benchmarking data sets. However, some points still have to be considered in the selection of the active compounds, decoys, and protein structures to obtain optimal benchmarking data sets.

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.