VPS13A knockdown impairs corticostriatal synaptic plasticity and locomotor behavior in a new mouse model of chorea-acanthocytosis.

Chorea-acanthocytosis (ChAc) is an inherited neurodegenerative movement disorder caused by VPS13A gene mutations leading to the absence of protein expression. The striatum is the most affected brain region in ChAc patients. However, the study of the VPS13A function in the brain has been poorly addressed. Here we generated a VPS13A knockdown (KD) model and aimed to elucidate the contribution of VPS13A to synaptic plasticity and neuronal communication in the corticostriatal circuit. First, we infected primary cortical neurons with miR30-shRNA against VPS13A and analyzed its effects on neuronal plasticity. VPS13A-KD neurons showed a higher degree of branching than controls, accompanied by decreased BDNF and PSD-95 levels, indicative of synaptic alterations. We then injected AAV-KD bilaterally in the frontal cortex and two different regions of the striatum of mice and analyzed the effects of VPS13A-KD on animal behavior and synaptic plasticity. VPS13A-KD mice showed modification of the locomotor behavior pattern, with increased exploratory behavior and hyperlocomotion. Corticostriatal dysfunction in VPS13A-KD mice was evidenced by impaired striatal long-term depression (LTD) after stimulation of cortical afferents, which was partially recovered by BDNF administration. VPS13A-KD did not lead to neuronal loss in the cortex or the striatum but induced a decrease in the neuronal release of CX3CL1 and triggered a microglial reaction, especially in the striatum. Notably, CX3CL1 administration partially restored the impaired corticostriatal LTD in VPS13A-KD mice. Our results unveil the involvement of VPS13A in neuronal connectivity modifying BDNF and CX3CL1 release. Moreover, the involvement of VPS13A in synaptic plasticity and motor behavior provides key information to further understand not only ChAc pathophysiology but also other neurological disorders.

Urinary metabolite quantitative trait loci in children and their interaction with dietary factors.

Human metabolism is influenced by genetic and environmental factors. Previous studies have identified over 23 loci associated with more than 26 urine metabolites levels in adults, which are known as urinary metabolite quantitative trait loci (metabQTLs). The aim of the present study is the identification for the first time of urinary metabQTLs in children and their interaction with dietary patterns. Association between genome-wide genotyping data and 44 urine metabolite levels measured by proton nuclear magnetic resonance spectroscopy was tested in 996 children from the Human Early Life Exposome project. Twelve statistically significant urine metabQTLs were identified, involving 11 unique loci and 10 different metabolites. Comparison with previous findings in adults revealed that six metabQTLs were already known, and one had been described in serum and three were involved the same locus as other reported metabQTLs but had different urinary metabolites. The remaining two metabQTLs represent novel urine metabolite-locus associations, which are reported for the first time in this study [single nucleotide polymorphism (SNP) rs12575496 for taurine, and the missense SNP rs2274870 for 3-hydroxyisobutyrate]. Moreover, it was found that urinary taurine levels were affected by the combined action of genetic variation and dietary patterns of meat intake as well as by the interaction of this SNP with beverage intake dietary patterns. Overall, we identified 12 urinary metabQTLs in children, including two novel associations. While a substantial part of the identified loci affected urinary metabolite levels both in children and in adults, the metabQTL for taurine seemed to be specific to children and interacted with dietary patterns.

Diagnostic accuracy of cerebrospinal fluid biomarkers in genetic prion diseases.

