Search

Background Hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM) are serious inherited heart diseases with various causative mutations identified. The full spectrum of causative mutations remains to be discovered, especially in understudied populations. Methods Here, we established the DOHA Registry and Biobank for cardiomyopathies in Qatar, followed by sequencing of 174 genes on 51 HCM and 53 DCM patients, and 31 relatives. Results In HCM, the analysis of 25 HCM‐associated genes showed that 20% of HCM cases had putative pathogenic variants for cardiomyopathy, mainly in sarcomere genes. Additional 49% of HCM cases had variants of uncertain significance, while 31% of HCM cases had likely benign variant(s) or had no variants identified within the analyzed HCM genes. In DCM, 56 putative DCM genes were analyzed. Eight percent of DCM cases had putative pathogenic variants for DCM, in the TTN gene while 70% of cases had variants of uncertain significance, in the analyzed DCM genes, that will need further pathogenicity assessment. Moreover, 22% of DCM cases remain unexplained, by having likely benign variant(s) or having no variants detected in any of the analyzed DCM genes. Conclusion We identified or replicated at least four recurrent variants among cardiomyopathy patients, which could be founder disease mutations in the Arabic population, including a frameshift variant (c.1371_1381dupTATCCAGTTAT) of unknown significance in the FKTN gene which seems to cause DCM in homozygosity, and HCM in heterozygosity. In vivo and/or in vitro functional validation need to be pursued in order to assess the pathogenicity of the identified variants., We established the DOHA Registry and Biobank for cardiomyopathies in Qatar, followed by sequencing of 174 genes on 51 HCM and 53 DCM patients, and 31 relatives. Twenty percent of HCM cases and eight% of DCM cases had pathogenic variants for cardiomyopathy, while 31% of HCM cases and 22% of DCM cases remain unexplained by having likely benign variant(s) or having no variants detected in any of the analysed genes.We replicated at least four recurrent variants among cardiomyopathy patients, which could be founder disease mutations in the Arabic population, including a frameshift variant (c.1371_1381dupTATCCAGTTAT) in the FKTN gene which seems to cause DCM in homozygosity and HCM, in heterozygosity.

Frontal fibrosing alopecia (FFA) is a recently described inflammatory and scarring type of hair loss affecting almost exclusively women. Despite a dramatic recent increase in incidence the aetiopathogenesis of FFA remains unknown. We undertake genome-wide association studies in females from a UK cohort, comprising 844 cases and 3,760 controls, a Spanish cohort of 172 cases and 385 controls, and perform statistical meta-analysis. We observe genome-wide significant association with FFA at four genomic loci: 2p22.2, 6p21.1, 8q24.22 and 15q2.1. Within the 6p21.1 locus, fine-mapping indicates that the association is driven by the HLA-B*07:02 allele. At 2p22.1, we implicate a putative causal missense variant in CYP1B1, encoding the homonymous xenobiotic- and hormone-processing enzyme. Transcriptomic analysis of affected scalp tissue highlights overrepresentation of transcripts encoding components of innate and adaptive immune response pathways. These findings provide insight into disease pathogenesis and characterise FFA as a genetically predisposed immuno-inflammatory disorder driven by HLA-B*07:02., Frontal fibrosing alopecia (FFA) features lichenoid cutaneous inflammation and scarring hair loss. Here, Tziotzios et al. identify four genetic loci associated with FFA by GWAS followed by Bayesian fine-mapping, co-localisation and HLA imputation which highlights HLA-B*07:02 as a risk factor.

Rationale: Ischemic stroke is among the leading causes of adult disability. Part of the variability in functional outcome after stroke has been attributed to genetic factors but no locus has been consistently associated with stroke outcome. Objective: Our aim was to identify genetic loci influencing the recovery process using accurate phenotyping to produce the largest GWAS (genome-wide association study) in ischemic stroke recovery to date. Methods and Results: A 12-cohort, 2-phase (discovery-replication and joint) meta-analysis of GWAS included anterior-territory and previously independent ischemic stroke cases. Functional outcome was recorded using 3-month modified Rankin Scale. Analyses were adjusted for confounders such as discharge National Institutes of Health Stroke Scale. A gene-based burden test was performed. The discovery phase (n=1225) was followed by open (n=2482) and stringent joint-analyses (n=1791). Those cohorts with modified Rankin Scale recorded at time points other than 3-month or incomplete data on previous functional status were excluded in the stringent analyses. Novel variants in PATJ (Pals1-associated tight junction) gene were associated with worse functional outcome at 3-month after stroke. The top variant was rs76221407 (G allele, β=0.40, P =1.70×10 −9 ). Conclusions: Our results identify a set of common variants in PATJ gene associated with 3-month functional outcome at genome-wide significance level. Future studies should examine the role of PATJ in stroke recovery and consider stringent phenotyping to enrich the information captured to unveil additional stroke outcome loci.

