High-dimensional immune profiling identifies a biomarker to monitor dimethyl fumarate response in multiple sclerosis.

Dimethyl fumarate (DMF) is an immunomodulatory treatment for multiple sclerosis (MS). Despite its wide clinical use, the mechanisms underlying clinical response are not understood. This study aimed to reveal immune markers of therapeutic response to DMF treatment in MS. For this purpose, we prospectively collected peripheral blood mononuclear cells (PBMCs) from a highly characterized cohort of 44 individuals with MS before and at 12 and 48 wk of DMF treatment. Single cells were profiled using high-dimensional mass cytometry. To capture the heterogeneity of different immune subsets, we adopted a bioinformatic multipanel approach that allowed cell population-cluster assignment of more than 50 different parameters, including lineage and activation markers as well as chemokine receptors and cytokines. Data were further analyzed in a semiunbiased fashion implementing a supervised representation learning approach to capture subtle longitudinal immune changes characteristic for therapy response. With this approach, we identified a population of memory T helper cells expressing high levels of neuroinflammatory cytokines (granulocyte-macrophage colony-stimulating factor [GM-CSF], interferon γ [IFNγ]) as well as CXCR3, whose abundance correlated with treatment response. Using spectral flow cytometry, we confirmed these findings in a second cohort of patients. Serum neurofilament light-chain levels confirmed the correlation of this immune cell signature with axonal damage. The identified cell population is expanded in peripheral blood under natalizumab treatment, substantiating a specific role in treatment response. We propose that depletion of GM-CSF-, IFNγ-, and CXCR3-expressing T helper cells is the main mechanism of action of DMF and allows monitoring of treatment response.

Meta-analysis of ACE inhibitor-induced angioedema identifies novel risk locus.

BACKGROUND: Angioedema is a rare but potentially life-threatening adverse drug reaction in patients receiving angiotensin-converting enzyme inhibitors (ACEis). Research suggests that susceptibility to ACEi-induced angioedema (ACEi-AE) involves both genetic and nongenetic risk factors. Genome- and exome-wide studies of ACEi-AE have identified the first genetic risk loci. However, understanding of the underlying pathophysiology remains limited. OBJECTIVE: We sought to identify further genetic factors of ACEi-AE to eventually gain a deeper understanding of its pathophysiology. METHODS: By combining data from 8 cohorts, a genome-wide association study meta-analysis was performed in more than 1000 European patients with ACEi-AE. Secondary bioinformatic analyses were conducted to fine-map associated loci, identify relevant genes and pathways, and assess the genetic overlap between ACEi-AE and other traits. Finally, an exploratory cross-ancestry analysis was performed to assess shared genetic factors in European and African-American patients with ACEi-AE. RESULTS: Three genome-wide significant risk loci were identified. One of these, located on chromosome 20q11.22, has not been implicated previously in ACEi-AE. Integrative secondary analyses highlighted previously reported genes (BDKRB2 [bradykinin receptor B2] and F5 [coagulation factor 5]) as well as biologically plausible novel candidate genes (PROCR [protein C receptor] and EDEM2 [endoplasmic reticulum degradation enhancing alpha-mannosidase like protein 2]). Lead variants at the risk loci were found with similar effect sizes and directions in an African-American cohort. CONCLUSIONS: The present results contributed to a deeper understanding of the pathophysiology of ACEi-AE by (1) providing further evidence for the involvement of bradykinin signaling and coagulation pathways and (2) suggesting, for the first time, the involvement of the fibrinolysis pathway in this adverse drug reaction. An exploratory cross-ancestry comparison implicated the relevance of the associated risk loci across diverse ancestries.

Relationships between neurotransmitter receptor densities and expression levels of their corresponding genes in the human hippocampus.