Genetic prion diseases are a rare and diverse group of fatal neurodegenerative disorders caused by pathogenic sequence variations in the prion protein gene, PRNP. Data on CSF biomarkers in patients with genetic prion diseases are limited and conflicting results have been reported for unclear reasons. Here, we aimed to analyse the diagnostic accuracy of CSF biomarkers currently used in prion clinical diagnosis in 302 symptomatic genetic prion disease cases from 11 prion diagnostic centres, encompassing a total of 36 different pathogenic sequence variations within the open reading frame of PRNP. CSF samples were assessed for the surrogate markers of neurodegeneration, 14-3-3 protein (14-3-3), total-tau protein (t-tau) and α-synuclein and for prion seeding activity through the real-time quaking-induced conversion assay. Biomarker results were compared with those obtained in healthy and neurological controls. For the most prevalent PRNP pathogenic sequence variations, biomarker accuracy and associations between biomarkers, demographic and genetic determinants were assessed. Additionally, the prognostic value of biomarkers for predicting total disease duration from symptom onset to death was investigated. High sensitivity of the four biomarkers was detected for genetic Creutzfeldt-Jakob disease associated with the E200K and V210I mutations, but low sensitivity was observed for mutations associated with Gerstmann-Sträussler-Scheinker syndrome and fatal familial insomnia. All biomarkers showed good to excellent specificity using the standard cut-offs often used for sporadic Creutzfeldt-Jakob disease. In genetic prion diseases related to octapeptide repeat insertions, the biomarker sensitivity correlated with the number of repeats. New genetic prion disease-specific cut-offs for 14-3-3, t-tau and α-synuclein were calculated. Disease duration in genetic Creutzfeldt-Jakob disease-E200K, Gerstmann-Sträussler-Scheinker-P102L and fatal familial insomnia was highly dependent on PRNP codon 129 MV polymorphism and was significantly associated with biomarker levels. In a large cohort of genetic prion diseases, the simultaneous analysis of CSF prion disease biomarkers allowed the determination of new mutation-specific cut-offs improving the discrimination of genetic prion disease cases and unveiled genetic prion disease-specific associations with disease duration.

Blood miRNA levels associated with ADHD traits in children across six European birth cohorts.

BACKGROUND: Attention-deficit/hyperactivity disorder (ADHD) is a prevalent and highly heritable neurodevelopmental disorder of major societal concern. Diagnosis can be challenging and there are large knowledge gaps regarding its etiology, though studies suggest an interplay of genetic and environmental factors involving epigenetic mechanisms. MicroRNAs (miRNAs) show promise as biomarkers of human pathology and novel therapies, and here we aimed to identify blood miRNAs associated with traits of ADHD as possible biomarker candidates and further explore their biological relevance. METHODS: Our study population consisted of 1126 children (aged 5-12 years, 46% female) from the Human Early Life Exposome study, a study spanning six ongoing population-based European birth cohorts. Expression profiles of miRNAs in whole blood samples were quantified by microarray and tested for association with ADHD-related measures of behavior and neuropsychological functions from questionnaires (Conner's Rating Scale and Child Behavior Checklist) and computer-based tests (the N-back task and Attention Network Test). RESULTS: We identified 29 miRNAs significantly associated (false discovery rate < .05) with the Conner's questionnaire-rated trait hyperactivity, 15 of which have been linked to ADHD in previous studies. Investigation into their biological relevance revealed involvement in several pathways related to neurodevelopment and function, as well as being linked with other neurodevelopmental or psychiatric disorders known to overlap with ADHD both in symptomology, genetic risk, and co-occurrence, such as autism spectrum disorder or schizophrenia. An additional three miRNAs were significantly associated with Conner's-rated inattention. No associations were found with questionnaire-rated total ADHD index or with computer-based tests. CONCLUSIONS: The large overlap of our hyperactivity-associated miRNAs with previous studies on ADHD is intriguing and warrant further investigation. Though this study should be considered explorative and preliminary, these findings contribute towards identifying a set of miRNAs for use as blood-based biomarkers to aid in earlier and easier ADHD diagnosis.

Phospho-RNA-Seq Highlights Specific Small RNA Profiles in Plasma Extracellular Vesicles.

Small RNAs (sRNAs) are bioactive molecules that can be detected in biofluids, reflecting physiological and pathological states. In plasma, sRNAs are found within extracellular vesicles (EVs) and in extravesicular compartments, offering potential sources of highly sensitive biomarkers. Deep sequencing strategies to profile sRNAs favor the detection of microRNAs (miRNAs), the best-known class of sRNAs. Phospho-RNA-seq, through the enzymatic treatment of sRNAs with T4 polynucleotide kinase (T4-PNK), has been recently developed to increase the detection of thousands of previously inaccessible RNAs. In this study, we investigated the value of phospho-RNA-seq on both the EVs and extravesicular plasma subfractions. Phospho-RNA-seq increased the proportion of sRNAs used for alignment and highlighted the diversity of the sRNA transcriptome. Unsupervised clustering analysis using sRNA counts matrices correctly classified the EVs and extravesicular samples only in the T4-PNK treated samples, indicating that phospho-RNA-seq stresses the features of sRNAs in each plasma subfraction. Furthermore, T4-PNK treatment emphasized specific miRNA variants differing in the 5'-end (5'-isomiRs) and certain types of tRNA fragments in each plasma fraction. Phospho-RNA-seq increased the number of tissue-specific messenger RNA (mRNA) fragments in the EVs compared with the extravesicular fraction, suggesting that phospho-RNA-seq favors the discovery of tissue-specific sRNAs in EVs. Overall, the present data emphasizes the value of phospho-RNA-seq in uncovering RNA-based biomarkers in EVs.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan.

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Efficient and flexible Integration of variant characteristics in rare variant association studies using integrated nested Laplace approximation.

