Identification of potential common genetic modifiers of neurofibromas: a genome-wide association study in 1333 patients with neurofibromatosis type 1.

BACKGROUND: Neurofibromatosis type 1 (NF1) is characterized by the highly variable and unpredictable development of benign peripheral nerve sheath tumours: cutaneous (cNFs), subcutaneous (scNFs) and plexiform (pNFs) neurofibromas. OBJECTIVES: To identify neurofibroma modifier genes, in order to develop a database of patients with NF1. METHODS: All patients were phenotypically evaluated by a medical practitioner using a standardized questionnaire and the causal NF1 variant identified. We enrolled 1333 patients with NF1 who were genotyped for > 7 million common variants. RESULTS: A genome-wide association case-only study identified a significant association with 9q21.33 in the pNF phenotype in the discovery cohort. Twelve, three and four regions suggestive of association at the P ≤ 1 × 10-6 threshold were identified for pNFs, cNFs and scNFs, respectively. Evidence of replication was observed for 4, 2 and 6 loci, including 168 candidate modifier protein-coding genes. Among the candidate modifier genes, some were implicated in the RAS-mitogen-activated protein kinase pathway, cell-cycle control and myelination. Using an original CRISPR/Cas9-based functional assay, we confirmed GAS1 and SPRED2 as pNF and scNF candidate modifiers, as their inactivation specifically affected NF1-mutant Schwann cell growth. CONCLUSIONS: Our study may shed new light on the pathogenesis of NF1-associated neurofibromas and will, hopefully, contribute to the development of personalized care for patients with this deleterious and life-threatening condition.

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers.

BACKGROUND: Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. OBJECTIVE: We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. METHODS: We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses. RESULTS: We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). CONCLUSIONS: We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Genome-wide association analysis in dilated cardiomyopathy reveals two new players in systolic heart failure on chromosomes 3p25.1 and 22q11.23.

AIMS: Our objective was to better understand the genetic bases of dilated cardiomyopathy (DCM), a leading cause of systolic heart failure. METHODS AND RESULTS: We conducted the largest genome-wide association study performed so far in DCM, with 2719 cases and 4440 controls in the discovery population. We identified and replicated two new DCM-associated loci on chromosome 3p25.1 [lead single-nucleotide polymorphism (SNP) rs62232870, P = 8.7 × 10-11 and 7.7 × 10-4 in the discovery and replication steps, respectively] and chromosome 22q11.23 (lead SNP rs7284877, P = 3.3 × 10-8 and 1.4 × 10-3 in the discovery and replication steps, respectively), while confirming two previously identified DCM loci on chromosomes 10 and 1, BAG3 and HSPB7. A genetic risk score constructed from the number of risk alleles at these four DCM loci revealed a 3-fold increased risk of DCM for individuals with 8 risk alleles compared to individuals with 5 risk alleles (median of the referral population). In silico annotation and functional 4C-sequencing analyses on iPSC-derived cardiomyocytes identify SLC6A6 as the most likely DCM gene at the 3p25.1 locus. This gene encodes a taurine transporter whose involvement in myocardial dysfunction and DCM is supported by numerous observations in humans and animals. At the 22q11.23 locus, in silico and data mining annotations, and to a lesser extent functional analysis, strongly suggest SMARCB1 as the candidate culprit gene. CONCLUSION: This study provides a better understanding of the genetic architecture of DCM and sheds light on novel biological pathways underlying heart failure.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium.

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy.

