Translational profiling of mouse dopaminoceptive neurons reveals region-specific gene expression, exon usage, and striatal prostaglandin E2 modulatory effects.

Forebrain dopamine-sensitive (dopaminoceptive) neurons play a key role in movement, action selection, motivation, and working memory. Their activity is altered in Parkinson's disease, addiction, schizophrenia, and other conditions, and drugs that stimulate or antagonize dopamine receptors have major therapeutic applications. Yet, similarities and differences between the various neuronal populations sensitive to dopamine have not been systematically explored. To characterize them, we compared translating mRNAs in the dorsal striatum and nucleus accumbens neurons expressing D1 or D2 dopamine receptor and prefrontal cortex neurons expressing D1 receptor. We identified genome-wide cortico-striatal, striatal D1/D2 and dorso/ventral differences in the translating mRNA and isoform landscapes, which characterize dopaminoceptive neuronal populations. Expression patterns and network analyses identified novel transcription factors with presumptive roles in these differences. Prostaglandin E2 (PGE2) was a candidate upstream regulator in the dorsal striatum. We pharmacologically explored this hypothesis and showed that misoprostol, a PGE2 receptor agonist, decreased the excitability of D2 striatal projection neurons in slices, and diminished their activity in vivo during novel environment exploration. We found that misoprostol also modulates mouse behavior including by facilitating reversal learning. Our study provides powerful resources for characterizing dopamine target neurons, new information about striatal gene expression patterns and regulation. It also reveals the unforeseen role of PGE2 in the striatum as a potential neuromodulator and an attractive therapeutic target.

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.

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.

Operant Training for Highly Palatable Food Alters Translating Messenger RNA in Nucleus Accumbens D2 Neurons and Reveals a Modulatory Role of Ncdn.

BACKGROUND: Highly palatable food triggers behavioral responses including strong motivation. These effects involve the reward system and dopamine neurons, which modulate neurons in the nucleus accumbens (NAc). The molecular mechanisms underlying the long-lasting effects of highly palatable food on feeding behavior are poorly understood. METHODS: We studied the effects of 2-week operant conditioning of mice with standard or isocaloric highly palatable food. We investigated the behavioral responses and dendritic spine modifications in the NAc. We compared the translating messenger RNA in NAc neurons identified by the type of dopamine receptors they express, depending on the kind of food and training. We tested the consequences of invalidation of an abundant downregulated gene, Ncdn. RESULTS: Operant conditioning for highly palatable food increased motivation for food even in well-fed mice. In wild-type mice, free choice between regular and highly palatable food increased weight compared with access to regular food only. Highly palatable food increased spine density in the NAc. In animals trained for highly palatable food, translating messenger RNAs were modified in NAc neurons expressing dopamine D2 receptors, mostly corresponding to striatal projection neurons, but not in neurons expressing D1 receptors. Knockout of Ncdn, an abundant downregulated gene, opposed the conditioning-induced changes in satiety-sensitive feeding behavior and apparent motivation for highly palatable food, suggesting that downregulation may be a compensatory mechanism. CONCLUSIONS: Our results emphasize the importance of messenger RNA alterations in D2 striatal projection neurons in the NAc in the behavioral consequences of highly palatable food conditioning and suggest a modulatory contribution of Ncdn downregulation.

Operant training for highly palatable food alters translating mRNA in nucleus accumbens D2 neurons and reveals a modulatory role of Neurochondrin.