Rare variants are thought to play an important role in the etiology of complex diseases and may explain a significant fraction of the missing heritability in genetic disease studies. Next-generation sequencing facilitates the association of rare variants in coding or regulatory regions with complex diseases in large cohorts at genome-wide scale. However, rare variant association studies (RVAS) still lack power when cohorts are small to medium-sized and if genetic variation explains a small fraction of phenotypic variance. Here we present a novel Bayesian rare variant Association Test using Integrated Nested Laplace Approximation (BATI). Unlike existing RVAS tests, BATI allows integration of individual or variant-specific features as covariates, while efficiently performing inference based on full model estimation. We demonstrate that BATI outperforms established RVAS methods on realistic, semi-synthetic whole-exome sequencing cohorts, especially when using meaningful biological context, such as functional annotation. We show that BATI achieves power above 70% in scenarios in which competing tests fail to identify risk genes, e.g. when risk variants in sum explain less than 0.5% of phenotypic variance. We have integrated BATI, together with five existing RVAS tests in the 'Rare Variant Genome Wide Association Study' (rvGWAS) framework for data analyzed by whole-exome or whole genome sequencing. rvGWAS supports rare variant association for genes or any other biological unit such as promoters, while allowing the analysis of essential functionalities like quality control or filtering. Applying rvGWAS to a Chronic Lymphocytic Leukemia study we identified eight candidate predisposition genes, including EHMT2 and COPS7A.

Huntington's disease brain-derived small RNAs recapitulate associated neuropathology in mice.

Progressive motor alterations and selective death of striatal medium spiny neurons (MSNs) are key pathological hallmarks of Huntington's disease (HD), a neurodegenerative condition caused by a CAG trinucleotide repeat expansion in the coding region of the huntingtin (HTT) gene. Most research has focused on the pathogenic effects of the resultant protein product(s); however, growing evidence indicates that expanded CAG repeats within mutant HTT mRNA and derived small CAG repeat RNAs (sCAG) participate in HD pathophysiology. The individual contribution of protein versus RNA toxicity to HD pathophysiology remains largely uncharacterized and the role of other classes of small RNAs (sRNA) that are strongly perturbed in HD is uncertain. Here, we demonstrate that sRNA produced in the putamen of HD patients (HD-sRNA-PT) are sufficient to induce HD pathology in vivo. Mice injected with HD-sRNA-PT show motor abnormalities, decreased levels of striatal HD-related proteins, disruption of the indirect pathway, and strong transcriptional abnormalities, paralleling human HD pathology. Importantly, we show that the specific blockage of sCAG mitigates HD-sRNA-PT neurotoxicity only to a limited extent. This observation prompted us to identify other sRNA species enriched in HD putamen with neurotoxic potential. We detected high levels of tRNA fragments (tRFs) in HD putamen, and we validated the neurotoxic potential of an Alanine derived tRF in vitro. These results highlight that HD-sRNA-PT are neurotoxic, and suggest that multiple sRNA species contribute to striatal dysfunction and general transcriptomic changes, favoring therapeutic strategies based on the blockage of sRNA-mediated toxicity.

Allele balance bias identifies systematic genotyping errors and false disease associations.

In recent years, next-generation sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors, and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the AB distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies. Availability: https://github.com/Francesc-Muyas/ABB.