BACKGROUND: Idiopathic membranous nephropathy is a major cause of the nephrotic syndrome in adults, but its etiologic basis is not fully understood. We investigated the genetic basis of biopsy-proven cases of idiopathic membranous nephropathy in a white population. METHODS: We performed independent genomewide association studies of single-nucleotide polymorphisms (SNPs) in patients with idiopathic membranous nephropathy from three populations of white ancestry (75 French, 146 Dutch, and 335 British patients). The patients were compared with racially matched control subjects; population stratification and quality controls were carried out according to standard criteria. Associations were calculated by means of a chi-square basic allele test; the threshold for significance was adjusted for multiple comparisons (with the Bonferroni method). RESULTS: In a joint analysis of data from the 556 patients studied (398 men), we identified significant alleles at two genomic loci associated with idiopathic membranous nephropathy. Chromosome 2q24 contains the gene encoding M-type phospholipase A(2) receptor (PLA(2)R1) (SNP rs4664308, P=8.6×10(-29)), previously shown to be the target of an autoimmune response. Chromosome 6p21 contains the gene encoding HLA complex class II HLA-DQ alpha chain 1 (HLA-DQA1) (SNP rs2187668, P=8.0×10(-93)). The association with HLA-DQA1 was significant in all three populations (P=1.8×10(-9), P=5.6×10(-27), and P=5.2×10(-36) in the French, Dutch, and British groups, respectively). The odds ratio for idiopathic membranous nephropathy with homozygosity for both risk alleles was 78.5 (95% confidence interval, 34.6 to 178.2). CONCLUSIONS: An HLA-DQA1 allele on chromosome 6p21 is most closely associated with idiopathic membranous nephropathy in persons of white ancestry. This allele may facilitate an autoimmune response against targets such as variants of PLA2R1. Our findings suggest a basis for understanding this disease and illuminate how adaptive immunity is regulated by HLA.

Risperidone response in patients with schizophrenia drives DNA methylation changes in immune and neuronal systems.

Background: The choice of efficient antipsychotic therapy for schizophrenia relies on a time-consuming trial-and-error approach, whereas the social and economic burdens of the disease call for faster alternatives. Material & methods: In a search for predictive biomarkers of antipsychotic response, blood methylomes of 28 patients were analyzed before and 4 weeks into risperidone therapy. Results: Several CpGs exhibiting response-specific temporal dynamics were identified in otherwise temporally stable methylomes and noticeable global response-related differences were observed between good and bad responders. These were associated with genes involved in immunity, neurotransmission and neuronal development. Polymorphisms in many of these genes were previously linked with schizophrenia etiology and antipsychotic response. Conclusion: Antipsychotic response seems to be shaped by both stable and medication-induced methylation differences.

The most common way to treat schizophrenia is antipsychotic medication. However, not all antipsychotics work for all patients. The only way to find a suitable antipsychotic is to prescribe one and wait, sometimes for months, to see if it works. Finding an alternative to this trial-and-error method would help reduce patient suffering and costs for healthcare systems. The idea is to look in the DNA of our blood cells for specific marks that can change in response to our lifestyle or health condition. These marks could help us predict how patients will react to the drug. In other words, they can serve as biomarkers of antipsychotic response. The current work examined the blood of schizophrenia patients before and 4 weeks after starting medication. The patients who did not respond well to the drug had different marks on the genes involved in immune defense and nervous system functioning. Some of these genes also play roles in the development of schizophrenia, whereas others can directly affect what happens to the drug in the patient's body. Although marks that predict how patients will react were not identified with certainty, valuable targets for future research were identified.

Layer myocardial strain is the most heritable echocardiographic trait.