BACKGROUND: Highly palatable food triggers behavioral alterations reminiscent of those induced by addictive drugs. These effects involve the reward system and dopamine neurons, which modulate neurons in the nucleus accumbens (NAc). The molecular mechanisms underlying the effects of highly palatable food on feeding behavior are poorly understood. METHODS: We studied the effects of 2-week operant conditioning of mice with standard or isocaloric highly palatable food. We investigated the behavioral effects and dendritic spine modifications in the NAc. We compared the translating mRNA in NAc neurons identified by the type of dopamine receptors they express, depending on the type of food and training. We tested the consequences of invalidation of an abundant downregulated gene, Ncdn (Neurochondrin). RESULTS: Operant conditioning for highly palatable food increases motivation for food even in well-fed mice. In control mice, free access to regular or highly palatable food results in increased weight as compared to regular food only. Highly palatable food increases spine density in the NAc. In animals trained for highly palatable food, translating mRNAs are modified in NAc dopamine D2-receptor-expressing neurons, mostly corresponding to striatal projection neurons, but not in those expressing D1-receptors. Knock-out of Ncdn, an abundant down-regulated gene, opposes the conditioning-induced changes in satiety-sensitive feeding behavior and apparent motivation for highly palatable food, suggesting down-regulation may be a compensatory mechanism. CONCLUSIONS: Our results emphasize the importance of mRNA alterations D2 striatal projection neurons in the NAc in the behavioral consequences of highly palatable food conditioning and suggest a modulatory contribution of Ncdn downregulation.

DNA methylation and hydroxymethylation characterize the identity of D1 and D2 striatal projection neurons.

Neuronal DNA modifications differ from those in other cells, including methylation outside CpG context and abundant 5-hydroxymethylation whose relevance for neuronal identities are unclear. Striatal projection neurons expressing D1 or D2 dopamine receptors allow addressing this question, as they share many characteristics but differ in their gene expression profiles, connections, and functional roles. We compare translating mRNAs and DNA modifications in these two populations. DNA methylation differences occur predominantly in large genomic clusters including differentially expressed genes, potentially important for D1 and D2 neurons. Decreased gene body methylation is associated with higher gene expression. Hydroxymethylation differences are more scattered and affect transcription factor binding sites, which can influence gene expression. We also find a strong genome-wide hydroxymethylation asymmetry between the two DNA strands, particularly pronounced at expressed genes and retrotransposons. These results identify novel properties of neuronal DNA modifications and unveil epigenetic characteristics of striatal projection neurons heterogeneity.

Multiple-cohort genetic association study reveals CXCR6 as a new chemokine receptor involved in long-term nonprogression to AIDS.

BACKGROUND: The compilation of previous genomewide association studies of AIDS shows a major polymorphism in the HCP5 gene associated with both control of the viral load and long-term nonprogression (LTNP) to AIDS. METHODS: To look for genetic variants that affect LTNP without necessary control of the viral load, we reanalyzed the genomewide data of the unique LTNP Genomics of Resistance to Immunodeficiency Virus (GRIV) cohort by excluding "elite controller" patients, who were controlling the viral load at very low levels (<100 copies/mL). RESULTS: The rs2234358 polymorphism in the CXCR6 gene was the strongest signal (P=2.5 x 10(-7); odds ratio, 1.85) obtained for the genomewide association study comparing the 186 GRIV LTNPs who were not elite controllers with 697 uninfected control subjects. This association was replicated in 3 additional independent European studies, reaching genomewide significance of P(combined)=9.7 x 10(-10). This association with LTNP is independent of the CCR2-CCR5 locus and the HCP5 polymorphisms. CONCLUSIONS: The statistical significance, the replication, and the magnitude of the association demonstrate that CXCR6 is likely involved in the molecular etiology of AIDS and, in particular, in LTNP, emphasizing the power of extreme-phenotype cohorts. CXCR6 is a chemokine receptor that is known as a minor coreceptor in human immunodeficiency virus type 1 infection but could participate in disease progression through its role as a mediator of inflammation.

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.

Enrichment or depletion of a GO category within a class of genes: which test?

MOTIVATION: A number of available program packages determine the significant enrichments and/or depletions of GO categories among a class of genes of interest. Whereas a correct formulation of the problem leads to a single exact null distribution, these GO tools use a large variety of statistical tests whose denominations often do not clarify the underlying P-value computations. SUMMARY: We review the different formulations of the problem and the tests they lead to: the binomial, chi2, equality of two probabilities, Fisher's exact and hypergeometric tests. We clarify the relationships existing between these tests, in particular the equivalence between the hypergeometric test and Fisher's exact test. We recall that the other tests are valid only for large samples, the test of equality of two probabilities and the chi2-test being equivalent. We discuss the appropriateness of one- and two-sided P-values, as well as some discreteness and conservatism issues. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

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 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.