Specific small-RNA signatures in the amygdala at premotor and motor stages of Parkinson's disease revealed by deep sequencing analysis.

MOTIVATION: Most computational tools for small non-coding RNAs (sRNA) sequencing data analysis focus in microRNAs (miRNAs), overlooking other types of sRNAs that show multi-mapping hits. Here, we have developed a pipeline to non-redundantly quantify all types of sRNAs, and extract patterns of expression in biologically defined groups. We have used our tool to characterize and profile sRNAs in post-mortem brain samples of control individuals and Parkinson's disease (PD) cases at early-premotor and late-symptomatic stages. RESULTS: Clusters of co-expressed sRNAs mapping onto tRNAs significantly separated premotor and motor cases from controls. A similar result was obtained using a matrix of miRNAs slightly varying in sequence (isomiRs). The present framework revealed sRNA alterations at premotor stages of PD, which might reflect initial pathogenic perturbations. This tool may be useful to discover sRNA expression patterns linked to different biological conditions. AVAILABILITY AND IMPLEMENTATION: The full code is available at http://github.com/lpantano/seqbuster CONTACT: lpantano@hsph.harvard.edu or eulalia.marti@crg.eu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia.

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.

Distinguishability and agreement with continuous data.

The analysis of concordance among repeated measures has received a huge amount of attention in the statistical literature leading to a range of different approaches. However, because all the approaches are able to assess the closeness among the readings taken on the same subject, the conclusions about the degree of concordance should be similar regardless the approach applied. Here, two indices to assess the concordance among continuous repeated measures, the intraclass correlation coefficient and the total deviation index, are applied and compared in two case examples. The first example concerns the repeatability of individual nutrient allocation strategy assessed by stable isotope analysis. The second example dealt with the assessment of the concordance of functional magnetic resonance imaging data that shows spatial correlation. The results differ depending upon the approach applied leading to contradictory conclusions about the degree of concordance. The reason behind these results is discussed reaching the conclusion that the total deviation index is just assessing agreement among repeated measurements, whereas the intraclass correlation coefficient assesses the concept of distinguishability among subjects that involves agreement among repeated measurements and spread of subjects at once. Therefore, the best way to select the right approach is to understand the right question behind the research hypothesis.

Common genetic variants associated with urinary phthalate levels in children: A genome-wide study.

INTRODUCTION: Phthalates, or dieters of phthalic acid, are a ubiquitous type of plasticizer used in a variety of common consumer and industrial products. They act as endocrine disruptors and are associated with increased risk for several diseases. Once in the body, phthalates are metabolized through partially known mechanisms, involving phase I and phase II enzymes. OBJECTIVE: In this study we aimed to identify common single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) associated with the metabolism of phthalate compounds in children through genome-wide association studies (GWAS). METHODS: The study used data from 1,044 children with European ancestry from the Human Early Life Exposome (HELIX) cohort. Ten phthalate metabolites were assessed in a two-void pooled urine collected at the mean age of 8 years. Six ratios between secondary and primary phthalate metabolites were calculated. Genome-wide genotyping was done with the Infinium Global Screening Array (GSA) and imputation with the Haplotype Reference Consortium (HRC) panel. PennCNV was used to estimate copy number variants (CNVs) and CNVRanger to identify consensus regions. GWAS of SNPs and CNVs were conducted using PLINK and SNPassoc, respectively. Subsequently, functional annotation of suggestive SNPs (p-value < 1E-05) was done with the FUMA web-tool. RESULTS: We identified four genome-wide significant (p-value < 5E-08) loci at chromosome (chr) 3 (FECHP1 for oxo-MiNP_oh-MiNP ratio), chr6 (SLC17A1 for MECPP_MEHHP ratio), chr9 (RAPGEF1 for MBzP), and chr10 (CYP2C9 for MECPP_MEHHP ratio). Moreover, 115 additional loci were found at suggestive significance (p-value < 1E-05). Two CNVs located at chr11 (MRGPRX1 for oh-MiNP and SLC35F2 for MEP) were also identified. Functional annotation pointed to genes involved in phase I and phase II detoxification, molecular transfer across membranes, and renal excretion. CONCLUSION: Through genome-wide screenings we identified known and novel loci implicated in phthalate metabolism in children. Genes annotated to these loci participate in detoxification, transmembrane transfer, and renal excretion.