AIMS: Myocardial deformation assessed by strain analysis represents a significant advancement in our assessment of cardiac mechanics. However, whether this variable is genetically heritable or whether all/most of its variability is related to environmental factors is currently unknown. We sought to determine the heritability of echocardiographically determined cardiac mechanics indices in a population setting. METHODS AND RESULTS: A total of 1357 initially healthy subjects (women 51.6%; 48.2 ± 14.1 years) were included in this study from 20-year follow-up after the fourth visit of the longitudinal familial STANISLAS cohort (Lorraine, France). Data were acquired using state-of-the-art cardiac ultrasound equipment, using acquisition and measurement protocols recommended by the EACVI (European Association of Cardiovascular Imaging)/ASE (American Society of Echocardiography)/Industry Task Force. Layer-specific global longitudinal strain (GLS) and global circumferential strain (full-wall, subendocardial, and subepicardial) and conventional structural and functional cardiac parameters and their potential heritability were assessed using restricted maximum likelihood analysis, with genetic relatedness matrix calculated from genome-wide association data. Indices of longitudinal/circumferential myocardial function and left ventricular (LV) ejection fraction had low heritability (ranging from 10% to 20%). Diastolic and standard LV function parameters had moderate heritability (ranging from 20% to 30%) except for end-systolic and end-diastolic volumes (30% and 45%, respectively). In contrast, global longitudinal subendocardial strain (GLSEndo)/global longitudinal subepicardial strain (GLSEpi) ratio had a high level of heritability (65%). Except for GLSEndo/GLSEpi ratio, a large percentage of variance remained unexplained (>50%). CONCLUSIONS: In our population cohort, GLSEndo/GLSEpi ratio had a high level of heritability, whereas other classical and mechanical LV function parameters did not. Given the increasing recognition of GLSEndo/GLSEpi ratio as an early/sensitive imaging biomarker of systolic dysfunction, our results suggest the possible existence of individual genetic predispositions to myocardial decline.

Association and performance of polygenic risk scores for breast cancer among French women presenting or not a familial predisposition to the disease.

BACKGROUND: Three partially overlapping breast cancer polygenic risk scores (PRS) comprising 77, 179 and 313 SNPs have been proposed for European-ancestry women by the Breast Cancer Association Consortium (BCAC) for improving risk prediction in the general population. However, the effect of these SNPs may vary from one country to another and within a country because of other factors. OBJECTIVE: To assess their associated risk and predictive performance in French women from (1) the CECILE population-based case-control study, (2) BRCA1 or BRCA2 (BRCA1/2) pathogenic variant (PV) carriers from the GEMO study, and (3) familial breast cancer cases with no BRCA1/2 PV and unrelated controls from the GENESIS study. RESULTS: All three PRS were associated with breast cancer in all studies, with odds ratios per standard deviation varying from 1.7 to 2.0 in CECILE and GENESIS, and hazard ratios varying from 1.1 to 1.4 in GEMO. The predictive performance of PRS313 in CECILE was similar to that reported in BCAC but lower than that in GENESIS (area under the receiver operating characteristic curve (AUC) = 0.67 and 0.75, respectively). PRS were less performant in BRCA2 and BRCA1 PV carriers (AUC = 0.58 and 0.54 respectively). CONCLUSION: Our results are in line with previous validation studies in the general population and in BRCA1/2 PV carriers. Additionally, we showed that PRS may be of clinical utility for women with a strong family history of breast cancer and no BRCA1/2 PV, and for those carrying a predicted PV in a moderate-risk gene like ATM, CHEK2 or PALB2.

Weak Association Between Genetic Markers of Hyperuricemia and Cardiorenal Outcomes: Insights From the STANISLAS Study Cohort With a 20-Year Follow-Up.

Background Hyperuricemia is associated with poor cardiovascular outcomes, although it is uncertain whether this relationship is causal in nature. This study aimed to: (1) assess the heritability of serum uric acid (SUA) levels, (2) conduct a genome-wide association study on SUA levels, and (3) investigate the association between certain single-nucleotide polymorphisms and target organ damage. Methods and Results The STANISLAS (Suivi Temporaire Annuel Non-Invasif de la Santé des Lorrains Assurés Sociaux) study cohort is a single-center longitudinal cohort recruited between 1993 and 1995 (visit 1), with a last visit (visit 4 [V4]) performed ≈20 years apart. Serum lipid profile, SUA, urinary albumin/creatinine ratio, estimated glomerular filtration rate, 24-hour ambulatory blood pressure monitoring, transthoracic echocardiography, pulse wave velocity, and genotyping for each participant were assessed at V4. A total of 1573 participants were included at V4, among whom 1417 had available SUA data at visit 1. Genome-wide association study results highlighted multiple single-nucleotide polymorphisms on the SLC2A9 gene linked to SUA levels. Carriers of the most associated mutated SLC2A9 allele (rs16890979) had significantly lower SUA levels. Although SUA level at V4 was highly associated with diabetes, prediabetes, higher body mass index, CRP (C-reactive protein) levels, estimated glomerular filtration rate variation (visit 1-V4), carotid intima-media thickness, and pulse wave velocity, rs16890979 was only associated with higher carotid intima-media thickness. Conclusions Our findings demonstrate that rs16890979, a genetic determinant of SUA levels located on the SLC2A9 gene, is associated with carotid intima-media thickness despite significant associations between SUA levels and several clinical outcomes, thereby lending support to the hypothesis of a link between SUA and cardiovascular disease.