Neurotransmitter receptors are key molecules in signal transmission, their alterations are associated with brain dysfunction. Relationships between receptors and their corresponding genes are poorly understood, especially in humans. We combined in vitro receptor autoradiography and RNA sequencing to quantify, in the same tissue samples (7 subjects), the densities of 14 receptors and expression levels of their corresponding 43 genes in the Cornu Ammonis (CA) and dentate gyrus (DG) of human hippocampus. Significant differences in receptor densities between both structures were found only for metabotropic receptors, whereas significant differences in RNA expression levels mostly pertained ionotropic receptors. Receptor fingerprints of CA and DG differ in shapes but have similar sizes; the opposite holds true for their "RNA fingerprints", which represent the expression levels of multiple genes in a single area. In addition, the correlation coefficients between receptor densities and corresponding gene expression levels vary widely and the mean correlation strength was weak-to-moderate. Our results suggest that receptor densities are not only controlled by corresponding RNA expression levels, but also by multiple regionally specific post-translational factors.

Integrative genomic analysis of pediatric T-cell lymphoblastic lymphoma reveals candidates of clinical significance.

T-cell lymphoblastic lymphoma (T-LBL) is a heterogeneous malignancy of lymphoblasts committed to T-cell lineage. The dismal outcomes (15%-30%) after T-LBL relapse warrant establishing risk-based treatment. To our knowledge, this study presents the first comprehensive, systematic, integrated, genome-wide analysis including relapsed cases that identifies molecular markers of prognostic relevance for T-LBL. NOTCH1 was identified as the putative driver for T-LBL. An activated NOTCH/PI3K-AKT signaling axis and alterations in cell cycle regulators constitute the core oncogenic program for T-LBL. Mutated KMT2D was identified as a prognostic marker. The cumulative incidence of relapse was 47% ± 17% in patients with KMT2D mutations, compared with 14% ± 3% in wild-type KMT2D. Structural analysis of the mutated domains of KMT2D revealed a plausible impact on structure and functional consequences. These findings provide new insights into the pathogenesis of T-LBL, including high translational potential. The ongoing LBL 2018 trial (www.clinicaltrials.gov #NCT04043494) allows for prospective validation and subsequent fine tuning of the stratification criteria for T-LBL risk groups to improve survival of pediatric patients.

Using polygenic scores and clinical data for bipolar disorder patient stratification and lithium response prediction: machine learning approach.

BACKGROUND: Response to lithium in patients with bipolar disorder is associated with clinical and transdiagnostic genetic factors. The predictive combination of these variables might help clinicians better predict which patients will respond to lithium treatment. AIMS: To use a combination of transdiagnostic genetic and clinical factors to predict lithium response in patients with bipolar disorder. METHOD: This study utilised genetic and clinical data (n = 1034) collected as part of the International Consortium on Lithium Genetics (ConLi+Gen) project. Polygenic risk scores (PRS) were computed for schizophrenia and major depressive disorder, and then combined with clinical variables using a cross-validated machine-learning regression approach. Unimodal, multimodal and genetically stratified models were trained and validated using ridge, elastic net and random forest regression on 692 patients with bipolar disorder from ten study sites using leave-site-out cross-validation. All models were then tested on an independent test set of 342 patients. The best performing models were then tested in a classification framework. RESULTS: The best performing linear model explained 5.1% (P = 0.0001) of variance in lithium response and was composed of clinical variables, PRS variables and interaction terms between them. The best performing non-linear model used only clinical variables and explained 8.1% (P = 0.0001) of variance in lithium response. A priori genomic stratification improved non-linear model performance to 13.7% (P = 0.0001) and improved the binary classification of lithium response. This model stratified patients based on their meta-polygenic loadings for major depressive disorder and schizophrenia and was then trained using clinical data. CONCLUSIONS: Using PRS to first stratify patients genetically and then train machine-learning models with clinical predictors led to large improvements in lithium response prediction. When used with other PRS and biological markers in the future this approach may help inform which patients are most likely to respond to lithium treatment.

Bipolar multiplex families have an increased burden of common risk variants for psychiatric disorders.

Multiplex families with a high prevalence of a psychiatric disorder are often examined to identify rare genetic variants with large effect sizes. In the present study, we analysed whether the risk for bipolar disorder (BD) in BD multiplex families is influenced by common genetic variants. Furthermore, we investigated whether this risk is conferred mainly by BD-specific risk variants or by variants also associated with the susceptibility to schizophrenia or major depression. In total, 395 individuals from 33 Andalusian BD multiplex families (166 BD, 78 major depressive disorder, 151 unaffected) as well as 438 subjects from an independent, BD case/control cohort (161 unrelated BD, 277 unrelated controls) were analysed. Polygenic risk scores (PRS) for BD, schizophrenia (SCZ), and major depression were calculated and compared between the cohorts. Both the familial BD cases and unaffected family members had higher PRS for all three psychiatric disorders than the independent controls, with BD and SCZ being significant after correction for multiple testing, suggesting a high baseline risk for several psychiatric disorders in the families. Moreover, familial BD cases showed significantly higher BD PRS than unaffected family members and unrelated BD cases. A plausible hypothesis is that, in multiplex families with a general increase in risk for psychiatric disease, BD development is attributable to a high burden of common variants that confer a specific risk for BD. The present analyses demonstrated that common genetic risk variants for psychiatric disorders are likely to contribute to the high incidence of affective psychiatric disorders in the multiplex families. However, the PRS explained only part of the observed phenotypic variance, and rare variants might have also contributed to disease development.

Genome-wide association study of panic disorder reveals genetic overlap with neuroticism and depression.

Panic disorder (PD) has a lifetime prevalence of 2-4% and heritability estimates of 40%. The contributory genetic variants remain largely unknown, with few and inconsistent loci having been reported. The present report describes the largest genome-wide association study (GWAS) of PD to date comprising genome-wide genotype data of 2248 clinically well-characterized PD patients and 7992 ethnically matched controls. The samples originated from four European countries (Denmark, Estonia, Germany, and Sweden). Standard GWAS quality control procedures were conducted on each individual dataset, and imputation was performed using the 1000 Genomes Project reference panel. A meta-analysis was then performed using the Ricopili pipeline. No genome-wide significant locus was identified. Leave-one-out analyses generated highly significant polygenic risk scores (PRS) (explained variance of up to 2.6%). Linkage disequilibrium (LD) score regression analysis of the GWAS data showed that the estimated heritability for PD was 28.0-34.2%. After correction for multiple testing, a significant genetic correlation was found between PD and major depressive disorder, depressive symptoms, and neuroticism. A total of 255 single-nucleotide polymorphisms (SNPs) with p < 1 × 10-4 were followed up in an independent sample of 2408 PD patients and 228,470 controls from Denmark, Iceland and the Netherlands. In the combined analysis, SNP rs144783209 showed the strongest association with PD (pcomb = 3.10 × 10-7). Sign tests revealed a significant enrichment of SNPs with a discovery p-value of <0.0001 in the combined follow up cohort (p = 0.048). The present integrative analysis represents a major step towards the elucidation of the genetic susceptibility to PD.

Association of polygenic score for major depression with response to lithium in patients with bipolar disorder.

Lithium is a first-line medication for bipolar disorder (BD), but only one in three patients respond optimally to the drug. Since evidence shows a strong clinical and genetic overlap between depression and bipolar disorder, we investigated whether a polygenic susceptibility to major depression is associated with response to lithium treatment in patients with BD. Weighted polygenic scores (PGSs) were computed for major depression (MD) at different GWAS p value thresholds using genetic data obtained from 2586 bipolar patients who received lithium treatment and took part in the Consortium on Lithium Genetics (ConLi+Gen) study. Summary statistics from genome-wide association studies in MD (135,458 cases and 344,901 controls) from the Psychiatric Genomics Consortium (PGC) were used for PGS weighting. Response to lithium treatment was defined by continuous scores and categorical outcome (responders versus non-responders) using measurements on the Alda scale. Associations between PGSs of MD and lithium treatment response were assessed using a linear and binary logistic regression modeling for the continuous and categorical outcomes, respectively. The analysis was performed for the entire cohort, and for European and Asian sub-samples. The PGSs for MD were significantly associated with lithium treatment response in multi-ethnic, European or Asian populations, at various p value thresholds. Bipolar patients with a low polygenic load for MD were more likely to respond well to lithium, compared to those patients with high polygenic load [lowest vs highest PGS quartiles, multi-ethnic sample: OR = 1.54 (95% CI: 1.18-2.01) and European sample: OR = 1.75 (95% CI: 1.30-2.36)]. While our analysis in the Asian sample found equivalent effect size in the same direction: OR = 1.71 (95% CI: 0.61-4.90), this was not statistically significant. Using PGS decile comparison, we found a similar trend of association between a high genetic loading for MD and lower response to lithium. Our findings underscore the genetic contribution to lithium response in BD and support the emerging concept of a lithium-responsive biotype in BD.

Pathway-Specific Genetic Risk for Alzheimer's Disease Differentiates Regional Patterns of Cortical Atrophy in Older Adults.

Brain aging is highly variable and represents a challenge to delimit aging from disease processes. Moreover, genetic factors may influence both aging and disease. Here we focused on this issue and investigated effects of multiple genetic loci previously identified to be associated with late-onset Alzheimer's disease (AD) on brain structure of older adults from a population sample. We calculated a genetic risk score (GRS) using genome-wide significant single-nucleotide polymorphisms from genome-wide association studies of AD and tested its effect on cortical thickness (CT). We observed a common pattern of cortical thinning (right inferior frontal, left posterior temporal, medial occipital cortex). To identify CT changes by specific biological processes, we subdivided the GRS effect according to AD-associated pathways and performed follow-up analyses. The common pattern from the main analysis was further differentiated by pathway-specific effects yielding a more bilateral pattern. Further findings were located in the superior parietal and mid/anterior cingulate regions representing 2 unique pathway-specific patterns. All patterns, except the superior parietal pattern, were influenced by apolipoprotein E. Our step-wise approach revealed atrophy patterns that partially resembled imaging findings in early stages of AD. Our study provides evidence that genetic burden for AD contributes to structural brain variability in normal aging.

No Reliable Association between Runs of Homozygosity and Schizophrenia in a Well-Powered Replication Study.

It is well known that inbreeding increases the risk of recessive monogenic diseases, but it is less certain whether it contributes to the etiology of complex diseases such as schizophrenia. One way to estimate the effects of inbreeding is to examine the association between disease diagnosis and genome-wide autozygosity estimated using runs of homozygosity (ROH) in genome-wide single nucleotide polymorphism arrays. Using data for schizophrenia from the Psychiatric Genomics Consortium (n = 21,868), Keller et al. (2012) estimated that the odds of developing schizophrenia increased by approximately 17% for every additional percent of the genome that is autozygous (ß = 16.1, CI(ß) = [6.93, 25.7], Z = 3.44, p = 0.0006). Here we describe replication results from 22 independent schizophrenia case-control datasets from the Psychiatric Genomics Consortium (n = 39,830). Using the same ROH calling thresholds and procedures as Keller et al. (2012), we were unable to replicate the significant association between ROH burden and schizophrenia in the independent PGC phase II data, although the effect was in the predicted direction, and the combined (original + replication) dataset yielded an attenuated but significant relationship between Froh and schizophrenia (ß = 4.86,CI(ß) = [0.90,8.83],Z = 2.40,p = 0.02). Since Keller et al. (2012), several studies reported inconsistent association of ROH burden with complex traits, particularly in case-control data. These conflicting results might suggest that the effects of autozygosity are confounded by various factors, such as socioeconomic status, education, urbanicity, and religiosity, which may be associated with both real inbreeding and the outcome measures of interest.

Copy number variation analysis and targeted NGS in 77 families with suspected Lynch syndrome reveals novel potential causative genes.

In many families with suspected Lynch syndrome (LS), no germline mutation in the causative mismatch repair (MMR) genes is detected during routine diagnostics. To identify novel causative genes for LS, the present study investigated 77 unrelated, mutation-negative patients with clinically suspected LS and a loss of MSH2 in tumor tissue. An analysis for genomic copy number variants (CNV) was performed, with subsequent next generation sequencing (NGS) of selected candidate genes in a subgroup of the cohort. Genomic DNA was genotyped using Illumina's HumanOmniExpress Bead Array. After quality control and filtering, 25 deletions and 16 duplications encompassing 73 genes were identified in 28 patients. No recurrent CNV was detected, and none of the CNVs affected the regulatory regions of MSH2. A total of 49 candidate genes from genomic regions implicated by the present CNV analysis and 30 known or assumed risk genes for colorectal cancer (CRC) were then sequenced in a subset of 38 patients using a customized NGS gene panel and Sanger sequencing. Single nucleotide variants were identified in 14 candidate genes from the CNV analysis. The most promising of these candidate genes were: (i) PRKCA, PRKDC, and MCM4, as a functional relation to MSH2 is predicted by network analysis, and (ii) CSMD1, as this is commonly mutated in CRC. Furthermore, six patients harbored POLE variants outside the exonuclease domain, suggesting that these might be implicated in hereditary CRC. Analyses in larger cohorts of suspected LS patients recruited via international collaborations are warranted to verify the present findings.

Common functional variants of the glutamatergic system in Autism spectrum disorder with high and low intellectual abilities.

The genetic architecture underlying Autism spectrum disorder (ASD) has been suggested to differ between individuals with lower (IQ ≤ 70; LIQ) and higher intellectual abilities (IQ > 70; HIQ). Among the identified pathomechanisms, the glutamatergic signalling pathway is of specific interest in ASD. We investigated 187 common functional variants of this neurotransmitter system for association with ASD and with symptom severity in two independent samples, a German (German-ALL: N = 583 families) and the Autism Genome Project cohort (AGP-ALL: N = 2001 families), split into HIQ, and LIQ subgroups. We did not identify any association withstanding correction for multiple testing. However, we report a replicated nominal significant under-transmission (OR < 0.79, p < 0.04) of the AKAP13 rs745191-T allele in both LIQ cohorts, but not in the much larger HIQ cohorts. At the phenotypic level, we nominally replicated associations of CAMK2A-rs2241694 with non-verbal communication in both combined LIQ and HIQ ASD cohorts. Variants PLD1-rs2124147 and ADCY1-rs2461127 were nominally associated with impaired non-verbal abilities and AKAP2-rs3739456 with repetitive behaviour in both LIQ cohorts. All four LIQ-associated genes are involved in G-protein coupled signal transduction, a downstream pathway of metabotropic glutamate receptor activation. We conclude that functional common variants of glutamatergic genes do not have a strong impact on ASD, but seem to moderately affect ASD risk and phenotypic expression. Since most of our nominally replicated hits were identified in the LIQ cohort, further investigation of the glutamatergic system in this subpopulation might be warranted.

Hypomorphic mutations in POLR3A are a frequent cause of sporadic and recessive spastic ataxia.

Despite extensive efforts, half of patients with rare movement disorders such as hereditary spastic paraplegias and cerebellar ataxias remain genetically unexplained, implicating novel genes and unrecognized mutations in known genes. Non-coding DNA variants are suspected to account for a substantial part of undiscovered causes of rare diseases. Here we identified mutations located deep in introns of POLR3A to be a frequent cause of hereditary spastic paraplegia and cerebellar ataxia. First, whole-exome sequencing findings in a recessive spastic ataxia family turned our attention to intronic variants in POLR3A, a gene previously associated with hypomyelinating leukodystrophy type 7. Next, we screened a cohort of hereditary spastic paraplegia and cerebellar ataxia cases (n = 618) for mutations in POLR3A and identified compound heterozygous POLR3A mutations in ∼3.1% of index cases. Interestingly, >80% of POLR3A mutation carriers presented the same deep-intronic mutation (c.1909+22G>A), which activates a cryptic splice site in a tissue and stage of development-specific manner and leads to a novel distinct and uniform phenotype. The phenotype is characterized by adolescent-onset progressive spastic ataxia with frequent occurrence of tremor, involvement of the central sensory tracts and dental problems (hypodontia, early onset of severe and aggressive periodontal disease). Instead of the typical hypomyelination magnetic resonance imaging pattern associated with classical POLR3A mutations, cases carrying c.1909+22G>A demonstrated hyperintensities along the superior cerebellar peduncles. These hyperintensities may represent the structural correlate to the cerebellar symptoms observed in these patients. The associated c.1909+22G>A variant was significantly enriched in 1139 cases with spastic ataxia-related phenotypes as compared to unrelated neurological and non-neurological phenotypes and healthy controls (P = 1.3 × 10-4). In this study we demonstrate that (i) autosomal-recessive mutations in POLR3A are a frequent cause of hereditary spastic ataxias, accounting for about 3% of hitherto genetically unclassified autosomal recessive and sporadic cases; and (ii) hypomyelination is frequently absent in POLR3A-related syndromes, especially when intronic mutations are present, and thus can no longer be considered as the unifying feature of POLR3A disease. Furthermore, our results demonstrate that substantial progress in revealing the causes of Mendelian diseases can be made by exploring the non-coding sequences of the human genome.

Genetic variants associated with response to lithium treatment in bipolar disorder: a genome-wide association study.

BACKGROUND: Lithium is a first-line treatment in bipolar disorder, but individual response is variable. Previous studies have suggested that lithium response is a heritable trait. However, no genetic markers of treatment response have been reproducibly identified. METHODS: Here, we report the results of a genome-wide association study of lithium response in 2563 patients collected by 22 participating sites from the International Consortium on Lithium Genetics (ConLiGen). Data from common single nucleotide polymorphisms (SNPs) were tested for association with categorical and continuous ratings of lithium response. Lithium response was measured using a well established scale (Alda scale). Genotyped SNPs were used to generate data at more than 6 million sites, using standard genomic imputation methods. Traits were regressed against genotype dosage. Results were combined across two batches by meta-analysis. FINDINGS: A single locus of four linked SNPs on chromosome 21 met genome-wide significance criteria for association with lithium response (rs79663003, p=1·37 × 10(-8); rs78015114, p=1·31 × 10(-8); rs74795342, p=3·31 × 10(-9); and rs75222709, p=3·50 × 10(-9)). In an independent, prospective study of 73 patients treated with lithium monotherapy for a period of up to 2 years, carriers of the response-associated alleles had a significantly lower rate of relapse than carriers of the alternate alleles (p=0·03268, hazard ratio 3·8, 95% CI 1·1-13·0). INTERPRETATION: The response-associated region contains two genes for long, non-coding RNAs (lncRNAs), AL157359.3 and AL157359.4. LncRNAs are increasingly appreciated as important regulators of gene expression, particularly in the CNS. Confirmed biomarkers of lithium response would constitute an important step forward in the clinical management of bipolar disorder. Further studies are needed to establish the biological context and potential clinical utility of these findings. FUNDING: Deutsche Forschungsgemeinschaft, National Institute of Mental Health Intramural Research Program.

Successful Replication of GWAS Hits for Multiple Sclerosis in 10,000 Germans Using the Exome Array.

Genome-wide association studies (GWAS) successfully identified various chromosomal regions to be associated with multiple sclerosis (MS). The primary aim of this study was to replicate reported associations from GWAS using an exome array in a large German study. German MS cases (n = 4,476) and German controls (n = 5,714) were genotyped using the Illumina HumanExome v1-Chip. Genotype calling was performed with the Illumina Genome Studio(TM) Genotyping Module, followed by zCall. Single-nucleotide polymorphisms (SNPs) in seven regions outside the human leukocyte antigen (HLA) region showed genome-wide significant associations with MS (P values < 5 × 10(-8) ). These associations have been reported previously. In addition, SNPs in three previously reported regions outside the HLA region yielded P values < 10(-5) . The effect of nine SNPs in the HLA region remained (P < 10(-5) ) after adjustment for other significant SNPs in the HLA region. All of these findings have been reported before or are driven by known risk loci. In summary, findings from previous GWAS for MS could be successfully replicated. We conclude that the regions identified in previous GWAS are also associated in the German population. This reassures the need for detailed investigations of the functional mechanisms underlying the replicated associations.

SUCLG2 identified as both a determinator of CSF Aß1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aß1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aß1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aß1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aß1-42.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy.

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

Runs of homozygosity and inbreeding in thyroid cancer.

BACKGROUND: Genome-wide association studies (GWASs) have identified several single-nucleotide polymorphisms (SNPs) influencing the risk of thyroid cancer (TC). Most cancer predisposition genes identified through GWASs function in a co-dominant manner, and studies have not found evidence for recessively functioning disease loci in TC. Our study examines whether homozygosity is associated with an increased risk of TC and searches for novel recessively acting disease loci. METHODS: Data from a previously conducted GWAS were used for the estimation of the proportion of phenotypic variance explained by all common SNPs, the detection of runs of homozygosity (ROH) and the determination of inbreeding to unravel their influence on TC. RESULTS: Inbreeding coefficients were significantly higher among cases than controls. Association on a SNP-by-SNP basis was controlled by using the false discovery rate at a level of q* < 0.05, with 34 SNPs representing true differences in homozygosity between cases and controls. The average size, the number and total length of ROHs per person were significantly higher in cases than in controls. A total of 16 recurrent ROHs of rather short length were identified although their association with TC risk was not significant at a genome-wide level. Several recurrent ROHs harbor genes associated with risk of TC. All of the ROHs showed significant evidence for natural selection (iHS, Fst, Fay and Wu's H). CONCLUSIONS: Our results support the existence of recessive alleles in TC susceptibility. Although regions of homozygosity were rather small, it might be possible that variants within these ROHs affect TC risk and may function in a recessive manner.

Association of the OPRM1 Variant rs1799971 (A118G) with Non-Specific Liability to Substance Dependence in a Collaborative de novo Meta-Analysis of European-Ancestry Cohorts.

The mu1 opioid receptor gene, OPRM1, has long been a high-priority candidate for human genetic studies of addiction. Because of its potential functional significance, the non-synonymous variant rs1799971 (A118G, Asn40Asp) in OPRM1 has been extensively studied, yet its role in addiction has remained unclear, with conflicting association findings. To resolve the question of what effect, if any, rs1799971 has on substance dependence risk, we conducted collaborative meta-analyses of 25 datasets with over 28,000 European-ancestry subjects. We investigated non-specific risk for "general" substance dependence, comparing cases dependent on any substance to controls who were non-dependent on all assessed substances. We also examined five specific substance dependence diagnoses: DSM-IV alcohol, opioid, cannabis, and cocaine dependence, and nicotine dependence defined by the proxy of heavy/light smoking (cigarettes-per-day >20 vs. ≤ 10). The G allele showed a modest protective effect on general substance dependence (OR = 0.90, 95% C.I. [0.83-0.97], p value = 0.0095, N = 16,908). We observed similar effects for each individual substance, although these were not statistically significant, likely because of reduced sample sizes. We conclude that rs1799971 contributes to mechanisms of addiction liability that are shared across different addictive substances. This project highlights the benefits of examining addictive behaviors collectively and the power of collaborative data sharing and meta-analyses.