SummaryEnvironmental exposures during early life play a critical role in life-course health, yet the molecular phenotypes underlying environmental effects on health are poorly understood. In the Human Early Life Exposome (HELIX) project, a multi-centre cohort of 1,301 mother-child pairs, we associated individual exposomes consisting of >100 chemical, physical and lifestyle exposures assessed in pregnancy and childhood, with multi-omics profiles (methylome, transcriptome, metabolome and proteins) in childhood. We identified 1,170 associations, 249 in pregnancy and 921 in childhood, which revealed potential biological responses and sources of exposure. The methylome best captures the persistent influence of pregnancy exposures, including maternal smoking; while childhood exposures were associated with features from all omics layers, revealing novel signatures for indoor air quality, essential trace elements, endocrine disruptors and weather conditions. This study provides a unique resource (https://helixomics.isglobal.org/) to guide future investigation on the biological effects of the early life exposome.

Eating disorders and substance use disorders frequently co-occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa (BN) and problem alcohol use (genetic correlation [rg], twin-based=0.23-0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome-wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge-eating, AN without binge-eating, and a BN factor score), and eight substance-use-related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder (MDD). Total sample sizes per phenotype ranged from ~2,400 to ~537,000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism-based genetic correlations between eating disorder and substance-use-related phenotypes. Significant positive genetic associations emerged between AUD and AN (rg=0.18; false discovery rate q=0.0006), cannabis initiation and AN (rg=0.23; qwith binge-eating (rg=0.27; q=0.0016). Conversely, significant negative genetic correlations were observed between three non-diagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge-eating (rgs=-0.19 to −0.23; qs

Background: The adverse health effects of early life exposure to tobacco smoking have been widely reported. In spite of this, the underlying molecular mechanisms of in utero and postnatal exposure to tobacco smoke are only partially understood. Here, we aimed to identify multi-layer molecular signatures associated with exposure to tobacco smoke in these two exposure windows. Methods: We investigated the associations of maternal smoking during pregnancy and childhood secondhand smoke (SHS) exposure with molecular features measured in 1203 European children (mean age 8.1 years) from the Human Early Life Exposome (HELIX) project. Molecular features, covering 4 layers, included blood DNA methylation and gene and miRNA transcription, plasma proteins, and sera and urinary metabolites. Results: Maternal smoking during pregnancy was associated with DNA methylation changes at 18 loci in child blood. DNA methylation at 5 of these loci was related to expression of the nearby genes. However, the expression of these genes themselves was only weakly associated with maternal smoking. Conversely, childhood SHS was not associated with blood DNA methylation or transcription patterns, but with reduced levels of several serum metabolites and with increased plasma PAI1 (plasminogen activator inhibitor-1), a protein that inhibits fibrinolysis. Some of the in utero and childhood smoking-related molecular marks showed dose-response trends, with stronger effects with higher dose or longer duration of the exposure. Conclusion: In this first study covering multi-layer molecular features, pregnancy and childhood exposure to tobacco smoke were associated with distinct molecular phenotypes in children. The persistent and dose-dependent changes in the methylome make CpGs good candidates to develop biomarkers of past exposure. Moreover, compared to methylation, the weak association of maternal smoking in pregnancy with gene expression suggests different reversal rates and a methylation-based memory to past exposures. Finally, certain metabolites and protein markers evidenced potential early biological effects of postnatal SHS, such as fibrinolysis. The study has received funding from the European Community’s Seventh Framework Programme (FP7/2007-206) under grant agreement no. 308333 (the HELIX project); and the H2020-EU.3.1.2. - Preventing Disease Programme under grant agreement no 874583 (the ATHLETE project). BiB received core infrastructure funding from the Wellcome Trust (WT101597MA) and a joint grant from the UK Medical Research Council (MRC) and Economic and Social Science Research Council (ESRC) (MR/N024397/1). INMA data collections were supported by grants from the Instituto de Salud Carlos III, CIBERESP, and the Generalitat de Catalunya-CIRIT. KANC was funded by the grant of the Lithuanian Agency for Science Innovation and Technology (6-04-2014_31V-66). The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. The RHEA project was financially supported by European projects and the Greek Ministry of Health (Program of Prevention of Obesity and Neurodevelopmental Disorders in Preschool Children, in Heraklion district, Crete, Greece: 2011–2014; “Rhea Plus”: Primary Prevention Program of Environmental Risk Factors for Reproductive Health, and Child Health: 2012–2015). The work was also supported by MICINN (MTM2015-68140-R) and Centro Nacional de Genotipado-CEGEN-PRB2-ISCIII. The CRG/UPF Proteomics Unit is part of the Spanish Infrastructure for Omics Technologies (ICTS OmicsTech), and it is a member of the ProteoRed PRB3 consortium which is supported by grant PT17/0019 of the PE I+D+i 2013-2016 from the Instituto de Salud Carlos III (ISCIII) and ERDF. We acknowledge support from the Spanish Ministry of Science and Innovation through the “Centro de Excelencia Severo Ochoa 2019-2023” Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program.

Autoinflammatory diseases (AIDs) were first described as clinical disorders characterized by recurrent episodes of seemingly unprovoked sterile inflammation. In the past few years, the identification of novel AIDs expanded their phenotypes toward more complex clinical pictures associating vasculopathy, autoimmunity, or immunodeficiency. Herein, we describe two unrelated patients suffering since the neonatal period from a complex disease mainly characterized by severe sterile inflammation, recurrent bacterial infections, and marked humoral immunodeficiency. Whole-exome sequencing detected a novel, de novo heterozygous PLCG2 variant in each patient (p.Ala708Pro and p.Leu845_Leu848del). A clear enhanced PLCγ2 activity for both variants was demonstrated by both ex vivo calcium responses of the patient's B cells to IgM stimulation and in vitro assessment of PLC activity. These data supported the autoinflammation and PLCγ2-associated antibody deficiency and immune dysregulation (APLAID) diagnosis in both patients. Immunological evaluation revealed a severe decrease of immunoglobulins and B cells, especially class-switched memory B cells, with normal T and NK cell counts. Analysis of bone marrow of one patient revealed a reduced immature B cell fraction compared with controls. Additional investigations showed that both PLCG2 variants activate the NLRP3-inflammasome through the alternative pathway instead of the canonical pathway. Collectively, the evidences here shown expand APLAID diversity toward more severe phenotypes than previously reported including dominantly inherited agammaglobulinemia, add novel data about its genetic basis, and implicate the alternative NLRP3-inflammasome activation pathway in the basis of sterile inflammation. This work has been partially funded by the following: CERCA Programme/Generalitat de Catalunya (JIA), SAF2015-68472-C2-1-R grant from the Spanish Ministry of Economy and Competitiveness co-financed by European Regional Development Fund (ERDF) (JIA), RTI2018-096824-B-C21 grant from the Spanish Ministry of Science, Innovation and Universities co-financed by ERDF (JIA), AC15/00027 grant from the Instituto de Salud Carlos III/Transnational Research Projects on Rare Diseases (JIA), PI14/00405 grant from the Instituto de Salud Carlos III co-financed by ERDF (RC), PI13/00174 grant from the Spanish Ministry of Economy and Competitiveness co-financed by ERDF (PP), SAF2017-88276-R grant from the Spanish Ministry of Economy, Industry and Competitiveness co-financed by ERDF (PP), ERC-2013-CoG project 614578 from the European Research Council (PP), 20859/PI/18 grant from Fundación Séneca (PP), SAF2015-68472-C2-2-R grant from the Spanish Ministry of Economy and Competitiveness co-financed by ERDF (FC) and RTI2018-096824-B-C22 grant from the Spanish Ministry of Science, Innovation and Universities co-financed by ERDF (FC). MK acknowledges support from CRUK (A16567) and MRC (P028160). H M-B is a Rio Hortega fellowship from Instituto de Salud Carlos III (CM14/00008) and DB acknowledges support from the UCL Impact Studentship.

The molecular basis of many health outcomes attributed to solar ultraviolet radiation (UVR) is unknown. We tested the hypothesis that they may originate from transcriptional changes in blood cells. This was determined by assessing the effect of fluorescent solar simulated radiation (FSSR) on the transcriptional profile of peripheral blood pre- and 6 h, 24 h and 48 h post-exposure in nine healthy volunteers. Expression of 20 genes was down-regulated and one was up-regulated at 6 h after FSSR. All recovered to baseline expression at 24 h or 48 h. These genes have been associated with immune regulation, cancer and blood pressure; health effects attributed to vitamin D via solar UVR exposure. Plasma 25-hydroxyvitamin D3 [25OHD3] levels increased over time after FSSR and were maximal at 48 h. The increase was more pronounced in participants with low basal 25OHD3 levels. Mediation analyses suggested that changes in gene expression due to FSSR were independent of 25OHD3 and blood cell subpopulations.

Annals of the rheumatic diseases 76(5), 906-913 (2017). doi:10.1136/annrheumdis-2016-210324, Published by BMJ Publ. Group, London

Summary Background Terrestrial ultraviolet (UV) radiation causes erythema, oxidative stress, DNA mutations and skin cancer. Skin can adapt to these adverse effects by DNA repair, apoptosis, keratinization and tanning. Objectives To investigate the transcriptional response to fluorescent solar‐simulated radiation (FSSR) in sun‐sensitive human skin in vivo. Methods Seven healthy male volunteers were exposed to 0, 3 and 6 standard erythemal doses (SED). Skin biopsies were taken at 6 h and 24 h after exposure. Gene and microRNA expression were quantified with next generation sequencing. A set of candidate genes was validated by quantitative polymerase chain reaction (qPCR); and wavelength dependence was examined in other volunteers through microarrays. Results The number of differentially expressed genes increased with FSSR dose and decreased between 6 and 24 h. Six hours after 6 SED, 4071 genes were differentially expressed, but only 16 genes were affected at 24 h after 3 SED. Genes for apoptosis and keratinization were prominent at 6 h, whereas inflammation and immunoregulation genes were predominant at 24 h. Validation by qPCR confirmed the altered expression of nine genes detected under all conditions; genes related to DNA repair and apoptosis; immunity and inflammation; pigmentation; and vitamin D synthesis. In general, candidate genes also responded to UVA1 (340–400 nm) and/or UVB (300 nm), but with variations in wavelength dependence and peak expression time. Only four microRNAs were differentially expressed by FSSR. Conclusions The UV radiation doses of this acute study are readily achieved daily during holidays in the sun, suggesting that the skin transcriptional profile of ‘typical’ holiday makers is markedly deregulated. What's already known about this topic? The skin's transcriptional profile underpins its adverse (i.e. inflammation) and adaptive molecular, cellular and clinical responses (i.e. tanning, hyperkeratosis) to solar ultraviolet radiation.Few studies have assessed microRNA and gene expression in vivo in humans, and there is a lack of information on dose, time and waveband effects. What does this study add? Acute doses of fluorescent solar‐simulated radiation (FSSR), of similar magnitude to those received daily in holiday situations, markedly altered the skin's transcriptional profiles.The number of differentially expressed genes was FSSR‐dose‐dependent, reached a peak at 6 h and returned to baseline at 24 h.The initial transcriptional response involved apoptosis and keratinization, followed by inflammation and immune modulation. In these conditions, microRNA expression was less affected than gene expression., Linked Comment: Hart. Br J Dermatol 2020; 182:1328–1329. Plain language summary available online Respond to this article

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness 1, affecting 0.9–4% of women and 0.3% of men 2–4, with twin-based heritability estimates of 50–60% 5. Mortality rates are higher than those in other psychiatric disorders 6, and outcomes are unacceptably poor 7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI) 8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

Birth weight (BW) variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. These associations have been proposed to reflect the lifelong consequences of an adverse intrauterine environment. In earlier work, we demonstrated that much of the negative correlation between BW and adult cardio-metabolic traits could instead be attributable to shared genetic effects. However, that work and other previous studies did not systematically distinguish the direct effects of an individual’s own genotype on BW and subsequent disease risk from indirect effects of their mother’s correlated genotype, mediated by the intrauterine environment. Here, we describe expanded genome-wide association analyses of own BW (n=321,223) and offspring BW (n=230,069 mothers), which identified 278 independent association signals influencing BW (214 novel). We used structural equation modelling to decompose the contributions of direct fetal and indirect maternal genetic influences on BW, implicating fetal- and maternal-specific mechanisms. We used Mendelian randomization to explore the causal relationships between factors influencing BW through fetal or maternal routes, for example, glycemic traits and blood pressure. Direct fetal genotype effects dominate the shared genetic contribution to the association between lower BW and higher type 2 diabetes risk, whereas the relationship between lower BW and higher later blood pressure (BP) is driven by a combination of indirect maternal and direct fetal genetic effects: indirect effects of maternal BP-raising genotypes act to reduce offspring BW, but only direct fetal genotype effects (once inherited) increase the offspring’s later BP. Instrumental variable analysis using maternal BW-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring BP. In successfully separating fetal from maternal genetic effects, this work represents an important advance in genetic studies of perinatal outcomes, and shows that the association between lower BW and higher adult BP is attributable to genetic effects, and not to intrauterine programming.

Juvenile idiopathic arthritis (JIA) is a complex rheumatic disease with both autoimmune and autoinflammatory components. Recently, familial cases of systemic-onset JIA have been attributed to mutations in LACC1/FAMIN. We describe three affected siblings from a Moroccan consanguineous family with an early-onset chronic, symmetric and erosive arthritis previously diagnosed as rheumatoid factor (RF)-negative polyarticular JIA. Autozygosity mapping identified four homozygous regions shared by all patients, located in chromosomes 3, 6 (n:2) and 13, containing over 330 genes. Subsequent whole exome sequencing identified two potential candidate variants within these regions (in FARS2 and LACC1/FAMIN). Genotyping of a cohort of healthy Moroccan individuals (n: 352) and bioinformatics analyses finally supported the frameshift c.128_129delGT mutation in the LACC1/FAMIN gene, leading to a truncated protein (p.Cys43Tyrfs*6), as the most probable causative gene defect. Additional targeted sequencing studies performed in patients with systemic-onset JIA (n:23) and RF-negative polyarticular JIA (n: 44) revealed no pathogenic LACC1/FAMIN mutations. Our findings support the homozygous genotype in the LACC1/FAMIN gene as the defect underlying the family here described with a recessively inherited severe inflammatory joint disease. Our evidences provide further support to the involvement of LACC1/FAMIN deficiency in different types of JIA in addition to the initially described systemic-onset JIA., The authors would like to thank the CRG Genomics Unit for assistance with whole exome sequencing. This work has been partially funded by: CERCA Programme/Generalitat de Catalunya (JIA, XE, SO), the PERIS program of the Generalitat de Catalunya grant SLT002/16/00310 (RR), the Spanish Ministry of Economy and Competitiveness co-financed by European Regional Development Fund (ERDF) grant SAF2015-68472-C2-1-R (JIA), the Instituto de Salud Carlos III/Transnational Research Projects on Rare Diseases (JIA) grant AC15/00027, the Spanish Society of Pediatric Rheumatology (JIA), the Secretaria d’Universitats i Recerca del Departament d’Economia grant 2009-SGR-1502 (XE) and the European Union Seventh Framework Programme (FP7/2007-2013) grant agreement no. 262055 (XE). We also acknowledge support of the Spanish Ministry of Economy and Competitiveness (MEIC) to the EMBL partnership, ‘Centro de Excelencia Severo Ochoa’.

Aim: Gestational diabetes mellitus (GDM) has been linked with adverse long-term health outcomes for the fetus and mother. These effects may be mediated by epigenetic modifications. Materials & methods: Genome-wide RNA sequencing was performed in placental tissue and maternal blood in six GDM and six non-GDM pregnancies. Promoter region DNA methylation was examined for selected genes and correlated with gene expression to examine an epigenetic modulator mechanism. Results: Reductions of mRNA expression and increases in promoter methylation were observed for G6PD in GDM women, and for genes encoding IGF-binding proteins in GDM-exposed placenta. Conclusion: GDM involves epigenetic attenuation of G6PD, which may lead to hyperglycemia and oxidative stress, and the IGF-axis, which may modulate fetal macrosomia.

The circulating transcriptome is a valuable source of cancer biomarkers, which, with the exception of microRNAs (miRNAs), remains relatively unexplored. To elucidate which RNAs are present in plasma from melanoma patients and which could be used to distinguish cancer patients from healthy individuals, we used next generation sequencing (NGS), and validation was carried out by qPCR and/or ddPCR. We identified 442 different microRNAs in samples, eleven of which were differentially expressed (p &lt, 0.05). Levels of miR-134-5p and miR-320a-3p were significantly down-regulated (p &lt, 0.001) in melanoma samples (n = 96) compared to healthy controls (n = 28). Differentially expressed protein-encoding mRNA 5&prime, fragments were enriched for the angiopoietin, p21-activated kinase (PAK), and EIF2 pathways. Levels of ATM1, AMFR, SOS1, and CD109 gene fragments were up-regulated (p &lt, 0.001) in melanoma samples (n = 144) compared to healthy controls (n = 41) (AUC = 0.825). Over 40% of mapped reads were YRNAs, a class of non-coding RNAs that to date has been little explored. Expression levels of RNY3P1, RNY4P1, and RNY4P25 were significantly higher in patients with stage 0 disease than either healthy controls or more advanced stage disease (p &lt, 0.001). In conclusion, we have identified a number of novel RNA biomarkers, which, most importantly, we validated in multi-center retrospective and prospective cohorts, suggesting potential diagnostic use of these RNA species.

Huntington's disease (HD) is a polyglutamine disorder caused by a CAG expansion in the Huntingtin (HTT) gene exon 1. This expansion encodes a mutant protein whose abnormal function is traditionally associated with HD pathogenesis; however, recent evidence has also linked HD pathogenesis to RNA stable hairpins formed by the mutant HTT expansion. Here, we have shown that a locked nucleic acid-modified antisense oligonucleotide complementary to the CAG repeat (LNA-CTG) preferentially binds to mutant HTT without affecting HTT mRNA or protein levels. LNA-CTGs produced rapid and sustained improvement of motor deficits in an R6/2 mouse HD model that was paralleled by persistent binding of LNA-CTG to the expanded HTT exon 1 transgene. Motor improvement was accompanied by a pronounced recovery in the levels of several striatal neuronal markers severely impaired in R6/2 mice. Furthermore, in R6/2 mice, LNA-CTG blocked several pathogenic mechanisms caused by expanded CAG RNA, including small RNA toxicity and decreased Rn45s expression levels. These results suggest that LNA-CTGs promote neuroprotection by blocking the detrimental activity of CAG repeats within HTT mRNA. The present data emphasize the relevance of expanded CAG RNA to HD pathogenesis, indicate that inhibition of HTT expression is not required to reverse motor deficits, and further suggest a therapeutic potential for LNA-CTG in polyglutamine disorders. This work was supported by the Spanish government through the Plan Nacional de I+D+I and cofunded by grants from the Instituto de Salud Carlos III (ISCIII) – Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER) (project PI11/02036, to EM, and PI13/01250, to EPN); the Spanish Ministerio de Economía y Competitividad (MINECO) (SAF2008-00357 and SAF2013-49108-R, to XE, and SAF2014-60551-R: iRPaD, to EM); and the Generalitat de Catalunya, Departament Economia i Coneixement, Secretaria Universitats i Recerca (AGAUR 2014 SGR-1138, to XE).

Scope: Changes in genetic variations affecting the taste receptor, type 2, member 38 (TAS2R38) may identify the interacting mechanism leading to obesity and potential associations with proteins partaking in innate immunity, such as surfactant protein D (SPD) and mannan-binding lectin (MBL). Methods and results: We evaluated haplotypes of the bitter-taste receptor TAS2R38 in an identification sample of 210 women in different weight conditions, including anorexia nervosa and obesity. The association with SPD and MBL was tested in an independent sample picturing general population (n = 534). The relationship with obesity was validated in an extended final sample of 1319 participants. In the sample comprised of women in extreme weight conditions, increased obesity was identified in AVI/AVI subjects (OR = 2.5 [1.06–6.11], p = 0.035). In the sample picturing general population, increased SPD and MBL concentrations were found in nonsmoking AVI carriers. In this cohort, smoking and obesity blunted associations between TAS2R38 haplotypes and SPD and MBL. In the extended sample, the association of AVI/AVI haplotypes with increased obesity was also identified (OR = 1.4 [0.99/1.85], p = 0.049), being more robust in subjects aged

Next-generation sequencing (NGS) has the capacity of carrier screening in gamete donation (GD) programs. We have developed and validated an NGS carrier-screening test (qCarrier test) that includes 200 genes associated with 368 disorders (277 autosomal recessive and 37 X-linked). Carrier screening is performed on oocyte donation candidates and the male partner of oocyte recipient. Carriers of X-linked conditions are excluded from the GD program, whereas donors are chosen who do not carry mutations for the same gene/disease as the recipients. The validation phase showed a high sensitivity (>99% sensitivity) detecting all single-nucleotide variants, 13 indels, and 25 copy-number variants included in the validation set. A total of 1,301 individuals were analysed with the qCarrier test, including 483 candidate oocyte donors and 635 receptor couples, 105 females receiving sperm donation, and 39 couples seeking pregnancy. We identified 56% of individuals who are carriers for at least one genetic condition and 1.7% of female donors who were excluded from the program due to a carrier state of X-linked conditions. Globally, 3% of a priori assigned donations had a high reproductive risk that could be minimized after testing. Genetic counselling at different stages is essential for helping to facilitate a successful and healthy pregnancy.

In 2013, the American College of Medical Genetics (ACMG) examined the issue of incidental findings in whole exome and whole genome sequencing, and introduced recommendations to search for, evaluate and report medically actionable variants in a set of 56 genes. At a debate held during the 2014 European Society for Human Genetics Conference (ESHG) in Milan, Italy, the first author of that paper presented this view in a debate session that did not end with a conclusive vote from the mainly European audience for or against reporting back actionable incidental findings. In this meeting report, we elaborate on the discussions held during a special meeting hosted at the ESHG in 2013 from posing the question 'How to reach a (European) consensus on reporting incidental findings and unclassified variants in diagnostic next generation sequencing'. We ask whether an European consensus exists on the reporting of incidental findings in genome diagnostics, and present a series of key issues that require discussion at both a national and European level in order to develop recommendations for handling incidental findings and unclassified variants in line with the legal and cultural particularities of individual European member states.

Microcephaly is a rare condition in which the occipitofrontal circumference in a child is more than two standard deviations below the mean of children of the same age and gender. It is mainly caused by genetic abnormalities that interfere with the growth of the cerebral cortex during early months of fetal development. We present a case of a 12 years old patient with microcephaly. To identify a possible genetic origin of the phenotype, we performed array CGH and exome sequencing in the patient. Exome sequencing revealed the presence of a de novo missense mutation in the TUBB5 gene (E401K). Mutations in the TUBB5 are mainly responsible for microcephaly but the clinical spectrum is wide, from patients with severe developmental delay, and the presence of different brain malformations, to patients with only slightly cognitive impairment and normal motor development. Our patient shows a milder phenotype than other patients carrying the same mutation. These differences in the clinical features suggest that other factors, presumably genetic or epigenetic, could be modulating clinical expressivity of TUBB5. It is therefore evident that more functional studies are needed to understand the pathology that underlies the clinical spectrum of tubulin associated disease states.

Anorexia nervosa (AN) is a complex and serious eating disorder, occurring in ~1% of individuals. Despite having the highest mortality rate of any psychiatric disorder, little is known about the aetiology of AN, and few effective treatments exist.Global efforts to collect large sample sizes of individuals with AN have been highly successful, and a recent study consequently identified the first genome-wide significant locus involved in AN. This result, coupled with other recent studies and epidemiological evidence, suggest that previous characterizations of AN as a purely psychiatric disorder are over-simplified. Rather, both neurological and metabolic pathways may also be involved.In order to elucidate more of the system-specific aetiology of AN, we applied transcriptomic imputation methods to 3,495 cases and 10,982 controls, collected by the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED). Transcriptomic Imputation (TI) methods approaches use machine-learning methods to impute tissue-specific gene expression from large genotype data using curated eQTL reference panels. These offer an exciting opportunity to compare gene associations across neurological and metabolic tissues. Here, we applied CommonMind Consortium (CMC) and GTEx-derived gene expression prediction models for 13 brain tissues and 12 tissues with potential metabolic involvement (adipose, adrenal gland, 2 colon, 3 esophagus, liver, pancreas, small intestine, spleen, stomach).We identified 35 significant gene-tissue associations within the large chromosome 12 region described in the recent PGC-ED GWAS. We applied forward stepwise conditional analyses and FINEMAP to associations within this locus to identify putatively causal signals. We identified four independently associated genes;RPS26, C12orf49, SUOX, andRDH16.We also identified two further genome-wide significant gene-tissue associations, both in brain tissues;REEP5, in the dorso-lateral pre-frontal cortex (DLPFC; p=8.52×10−07), andCUL3, in the caudate basal ganglia (p=1.8×10−06). These genes are significantly enriched for associations with anthropometric phenotypes in the UK BioBank, as well as multiple psychiatric, addiction, and appetite/satiety pathways. Our results support a model of AN risk influenced by both metabolic and psychiatric factors.

International audience; Mitral valve prolapse (MVP) is a common cardiac valve disease that affects nearly 1 in 40 individuals. It can manifest as mitral regurgitation and is the leading indication for mitral valve surgery. Despite a clear heritable component, the genetic aetiology leading to non-syndromic MVP has remained elusive. Four affected individuals from a large multigenerational family segregating non-syndromic MVP underwent capture sequencing of the linked interval on chromosome 11. We report a missense mutation in the DCHS1 gene, the human homologue of the Drosophila cell polarity gene dachsous (ds), that segregates with MVP in the family. Morpholino knockdown of the zebrafish homologue dachsous1b resulted in a cardiac atrioventricular canal defect that could be rescued by wild-type human DCHS1, but not by DCHS1 messenger RNA with the familial mutation. Further genetic studies identified two additional families in which a second deleterious DCHS1 mutation segregates with MVP. Both DCHS1 mutations reduce protein stability as demonstrated in zebrafish, cultured cells and, notably, in mitral valve interstitial cells (MVICs) obtained during mitral valve repair surgery of a proband. Dchs1(+/-) mice had prolapse of thickened mitral leaflets, which could be traced back to developmental errors in valve morphogenesis. DCHS1 deficiency in MVP patient MVICs, as well as in Dchs1(+/-) mouse MVICs, result in altered migration and cellular patterning, supporting these processes as aetiological underpinnings for the disease. Understanding the role of DCHS1 in mitral valve development and MVP pathogenesis holds potential for therapeutic insights for this very common disease.

Many individuals with Parkinson's disease (PD) develop cognitive deficits, and a phenotypic and molecular overlap between neurodegenerative diseases exists. We investigated the contribution of rare variants in seven genes of known relevance to dementias (β-amyloid precursor protein (APP), PSEN1/2, MAPT (microtubule-associated protein tau), fused in sarcoma (FUS), granulin (GRN) and TAR DNA-binding protein 43 (TDP-43)) to PD and PD plus dementia (PD+D) in a discovery sample of 376 individuals with PD and followed by the genotyping of 25 out of the 27 identified variants with a minor allele frequency5% in 975 individuals with PD, 93 cases with Lewy body disease on neuropathological examination, 613 individuals with Alzheimer's disease (AD), 182 cases with frontotemporal dementia and 1014 general population controls. Variants identified in APP were functionally followed up by Aβ mass spectrometry in transiently transfected HEK293 cells. PD+D cases harbored more rare variants across all the seven genes than PD individuals without dementia, and rare variants in APP were more common in PD cases overall than in either the AD cases or controls. When additional controls from publically available databases were added, one rare variant in APP (c.1795GA(p.(E599K))) was significantly associated with the PD phenotype but was not found in either the PD cases or controls of an independent replication sample. One of the identified rare variants (c.2125GA (p.(G709S))) shifted the Aβ spectrum from Aβ40 to Aβ39 and Aβ37. Although the precise mechanism remains to be elucidated, our data suggest a possible role for APP in modifying the PD phenotype as well as a general contribution of genetic factors to the development of dementia in individuals with PD.

Tumors are composed of an evolving population of cells subjected to tissue-specific selection, which fuels tumor heterogeneity and ultimately complicates cancer driver gene identification. Here, we integrate cancer cell fraction, population recurrence, and functional impact of somatic mutations as signatures of selection into a Bayesian inference model for driver prediction. In an in-depth benchmark, we demonstrate that our model, cDriver, outperforms competing methods when analyzing solid tumors, hematological malignancies, and pan-cancer datasets. Applying cDriver to exome sequencing data of 21 cancer types from 6,870 individuals revealed 98 unreported tumor type-driver gene connections. These novel connections are highly enriched for chromatin-modifying proteins, hinting at a universal role of chromatin regulation in cancer etiology. Although infrequently mutated as single genes, we show that chromatin modifiers are altered in a large fraction of cancer patients. In summary, we demonstrate that integration of evolutionary signatures is key for identifying mutational driver genes, thereby facilitating the discovery of novel therapeutic targets for cancer treatment.

Tumors are composed of an evolving population of cells subjected to tissue-specific selection, which fuels tumor heterogeneity and ultimately complicates cancer driver gene identification. Here, we integrate cancer cell fraction, population recurrence, and functional impact of somatic mutations as signatures of selection into a Bayesian model for driver prediction. We demonstrate that our model, cDriver, outperforms competing methods when analyzing solid tumors, hematological malignancies, and pan-cancer datasets. Applying cDriver to exome sequencing data of 21 cancer types from 6,870 individuals revealed 98 unreported tumor type-driver gene connections. These novel connections are highly enriched for chromatin-modifying proteins, hinting at a universal role of chromatin regulation in cancer etiology. Although infrequently mutated as single genes, we show that chromatin modifiers are altered in a large fraction of cancer patients. In summary, we demonstrate that integration of evolutionary signatures is key for identifying mutational driver genes, thereby facilitating the discovery of novel therapeutic targets for cancer treatment. We acknowledge support of the Spanish Ministry of Economy and Competitiveness, 'Centro de Excelencia Severo Ochoa 2013-2017'. We acknowledge the support of the CERCA Programme/Generalitat de Catalunya. This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 635290. Luis Zapata has been supported by the International PhD scholarship program of La Caixa at CRG.

The term cystic fibrosis (CF)-like disease is used to describe patients with a borderline sweat test and suggestive CF clinical features but without two CFTR(cystic fibrosis transmembrane conductance regulator) mutations. We have performed the extensive molecular analysis of four candidate genes (SCNN1A, SCNN1B, SCNN1G and SERPINA1) in a cohort of 10 uncharacterized patients with CF and CF-like disease. We have used whole-exome sequencing to characterize mutations in the CFTR gene and these four candidate genes. CFTR molecular analysis allowed a complete characterization of three of four CF patients. Candidate variants in SCNN1A, SCNN1B, SCNN1G and SERPINA1 in six patients with CF-like phenotypes were confirmed by Sanger sequencing and were further supported by in silico predictive analysis, pedigree studies, sweat test in other family members, and analysis in CF patients and healthy subjects. Our results suggest that CF-like disease probably results from complex genotypes in several genes in an oligogenic form, with rare variants interacting with environmental factors.

Fragile X-associated tremor/ataxia syndrome (FXTAS) is a late-onset neurodegenerative disorder associated with FMR1 gene premutation alleles (55-200 CGG repeats). Fragile X-associated tremor/ataxia syndrome clinical core features include action tremor, gait ataxia, cognitive deficits progressing to dementia, and frequently parkinsonism. Although the pathogenic molecular mechanism of FXTAS is not completely understood, the restriction of the phenotype to the FMR1 premutation range has given rise to a model based on a RNA toxic gain-of-function. Since the identification of the first microRNAs (miRNAs) and their role in normal development, several studies have associated them with neurodegenerative diseases such as Parkinson, Alzheimer and Huntington diseases, suggesting that they play a key role in brain development, as well as in its morphogenesis. Herein, we present the characterization of miRNA expression profiles in FXTAS male patients using deep sequencing-based technologies and microarray technology. Deep sequencing analysis evidenced 83 miRNAs that were significantly deregulated whereas microarray analysis showed 31. When comparing these results, 14 miRNAs were found deregulated in FXTAS patients. MiR-424 and miR-574-3p showed significant fold change adjusted P-values in both platforms in FXTAS patients. MiR-424 has been founded substantially and specifically enriched in human cerebral cortical white matter of Alzheimer disease patients, which, together with cerebral atrophy, is a prominent imaging finding in individuals with FXTAS. The study provides the first systematic evidence of differential miRNA expression changes in FXTAS blood samples. Although further studies are necessary to better characterize the miRNA function in FXTAS disorder, our results suggest that they might contribute to its pathogenesis.

Background Here we have developed a novel and much more efficient strategy for the complete molecular characterisation of the cystic fibrosis (CF) transmembrane regulator (CFTR) gene, based on multiplexed targeted resequencing. We have tested this approach in a cohort of 92 samples with previously characterised CFTR mutations and polymorphisms. Methods After enrichment of the pooled barcoded DNA libraries with a custom NimbleGen SeqCap EZ Choice array (Roche) and sequencing with a HiSeq2000 (Illumina) sequencer, we applied several bioinformatics tools to call mutations and polymorphisms in CFTR. Results The combination of several bioinformatics tools allowed us to detect all known pathogenic variants (point mutations, short insertions/deletions, and large genomic rearrangements) and polymorphisms (including the poly-T and poly-thymidine-guanine polymorphic tracts) in the 92 samples. In addition, we report the precise characterisation of the breakpoints of seven genomic rearrangements in CFTR, including those of a novel deletion of exon 22 and a complex 85 kb inversion which includes two large deletions affecting exons 4–8 and 12–21, respectively. Conclusions This work is a proof-of-principle that targeted resequencing is an accurate and cost-effective approach for the genetic testing of CF and CFTR-related disorders (ie, male infertility) amenable to the routine clinical practice, and ready to substitute classical molecular methods in medical genetics.

Fibromyalgia is a chronic condition characterized by widespread pain, fatigue, non-restorative sleep and cognitive difficulties that affects 2–4% of the general population. Recently a possible relationship between the FMR1 premutation and fibromyalgia has been pointed out. In attempt to gather more data we screened for the FMR1 CGG expansion 700 DNA samples from unrelated fibromyalgia patients. This data might be useful for evaluating the incorporation of this test in rheumatologic procedures for women with fibromyalgia. The observed frequency of FMR1 premutation carriers (3 of 700, 0.4%) is not significantly different from the estimated rate in the general female population (1/250–1/400) (P = 0.539, P = 0.716). Clinical examination of the FMR1 premutation carriers identified revealed that all of them had important neurological symptoms with regard to muscular symptoms, neurocognitive alterations and neurovegetative impairments. With regard to other clinical aspects of the disease the cases apparently did not differ from the average fibromyalgia patients. On the basis of our results an FMR1 screening among fibromyalgia female patients would not be recommended. However it would be worthwhile to further evaluate the different clinical presentations that fibromyalgia patients might present based on their FMR1 premutation carrier status.

Genetic diagnostics of phenylketonuria (PKU) and tetrahydrobiopterin (BH4) deficient hyperphenylalaninemia (BH4DH) rely on methods that scan for known mutations or on laborious molecular tools that use Sanger sequencing. We have implemented a novel and much more efficient strategy based on high-throughput multiplex-targeted resequencing of four genes (PAH, GCH1, PTS, and QDPR) that, when affected by loss-of-function mutations, cause PKU and BH4DH. We have validated this approach in a cohort of 95 samples with the previously known PAH, GCH1, PTS, and QDPR mutations and one control sample. Pooled barcoded DNA libraries were enriched using a custom NimbleGen SeqCap EZ Choice array and sequenced using a HiSeq2000 sequencer. The combination of several robust bioinformatics tools allowed us to detect all known pathogenic mutations (point mutations, short insertions/deletions, and large genomic rearrangements) in the 95 samples, without detecting spurious calls in these genes in the control sample. We then used the same capture assay in a discovery cohort of 11 uncharacterized HPA patients using a MiSeq sequencer. In addition, we report the precise characterization of the breakpoints of four genomic rearrangements in PAH, including a novel deletion of 899 bp in intron 3. Our study is a proof-of-principle that high-throughput-targeted resequencing is ready to substitute classical molecular methods to perform differential genetic diagnosis of hyperphenylalaninemias, allowing the establishment of specifically tailored treatments a few days after birth.

Symptoms caused by bacterial, viral and malarial infections usually overlap and aetiologic diagnosis is difficult. Patient management in low-resource countries with limited laboratory services has been based predominantly on clinical evaluation and syndromic approaches. However, such clinical assessment has limited accuracy both for identifying the likely aetiological cause and for the early recognition of patients who will progress to serious or fatal disease. Plasma-detectable biomarkers that rapidly and accurately diagnose severe infectious diseases could reduce morbidity and decrease the unnecessary use of usually scarce therapeutic drugs. The discovery of microRNAs (miRNAs) has opened exciting new avenues to identify blood biomarkers of organ-specific injury. This review assesses current knowledge on the relationship between malaria disease and miRNAs, and evaluates how future research might lead to the use of these small molecules for identifying patients with severe malaria disease and facilitate treatment decisions. This work was supported by ISCIII-Subdirección General de Evaluación Plan Nacional I+D+I 2013-2016 (grant PI13/01478 cofunded by the Fondo Europeo de Desarrollo Regional [FEDER] and CES10/021-I3SNS to AM), the Departament d’Universitats i Recerca de la Generalitat de Catalunya (Agencia de Gestión de Ayudas Universitarias y de Investigación (AGAUR); grant 2014SGR263) and the CERCA Institutes Integration Program (SUMA 2013; Secretaría de Universidades e Investigación del Departamento de Economía y Conocimiento de la Generalidad de Cataluña, AGAUR and Fundación Institució dels Centres de Recerca de Catalunya). QB has a fellowship from the programme Miguel Servet of the ISCIII (Plan Nacional de I+D+I 2008-2011, grant number: CP11/00269). XE is supported by the Generalitat de Catalunya AGAUR 2014 SGR-1138 and the Spanish Ministry of Economy and Competitiveness, ‘Centro de Excelencia Severo Ochoa 2013-2017′, SEV-2012-0208.

(1) The objective of this study is to analyze differences in smell–taste capacity between females in extreme weight/eating conditions (EWC) and (2) to explore the interaction between smell/taste capacity, gastric hormones, eating behavior and body mass index (BMI). The sample comprised 239 females in EWC [64 Anorexia nervosa (AN) and 80 age-matched healthy-weight controls, and 59 obese and 36 age-matched healthy-weight controls]. Smell and taste assessments were performed through “Sniffin’ Sticks” and “Taste Strips,” respectively. The assessment measures included the eating disorders inventory-2, the symptom check list 90-revised, and The Dutch Eating Behavior Questionnaire, as well as peptides from the gastrointestinal tract [Ghrelin, peptide YY, cholecystokinin]. Smell capacity was differentially associated across EWC groups. Smell was clearly impaired in obese participants and increased in AN (hyposmia in Obesity was 54.3 and 6.4 % in AN), but taste capacity did not vary across EWC. Ghrelin levels were significantly decreased in obese subjects and were related to smell impairment. EWC individuals showed a distinct smell profile and circulating ghrelin levels compared to controls. Smell capacity and ghrelin may act as moderators of emotional eating and BMI This manuscript was supported by grants from Instituto de Salud Carlos III (Fondo de Investigación Sanitaria FIS PI11/210, FIS PI14/290 and CIBERobn) and Fondo Europeo de Desarrollo Regional (FEDER). CIBERObn is an initiative of ISCIII. Jose C. Fernández-García is recipient of a ‘Rio Hortega’ contract from ‘Instituto de Salud Carlos III,’ Madrid, Spain (CM12/00059). L. Garrido-Sánchez is supported by a fellowship from the Fondo de Investigación Sanitaria (FIS) “Miguel Servet I” MS13/00188-CP13/00188.

Common single-nucleotide polymorphisms (SNPs) account for a large proportion of the heritability of obsessive-compulsive disorder (OCD). Co-ocurrence of OCD and schizophrenia is commoner than expected based on their respective prevalences, complicating the clinical management of patients. This study addresses two main objectives: to identify particular genes associated with OCD by SNP-based and gene-based tests; and to test the existence of a polygenic risk shared with schizophrenia. The primary analysis was an exon-focused genome-wide association study of 370 OCD cases and 443 controls from Spain. A polygenic risk model based on the Psychiatric Genetics Consortium schizophrenia data set (PGC-SCZ2) was tested in our OCD data. A polygenic risk model based on our OCD data was tested on previous data of schizophrenia from our group. The most significant association at the gene-based test was found at DNM3 (P=7.9 × 10(-5)), a gene involved in synaptic vesicle endocytosis. The polygenic risk model from PGC-SCZ2 data was strongly associated with disease status in our OCD sample, reaching its most significant value after removal of the major histocompatibility complex region (lowest P=2.3 × 10(-6), explaining 3.7% of the variance). The shared polygenic risk was confirmed in our schizophrenia data. In conclusion, DNM3 may be involved in risk to OCD. The shared polygenic risk between schizophrenia and OCD may be partially responsible for the frequent comorbidity of both disorders, explaining epidemiological data on cross-disorder risk. This common etiology may have clinical implications. The study was supported by Fundación María José Jove, Fondo Europeo de Desarrollo Regional (FEDER), Xunta de Galicia; and by grants from the Instituto de Salud Carlos III FIS PI13/01136 to AC, CP11/00163 to JC, PI13/00918 to JMM, PI14/00413 to PA; the Generalitat de Catalunya AGAUR 2014 SGR-1138; the Spanish MINIECO SAF2013-49108- R Plan Estatal; and the European Commission 7th Framework Program, Project N. 262055 (ESGI) to XE. LD is supported by a Severo Ochoa fellowship of the Spanish MINIECO.

BACKGROUND: A brain-derived neurotrophic factor (BDNF) prodomain single-nucleotide polymorphism resulting in a valine to methionine substitution (Val66Met) has been associated with depression-related phenotypes and brain alterations involving regions consistently associated with major depressive disorder (MDD). The aim of our study was to evaluate the association of regional gray matter (GM) volume within the hippocampus and other unpredicted regions at the whole-brain level with the BDNF Val66Met polymorphism in MDD patients with melancholic features and their impact on treatment outcome. METHODS: A sample of 37 MDD inpatients was assessed with three-dimensional magnetic resonance imaging (1.5-T scanner). GM volume was analyzed with voxel-based morphometry (VBM) using Statistical Parametric Mapping (SPM5). The BDNF Val66Met variant was genotyped using SNPlex technology. MDD patients were classified according to genotype distribution under a dominant model of inheritance and thus comparing Val66 homozygotes (n = 22) versus Met66 carriers (n = 15). RESULTS: A significant GM volume reduction in the left hippocampus was observed in Met66 carriers. Conversely, in the same group, a volume increase in the right orbitofrontal cortex was detected. Moreover, a significant negative correlation between left hippocampal volume and days to remission was found in Val66 homozygotes, whereas right orbitofrontal volume was inversely correlated to days to remission in Met66 carriers. CONCLUSIONS: Our results suggest that the Val66Met BDNF variant may have a differential impact on the brain structure of melancholic patients with possible treatment outcome implications. Language: en

A nested case-control association study was designed to investigate the influence of maternal and fetal copy number variants (CNVs) on reproductive outcomes. Genotypes of ten CNVs encompassing GST and CYP genes were assessed. Significant associations were only found for child CNV genotypes. In particular, the child GSTM1 insertion allele was associated with prematurity protection (odds ratio, 95% CI: 0.67, 0.51-0.89; P < 0.01), whereas the child GSTT2B insertion allele was associated with an increased risk of being small for gestational age (odds ratio, 95% CI: 1.33, 1.07-1.67; P = 0.01). The study highlights the role of the fetal genome in prenatal development and also the need to analyse CNVs in a systematic manner.

Objective: The common single nucleotide polymorphism (SNP) rs9939609 in the fat mass and obesity-associated gene (FTO) is associated with obesity. As genetic variants associated with weight regulation might also be implicated in the etiology of eating disorders, we evaluated whether SNP rs9939609 is associated with bulimia nervosa (BN) and anorexia nervosa (AN). Methods: Association of rs9939609 with BN and AN was assessed in 689 patients with AN, 477 patients with BN, 984 healthy non-population-based controls, and 3,951 population-based controls (KORA-S4). Based on the familial and premorbid occurrence of obesity in patients with BN, we hypothesized an association of the obesity risk A-allele with BN. Results: In accordance with our hypothesis, we observed evidence for association of the rs9939609 A-allele with BN when compared to the non-population-based controls (unadjusted odds ratio (OR) = 1.142, one-sided 95% confidence interval (CI) 1.001–∞; one-sided p = 0.049) and a trend in the population-based controls (OR = 1.124, one-sided 95% CI 0.932–∞; one-sided p = 0.056). Interestingly, compared to both control groups, we further detected a nominal association of the rs9939609 A-allele to AN (OR = 1.181, 95% CI 1.027–1.359, two-sided p = 0.020 or OR = 1.673, 95% CI 1.101–2.541, two-sided p = 0.015,). Conclusion: Our data suggest that the obesity-predisposing FTO allele might be relevant in both AN and BN.

Epigenetic changes may play a role in the occurrence of asthma-related phenotypes.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.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.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.DNA methylation of ALOX12 may be an epigenetic biomarker for the risk of asthma-related phenotypes.

Arg/Arg homozygotes for the Gly16Arg polymorphism in the β₂-adrenoreceptor gene (ADRB2) have a reduced response to short-acting β₂-agonists but no effect has been associated with long-acting β₂-agonists (LABAs). We selected 604 subjects with current asthma from the European Community Respiratory Health Study to evaluate whether asthma control and lung function decline were associated with Gly16Arg polymorphism, and to test whether LABA or inhaled corticosteroid (ICS) use modified these effects. There was an increased risk of noncontrolled asthma (OR 1.33, 95% CI 1.01-1.75; p = 0.046) for each Arg allele. Among nonusers of ICS, the odds ratio of noncontrolled asthma among Arg/Arg versus Gly/Gly subjects was 2.73 (95% CI 1.28-5.82; p = 0.009). No increased risk of noncontrolled asthma associated with the Arg allele was observed among ICS and/or LABA users. For each Arg allele, a mean ± se decrease in decline in forced expiratory volume in 1 s of 7.7 ± 2.5 mL·yr⁻¹ was found (p-value for trend 0.003), irrespective of ICS or LABA use. Arg/Arg subjects had an increased risk of bronchial hyperresponsiveness (BHR) versus Gly/Gly subjects, with an odds ratio of 2.51 (95% CI 1.12-5.63; p = 0.025) if they did not use ICS. The Arg allele was associated with poorer asthma control, a steeper lung function decline and BHR. Absence of genotypic effects on asthma control among ICS users may be due to reversed β₂-adrenoreceptor desensitisation.

Disponible desde www.ncbi.nlm.nih.gov. -- Autores: Muiños-Gimeno, Margarita ... [et al.], [Background] The involvement of microRNAs (miRNAs) in neuronal differentiation and synaptic plasticity suggests a role for miRNAs in psychiatric disorders; association analyses and functional approaches were used to evaluate the implication of miRNAs in the susceptibility for panic disorder., [Methods] Case-control studies for 712 single-nucleotide polymorphisms (SNPs) tagging 325 human miRNA regions were performed in 203 Spanish patients with panic disorder and 341 control subjects. A sample of 321 anxiety patients and 642 control subjects from Finland and 102 panic disorder patients and 829 control subjects from Estonia was used as a replica. Reporter-gene assays and miRNA overexpression experiments in neuroblastoma cells were used to functionally evaluate the spectrum of genes regulated by the associated miRNAs., [Results] Two SNPs associated with panic disorder: rs6502892 tagging miR-22 (p < .0002), and rs11763020 tagging miR-339 (p < .00008). Other SNPs tagging miR-138-2, miR-488, miR-491, and miR-148a regions associated with different panic disorder phenotypes. Replication in the north-European sample supported several of these associations, although they did not pass correction for multiple testing. Functional studies revealed that miR-138-2, miR-148a, and miR-488 repress (30%–60%) several candidate genes for panic disorder—GABRA6, CCKBR and POMC, respectively—and that miR-22 regulates four other candidate genes: BDNF, HTR2C, MAOA, and RGS2. Transcriptome analysis of neuroblastoma cells transfected with miR-22 and miR-488 showed altered expression of a subset of predicted target genes for these miRNAs and of genes that might be affecting physiological pathways related to anxiety., [Conclusions] This work represents the first report of a possible implication of miRNAs in the etiology of panic disorder., Este trabajo fue financiado por el "Instituto Carlos III y el Fondo de Investigaciones Sanitarias" (nº CIBER-CB06/02/0058, CIBER-SAM, FIS/ISCIII: P1052565, ISCIII: GO3/184, FI05/0006 para Margarita Muiños-Gimeno), la "Fundación Marató-TV3" (nº 014 331), el "Departament d'Universitats Innovació i Empresa, Generalitat de Catalunya" (nº 2005SGR00008 y 2009SGR1435), el Ministerio de Ciencia e Innovación (nº SAF2008-00357), Sexto Programa Marco de la Unión Europea Proyecto Integrado SIROCCO (nº LSHG-CT-2006-037900). Andres Metspalu fue apoyado por el Séptimo Programa Marco Unión de la Europea (nº ENGAGE 201413, ECOGENE # 205419 EBC, y OPENGENE # 245536), la Unión Europea a través del Fondo Europeo de Desarrollo Regional (Centro de Excelencia en Genómica) y el Gobierno de Estonia (nº SF0180142s08). El Centro Español Nacional de Genotipado (CeGen) fue apoyado por "Genoma España". Academia de Finlandia, Yrjö and Tuulikki Ilvonen Foundation y la Fundación Sigrid Juselius. Fundación Científica de Estonia (nº 7034).

Objective: Prior studies investigating the association between APOE alleles epsilon 2/epsilon 4 and risk of intracerebral hemorrhage (ICH) have been inconsistent and limited to small sample sizes, and did not account for confounding by population stratification or determine which genetic risk model was best applied. Methods: We performed a large-scale genetic association study of 2189 ICH cases and 4041 controls from 7 cohorts, which were analyzed using additive models for epsilon 2 and epsilon 4. Results were subsequently meta-analyzed using a random effects model. A proportion of the individuals (322 cases, 357 controls) had available genome-wide data to adjust for population stratification. Results: Alleles epsilon 2 and epsilon 4 were associated with lobar ICH at genome-wide significance levels (odds ratio [OR] = 1.82, 95% confidence interval [CI] = 1.50-2.23, p = 6.6 x 10(-10); and OR = 2.20, 95%CI = 1.85-2.63, p = 2.4 x 10(-11), respectively). Restriction of analysis to definite/probable cerebral amyloid angiopathy ICH uncovered a stronger effect. Allele epsilon 4 was also associated with increased risk for deep ICH (OR = 1.21, 95% CI = 1.08-1.36, p = 2.6 x 10(-4)). Risk prediction evaluation identified the additive model as best for describing the effect of APOE genotypes. Interpretation: APOE epsilon 2 and epsilon 4 are independent risk factors for lobar ICH, consistent with their known associations with amyloid biology. In addition, we present preliminary findings on a novel association between APOE epsilon 4 and deep ICH. Finally, we demonstrate that an additive model for these APOE variants is superior to other forms of genetic risk modeling previously applied. ANN NEUROL 2010;68:934-943

Background Copy number variations (CNV) have become an important source of human genome variability noteworthy to consider when studying genetic susceptibility to complex diseases. As recent studies have found evidences for the potential involvement of CNVs in psychiatric disorders, we have studied the dosage effect of structural genome variants as a possible susceptibility factor for different psychiatric disorders in a candidate gene approach. Methods After selection of 68 psychiatric disorders’ candidate genes overlapping with CNVs, MLPA assays were designed to determine changes in copy number of these genes. The studied sample consisted of 724 patients with psychiatric disorders (accounting for anxiety disorders, mood disorders, eating disorders and schizophrenia) and 341 control individuals. Results CNVs were detected in 30 out of the 68 genes screened, indicating that a considerable proportion of neuronal pathways genes contain CNVs. When testing the overall burden of rare structural genomic variants in the different psychiatric disorders compared to control individuals, there was no statistically significant difference in the total amount of gains and losses. However, 14 out of the 30 changes were only found in psychiatric disorder patients but not in control individuals. These genes include GRM7, previously associated to major depression disorder and bipolar disorder, SLC6A13, in anxiety disorders, and S100B, SSTR5 and COMT in schizophrenia. Conclusions Although we have not been able to found a clear association between the studied CNVs and psychiatric disorders, the rare variants found only within the patients could account for a step further towards understanding the pathophysiology of psychiatric disorders.

Background The G protein-coupled receptor neuropeptide S receptor 1 (NPSR1) and its ligand neuropeptide S (NPS) form a signaling system mainly implicated in susceptibility to asthma and inflammatory disorders in humans and regulation of anxiety and arousal in rodents. We addressed here the role of NPS and NPSR1 as susceptibility genes for human anxiety disorders. Methods We performed comprehensive association analysis of genetic variants in NPS and NPSR1 in three independent study samples. We first studied a population-based sample (Health 2000, Finland) of 321 anxiety disorder patients and 1317 control subjects and subsequently a Spanish clinical panic disorder sample consisting of 188 cases and 315 control subjects. In addition, we examined a birth cohort of 2020 children (Barn Allergi Miljo Stockholm Epidemiologi [BAMSE], Sweden). We then tested whether alleles of the most significantly associated single nucleotide polymorphisms alter DNA-protein complex formation in electrophoretic mobility shift assays. Finally, we compared acute stress responses on the gene expression level in wild-type and Npsr1 −/− mice. Results We confirmed previously observed epidemiological association between anxiety and asthma in two population-based cohorts. Single nucleotide polymorphisms within NPS and NPSR1 associated with panic disorder diagnosis in the Finnish and Spanish samples and with parent-reported anxiety/depression in the BAMSE sample. Moreover, some of the implicated single nucleotide polymorphisms potentially affect transcription factor binding. Expression of neurotrophin-3, a neurotrophic factor connected to stress and panic reaction, was significantly downregulated in brain regions of stressed Npsr1 −/− mice, whereas interleukin-1 beta, an active stress-related immunotransmitter, was upregulated. Conclusions Our results suggest that NPS-NPSR1 signaling is likely involved in anxiety.

Objective The risk of rheumatoid arthritis (RA) is increased in the offspring of individuals affected with various autoimmune disorders, including psoriasis. Recently, the deletion of 2 genes from the late cornified envelope (LCE) gene cluster, LCE3C and LCE3B, has been associated with psoriasis in several populations. The purpose of this study was to assess whether this polymorphic gene deletion could also be involved in susceptibility to RA. Methods We tested for association between the LCE3C_LCE3B copy number variant and a single-nucleotide polymorphism in strong linkage disequilibrium with this variant (rs4112788) and RA in 2 independent case–control data sets (197 and 400 samples from patients with RA, respectively, and 411 and 567 samples from control subjects, respectively), collected at 4 Spanish hospitals. All samples were directly typed for presence of the LCE3C_LCE3B deletion (LCE3C_LCE3B-del) by polymerase chain reaction, and association analysis was performed using the SNPassoc R package. Results An association of homozygosity for the LCE3C_LCE3B-del and rs4112788 C allele with the risk of RA was observed in the first data set and was replicated in an independent case–control set. A combined analysis showed an overall P value of 0.0012 (odds ratio [OR] 1.45, 95% confidence interval [95% CI] 1.16–1.81) for association of the LCE3C_LCE3B-del. When the analysis was stratified for serologic data, we observed association in anti–cyclic citrullinated peptide (anti-CCP)–positive patients (P = 0.012, OR 1.51 [95% CI 1.09–2.13]) but not in anti-CCP–negative patients. Conclusion We have identified an association between the LCE3C_LCE3B-del and RA, and we have verified a pleiotropic effect of a common genetic risk factor (LCE3C_LCE3B-del) for autoimmune diseases that is involved in both psoriasis and RA.

Objective: We have assessed whether CFTR gene has a major impact on chronic pancreatitis (CP) pathogenesis than that provided by the CFTR mutations. For this aim, we have evaluated clinical parameters, CFTR mutations, and 3 potential regulatory CFTR variants (coding single-nucleotide polymorphisms): c.1540A>G, c.2694T>G, and c.4521G>A. Methods: CFTR gene analysis was performed in a cohort of 136 CP patients and 93 controls from Spanish population using current scanning techniques (single-strand conformation polymorphism/heteroduplex, denaturing gradient gel electrophoresis, and denaturing high-performance liquid chromatography) and direct sequencing. Results: A higher frequency of CFTR mutations were observed in patients (39%) than in controls (15%; P A (P-corrected = 0.036). Furthermore, the multi-marker analysis revealed that the 1540A; 2694G; 4521A (AGA) haplotype was more prevalent in CP than controls (P-corrected = 0.042). Remarkably, this association was unrelated to CF-causing mutations (P = 0.006). Conclusions: Our results corroborate the higher susceptibility of CF carriers to CP and, furthermore, suggest that the AGA haplotype could contribute to an increased risk in the development of CP irrespective of other CF-causing mutations.

Disruptions in circadian rhythms have been described in mood disorders (MD), but the involvement of genetic variation in genes pertaining to the molecular circadian machinery in the susceptibility to MD has not been conclusively determined. We examined 209 single-nucleotide polymorphisms (SNPs) covering 19 circadian genes (ADCYAP1, ARNTL, ARNTL2, BHLHB2, BHLHB3, CLOCK, CRY1, CRY2, CSNK1E, DBP, NPAS2, NR1D1, PER1, PER2, PER3, RORA, TIMELESS, VIP, and VIPR2) in a sample of 534 MD patients (335 with unipolar major mood depression (MDD) and 199 with bipolar disorder (BD)) and 440 community-based screened controls. Nominally, statistically significant associations were found in 15 circadian genes. The gene-wide test, corrected for the number of SNPs analyzed in each gene, identified significant associations in CRY1 (rs2287161), NPAS2 (rs11123857), and VIPR2 (rs885861) genes with the combined MD sample. In the MDD subsample, the same SNPs in CRY1 and NPAS2 of the combined sample remained associated, whereas in the BD subsample CLOCK (rs10462028) and VIP (rs17083008) were specifically associated. The association with an SNP located 3′ near CRY1 gene in MDD remained statistically significant after permutation correction at experiment level (p=0.007). Significant additive effects were found between the SNPs that were statistically significant at the gene-wide level. We also found evidence of associations between two-marker haplotypes in CRY1 and NPAS2 genes and MD. Our data support the contribution of the circadian system to the genetic susceptibility to MD and suggest that different circadian genes may have specific effects on MD polarity.

Background Attention-deficit/hyperactivity disorder (ADHD) is a childhood-onset neuropsychiatric disease that persists into adulthood in at least 30% of patients. There is evidence suggesting that abnormal left-right brain asymmetries in ADHD patients may be involved in a variety of ADHD-related cognitive processes, including sustained attention, working memory, response inhibition and planning. Although mechanisms underlying cerebral lateralization are unknown, left-right cortical asymmetry has been associated with transcriptional asymmetry at embryonic stages and several genes differentially expressed between hemispheres have been identified. Methods We selected six functional candidate genes showing at least 1.9-fold differential expression between hemispheres ( BAIAP2 , DAPPER1 , LMO4 , NEUROD6 , ATP2B3 , and ID2 ) and performed a case-control association study in an initial Spanish sample of 587 ADHD patients (270 adults and 317 children) and 587 control subjects. Results The single- and multiple-marker analysis provided evidence for a contribution of BAIAP2 to adulthood ADHD ( p = .0026 and p = .0016, respectively). We thus tested BAIAP2 for replication in two independent adult samples from Germany (639 ADHD patients and 612 control subjects) and Norway (417 ADHD cases and 469 control subjects). While no significant results were observed in the Norwegian sample, we replicated the initial association between BAIAP2 and adulthood ADHD in the German population ( p = .0062). Conclusions Our results support the participation of BAIAP2 in the continuity of ADHD across life span, at least in some of the populations analyzed, and suggest that genetic factors potentially influencing abnormal cerebral lateralization may be involved in this disorder.

Background Several clinical and genetic studies have focused on the role of brain-derived neurotrophic factor (BDNF) in the pathophysiology of various mental disorders. Recent lines of evidence regarding the network hypothesis of treatment outcome point towards the involvement of BDNF variants in the pharmacologic response in mood disorders (MD). Furthermore, there is strong evidence of a role for the serotonergic system in the pathophysiology and treatment of OCD, and upregulation of BDNF has been observed with various classes of antidepressants, including selective serotonin reuptake inhibitors (SSRI). Thus, we hypothesized that the BDNF gene might also be associated with treatment outcome in OCD. Methods We performed a single-marker and haplotype association study of eight tag single nucleotide polymorphisms in the BDNF genomic region and related this to pharmacologic response in a sample of 131 OCD patients. Results We found an association for a haplotype containing two single nucleotide polymorphisms that have previously been reported to be associated with treatment outcome in MD (rs908867 and rs1491850). Conclusions Our results support the hypothesis that the BDNF gene is involved in the response to psychopharmacologic treatment even though these preliminary findings await replication in a follow-up sample.

Objective As a key player in modulating both human physiological and behavioural functions including anxiety, perception and in particular appetite, serotonin (5-hydroxytryptamine, 5-HT) is likely to be involved in the aetiology of eating disorders. Studies showing serotonin receptor type 3 (5-HT 3) receptors to mediate food intake depression (anorexic response) have triggered our interest in investigating the putative role of variants in the 5-HT 3 receptor genes, HTR3A and HTR3B, in the susceptibility to anorexia nervosa (AN) and bulimia nervosa (BN). Methods Two hundred and sixty-five patients with AN and 91 patients with BN as well as 191 healthy controls served as a pilot study group for mutational analysis by direct sequencing. Variants showing a significant association were subsequently genotyped in an independent Spanish cohort of 78 patients with AN and 119 patients with BN as well as 331 healthy controls for replication purposes. Results In the pilot study, we found the coding HTR3B variant, p.Y129S, (rs1176744, P= 0.004, odds ratio = 2.06) to be associated with the restrictive subtype of AN. The association was confirmed in the Spanish study group (P=0.034 , odds ratio=2.26). Conclusion Our study provides first evidence for an involvement of 5-HT 3 variants in the aetiopathology of eating disorders in humans.

A fundamental difficulty in human genetics research is the identification of the spectrum of genetic variants that contribute to the susceptibility to common/complex disorders. We tested here the hypothesis that functional genetic variants may confer susceptibility to several related common disorders. We analyzed five main psychiatric diagnostic categories (substance-abuse, anxiety, eating, psychotic, and mood disorders) and two different control groups, representing a total of 3,214 samples, for 748 promoter and non-synonymous single nucleotide polymorphisms (SNPs) at 306 genes involved in neurotransmission and/or neurodevelopment. We identified strong associations to individual disorders, such as growth hormone releasing hormone (GHRH) with anxiety disorders, prolactin regulatory element (PREB) with eating disorders, ionotropic kainate glutamate receptor 5 (GRIK5) with bipolar disorder and several SNPs associated to several disorders, that may represent individual and related disease susceptibility factors. Remarkably, a functional SNP, rs945032, located in the promoter region of the bradykinin receptor B2 gene (BDKRB2) was associated to three disorders (panic disorder, substance abuse, and bipolar disorder), and two additional BDKRB2 SNPs to obsessive-compulsive disorder and major depression, providing evidence for common variants of susceptibility to several related psychiatric disorders. The association of BDKRB2 (odd ratios between 1.65 and 3.06) to several psychiatric disorders supports the view that a common genetic variant could confer susceptibility to clinically related phenotypes, and defines a new functional hint in the pathophysiology of psychiatric diseases.

Genetic and functional data indicate that variation in the expression of the neurotrophin-3 receptor gene (NTRK3) may have an impact on neuronal plasticity, suggesting a role for NTRK3 in the pathophysiology of anxiety disorders. MicroRNA (miRNA) posttranscriptional gene regulators act by base-pairing to specific sequence sites, usually at the 3'UTR of the target mRNA. Variants at these sites might result in gene expression changes contributing to disease susceptibility. We investigated genetic variation in two different isoforms of NTRK3 as candidate susceptibility factors for anxiety by resequencing their 3'UTRs in patients with panic disorder (PD), obsessive-compulsive disorder (OCD), and in controls. We have found the C allele of rs28521337, located in a functional target site for miR-485-3p in the truncated isoform of NTRK3, to be significantly associated with the hoarding phenotype of OCD. We have also identified two new rare variants in the 3'UTR of NTRK3, ss102661458 and ss102661460, each present only in one chromosome of a patient with PD. The ss102661458 variant is located in a functional target site for miR-765, and the ss102661460 in functional target sites for two miRNAs, miR-509 and miR-128, the latter being a brain-enriched miRNA involved in neuronal differentiation and synaptic processing. Interestingly, these two variants significantly alter the miRNA-mediated regulation of NTRK3, resulting in recovery of gene expression. These data implicate miRNAs as key posttranscriptional regulators of NTRK3 and provide a framework for allele-specific miRNA regulation of NTRK3 in anxiety disorders. This work was supported by the “Instituto Carlos III and Fondo de Investigaciones Sanitarias” [CIBER-CB06/02/0058, CIBER-SAM, FIS/ISCIII:P1052565, ISCIII:GO3/184], the “Fundació la Marató-TV3” [014331], the “Departament d’Universitats Innovació i Empresa, Generalitat de Catalunya” [2005SGR00008] and the European Union Sixth Framework Programme Integrated Project SIROCCO [Grant LSHG-CT-2006-037900]. The Spanish National Genotyping Center (CeGen) is supported by “Genoma España” and the “Ministerio de Educación y Ciencia” (Spanish Government). M.M is a recipient of a FIS fellowship [FI05/0006]. Y.E was supported by the “Ramón y Cajal” Program [Spanish Ministry of Science and Education]. We would like to thank J.M. Mercader for his help in technical issues and suggestions as well as B. Cormand, A. Macaya, R. Corominas and E. Cuenca (Hospital Universitari Vall d'Hebron, Barcelona) for kindly proving us with control samples.