Impact of natriuretic peptide polymorphisms on diastolic and metabolic function in a populational cohort: insights from the STANISLAS cohort.

AIMS: Elevated brain natriuretic peptide (BNP) and the N-terminal fragment of its pro-hormone (NT-proBNP) have become established biomarkers for heart failure and are associated with cardiovascular morbidity and mortality. Investigating sources of inter-individual heterogeneity, particularly genetic factors, could help better identify patients at risk of future cardiovascular disease. The aim of this study was to estimate the heritability of circulating NT-proBNP levels, to perform a genome-wide association study (GWAS) and gene-candidate analysis focused on NPPB-NPPA genes on these levels, and to examine their association with cardiovascular or metabolic outcomes. METHODS AND RESULTS: A total of 1555 individuals from the STANISLAS study were included. The heritability of circulating NT-proBNP levels was estimated at 15%, with seven single nucleotide polymorphisms (SNPs) reaching the significant threshold in the GWAS. All above SNPs were located on the same gene cluster constituted of MTHFR, CLCN6, NPPA, NPPB, and C1orf167. NPPA gene expression was also associated with NT-proBNP levels. Moreover, six other SNPs from NPPA-NPPB genes were associated with diastolic function (lateral e' on echocardiography) and metabolic features (glycated haemoglobin). CONCLUSIONS: The heritability of natriuretic peptides appears relatively low (15%) and mainly based on the same gene cluster constituted of MTHFR, CLCN6, NPPA, NPPB, and C1orf167. Natriuretic peptide polymorphisms are associated with natriuretic peptide levels and diastolic function. These results suggest that natriuretic peptide polymorphisms may have an impact in the early stages of cardiovascular and metabolic disease.

Fatty acid desaturase genetic variations and dietary omega-3 fatty acid intake associate with arterial stiffness.

Aims: Long-chain polyunsaturated fatty acids (PUFAs) generate diverse bioactive lipid mediators, which tightly regulate vascular inflammation. The effects of omega-3 PUFA supplementation in cardiovascular prevention however remain controversial. In addition to direct dietary intake, fatty acid desaturases (FADS) determine PUFA levels. Increased arterial stiffness represents an independent predictor of mortality and cardiovascular events. The aim of the present study was to determine the association of PUFA intake, FADS1 genotype, and FADS expression with arterial stiffness. Methods and results: A cross-sectional population-based cohort study of 1464 participants without overt cardiovascular disease was conducted. Dietary intake was assessed using a food frequency questionnaire. Arterial stiffness was assessed by carotid-femoral pulse wave velocity (cfPWV), and the FADS1 locus variant was determined. Blood cell transcriptomics was performed in a subset of 410 individuals. Pulse wave velocity was significantly associated with the FADS1 locus variant. Differential associations between PWV and omega-3 PUFA intake were observed depending on the FADS1 genotype. High omega-3 PUFA intake attenuated the FADS1 genotype-dependent associations. Carriers of the minor FADS1 locus variant exhibited increased expression of FADS2, which is associated with PWV. Conclusion: Taken together, these findings point to FADS1 genotype-dependent associations of omega-3 PUFA intake on subclinical cardiovascular disease. These findings may have implications for identifying responders and non-responders to omega-3 PUFA supplementation and open up for personalized dietary counselling in cardiovascular prevention.

Accelerated genome sequencing with controlled costs for infants in intensive care units: a feasibility study in a French hospital network.

Obtaining a rapid etiological diagnosis for infants with early-onset rare diseases remains a major challenge. These diseases often have a severe presentation and unknown prognosis, and the genetic causes are very heterogeneous. In a French hospital network, we assessed the feasibility of performing accelerated trio-genome sequencing (GS) with limited additional costs by integrating urgent requests into the routine workflow. In addition to evaluating our capacity for such an approach, this prospective multicentre pilot study was designed to identify pitfalls encountered during its implementation. Over 14 months, we included newborns and infants hospitalized in neonatal or paediatric intensive care units with probable genetic disease and in urgent need for etiological diagnosis to guide medical care. The duration of each step and the pitfalls were recorded. We analysed any deviation from the planned schedule and identified obstacles. Trio-GS was performed for 37 individuals, leading to a molecular diagnosis in 18/37 (49%), and 21/37 (57%) after reanalysis. Corrective measures and protocol adaptations resulted in a median duration of 42 days from blood sampling to report. Accelerated trio-GS is undeniably valuable for individuals in an urgent care context. Such a circuit should coexist with a rapid or ultra-rapid circuit, which, although more expensive, can be used in particularly urgent cases. The drop in GS costs should result in its generalized use for diagnostic purposes and lead to a reduction of the costs of rapid GS.

Heritability of a resting heart rate in a 20-year follow-up family cohort with GWAS data: Insights from the STANISLAS cohort.

BACKGROUND: The association between resting heart rate (HR) and cardiovascular outcomes, especially heart failure, is now well established. However, whether HR is mainly an integrated marker of risk associated with other features, or rather a genetic origin risk marker, is still a matter for debate. Previous studies reported a heritability ranging from 14% to 65%. DESIGN: We assessed HR heritability in the STANISLAS family-study, based on the data of four visits performed over a 20-year period, and adjusted for most known confounding effects. METHODS: These analyses were conducted using a linear mixed model, adjusted on age, sex, tea or coffee consumption, beta-blocker use, physical activity, tobacco use, and alcohol consumption to estimate the variance captured by additive genetic effects, via average information restricted maximum likelihood analysis, with both self-reported pedigree and genetic relatedness matrix (GRM) calculated from genome-wide association study data. RESULTS: Based on the data of all visits, the HR heritability (h2) estimate was 23.2% with GRM and 24.5% with pedigree. However, we found a large heterogeneity of HR heritability estimations when restricting the analysis to each of the four visits (h2 from 19% to 39% using pedigree, and from 14% to 32% using GRM). Moreover, only a little part of variance was explained by the common household effect (<5%), and half of the variance remained unexplained. CONCLUSION: Using a comprehensive analysis based on a family cohort, including the data of multiple visits and GRM, we found that HR variability is about 25% from genetic origin, 25% from repeated measures and 50% remains unexplained.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians.

Differentiated thyroid carcinoma (DTC) incidence is characterized by wide ethnic and geographic variations, with high incidence rates observed in Oceanian populations. Genome-wide association studies (GWAS) identified mainly four DTC susceptibility loci at 9q22.33, 14q13.3, 2q35 and 8p12. Here we performed fine-mapping of the 2q35 and 8p12 loci in the population of the EPITHYR consortium that includes Europeans, Melanesians and Polynesians to identify likely causal variants for DTC risk. We conducted a colocalization analysis using eQTLs data to determine the SNPs with the highest probability of causality. At 2q35, we highlighted rs16857609 located in DIRC3. This SNP has a high probability of causality in the three populations, and a significant association in Europeans (OR = 1.4, p = 1.9 x 10-10). It is also associated with expression of DIRC3 and of the nearby gene IGFBP5 in thyroid tumour cells. At 8p12, we identified rs7844425 which was significantly associated with DTC in Europeans (OR = 1.32, p = 7.6 x 10-8) and rs2439304, which was highlighted by the colocalization analysis but only moderately associated with DTC in our dataset (OR = 1.2, p = 0.001). These SNPs are linked to the expression of NRG1 in thyroid tissue. Hence, our study identified novel variants at 2q35 and 8p12 to be prioritized for further functional studies.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma.

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Gene expression and response prediction to amisulpride in the OPTiMiSE first episode psychoses.

A fundamental shortcoming in the current treatment of schizophrenia is the lack of valid criteria to predict who will respond to antipsychotic treatment. The identification of blood-based biological markers of the therapeutic response would enable clinicians to identify the subgroup of patients in whom conventional antipsychotic treatment is ineffective and offer alternative treatments. As part of the Optimisation of Treatment and Management of Schizophrenia in Europe (OPTiMiSE) programme, we conducted an RNA-Seq analysis on 188 subjects with first episode psychosis, all of whom were subsequently treated with amisulpride for 4 weeks. We compared gene expression on total RNA from patients' blood before and after treatment and identified 32 genes for which the expression changed after treatment in good responders only. These findings were replicated in an independent sample of 24 patients with first episode psychosis. Six genes showed a significant difference in expression level between good and poor responders before starting treatment, allowing to predict treatment outcome with a predictive value of 93.8% when combined with clinical features. Collectively, these findings identified new mechanisms to explain symptom improvement after amisulpride medication and highlight the potential of combining gene expression profiling with clinical data to predict treatment response in first episode psychoses.

PCSK9 Protein and rs562556 Polymorphism Are Associated With Arterial Plaques in Healthy Middle-Aged Population: The STANISLAS Cohort.

Background PCSK9 (Proprotein convertase subtilisin/kexin type 9) binds low-density lipoprotein receptor, preventing its recycling. PCSK9 is a risk predictor and a biotarget in atherosclerosis. The PCSK9-rs562556 variant has been reported as a gain-of-function mutation. The aim of this study was to determine whether the PCSK9-low-density lipoprotein receptor-rs562556 axis is associated with carotid artery plaques between 2 visits separated by almost 20 years in a longitudinal population cohort. Methods and Results The STANISLAS (Suivi Temporaire Annuel Non-Invasif de la Santé des Lorrains Assurés Sociaux) cohort is a longitudinal familial cohort from the Lorraine region of France. Participants attending 2 visits (visit 1 and visit 4) separated by 18.5 years (mean) were included (n=997). Carotid artery plaques were determined with standardized vascular echography. The mean age of the adult population at visit 1 was 42±5 years. At visit 4, 203 (20.4%) participants had arterial plaques. Participants who developed arterial plaques were older (42.7±5.4 versus 41.7±4.7 years), more often male (60% versus 49%), smokers (29% versus 18%), with diabetes mellitus (6% versus 3%), and higher cholesterol levels (low-density lipoprotein cholesterol, 1.6±0.4 versus 1.5±0.3 g/L) (all P<0.05). The independent factors associated with arterial plaques were age, smoking, and low-density lipoprotein cholesterol. Higher PCSK9 levels were associated with arterial plaques on top of the clinical model (odds ratio, 2.14; 95% CI,= 1.28-3.58); the missense mutation coding the single-nucleotide polymorphism rs562556 was associated with both higher PCSK9 concentration and incident carotid arterial plaques. Conclusions Higher PCSK9 concentration was associated with the development of arterial plaques almost 20 years in advance in a healthy middle-aged population. Mutations of the single-nucleotide polymorphism rs562556 associated with both PCSK9 levels and arterial plaques reinforce the potential causality of our findings. PCSK9 inhibitors could be useful for primary cardiovascular prevention.

Sex specific associations in genome wide association analysis of renal cell carcinoma.

Renal cell carcinoma (RCC) has an undisputed genetic component and a stable 2:1 male to female sex ratio in its incidence across populations, suggesting possible sexual dimorphism in its genetic susceptibility. We conducted the first sex-specific genome-wide association analysis of RCC for men (3227 cases, 4916 controls) and women (1992 cases, 3095 controls) of European ancestry from two RCC genome-wide scans and replicated the top findings using an additional series of men (2261 cases, 5852 controls) and women (1399 cases, 1575 controls) from two independent cohorts of European origin. Our study confirmed sex-specific associations for two known RCC risk loci at 14q24.2 (DPF3) and 2p21(EPAS1). We also identified two additional suggestive male-specific loci at 6q24.3 (SAMD5, male odds ratio (ORmale) = 0.83 [95% CI = 0.78-0.89], Pmale = 1.71 × 10-8 compared with female odds ratio (ORfemale) = 0.98 [95% CI = 0.90-1.07], Pfemale = 0.68) and 12q23.3 (intergenic, ORmale = 0.75 [95% CI = 0.68-0.83], Pmale = 1.59 × 10-8 compared with ORfemale = 0.93 [95% CI = 0.82-1.06], Pfemale = 0.21) that attained genome-wide significance in the joint meta-analysis. Herein, we provide evidence of sex-specific associations in RCC genetic susceptibility and advocate the necessity of larger genetic and genomic studies to unravel the endogenous causes of sex bias in sexually dimorphic traits and diseases like RCC.

Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples.

Next-generation sequencing (NGS) studies are becoming routinely used for the detection of novel and clinically actionable DNA variants at a pangenomic scale. Such analyses are now used in the clinical practice to enable precision medicine. Formalin-fixed paraffin-embedded (FFPE) tissues are still one of the most abundant source of cancer clinical specimen, unfortunately this method of preparation is known to degrade DNA and therefore compromise subsequent analysis. Some studies have reported that variant detection can be performed on FFPE samples sequenced with NGS techniques, but few or none have done an in-depth coverage analysis and compared the influence of different state-of-the-art FFPE DNA extraction kits on the quality of the variant calling. Here, we generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell™ RSC DNA FFPE Kit from Promega). We determined the rate of concordance of called variants between matched FF and FFPE samples on all common variants (representing at least 86% of the total number of variants for SNVs). The concordance rate is very high between all matched FF / FFPE pairs, with equivalent values for the three kits we analyzed. On the other hand, when looking at the difference between the total number of variants in FF and FFPE, we find a significant variation for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but statistically significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