MicroRNA profiling of Parkinson's disease brains identifies early downregulation of miR-34b/c which modulate mitochondrial function.

MicroRNAs (miRNAs) are post-transcriptional gene expression regulators, playing key roles in neuronal development, plasticity and disease. Parkinson's disease (PD) is the second most common neurodegenerative disorder, characterized by the presence of protein inclusions or Lewy bodies and a progressive loss of dopaminergic neurons in the midbrain. Here, we have evaluated miRNA expression deregulation in PD brain samples. MiRNA expression profiling revealed decreased expression of miR-34b and miR-34c in brain areas with variable neuropathological affectation at clinical (motor) stages (Braak stages 4 and 5) of the disease, including the amygdala, frontal cortex, substantia nigra and cerebellum. Furthermore, misregulation of miR-34b/c was detected in pre-motor stages (stages 1-3) of the disease, and thus in cases that did not receive any PD-related treatment during life. Depletion of miR-34b or miR-34c in differentiated SH-SY5Y dopaminergic neuronal cells resulted in a moderate reduction in cell viability that was accompanied by altered mitochondrial function and dynamics, oxidative stress and reduction in total cellular adenosin triphosphate content. MiR-34b/c downregulation was coupled to a decrease in the expression of DJ1 and Parkin, two proteins associated to familial forms of PD that also have a role in idiopathic cases. Accordingly, DJ1 and Parkin expression was reduced in PD brain samples displaying strong miR-34b/c downregulation. We propose that early deregulation of miR-34b/c in PD triggers downstream transcriptome alterations underlying mitochondrial dysfunction and oxidative stress, which ultimately compromise cell viability. A better understanding of the cellular pathways controlling and/or controlled by miR-34b/c should allow identification of targets for development of therapeutic approaches.

Upregulation of a small vault RNA (svtRNA2-1a) is an early event in Parkinson disease and induces neuronal dysfunction.

MicroRNAs (miRNAs) and other small non-coding RNAs (sncRNAs) are post-transcriptional regulators of gene expression, playing key roles in neuronal development, plasticity, and disease. Transcriptome deregulation caused by miRNA dysfunction has been associated to neurodegenerative diseases. Parkinson disease (PD) is the second most common neurodegenerative disease showing deregulation of the coding and small non-coding transcriptome. On profiling sncRNA in PD brain areas differently affected, we found that upregulation of a small vault RNA (svtRNA2-1a) is widespread in PD brains, occurring early in the course of the disease (at pre-motor stages). SvtRNA2-1a biogenesis was dependent on Dicer activity on its precursor (vtRNA2-1) but independent of Drosha endonuclease, unlike the canonical miRNAs. Although endogenous svtRNA2-1a was enriched in Ago-2 immunoprecipitates in differentiated SH-SY5Y neuronal cells, overexpression of svtRNA2-1a induced subtle transcriptomic changes, suggesting that gene expression regulation may involve other mechanisms than mRNA decay only. Function enrichment analysis of the genes deregulated by svtRNA2-1a overexpression or svtRNA2-1a predicted targets identified pathways related to nervous system development and cell type specification. The expression pattern of svtRNA2-1a during development and aging of the human brain and the detrimental consequences of a svtRNA2-1a mimic overexpression in neuronal cells further indicate that low svtRNA2-1a levels may be important for the maintenance of neurons. Our results suggest that early svtRNA2-1a upregulation in PD may contribute to perturbations of gene expression networks, underlying metabolic impairment and cell dysfunction. A better understanding of the pathways regulated by svtRNA2-a, and also the mechanisms regulating its expression should facilitate the identification of new targets for therapeutic approaches in PD.

DNA hypomethylation at ALOX12 is associated with persistent wheezing in childhood.

RATIONALE: Epigenetic changes may play a role in the occurrence of asthma-related phenotypes. OBJECTIVES: To identify epigenetic marks in terms of DNA methylation of asthma-related phenotypes in childhood, and to assess the effect of prenatal exposures and genetic variation on these epigenetic marks. METHODS: Data came from two cohorts embedded in the Infancia y Medio Ambiente (INMA) PROJECT: Menorca (n = 122) and Sabadell (n = 236). Wheezing phenotypes were defined at age 4-6 years. Cytosine-guanine (CpG) dinucleotide site DNA methylation differences associated with wheezing phenotypes were screened in children of the Menorca study using the Illumina GoldenGate Panel I. Findings were validated and replicated using pyrosequencing. Information on maternal smoking and folate supplement use was obtained through questionnaires. Dichlorodiphenyldichloroethylene was measured in cord blood or maternal serum. Genotypes were extracted from genome-wide data. MEASUREMENT AND MAIN RESULTS: Screening identified lower DNA methylation at a CpG site in the arachidonate 12-lipoxygenase (ALOX12) gene in children having persistent wheezing compared with those never wheezed (P = 0.003). DNA hypomethylation at ALOX12 loci was associated with higher risk of persistent wheezing in the Menorca study (odds ratio per 1% methylation decrease, 1.13; 95% confidence interval, 0.99-1.29; P = 0.077) and in the Sabadell study (odds ratio, 1.16; 95% confidence interval, 1.03-1.37; P = 0.017). Higher levels of prenatal dichlorodiphenyldichloroethylene were associated with DNA hypomethylation of ALOX12 in the Menorca study (P = 0.033), but not in the Sabadell study (P = 0.377). ALOX12 DNA methylation was strongly determined by underlying genetic polymorphisms. CONCLUSIONS: DNA methylation of ALOX12 may be an epigenetic biomarker for the risk of asthma-related phenotypes.

Genome-Wide Association Study of Blood Mercury in European Pregnant Women and Children.

Mercury has high industrial utility and is present in many products, and environmental contamination and occupational exposure are widespread. There are numerous biological systems involved in the absorption, metabolism, and excretion of Hg, and it is possible that some systems may be impacted by genetic variation. If so, genotype may affect tissue concentrations of Hg and subsequent toxic effects. Genome-wide association testing was performed on blood Hg samples from pregnant women of the Avon Longitudinal Study of Parents and Children (n = 2893) and children of the Human Early Life Exposome (n = 1042). Directly-genotyped single-nucleotide polymorphisms (SNPs) were imputed to the Haplotype Reference Consortium r1.1 panel of whole genotypes and modelled againstlog-transformed Hg. Heritability was estimated using linkage disequilibrium score regression. The heritability of Hg was estimated as 24.0% (95% CI: 16.9% to 46.4%) in pregnant women, but could not be determined in children. There were 16 SNPs associated with Hg in pregnant women above a suggestive p-value threshold (p < 1 × 10-5), and 21 for children. However, no SNP passed this threshold in both studies, and none were genome-wide significant (p < 5 × 10-8). SNP-Hg associations were highly discordant between women and children, and this may reflect differences in metabolism, a gene-age interaction, or dose-response effects. Several suggestive variants had plausible links to Hg metabolism, such as rs146099921 in metal transporter SLC39A14, and two variants (rs28618224, rs7154700) in potassium voltage-gated channel genes. The findings would benefit from external validation, as suggestive results may contain both true associations and false positives.

Potentially causal associations between placental DNA methylation and schizophrenia and other neuropsychiatric disorders.

Increasing evidence supports the role of placenta in neurodevelopment and potentially, in the later onset of neuropsychiatric disorders. Recently, methylation quantitative trait loci (mQTL) and interaction QTL (iQTL) maps have proven useful to understand SNP-genome wide association study (GWAS) relationships, otherwise missed by conventional expression QTLs. In this context, we propose that part of the genetic predisposition to complex neuropsychiatric disorders acts through placental DNA methylation (DNAm). We constructed the first public placental cis-mQTL database including nearly eight million mQTLs calculated in 368 fetal placenta DNA samples from the INMA project, ran cell type- and gestational age-imQTL models and combined those data with the summary statistics of the largest GWAS on 10 neuropsychiatric disorders using Summary-based Mendelian Randomization (SMR) and colocalization. Finally, we evaluated the influence of the DNAm sites identified on placental gene expression in the RICHS cohort. We found that placental cis-mQTLs are highly enriched in placenta-specific active chromatin regions, and useful to map the etiology of neuropsychiatric disorders at prenatal stages. Specifically, part of the genetic burden for schizophrenia, bipolar disorder and major depressive disorder confers risk through placental DNAm. The potential causality of several of the observed associations is reinforced by secondary association signals identified in conditional analyses, regional pleiotropic methylation signals associated to the same disorder, and cell type-imQTLs, additionally associated to the expression levels of relevant immune genes in placenta. In conclusion, the genetic risk of several neuropsychiatric disorders could operate, at least in part, through DNAm and associated gene expression in placenta.