New chimeric RNAs in acute myeloid leukemia.

Background: High-throughput next generation sequencing (NGS) technologies enable the detection of biomarkers used for tumor classification, disease monitoring and cancer therapy. Whole-transcriptome analysis using RNA-seq is important, not only as a means of understanding the mechanisms responsible for complex diseases but also to efficiently identify novel genes/exons, splice isoforms, RNA editing, allele-specific mutations, differential gene expression and fusion-transcripts or chimeric RNA (chRNA). Methods: We used Crac, a tool that uses genomic locations and local coverage to classify biological events and directly infer splice and chimeric junctions within a single read. Crac's algorithm extracts transcriptional chimeric events irrespective of annotation with a high sensitivity, and CracTools was used to aggregate, annotate and filter the chRNA reads. The selected chRNA candidates were validated by real time PCR and sequencing.  In order to check the tumor specific expression of chRNA, we analyzed a publicly available dataset using a new tag search approach. Results:  We present data related to acute myeloid leukemia (AML) RNA-seq analysis. We highlight novel biological cases of chRNA, in addition to previously well characterized leukemia chRNA. We have identified and validated 17 chRNAs among 3 AML patients: 10 from an AML patient with a translocation between chromosomes 15 and 17 (AML-t(15;17), 4  from patient with normal karyotype (AML-NK) 3 from a patient with chromosomal 16 inversion (AML-inv16). The new fusion transcripts can be classified into four groups according to the exon organization. Conclusions:  All groups suggest complex but distinct synthesis mechanisms involving either collinear exons of different genes, non-collinear exons, or exons of different chromosomes. Finally, we check tumor-specific expression in a larger RNA-seq AML cohort and identify new AML biomarkers that could improve diagnosis and prognosis of AML.