Genetic variants and abnormal processing of pre-miR-182, a circadian clock modulator, in major depression patients with late insomnia.

Previous studies in mice have reported five different microRNAs (miRNAs; miR-219-1/132/183/96/182) to be modulators of the endogenous circadian clock and have presented experimental evidence for some of the genes involved in the molecular clock machinery as target sites. Moreover, disruption of circadian rhythms has long been implicated in the pathophysiology of major depression (MD). We investigated these miRNAs and some of their target sites at the sequence and functional levels as possible predisposing factors for susceptibility to MD and related chronobiological subphenotypes. Mutational screening was performed in a sample of 359 MD patients and 341 control individuals. We found a significant association between the T allele of the rs76481776 polymorphism in the pre-miR-182 and late insomnia in MD patients. Previous studies have reported an association between insomnia and CLOCK gene, a predicted miR-182 target site. A significant overexpression of miR-182 was detected by quantitative real-time polymerase chain reaction in cells transfected with the mutated form of the pre-miR-182 when compared with wild-type form. Moreover, a significant reduction in luciferase activity of plasmids with 3' UTR of ADCY6, CLOCK and DSIP genes was shown when transfecting cells with the mutated form of pre-miR-182 compared with cells that did not express miR-182. These data indicate that abnormal processing of pre-miR-182 in patients carrying the T allele of the rs76481776 polymorphism may contribute to the dysregulation of circadian rhythms in MD patients with insomnia, which could influence expression levels of the mature form of miR-182 and might increase downregulation in some of its target genes.

Association between the NMDA glutamate receptor GRIN2B gene and obsessive-compulsive disorder.

BACKGROUND: Recent data from neuroimaging, genetic and clinical trials and animal models suggest a role for altered glutamatergic neuro transmission in the pathogenesis of obsessive-compulsive disorder (OCD). The aim of this study was to investigate whether variants in the GRIN2B gene, the gene encoding the NR2 subunit of the N-methyl-D-aspartate (NMDA) glutamate receptor, may contribute to genetic susceptibility to OCD or to different OCD subphenotypes. METHODS: Between 2003 and 2008, we performed a case-control association study in which we genotyped 10 tag single-nucleotide polymorphisms (SNPs) in the 3' untranslated region (3' UTR) of GRIN2B. We performed SNP association and haplotype analysis considering the OCD diagnosis and different OCD subphenotypes: early-onset OCD, comorbid tic disorders and OCD clinical symptom dimensions. RESULTS: We enrolled 225 patients with OCD and 279 controls recruited from the OCD Clinic at Bellvitge Hospital (Barcelona, Spain). No significant difference in the distribution of alleles or genotypes was detected between patients with OCD and controls. Nonetheless, on analyzing OCD subphenotypes, the rs1805476 SNP in male patients (95% confidence interval [CI] 1.37-4.22, p = 0.002) and a 4-SNP haplotype in the whole sample (rs1805476, rs1805501, rs1805502 and rs1805477; odds ratio 1.92, 95% CI 1.22-3.01; permutation p = 0.023) were significantly associated with the presence of contamination obsessions and cleaning compulsions. LIMITATIONS: Study limitations included the risk of population stratification associated with the case-control design, use of psychiatrically unscreened blood donors as the control group, reduced sample size of participants with certain OCD subphenotypes and tested polymorphisms limited to 3' UTR and exon 13 of GRIN2B. CONCLUSION: Our results converge with recent data suggesting a possible contribution of glutamatergic variants to the genetic vulnerability to OCD or at least to certain OCD manifestations. The dissection of OCD into more homogeneous subphenotypes may constitute a useful tool to disentangle the complex genetic basis of the disorder.