Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies

Apolipoprotein B gene expression and regulation in relation to Alzheimer's disease pathophysiology

Apolipoprotein B (APOB), a receptor-binding protein present in cholesterol-rich lipoproteins, has been implicated in Alzheimer's disease (AD). High levels of APOB-containing low-density lipoproteins (LDL) are linked to the pathogenesis of both early-onset familial and late onset sporadic AD. Rare coding mutations in the APOB gene are associated with familial AD, suggesting a role for APOB-bound lipoproteins in the central nervous system. This research explores APOB gene regulation across the AD spectrum using four cohorts: BRAINEAC (elderly control brains), DBCBB (controls, AD brains), ROSMAP (controls, MCI, AD brains), and ADNI (control, MCI, AD clinical subjects). APOB protein levels, measured via mass spectrometry and ELISA, positively correlated with AD pathology indices and cognition, while APOB mRNA levels showed negative correlations. Brain APOB protein levels also correlated with cortical Aβ levels. A common coding variant in the APOB gene locus affected its expression but didn't impact AD risk or brain cholesterol concentrations, except for 24-S-hydroxycholesterol. Polymorphisms in the CYP27A1 gene, notably rs4674344, were associated with APOB protein levels. A negative correlation was observed between brain APOB gene expression and AD biomarker levels. CSF APOB correlated with Tau pathology in presymptomatic subjects, while cortical APOB was strongly associated with cortical Aβ deposition in late-stage AD. The study discusses the potential link between blood-brain barrier dysfunction and AD symptoms in relation to APOB neurobiology. Overall, APOB's involvement in lipoprotein metabolism appears to influence AD pathology across different stages of the disease.

PTPRS is a novel marker for early Tau pathology and synaptic integrity in Alzheimer's disease

We examined the role of protein tyrosine phosphatase receptor sigma (PTPRS) in the context of Alzheimer's disease and synaptic integrity. Publicly available datasets (BRAINEAC, ROSMAP, ADC1) and a cohort of asymptomatic but "at risk" individuals (PREVENT-AD) were used to explore the relationship between PTPRS and various Alzheimer's disease biomarkers. We identified that PTPRS rs10415488 variant C shows features of neuroprotection against early Tau pathology and synaptic degeneration in Alzheimer's disease. This single nucleotide polymorphism correlated with higher PTPRS transcript abundance and lower p(181)Tau and GAP-43 levels in the CSF. In the brain, PTPRS protein abundance was significantly correlated with the quantity of two markers of synaptic integrity: SNAP25 and SYT-1. We also found the presence of sexual dimorphism for PTPRS, with higher CSF concentrations in males than females. Male carriers for variant C were found to have a 10-month delay in the onset of AD. We thus conclude that PTPRS acts as a neuroprotective receptor in Alzheimer's disease. Its protective effect is most important in males, in whom it postpones the age of onset of the disease.

PTPRS is a novel marker for early tau pathology and synaptic integrity in Alzheimer's disease

We examined the role of protein tyrosine phosphatase receptor sigma (PTPRS) in the context of Alzheimer's disease and synaptic integrity. Publicly available datasets (BRAINEAC, ROSMAP, ADC1) and a cohort of asymptomatic but "at risk" individuals (PREVENT-AD) were used to explore the relationship between PTPRS and various Alzheimer's disease biomarkers. We identified that PTPRS rs10415488 variant C shows features of neuroprotection against early tau pathology and synaptic degeneration in Alzheimer's disease. This single nucleotide polymorphism correlated with higher PTPRS transcript abundance and lower P-tau181 and GAP-43 levels in the CSF. In the brain, PTPRS protein abundance was significantly correlated with the quantity of two markers of synaptic integrity: SNAP25 and SYT-1. We also found the presence of sexual dimorphism for PTPRS, with higher CSF concentrations in males than females. Male carriers for variant C were found to have a 10-month delay in the onset of AD. We thus conclude that PTPRS acts as a neuroprotective receptor in Alzheimer's disease. Its protective effect is most important in males, in whom it postpones the age of onset of the disease.

A miR-137-Related Biological Pathway of Risk for Schizophrenia Is Associated With Human Brain Emotion Processing

BACKGROUND

miR-137 is a microRNA involved in brain development, regulating neurogenesis and neuronal maturation. Genome-wide association studies have implicated miR-137 in schizophrenia risk but do not explain its involvement in brain function and underlying biology. Polygenic risk for schizophrenia mediated by miR-137 targets is associated with working memory, although other evidence points to emotion processing. We characterized the functional brain correlates of miR-137 target genes associated with schizophrenia while disentangling previously reported associations of miR-137 targets with working memory and emotion processing.

METHODS

Using RNA sequencing data from postmortem prefrontal cortex (N = 522), we identified a coexpression gene set enriched for miR-137 targets and schizophrenia risk genes. We validated the relationship of this set to miR-137 in vitro by manipulating miR-137 expression in neuroblastoma cells. We translated this gene set into polygenic scores of coexpression prediction and associated them with functional magnetic resonance imaging activation in healthy volunteers (n1 = 214; n2 = 136; n3 = 2075; n4 = 1800) and with short-term treatment response in patients with schizophrenia (N = 427).

RESULTS

In 4652 human participants, we found that 1) schizophrenia risk genes were coexpressed in a biologically validated set enriched for miR-137 targets; 2) increased expression of miR-137 target risk genes was mediated by low prefrontal miR-137 expression; 3) alleles that predict greater gene set coexpression were associated with greater prefrontal activation during emotion processing in 3 independent healthy cohorts (n1, n2, n3) in interaction with age (n4); and 4) these alleles predicted less improvement in negative symptoms following antipsychotic treatment in patients with schizophrenia.

CONCLUSIONS

The functional translation of miR-137 target gene expression linked with schizophrenia involves the neural substrates of emotion processing.

Bi-allelic genetic variants in the translational GTPases GTPBP1 and GTPBP2 cause a distinct identical neurodevelopmental syndrome

The homologous genes GTPBP1 and GTPBP2 encode GTP-binding proteins 1 and 2, which are involved in ribosomal homeostasis. Pathogenic variants in GTPBP2 were recently shown to be an ultra-rare cause of neurodegenerative or neurodevelopmental disorders (NDDs). Until now, no human phenotype has been linked to GTPBP1. Here, we describe individuals carrying bi-allelic GTPBP1 variants that display an identical phenotype with GTPBP2 and characterize the overall spectrum of GTP-binding protein (1/2)-related disorders. In this study, 20 individuals from 16 families with distinct NDDs and syndromic facial features were investigated by whole-exome (WES) or whole-genome (WGS) sequencing. To assess the functional impact of the identified genetic variants, semi-quantitative PCR, western blot, and ribosome profiling assays were performed in fibroblasts from affected individuals. We also investigated the effect of reducing expression of CG2017, an ortholog of human GTPBP1/2, in the fruit fly Drosophila melanogaster. Individuals with bi-allelic GTPBP1 or GTPBP2 variants presented with microcephaly, profound neurodevelopmental impairment, pathognomonic craniofacial features, and ectodermal defects. Abnormal vision and/or hearing, progressive spasticity, choreoathetoid movements, refractory epilepsy, and brain atrophy were part of the core phenotype of this syndrome. Cell line studies identified a loss-of-function (LoF) impact of the disease-associated variants but no significant abnormalities on ribosome profiling. Reduced expression of CG2017 isoforms was associated with locomotor impairment in Drosophila. In conclusion, bi-allelic GTPBP1 and GTPBP2 LoF variants cause an identical, distinct neurodevelopmental syndrome. Mutant CG2017 knockout flies display motor impairment, highlighting the conserved role for GTP-binding proteins in CNS development across species.

Genetic topography and cortical cell loss in Huntington's disease link development and neurodegeneration

Cortical cell loss is a core feature of Huntington's disease (HD), beginning many years before clinical motor diagnosis, during the premanifest stage. However, it is unclear how genetic topography relates to cortical cell loss. Here, we explore the biological processes and cell types underlying this relationship and validate these using cell-specific post-mortem data. Eighty premanifest participants on average 15 years from disease onset and 71 controls were included. Using volumetric and diffusion MRI we extracted HD-specific whole brain maps where lower grey matter volume and higher grey matter mean diffusivity, relative to controls, were used as proxies of cortical cell loss. These maps were combined with gene expression data from the Allen Human Brain Atlas (AHBA) to investigate the biological processes relating genetic topography and cortical cell loss. Cortical cell loss was positively correlated with the expression of developmental genes (i.e. higher expression correlated with greater atrophy and increased diffusivity) and negatively correlated with the expression of synaptic and metabolic genes that have been implicated in neurodegeneration. These findings were consistent for diffusion MRI and volumetric HD-specific brain maps. As wild-type huntingtin is known to play a role in neurodevelopment, we explored the association between wild-type huntingtin (HTT) expression and developmental gene expression across the AHBA. Co-expression network analyses in 134 human brains free of neurodegenerative disorders were also performed. HTT expression was correlated with the expression of genes involved in neurodevelopment while co-expression network analyses also revealed that HTT expression was associated with developmental biological processes. Expression weighted cell-type enrichment (EWCE) analyses were used to explore which specific cell types were associated with HD cortical cell loss and these associations were validated using cell specific single nucleus RNAseq (snRNAseq) data from post-mortem HD brains. The developmental transcriptomic profile of cortical cell loss in preHD was enriched in astrocytes and endothelial cells, while the neurodegenerative transcriptomic profile was enriched for neuronal and microglial cells. Astrocyte-specific genes differentially expressed in HD post-mortem brains relative to controls using snRNAseq were enriched in the developmental transcriptomic profile, while neuronal and microglial-specific genes were enriched in the neurodegenerative transcriptomic profile. Our findings suggest that cortical cell loss in preHD may arise from dual pathological processes, emerging as a consequence of neurodevelopmental changes, at the beginning of life, followed by neurodegeneration in adulthood, targeting areas with reduced expression of synaptic and metabolic genes. These events result in age-related cell death across multiple brain cell types.

Atlas of the aging mouse brain reveals white matter as vulnerable foci

Aging is the key risk factor for cognitive decline, yet the molecular changes underlying brain aging remain poorly understood. Here, we conducted spatiotemporal RNA sequencing of the mouse brain, profiling 1,076 samples from 15 regions across 7 ages and 2 rejuvenation interventions. Our analysis identified a brain-wide gene signature of aging in glial cells, which exhibited spatially defined changes in magnitude. By integrating spatial and single-nucleus transcriptomics, we found that glial aging was particularly accelerated in white matter compared with cortical regions, whereas specialized neuronal populations showed region-specific expression changes. Rejuvenation interventions, including young plasma injection and dietary restriction, exhibited distinct effects on gene expression in specific brain regions. Furthermore, we discovered differential gene expression patterns associated with three human neurodegenerative diseases, highlighting the importance of regional aging as a potential modulator of disease. Our findings identify molecular foci of brain aging, providing a foundation to target age-related cognitive decline.

Human hippocampal astrocytes: Computational dissection of their transcriptome, sexual differences and exosomes across ageing and mild-cognitive impairment

The role of astrocytes in Alzheimer's disease is often disregarded. Hence, characterization of astrocytes along their early evolution toward Alzheimer would be greatly beneficial. However, due to their exquisite responsiveness, in vivo studies are difficult. So public microarray data of hippocampal homogenates from (healthy) young, (healthy) elder and elder with mild cognitive impairment (MCI) were subjected to re-analysis by a multi-step computational pipeline. Ontologies and pathway analyses were compared after determining the differential genes that, belonging to astrocytes, have splice forms. Likewise, the subset of molecules exportable to exosomes was also determined. The results showed that astrocyte's phenotypes changed significantly. While already 'activated' astrocytes were found in the younger group, major changes occurred during ageing (increased vascular remodelling and response to mechanical stimulus, diminished long-term potentiation and increased long-term depression). MCI's astrocytes showed some 'rejuvenated' features, but their sensitivity to shear stress was markedly lost. Importantly, most of the changes showed to be sex biassed. Men's astrocytes are enriched in a type 'endfeet-astrocytome', whereas women's astrocytes appear close to the 'scar-forming' type (prone to endothelial dysfunction, hypercholesterolemia, loss of glutamatergic synapses, Ca+2 dysregulation, hypoxia, oxidative stress and 'pro-coagulant' phenotype). In conclusion, the computational dissection of the networks based on the hippocampal gene isoforms provides a relevant proxy to in vivo astrocytes, also revealing the occurrence of sexual differences. Analyses of the astrocytic exosomes did not provide an acceptable approximation to the overall functioning of astrocytes in the hippocampus, probably due to the selective cellular mechanisms which charge the cargo molecules.

HnRNP Pathologies in Frontotemporal Lobar Degeneration

Frontotemporal dementia (FTD) is the second most common form of young-onset (<65 years) dementia. Clinically, it primarily manifests as a disorder of behavioural, executive, and/or language functions. Pathologically, frontotemporal lobar degeneration (FTLD) is the predominant cause of FTD. FTLD is a proteinopathy, and the main pathological proteins identified so far are tau, TAR DNA-binding protein 43 (TDP-43), and fused in sarcoma (FUS). As TDP-43 and FUS are members of the heterogeneous ribonucleic acid protein (hnRNP) family, many studies in recent years have expanded the research on the relationship between other hnRNPs and FTLD pathology. Indeed, these studies provide evidence for an association between hnRNP abnormalities and FTLD. In particular, several studies have shown that multiple hnRNPs may exhibit nuclear depletion and cytoplasmic mislocalisation within neurons in FTLD cases. However, due to the diversity and complex association of hnRNPs, most studies are still at the stage of histological discovery of different hnRNP abnormalities in FTLD. We herein review the latest studies relating hnRNPs to FTLD. Together, these studies outline an important role of multiple hnRNPs in the pathogenesis of FTLD and suggest that future research into FTLD should include the whole spectrum of this protein family.

Age-related telomere attrition in the human putamen

Age is a major risk factor for neurodegenerative diseases. Shortening of leucocyte telomeres with advancing age, arguably a measure of "biological" age, is a known phenomenon and epidemiologically correlated with age-related disease. The main mechanism of telomere shortening is cell division, rendering telomere length in post-mitotic cells presumably stable. Longitudinal measurement of human brain telomere length is not feasible, and cross-sectional cortical brain samples so far indicated no attrition with age. Hence, age-related changes in telomere length in the brain and the association between telomere length and neurodegenerative diseases remain unknown. Here, we demonstrate that mean telomere length in the putamen, a part of the basal ganglia, physiologically shortens with age, like leukocyte telomeres. This was achieved by using matched brain and leukocyte-rich spleen samples from 98 post-mortem healthy human donors. Using spleen telomeres as a reference, we further found that mean telomere length was brain region-specific, as telomeres in the putamen were significantly shorter than in the cerebellum. Expression analyses of genes involved in telomere length regulation and oxidative phosphorylation revealed that both region- and age-dependent expression pattern corresponded with region-dependent telomere length dynamics. Collectively, our results indicate that mean telomere length in the human putamen physiologically shortens with advancing age and that both local and temporal gene expression dynamics correlate with this, pointing at a potential mechanism for the selective, age-related vulnerability of the nigro-striatal network.

Examining the biological mechanisms of human mental disorders resulting from gene-environment interdependence using novel functional genomic approaches

We explore how functional genomics approaches that integrate datasets from human and non-human model systems can improve our understanding of the effect of gene-environment interplay on the risk for mental disorders. We start by briefly defining the G-E paradigm and its challenges and then discuss the different levels of regulation of gene expression and the corresponding data existing in humans (genome wide genotyping, transcriptomics, DNA methylation, chromatin modifications, chromosome conformational changes, non-coding RNAs, proteomics and metabolomics), discussing novel approaches to the application of these data in the study of the origins of mental health. Finally, we discuss the multilevel integration of diverse types of data. Advance in the use of functional genomics in the context of a G-E perspective improves the detection of vulnerabilities, informing the development of preventive and therapeutic interventions.

Gene variations of glutamate metabolism pathway and epilepsy

Background

Epilepsy is a paroxysmal disorder of the brain, caused by an imbalance of neuronal excitation and inhibition. Glutamate is the most important excitatory neurotransmitter in the brain and plays an important role in epileptogenesis. Mutations in genes at any step/component of the glutamate metabolic pathway may lead to the development of epilepsy or epileptic encephalopathy.

Methods

Clinical history of 3 epilepsy patients with genetic variations of the glutamate metabolism pathway was collected. Electroencephalogram recording and magnetic resonance imaging  were performed in each patient. We also reviewed recent literature for a variety of the genetic variations involved in epilepsy.

Results

Case 1 was a SLC1A2 mutation-carrier diagnosed with developmental and epileptic encephalopathy (DEE) 41, whose seizures decreased after start of the ketogenic diet. Case 2 carried a GRIN2A gene mutation and was seizure-free for three years after taking levetiracetam and vitamin B6. Case 3 was a GRIN2B mutation-carrier  diagnosed with DEE 27, who seizures diminished after taking oxcarbazepine.

Conclusions

Preclinical and clinical evidence supports the therapeutic potential of glutamatergic signaling-targeting treatments for epilepsy. More studies are needed to discover novel DEE-related genetic mutations in the glutamate metabolic pathway.

Differential Gene Expression in Sporadic and Genetic Forms of Alzheimer's Disease and Frontotemporal Dementia in Brain Tissue and Lymphoblastoid Cell Lines

Sporadic early-onset Alzheimer's disease (EOAD) and autosomal dominant Alzheimer's disease (ADAD) provide the opportunity to investigate the physiopathological mechanisms in the absence of aging, present in late-onset forms. Frontotemporal dementia (FTD) causes early-onset dementia associated to tau or TDP43 protein deposits. A 15% of FTD cases are caused by mutations in C9orf72, GRN, or MAPT genes. Lymphoblastoid cell lines (LCLs) have been proposed as an alternative to brain tissue for studying earlier phases of neurodegenerative diseases. The aim of this study is to investigate the expression profile in EOAD, ADAD, and sporadic and genetic FTD (sFTD and gFTD, respectively), using brain tissue and LCLs. Sixty subjects of the following groups were included: EOAD, ADAD, sFTD, gFTD, and controls. Gene expression was analyzed with Clariom D microarray (Affymetrix). Brain tissue pairwise comparisons revealed six common differentially expressed genes (DEG) for all the patients' groups compared with controls: RGS20, WIF1, HSPB1, EMP3, S100A11 and GFAP. Common up-regulated biological pathways were identified both in brain and LCLs (including inflammation and glial cell differentiation), while down-regulated pathways were detected mainly in brain tissue (including synaptic signaling, metabolism and mitochondrial dysfunction). CD163, ADAMTS9 and LIN7A gene expression disruption was validated by qPCR in brain tissue and NrCAM in LCLs in their respective group comparisons. In conclusion, our study highlights neuroinflammation, metabolism and synaptic signaling disturbances as common altered pathways in different AD and FTD forms. The use of LCLs might be appropriate for studying early immune system and inflammation, and some neural features in neurodegenerative dementias.

Mesial Temporal Lobe Epilepsy (MTLE) Drug-Refractoriness Is Associated With P2X7 Receptors Overexpression in the Human Hippocampus and Temporal Neocortex and May Be Predicted by Low Circulating Levels of miR-22

Objective: ATP-gated ionotropic P2X7 receptors (P2X7R) actively participate in epilepsy and other neurological disorders. Neocortical nerve terminals of patients with Mesial Temporal Lobe Epilepsy with Hippocampal Sclerosis (MTLE-HS) express higher P2X7R amounts. Overexpression of P2X7R bolsters ATP signals during seizures resulting in glial cell activation, cytokines production, and GABAergic rundown with unrestrained glutamatergic excitation. In a mouse model of status epilepticus, increased expression of P2X7R has been associated with the down-modulation of the non-coding micro RNA, miR-22. MiR levels are stable in biological fluids and normally reflect remote tissue production making them ideal disease biomarkers. Here, we compared P2X7R and miR-22 expression in epileptic brains and in the serum of patients with MTLE-HS, respectively.

Methods: Quantitative RT-PCR was used to evaluate the expression of P2X7R in the hippocampus and anterior temporal lobe of 23 patients with MTLE-HS and 10 cadaveric controls. Confocal microscopy and Western blot analysis were performed to assess P2X7R protein amounts. MiR-22 expression was evaluated in cell-free sera of 40 MTLE-HS patients and 48 healthy controls.

Results: Nerve terminals of the hippocampus and neocortical temporal lobe of MTLE-HS patients overexpress (p &lt; 0.05) an 85 kDa P2X7R protein whereas the normally occurring 67 kDa receptor protein dominates in the brain of the cadaveric controls. Contrariwise, miR-22 serum levels are diminished (p &lt; 0.001) in MTLE-HS patients compared to age-matched control blood donors, a situation that is more evident in patients requiring multiple (&gt;3) anti-epileptic drug (AED) regimens.

Conclusion: Data show that there is an inverse relationship between miR-22 serum levels and P2X7R expression in the hippocampus and neocortex of MTLE-HS patients, which implies that measuring serum miR-22 may be a clinical surrogate of P2X7R brain expression in the MTLE-HS. Moreover, the high area under the ROC curve (0.777; 95% CI 0.629–0.925; p = 0.001) suggests that low miR-22 serum levels may be a sensitive predictor of poor response to AEDs among MTLE-HS patients. Results also anticipate that targeting the miR-22/P2X7R axis may be a good strategy to develop newer AEDs.

A novel missense variant in the LMNB2 gene causes progressive myoclonus epilepsy

Progressive myoclonus epilepsies (PMEs) are a group of disorders embracing myoclonus, seizures, and neurological dysfunctions. Because of the genetic and clinical heterogeneity, a large proportion of PMEs cases have remained molecularly undiagnosed. The present study aimed to determine the underlying genetic factors that contribute to the PME phenotype in an Iranian female patient. We describe a consanguineous Iranian family with autosomal recessive PME that had remained undiagnosed despite extensive genetic and pathological tests. After performing neuroimaging and clinical examinations, due to heterogeneity of PMEs, the proband was subjected to paired-end whole-exome sequencing and the candidate variant was confirmed by Sanger sequencing. Various in-silico tools were also used to predict the pathogenicity of the variant. In this study, we identified a novel homozygous missense variant (NM_032737.4:c.472C > T; p.(Arg158Trp)) in the LMNB2 gene (OMIM: 150341) as the most likely disease-causing variant. Neuroimaging revealed a progressive significant generalized atrophy in the cerebral and cerebellum without significant white matter signal changes. Video-electroencephalography monitoring showed a generalized pattern of high-voltage sharp waves in addition to multifocal spikes and waves compatible with mixed type seizures and epileptic encephalopathic pattern. Herein, we introduce the second case of PME caused by a novel variant in the LMNB2 gene. This study also underscores the potentiality of next-generation sequencing in the genetic diagnosis of patients with neurologic diseases with an unknown cause.

Multi-omics signatures of alcohol use disorder in the dorsal and ventral striatum

Alcohol Use Disorder (AUD) is a major contributor to global mortality and morbidity. Postmortem human brain tissue enables the investigation of molecular mechanisms of AUD in the neurocircuitry of addiction. We aimed to identify differentially expressed (DE) genes in the ventral and dorsal striatum between individuals with AUD and controls, and to integrate the results with findings from genome- and epigenome-wide association studies (GWAS/EWAS) to identify functionally relevant molecular mechanisms of AUD. DNA-methylation and gene expression (RNA-seq) data was generated from postmortem brain samples of 48 individuals with AUD and 51 controls from the ventral striatum (VS) and the dorsal striatal regions caudate nucleus (CN) and putamen (PUT). We identified DE genes using DESeq2, performed gene-set enrichment analysis (GSEA), and tested enrichment of DE genes in results of GWASs using MAGMA. Weighted correlation network analysis (WGCNA) was performed for DNA-methylation and gene expression data and gene overlap was tested. Differential gene expression was observed in the dorsal (FDR < 0.05), but not the ventral striatum of AUD cases. In the VS, DE genes at FDR < 0.25 were overrepresented in a recent GWAS of problematic alcohol use. The ARHGEF15 gene was upregulated in all three brain regions. GSEA in CN and VS pointed towards cell-structure associated GO-terms and in PUT towards immune pathways. The WGCNA modules most strongly associated with AUD showed strong enrichment for immune response and inflammation pathways. Our integrated analysis of multi-omics data sets provides further evidence for the importance of immune- and inflammation-related processes in AUD.

Functional characterisation of the amyotrophic lateral sclerosis risk locus GPX3/TNIP1

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a complex, late-onset, neurodegenerative disease with a genetic contribution to disease liability. Genome-wide association studies (GWAS) have identified ten risk loci to date, including the TNIP1/GPX3 locus on chromosome five. Given association analysis data alone cannot determine the most plausible risk gene for this locus, we undertook a comprehensive suite of in silico, in vivo and in vitro studies to address this.

METHODS

The Functional Mapping and Annotation (FUMA) pipeline and five tools (conditional and joint analysis (GCTA-COJO), Stratified Linkage Disequilibrium Score Regression (S-LDSC), Polygenic Priority Scoring (PoPS), Summary-based Mendelian Randomisation (SMR-HEIDI) and transcriptome-wide association study (TWAS) analyses) were used to perform bioinformatic integration of GWAS data (Ncases = 20,806, Ncontrols = 59,804) with 'omics reference datasets including the blood (eQTLgen consortium N = 31,684) and brain (N = 2581). This was followed up by specific expression studies in ALS case-control cohorts (microarray Ntotal = 942, protein Ntotal = 300) and gene knockdown (KD) studies of human neuronal iPSC cells and zebrafish-morpholinos (MO).

RESULTS

SMR analyses implicated both TNIP1 and GPX3 (p < 1.15 × 10-6), but there was no simple SNP/expression relationship. Integrating multiple datasets using PoPS supported GPX3 but not TNIP1. In vivo expression analyses from blood in ALS cases identified that lower GPX3 expression correlated with a more progressed disease (ALS functional rating score, p = 5.5 × 10-3, adjusted R2 = 0.042, Beffect = 27.4 ± 13.3 ng/ml/ALSFRS unit) with microarray and protein data suggesting lower expression with risk allele (recessive model p = 0.06, p = 0.02 respectively). Validation in vivo indicated gpx3 KD caused significant motor deficits in zebrafish-MO (mean difference vs. control ± 95% CI, vs. control, swim distance = 112 ± 28 mm, time = 1.29 ± 0.59 s, speed = 32.0 ± 2.53 mm/s, respectively, p for all < 0.0001), which were rescued with gpx3 expression, with no phenotype identified with tnip1 KD or gpx3 overexpression.

CONCLUSIONS

These results support GPX3 as a lead ALS risk gene in this locus, with more data needed to confirm/reject a role for TNIP1. This has implications for understanding disease mechanisms (GPX3 acts in the same pathway as SOD1, a well-established ALS-associated gene) and identifying new therapeutic approaches. Few previous examples of in-depth investigations of risk loci in ALS exist and a similar approach could be applied to investigate future expected GWAS findings.

Simultaneous Isolation of High-Quality RNA and DNA From Postmortem Human Central Nervous System Tissues for Omics Studies

Multi-omics approaches are increasingly being adopted to understand the complex networks underlying disease. The coisolation of high-quality nucleotides from affected tissues is paramount for the parallel analysis of transcriptomic, genomic, and epigenomic data sets. Although nucleotides extracted from postmortem central nervous system (CNS) tissue are widely used in the study of neurodegenerative disease, assessment of methods for the simultaneous isolation of DNA and RNA is limited. Herein, we describe a strategy for the isolation of high-quality DNA and RNA from postmortem human tissue from 7 CNS regions. Motor cortex, frontal cortex, hippocampus, occipital cortex, anterior cingulate cortex, cerebellum, and spinal cord tissues were obtained from 22 individuals diagnosed with motor neuron disease (MND) and 13 neurologically normal controls (n = 245 tissues). We demonstrated that the Qiagen AllPrep DNA/RNA kit consistently isolated DNA and RNA of high yield and quality from all 6 brain regions. Importantly, phenol-chloroform-based extraction was required to isolate high-yield RNA from spinal cord. RNA sequencing using RNA extracted from 6 CNS regions (n = 60) generated high-quality transcriptomes. Hierarchical clustering of data from motor cortex, using an MND susceptibility gene panel and marker genes of disease-associated microglia, demonstrated that MND-specific gene expression signatures could be detected in the transcriptome data.

Transcriptomic Changes Highly Similar to Alzheimer's Disease Are Observed in a Subpopulation of Individuals During Normal Brain Aging

Aging is a major risk factor for late-onset Alzheimer’s disease (LOAD). How aging contributes to the development of LOAD remains elusive. In this study, we examined multiple large-scale transcriptomic datasets from both normal aging and LOAD brains to understand the molecular interconnection between aging and LOAD. We found that shared gene expression changes between aging and LOAD are mostly seen in the hippocampal and several cortical regions. In the hippocampus, the expression of phosphoprotein, alternative splicing and cytoskeleton genes are commonly changed in both aging and AD, while synapse, ion transport, and synaptic vesicle genes are commonly down-regulated. Aging-specific changes are associated with acetylation and methylation, while LOAD-specific changes are more related to glycoprotein (both up- and down-regulations), inflammatory response (up-regulation), myelin sheath and lipoprotein (down-regulation). We also found that normal aging brain transcriptomes from relatively young donors (45–70 years old) clustered into several subgroups and some subgroups showed gene expression changes highly similar to those seen in LOAD brains. Using brain transcriptomic datasets from another cohort of older individuals (>70 years), we found that samples from cognitively normal older individuals clustered with the “healthy aging” subgroup while AD samples mainly clustered with the “AD similar” subgroups. This may imply that individuals in the healthy aging subgroup will likely remain cognitively normal when they become older and vice versa. In summary, our results suggest that on the transcriptome level, aging and LOAD have strong interconnections in some brain regions in a subpopulation of cognitively normal aging individuals. This supports the theory that the initiation of LOAD occurs decades earlier than the manifestation of clinical phenotype and it may be essential to closely study the “normal brain aging” to identify the very early molecular events that may lead to LOAD development.

Multi-omic insights into Parkinson's Disease: From genetic associations to functional mechanisms

Genome-Wide Association Studies (GWAS) have elucidated the genetic components of Parkinson's Disease (PD). However, because the vast majority of GWAS association signals fall within non-coding regions, translating these results into an interpretable, mechanistic understanding of the disease etiology remains a major challenge in the field. In this review, we provide an overview of the approaches to prioritize putative causal variants and genes as well as summarise the primary findings of previous studies. We then discuss recent efforts to integrate multi-omics data to identify likely pathogenic cell types and biological pathways implicated in PD pathogenesis. We have compiled full summary statistics of cell-type, tissue, and phentoype enrichment analyses from multiple studies of PD GWAS and provided them in a standardized format as a resource for the research community (https://github.com/RajLabMSSM/PD_omics_review). Finally, we discuss the experimental, computational, and conceptual advances that will be necessary to fully elucidate the effects of functional variants and genes on cellular dysregulation and disease risk.

HnRNP K mislocalisation is a novel protein pathology of frontotemporal lobar degeneration and ageing and leads to cryptic splicing

Heterogeneous nuclear ribonucleoproteins (HnRNPs) are a group of ubiquitously expressed RNA-binding proteins implicated in the regulation of all aspects of nucleic acid metabolism. HnRNP K is a member of this highly versatile hnRNP family. Pathological redistribution of hnRNP K to the cytoplasm has been linked to the pathogenesis of several malignancies but, until now, has been underexplored in the context of neurodegenerative disease. Here we show hnRNP K mislocalisation in pyramidal neurons of the frontal cortex to be a novel neuropathological feature that is associated with both frontotemporal lobar degeneration and ageing. HnRNP K mislocalisation is mutually exclusive to TDP-43 and tau pathological inclusions in neurons and was not observed to colocalise with mitochondrial, autophagosomal or stress granule markers. De-repression of cryptic exons in RNA targets following TDP-43 nuclear depletion is an emerging mechanism of potential neurotoxicity in frontotemporal lobar degeneration and the mechanistically overlapping disorder amyotrophic lateral sclerosis. We silenced hnRNP K in neuronal cells to identify the transcriptomic consequences of hnRNP K nuclear depletion. Intriguingly, by performing RNA-seq analysis we find that depletion of hnRNP K induces 101 novel cryptic exon events. We validated cryptic exon inclusion in an SH-SY5Y hnRNP K knockdown and in FTLD brain exhibiting hnRNP K nuclear depletion. We, therefore, present evidence for hnRNP K mislocalisation to be associated with FTLD and for this to induce widespread changes in splicing.

Physiological and Pathological Ageing of Astrocytes in the Human Brain

Ageing is the greatest risk factor for dementia, although physiological ageing by itself does not lead to cognitive decline. In addition to ageing, APOE ε4 is genetically the strongest risk factor for Alzheimer's disease and is highly expressed in astrocytes. There are indications that human astrocytes change with age and upon expression of APOE4. As these glial cells maintain water and ion homeostasis in the brain and regulate neuronal transmission, it is likely that age- and APOE4-related changes in astrocytes have a major impact on brain functioning and play a role in age-related diseases. In this review, we will discuss the molecular and morphological changes of human astrocytes in ageing and the contribution of APOE4. We conclude this review with a discussion on technical issues, innovations, and future perspectives on how to gain more knowledge on astrocytes in the human ageing brain.

Prediction of Alzheimer's disease using multi-variants from a Chinese genome-wide association study

Previous genome-wide association studies have identified dozens of susceptibility loci for sporadic Alzheimer’s disease, but few of these loci have been validated in longitudinal cohorts. Establishing predictive models of Alzheimer’s disease based on these novel variants is clinically important for verifying whether they have pathological functions and provide a useful tool for screening of disease risk. In the current study, we performed a two-stage genome-wide association study of 3913 patients with Alzheimer’s disease and 7593 controls and identified four novel variants (rs3777215, rs6859823, rs234434, and rs2255835; P_combined = 3.07 × 10⁻¹⁹, 2.49 × 10⁻²³, 1.35 × 10⁻⁶⁷, and 4.81 × 10⁻⁹, respectively) as well as nine variants in the apolipoprotein E region with genome-wide significance (P < 5.0 × 10⁻⁸). Literature mining suggested that these novel single nucleotide polymorphisms are related to amyloid precursor protein transport and metabolism, antioxidation, and neurogenesis. Based on their possible roles in the development of Alzheimer’s disease, we used different combinations of these variants and the apolipoprotein E status and successively built 11 predictive models. The predictive models include relatively few single nucleotide polymorphisms useful for clinical practice, in which the maximum number was 13 and the minimum was only four. These predictive models were all significant and their peak of area under the curve reached 0.73 both in the first and second stages. Finally, these models were validated using a separate longitudinal cohort of 5474 individuals. The results showed that individuals carrying risk variants included in the models had a shorter latency and higher incidence of Alzheimer’s disease, suggesting that our models can predict Alzheimer’s disease onset in a population with genetic susceptibility. The effectiveness of the models for predicting Alzheimer’s disease onset confirmed the contributions of these identified variants to disease pathogenesis. In conclusion, this is the first study to validate genome-wide association study-based predictive models for evaluating the risk of Alzheimer’s disease onset in a large Chinese population. The clinical application of these models will be beneficial for individuals harbouring these risk variants, and particularly for young individuals seeking genetic consultation.

Implication of specific retinal cell-type involvement and gene expression changes in AMD progression using integrative analysis of single-cell and bulk RNA-seq profiling

Age-related macular degeneration (AMD) is a blinding eye disease with no unifying theme for its etiology. We used single-cell RNA sequencing to analyze the transcriptomes of ~ 93,000 cells from the macula and peripheral retina from two adult human donors and bulk RNA sequencing from fifteen adult human donors with and without AMD. Analysis of our single-cell data identified 267 cell-type-specific genes. Comparison of macula and peripheral retinal regions found no cell-type differences but did identify 50 differentially expressed genes (DEGs) with about 1/3 expressed in cones. Integration of our single-cell data with bulk RNA sequencing data from normal and AMD donors showed compositional changes more pronounced in macula in rods, microglia, endothelium, Müller glia, and astrocytes in the transition from normal to advanced AMD. KEGG pathway analysis of our normal vs. advanced AMD eyes identified enrichment in complement and coagulation pathways, antigen presentation, tissue remodeling, and signaling pathways including PI3K-Akt, NOD-like, Toll-like, and Rap1. These results showcase the use of single-cell RNA sequencing to infer cell-type compositional and cell-type-specific gene expression changes in intact bulk tissue and provide a foundation for investigating molecular mechanisms of retinal disease that lead to new therapeutic targets.

A variant near DHCR24 associates with microstructural properties of white matter and peripheral lipid metabolism in adolescents

Visceral adiposity has been associated with altered microstructural properties of white matter in adolescents. Previous evidence suggests that circulating phospholipid PC(16:0/2:0) may mediate this association. To investigate the underlying biology, we performed a genome-wide association study (GWAS) of the shared variance of visceral fat, PC(16:0/2:0), and white matter microstructure in 872 adolescents from the Saguenay Youth Study. We further studied the metabolomic profile of the GWAS-lead variant in 931 adolescents. Visceral fat and white matter microstructure were assessed with magnetic resonance imaging. Circulating metabolites were quantified with serum lipidomics and metabolomics. We identified a genome-wide significant association near DHCR24 (Seladin-1) encoding a cholesterol-synthesizing enzyme (rs588709, p = 3.6 × 10-8); rs588709 was also associated nominally with each of the three traits (white matter microstructure: p = 2.1 × 10-6, PC(16:0/2:0): p = 0.005, visceral fat: p = 0.010). We found that the metabolic profile associated with rs588709 resembled that of a TM6SF2 variant impacting very low-density lipoprotein (VLDL) secretion and was only partially similar to that of a HMGCR variant. This suggests that the effect of rs588709 on VLDL lipids may arise due to altered phospholipid rather than cholesterol metabolism. The rs588709 was also nominally associated with circulating concentrations of omega-3 fatty acids in interaction with visceral fat and PC(16:0/2:0), and these fatty acid measures showed robust associations with white matter microstructure. Overall, the present study provides evidence that the DHCR24 locus may link peripheral metabolism to brain microstructure, an association with implications for cognitive impairment.

Transcriptome analysis of human dorsal striatum implicates attenuated canonical WNT signaling in neuroinflammation and in age-related impairment of striatal neurogenesis and synaptic plasticity

BACKGROUND

Motor and cognitive decline as part of the normal aging process is linked to alterations in synaptic plasticity and reduction of adult neurogenesis in the dorsal striatum. Neuroinflammation, particularly in the form of microglial activation, is suggested to contribute to these age-associated changes.

OBJECTIVE AND METHODS

To explore the molecular basis of alterations in striatal function during aging we analyzed RNA-Seq data for 117 postmortem human dorsal caudate samples and 97 putamen samples acquired through GTEx.

RESULTS

Increased expression of neuroinflammatory transcripts including TREM2, MHC II molecules HLA-DMB, HLA-DQA2, HLA-DPA1, HLA-DPB1, HLA-DMA and HLA-DRA, complement genes C1QA, C1QB, CIQC and C3AR1, and MHCI molecules HLA-B and HLA-F was identified. We also identified down-regulation of transcripts involved in neurogenesis, synaptogenesis, and synaptic pruning, including DCX, CX3CL1, and CD200, and the canonical WNTs WNT7A, WNT7B, and WNT8A. The canonical WNT signaling pathway has previously been shown to mediate adult neurogenesis and synapse formation and growth. Recent findings also highlight the link between WNT/β-catenin signaling and inflammation pathways.

CONCLUSIONS

These findings suggest that age-dependent attenuation of canonical WNT signaling plays a pivotal role in regulating striatal plasticity during aging. Dysregulation of WNT/β-catenin signaling via astrocyte-microglial interactions is suggested to be a novel mechanism that drives the decline of striatal neurogenesis and altered synaptic connectivity and plasticity, leading to a subsequent decrease in motor and cognitive performance with age. These findings may aid in the development of therapies targeting WNT/β-catenin signaling to combat cognitive and motor impairments associated with aging.

A pH-eQTL Interaction at the RIT2-SYT4 Parkinson's Disease Risk Locus in the Substantia Nigra

Parkinson's disease causes severe motor and cognitive disabilities that result from the progressive loss of dopamine neurons in the substantia nigra. The rs12456492 variant in the RIT2 gene has been repeatedly associated with increased risk for Parkinson's disease. From a transcriptomic perspective, a meta-analysis found that RIT2 gene expression is correlated with pH in the human brain. To assess these pH associations in relation to Parkinson's disease risk, we examined the two datasets that assayed rs12456492, gene expression, and pH in the postmortem human brain. Using the BrainEAC dataset, we replicate the positive correlation between RIT2 gene expression and pH in the human brain (n = 100). Furthermore, we found that the relationship between expression and pH is influenced by rs12456492. When tested across ten brain regions, this interaction is specifically found in the substantia nigra. A similar association was found for the co-localized SYT4 gene. In addition, SYT4 associations are stronger in a combined model with both genes, and the SYT4 interaction appears to be specific to males. In the Genotype-Tissue Expression (GTEx) dataset, the pH associations involving rs12456492 and expression of either SYT4 and RIT2 were not seen. This null finding may be due to the short postmortem intervals of the GTEx tissue samples. In the BrainEAC data, we tested the effect of postmortem interval and only observed the interactions in samples with the longer intervals. These previously unknown associations suggest novel roles for rs12456492, RIT2, and SYT4 in the regulation and response to pH in the substantia nigra.

Cholinergic receptor nicotinic beta 3 subunit polymorphisms and smoking in male Chinese patients with schizophrenia

OBJECTIVE

The aim of this study was to examine if cholinergic receptor nicotinic beta 3 subunit (CHRNB3) was a common genetic basis for both nicotine dependence and schizophrenia.

METHODS

Two CHRNB3 promoter single nucleotide polymorphisms (SNPs) were genotyped in 773 patients with schizophrenia and 302 healthy volunteers. Associations between smoking, schizophrenia, smoking+schizophrenia and CHRNB3 were analyzed. The mRNA expression of CHRNB3 in human brains was examined, and the expression correlations between CHRNB3 and dopaminergic and GABAergic receptor genes were evaluated.

RESULTS

The association between CHRNB3 and smoking was significant in the total sample, less significant in the smoking with schizophrenia, and suggestive in the smoking without schizophrenia. CHRNB3 had significant mRNA expression that was correlated with dopaminergic or GABAergic receptor expression in human brains. The two CHRNB3 SNPs had significant cis-acting regulatory effects on CHRNB3 mRNA expression.

CONCLUSIONS

Risk for smoking behavior was associated with CHRNB3. CHRNB3 mRNA is abundant in human brain and could play important role in the pathogenesis of smoking behavior.

Practical Euthanasia Method for Common Sea Stars (Asterias rubens) That Allows for High-Quality RNA Sampling

Sea stars in research are often lethally sampled without available methodology to render them insensible prior to sampling due to concerns over sufficient sample quality for applied molecular techniques. The objectives of this study were to describe an inexpensive and effective two-step euthanasia method for adult common sea stars (Asterias rubens) and to demonstrate that high-quality RNA samples for further use in downstream molecular analyses can be obtained from pyloric ceca of MgCl₂-immersed sea stars. Adult common sea stars (n = 15) were immersed in a 75 g/L magnesium chloride solution until they were no longer reactive to having their oral surface tapped with forceps (mean: 4 min, range 2-7 min), left immersed for an additional minute, and then sampled with sharp scissors. RNA from pyloric ceca (n = 10) was isolated using a liquid-liquid method, then samples were treated with DNase and analyzed for evaluation of RNA integrity number (RIN) for assessment of the quantity and purity of intact RNA. Aversive reactions to magnesium chloride solution were not observed and no sea stars regained spontaneous movement or reacted to sampling. The calculated RIN ranged from 7.3-9.8, demonstrating that the combination of animal welfare via the use of anesthesia and sampling for advanced molecular techniques is possible using this low-cost technique.

MIR-NATs repress MAPT translation and aid proteostasis in neurodegeneration

The human genome expresses thousands of natural antisense transcripts (NAT) that can regulate epigenetic state, transcription, RNA stability or translation of their overlapping genes1,2. Here we describe MAPT-AS1, a brain-enriched NAT that is conserved in primates and contains an embedded mammalian-wide interspersed repeat (MIR), which represses tau translation by competing for ribosomal RNA pairing with the MAPT mRNA internal ribosome entry site3. MAPT encodes tau, a neuronal intrinsically disordered protein (IDP) that stabilizes axonal microtubules. Hyperphosphorylated, aggregation-prone tau forms the hallmark inclusions of tauopathies4. Mutations in MAPT cause familial frontotemporal dementia, and common variations forming the MAPT H1 haplotype are a significant risk factor in many tauopathies5 and Parkinson's disease. Notably, expression of MAPT-AS1 or minimal essential sequences from MAPT-AS1 (including MIR) reduces-whereas silencing MAPT-AS1 expression increases-neuronal tau levels, and correlate with tau pathology in human brain. Moreover, we identified many additional NATs with embedded MIRs (MIR-NATs), which are overrepresented at coding genes linked to neurodegeneration and/or encoding IDPs, and confirmed MIR-NAT-mediated translational control of one such gene, PLCG1. These results demonstrate a key role for MAPT-AS1 in tauopathies and reveal a potentially broad contribution of MIR-NATs to the tightly controlled translation of IDPs6, with particular relevance for proteostasis in neurodegeneration.

Advances in protein-protein interaction network analysis for Parkinson's disease

Protein-protein interactions (PPIs) are a key component of the subcellular molecular networks which enable cells to function. Due to their importance in homeostasis, alterations to the networks can be detrimental, leading to cellular dysfunction and ultimately disease states. Parkinson's disease (PD) is a progressive neurodegenerative condition with multifactorial aetiology, spanning genetic variation and environmental modifiers. At a molecular and systems level, the characterisation of PD is the focus of extensive research, largely due to an unmet need for disease modifying therapies. PPI network analysis approaches are a valuable strategy to accelerate our understanding of the molecular crosstalk and biological processes underlying PD pathogenesis, especially due to the complex nature of this disease. In this review, we describe the utility of PPI network approaches in modelling complex systems, focusing on previous work in PD research. We discuss four principal strategies for using PPI network approaches: to infer PD related cellular functions, pathways and novel genes; to support genomics studies; to study the interactome of single PD related genes; and to compare the molecular basis of PD to other neurodegenerative disorders. This is an evolving area of research which is likely to further expand as omics data generation and availability increase. These approaches complement and bridge-the-gap between genetics and functional research to inform future investigations. In this review we outline several limitations that require consideration, acknowledging that ongoing challenges in this field continue to be addressed and the refinement of these approaches will facilitate further advances using PPI network analysis for understanding complex diseases.

The Spectrum of PRRT2-Associated Disorders: Update on Clinical Features and Pathophysiology

Mutations in the PRRT2 (proline-rich transmembrane protein 2) gene have been identified as the main cause of an expanding spectrum of disorders, including paroxysmal kinesigenic dyskinesia and benign familial infantile epilepsy, which places this gene at the border between epilepsy and movement disorders. The clinical spectrum has largely expanded to include episodic ataxia, hemiplegic migraine, and complex neurodevelopmental disorders in cases with biallelic mutations. Prior to the discovery of PRRT2 as the causative gene for this spectrum of disorders, the sensitivity of paroxysmal kinesigenic dyskinesia to anticonvulsant drugs regulating ion channel function as well as the co-occurrence of epilepsy in some patients or families fostered the hypothesis this could represent a channelopathy. However, recent evidence implicates PRRT2 in synapse functioning, which disproves the "channel hypothesis" (although PRRT2 modulates ion channels at the presynaptic level), and justifies the classification of these conditions as synaptopathies, an emerging rubric of brain disorders. This review aims to provide an update of the clinical and pathophysiologic features of PRRT2-associated disorders.

Type 2 diabetes mellitus-associated transcriptome alterations in cortical neurones and associated neurovascular unit cells in the ageing brain

Type 2 diabetes mellitus (T2D), characterised by peripheral insulin resistance, is a risk factor for dementia. In addition to its contribution to small and large vessel disease, T2D may directly damage cells of the brain neurovascular unit. In this study, we investigated the transcriptomic changes in cortical neurones, and associated astrocytes and endothelial cells of the neurovascular unit, in the ageing brain. Neurone, astrocyte, and endothelial cell-enriched mRNA, obtained by immuno-laser capture microdissection of temporal cortex (Brodmann area 21/22) from 6 cases with self-reported T2D in the Cognitive Function and Ageing Study neuropathology cohort, and an equal number of age and sex-matched controls, was assessed by microarray analysis. Integrated Molecular Pathway Level Analysis was performed using the Kyoto Encyclopaedia of Genes and Genomes database on significantly differentially expressed genes, defined as P < 0.05 and fold-change ± 1.2. Hub genes identified from Weighted Gene Co-expression Network Analysis were validated in neurones using the NanoString nCounter platform. The expression and cellular localisation of proteins encoded by selected candidate genes were confirmed by immunohistochemistry. 912, 2202, and 1227 genes were significantly differentially expressed between cases with self-reported T2D and controls in neurones, astrocytes, and endothelial cells respectively. Changes in cortical neurones included alterations in insulin and other signalling pathways, cell cycle, cellular senescence, inflammatory mediators, and components of the mitochondrial respiratory electron transport chain. Impaired insulin signalling was shared by neurovascular unit cells with, additionally, apoptotic pathway changes in astrocytes and dysregulation of advanced glycation end-product signalling in endothelial cells. Transcriptomic analysis identified changes in key cellular pathways associated with T2D that may contribute to neuronal damage and dysfunction. These effects on brain cells potentially contribute to a diabetic dementia, and may provide novel approaches for therapeutic intervention.

Cellular correlates of cortical thinning throughout the lifespan

Cortical thinning occurs throughout the entire life and extends to late-life neurodegeneration, yet the neurobiological substrates are poorly understood. Here, we used a virtual-histology technique and gene expression data from the Allen Human Brain Atlas to compare the regional profiles of longitudinal cortical thinning through life (4004 magnetic resonance images [MRIs]) with those of gene expression for several neuronal and non-neuronal cell types. The results were replicated in three independent datasets. We found that inter-regional profiles of cortical thinning related to expression profiles for marker genes of CA1 pyramidal cells, astrocytes and, microglia during development and in aging. During the two stages of life, the relationships went in opposite directions: greater gene expression related to less thinning in development and vice versa in aging. The association between cortical thinning and cell-specific gene expression was also present in mild cognitive impairment and Alzheimer's Disease. These findings suggest a role of astrocytes and microglia in promoting and supporting neuronal growth and dendritic structures through life that affects cortical thickness during development, aging, and neurodegeneration. Overall, the findings contribute to our understanding of the neurobiology underlying variations in MRI-derived estimates of cortical thinning through life and late-life disease.

Parkinson's Disease Master Regulators on Substantia Nigra and Frontal Cortex and Their Use for Drug Repositioning

Parkinson's disease (PD) is among the most prevalent neurodegenerative diseases. Available evidences support the view of PD as a complex disease, being the outcome of interactions between genetic and environmental factors. In face of diagnosis and therapy challenges, and the elusive PD etiology, the use of alternative methodological approaches for the elucidation of the disease pathophysiological mechanisms and proposal of novel potential therapeutic interventions has become increasingly necessary. In the present study, we first reconstructed the transcriptional regulatory networks (TN), centered on transcription factors (TF), of two brain regions affected in PD, the substantia nigra pars compacta (SNc) and the frontal cortex (FCtx). Then, we used case-control studies data from these regions to identify TFs working as master regulators (MR) of the disease, based on region-specific TNs. Twenty-nine regulatory units enriched with differentially expressed genes were identified for the SNc, and twenty for the FCtx, all of which were considered MR candidates for PD. Three consensus MR candidates were found for SNc and FCtx, namely ATF2, SLC30A9, and ZFP69B. In order to search for novel potential therapeutic interventions, we used these consensus MR candidate signatures as input to the Connectivity Map (CMap), a computational drug repositioning webtool. This analysis resulted in the identification of four drugs that reverse the expression pattern of all three MR consensus simultaneously, benperidol, harmaline, tubocurarine chloride, and vorinostat, thus suggested as novel potential PD therapeutic interventions.

Replicated risk CACNA1C variants for major psychiatric disorders may serve as potential therapeutic targets for the shared depressive endophenotype

Genome-wide association studies (GWASs) have reported numerous associations between risk variants and major psychiatric disorders (MPDs) including schizophrenia (SCZ), bipolar disorder (BPD), major depressive disorder (MDD) and others. We reviewed all of the published GWASs, and extracted the genome-wide significant (p<10^-6) and replicated associations between risk SNPs and MPDs. We found the associations of 6 variants located in 6 genes, including L type voltage-gated calcium channel (LTCCs) subunit alpha1 C gene (CACNA1C), that were genome-wide significant (2.0×10 ^-8 ≤p≤1.0×10 ^-6 ) and replicated at single-point level across at least two GWASs. Among them, the associations between MPDs and rs1006737 within CACNA1C are most robust. Thus, as a next step, the expression of the replicated risk genes in human hippocampus was analyzed. We found CACNA1C had significant mRNA expression in human hippocampus in two independent cohorts. Finally, we tried to elucidate the roles of venlafaxine and ω-3 PUFAs in the mRNA expression regulation of the replicated risk genes in hippocampus. We used cDNA chip-based microarray profiling to explore the transcriptome-wide mRNA expression regulation by ω-3 PUFAs (0.72/kg/d) and venlafaxine (0.25/kg/d) treatment in chronic mild stress (CMS) rats. ω-3 PUFAs and venlafaxine treatment elicited significant CACNA1C up-regulation. We concluded that CACNA1C might confer the genetic vulnerability to the shared depressive symptoms across MPDs and CACNA1C might be the therapeutic target for depressive endophenotype as well.

Meta-analysis of up to 622,409 individuals identifies 40 novel smoking behaviour associated genetic loci

Smoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been associated with smoking behaviour-related traits. We tested up to 235,116 single nucleotide variants (SNVs) on the exome-array for association with smoking initiation, cigarettes per day, pack-years, and smoking cessation in a fixed effects meta-analysis of up to 61 studies (up to 346,813 participants). In a subset of 112,811 participants, a further one million SNVs were also genotyped and tested for association with the four smoking behaviour traits. SNV-trait associations with P < 5 × 10-8 in either analysis were taken forward for replication in up to 275,596 independent participants from UK Biobank. Lastly, a meta-analysis of the discovery and replication studies was performed. Sixteen SNVs were associated with at least one of the smoking behaviour traits (P < 5 × 10-8) in the discovery samples. Ten novel SNVs, including rs12616219 near TMEM182, were followed-up and five of them (rs462779 in REV3L, rs12780116 in CNNM2, rs1190736 in GPR101, rs11539157 in PJA1, and rs12616219 near TMEM182) replicated at a Bonferroni significance threshold (P < 4.5 × 10-3) with consistent direction of effect. A further 35 SNVs were associated with smoking behaviour traits in the discovery plus replication meta-analysis (up to 622,409 participants) including a rare SNV, rs150493199, in CCDC141 and two low-frequency SNVs in CEP350 and HDGFRP2. Functional follow-up implied that decreased expression of REV3L may lower the probability of smoking initiation. The novel loci will facilitate understanding the genetic aetiology of smoking behaviour and may lead to the identification of potential drug targets for smoking prevention and/or cessation.

Comparison Between Expression Microarrays and RNA-Sequencing Using UKBEC Dataset Identified a trans-eQTL Associated with MPZ Gene in Substantia Nigra

In recent years, the advantages of RNA-sequencing (RNA-Seq) have made it the platform of choice for measuring gene expression over traditional microarrays. However, RNA-Seq comes with bioinformatical challenges and higher computational costs. Therefore, this study set out to assess whether the increased depth of transcriptomic information facilitated by RNA-Seq is worth the increased computation over microarrays, specifically at three levels: absolute expression levels, differentially expressed genes identification, and expression QTL (eQTL) mapping in regions of the human brain. Using the United Kingdom Brain Expression Consortium (UKBEC) dataset, there is high agreement of gene expression levels measured by microarrays and RNA-seq when quantifying absolute expression levels and when identifying differentially expressed genes. These findings suggest that depending on the aims of a study, the relative ease of working with microarray data may outweigh the computational time and costs of RNA-Seq pipelines. On the other, there was low agreement when mapping eQTLs. However, a number of eQTLs associated with genes that play important roles in the brain were found in both platforms. For example, a trans-eQTL was mapped that is associated with the MPZ gene in the substantia nigra. These eQTLs that we have highlighted are extremely promising candidates that merit further investigation.

Alzheimer's disease-related dysregulation of mRNA translation causes key pathological features with ageing

Alzheimer's disease (AD) is characterised by Aβ and tau pathology as well as synaptic degeneration, which correlates best with cognitive impairment. Previous work suggested that this pathological complexity may result from changes in mRNA translation. Here, we studied whether mRNA translation and its underlying signalling are altered in an early model of AD, and whether modelling this deficiency in mice causes pathological features with ageing. Using an unbiased screen, we show that exposure of primary neurons to nanomolar amounts of Aβ increases FMRP-regulated protein synthesis. This selective regulation of mRNA translation is dependent on a signalling cascade involving MAPK-interacting kinase 1 (Mnk1) and the eukaryotic initiation factor 4E (eIF4E), and ultimately results in reduction of CYFIP2, an FMRP-binding protein. Modelling this CYFIP2 reduction in mice, we find age-dependent Aβ accumulation in the thalamus, development of tau pathology in entorhinal cortex and hippocampus, as well as gliosis and synapse loss in the hippocampus, together with deficits in memory formation. Therefore, we conclude that early stages of AD involve increased translation of specific CYFIP2/FMRP-regulated transcripts. Since reducing endogenous CYFIP2 expression is sufficient to cause key features of AD with ageing in mice, we suggest that prolonged activation of this pathway is a primary step toward AD pathology, highlighting a novel direction for therapeutic targeting.

Genome-wide association study of cerebral small vessel disease reveals established and novel loci

Intracerebral haemorrhage and small vessel ischaemic stroke (SVS) are the most acute manifestations of cerebral small vessel disease, with no established preventive approaches beyond hypertension management. Combined genome-wide association study (GWAS) of these two correlated diseases may improve statistical power to detect novel genetic factors for cerebral small vessel disease, elucidating underlying disease mechanisms that may form the basis for future treatments. Because intracerebral haemorrhage location is an adequate surrogate for distinct histopathological variants of cerebral small vessel disease (lobar for cerebral amyloid angiopathy and non-lobar for arteriolosclerosis), we performed GWAS of intracerebral haemorrhage by location in 1813 subjects (755 lobar and 1005 non-lobar) and 1711 stroke-free control subjects. Intracerebral haemorrhage GWAS results by location were meta-analysed with GWAS results for SVS from MEGASTROKE, using 'Multi-Trait Analysis of GWAS' (MTAG) to integrate summary data across traits and generate combined effect estimates. After combining intracerebral haemorrhage and SVS datasets, our sample size included 241 024 participants (6255 intracerebral haemorrhage or SVS cases and 233 058 control subjects). Genome-wide significant associations were observed for non-lobar intracerebral haemorrhage enhanced by SVS with rs2758605 [MTAG P-value (P) = 2.6 × 10-8] at 1q22; rs72932727 (P = 1.7 × 10-8) at 2q33; and rs9515201 (P = 5.3 × 10-10) at 13q34. In the GTEx gene expression library, rs2758605 (1q22), rs72932727 (2q33) and rs9515201 (13q34) are significant cis-eQTLs for PMF1 (P = 1 × 10-4 in tibial nerve), NBEAL1, FAM117B and CARF (P < 2.1 × 10-7 in arteries) and COL4A2 and COL4A1 (P < 0.01 in brain putamen), respectively. Leveraging S-PrediXcan for gene-based association testing with the predicted expression models in tissues related with nerve, artery, and non-lobar brain, we found that experiment-wide significant (P < 8.5 × 10-7) associations at three genes at 2q33 including NBEAL1, FAM117B and WDR12 and genome-wide significant associations at two genes including ICA1L at 2q33 and ZCCHC14 at 16q24. Brain cell-type specific expression profiling libraries reveal that SEMA4A, SLC25A44 and PMF1 at 1q22 and COL4A1 and COL4A2 at 13q34 were mainly expressed in endothelial cells, while the genes at 2q33 (FAM117B, CARF and NBEAL1) were expressed in various cell types including astrocytes, oligodendrocytes and neurons. Our cross-phenotype genetic study of intracerebral haemorrhage and SVS demonstrates novel genome-wide associations for non-lobar intracerebral haemorrhage at 2q33 and 13q34. Our replication of the 1q22 locus previous seen in traditional GWAS of intracerebral haemorrhage, as well as the rediscovery of 13q34, which had previously been reported in candidate gene studies with other cerebral small vessel disease-related traits strengthens the credibility of applying this novel genome-wide approach across intracerebral haemorrhage and SVS.

Genetic Variation of a DRD2 Co-expression Network is Associated with Changes in Prefrontal Function After D2 Receptors Stimulation

Dopamine D2 receptors (D2Rs) contribute to the inverted U-shaped relationship between dopamine signaling and prefrontal function. Genetic networks from post-mortem human brain revealed 84 partner genes co-expressed with DRD2. Moreover, eight functional single nucleotide polymorphisms combined into a polygenic co-expression index (PCI) predicted co-expression of this DRD2 network and were associated with prefrontal function in humans. Here, we investigated the non-linear association of the PCI with behavioral and Working Memory (WM) related brain response to pharmacological D2Rs stimulation. Fifty healthy volunteers took part in a double-blind, placebo-controlled, functional MRI (fMRI) study with bromocriptine and performed the N-Back task. The PCI by drug interaction was significant on both WM behavioral scores (P = 0.046) and related prefrontal activity (all corrected P < 0.05) using a polynomial PCI model. Non-linear responses under placebo were reversed by bromocriptine administration. fMRI results on placebo were replicated in an independent sample of 50 participants who did not receive drug administration (P = 0.034). These results match earlier evidence in non-human primates and confirm the physiological relevance of this DRD2 co-expression network. Results show that in healthy subjects, different alleles evaluated as an ensemble are associated with non-linear prefrontal responses. Therefore, brain response to a dopaminergic drug may depend on a complex system of allelic patterns associated with DRD2 co-expression.

Investigating the Causal Effect of Brain Expression of CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10 Genes on Multiple Sclerosis: A Two-Sample Mendelian Randomization Approach

Multiple Sclerosis (MS) exhibits considerable heterogeneity in phenotypic expression, course, prognosis and response to therapy. This suggests this disease involves multiple, as yet poorly understood, causal mechanisms. In this work we assessed the possible causal link between gene expression level of five selected genes related to the pro-inflammatory NF-κB signaling pathway (i.e., CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10) in ten different brain tissues (i.e., cerebellum, frontal cortex, hippocampus, medulla, occipital cortex, putamen, substantia nigra, thalamus, temporal cortex and intralobular white matter) and MS. We adopted a two-stage Mendelian Randomization (MR) approach for the estimation of the causal effects of interest, based on summary-level data from 20 multiplex Sardinian families and data provided by the United Kingdom Brain Expression Consortium (UKBEC). Through Radial-MR and Cochrane's Q statistics we identified and removed genetic variants which are most likely to be invalid instruments. To estimate the total causal effect, univariable MR was carried out separately for each gene and brain region. We used Inverse-Variance Weighted estimator (IVW) as main analysis and MR-Egger Regression estimator (MR-ER) and Weighted Median Estimator (WME) as sensitivity analysis. As these genes belong to the same pathway and thus they can be closely related, we also estimated their direct causal effects by applying IVW and MR-ER within a multivariable MR (MVMR) approach using set of genetic instruments specific and common (composite) to each multiple exposures represented by the expression of the candidate genes. Univariate MR analysis showed a significant positive total causal effect for CCL2 and NFKB1 respectively in medulla and cerebellum. MVMR showed a direct positive causal effect for NFKB1 and TNFRSF1A, and a direct negative causal effect for CCL2 in cerebellum; while in medulla we observed a direct positive causal effect for CCL2. Since in general we observed a different magnitude for the gene specific causal effect we hypothesize that in cerebellum and medulla the effect of each gene expression is direct but also mediated by the others. These results confirm the importance of the involvement of NF-κB signaling pathway in brain tissue for the development of the disease and improve our understanding in the pathogenesis of MS.

Integrative analyses indicate an association between ITIH3 polymorphisms with autism spectrum disorder

It is challenge to pinpoint the functional variants among numerous genetic variants. Investigating the spatial dynamics of the human brain transcriptome for genes and exploring the expression quantitative trait loci data may provide the potential direction to identify the functional variants among autism spectrum disorders (ASD) patients. In order to explore the association of ITIH3 with ASD, the present study included three components: identifying the spatial-temporal expression of ITIH3 in the developing human brain using the expression data from the Allen Institute for Brain Science; examining the cis-acting regulatory effect of SNPs on the ITIH3 expression using UK Brain Expression Consortium database; validating the effect of identified SNPs using a case-control study with samples of 602 cases and 604 controls. The public expression data showed that ITIH3 may have a role in the development of human brain and suggested a cis-eQTL effect for rs2535629 and rs3617 on ITIH3 in the hippocampus. Genetic analysis of the above two SNPs suggested that the over-dominant model of rs2535629 was significantly associated with decreased risk of ASD. Convergent lines of evidence supported ITIH3 rs25352629 as a susceptibility variant for ASD.

Corticosteroids and Regional Variations in Thickness of the Human Cerebral Cortex across the Lifespan

Exposures to life stressors accumulate across the lifespan, with possible impact on brain health. Little is known, however, about the mechanisms mediating age-related changes in brain structure. We use a lifespan sample of participants (n = 21 251; 4-97 years) to investigate the relationship between the thickness of cerebral cortex and the expression of the glucocorticoid- and the mineralocorticoid-receptor genes (NR3C1 and NR3C2, respectively), obtained from the Allen Human Brain Atlas. In all participants, cortical thickness correlated negatively with the expression of both NR3C1 and NR3C2 across 34 cortical regions. The magnitude of this correlation varied across the lifespan. From childhood through early adulthood, the profile similarity (between NR3C1/NR3C2 expression and thickness) increased with age. Conversely, both profile similarities decreased with age in late life. These variations do not reflect age-related changes in NR3C1 and NR3C2 expression, as observed in 5 databases of gene expression in the human cerebral cortex (502 donors). Based on the co-expression of NR3C1 (and NR3C2) with genes specific to neural cell types, we determine the potential involvement of microglia, astrocytes, and CA1 pyramidal cells in mediating the relationship between corticosteroid exposure and cortical thickness. Therefore, corticosteroids may influence brain structure to a variable degree throughout life.

Transcriptomic signatures of brain regional vulnerability to Parkinson's disease

The molecular mechanisms underlying caudal-to-rostral progression of Lewy body pathology in Parkinson's disease remain poorly understood. Here, we identified transcriptomic signatures across brain regions involved in Braak Lewy body stages in non-neurological adults from the Allen Human Brain Atlas. Among the genes that are indicative of regional vulnerability, we found known genetic risk factors for Parkinson's disease: SCARB2, ELOVL7, SH3GL2, SNCA, BAP1, and ZNF184. Results were confirmed in two datasets of non-neurological subjects, while in two datasets of Parkinson's disease patients we found altered expression patterns. Co-expression analysis across vulnerable regions identified a module enriched for genes associated with dopamine synthesis and microglia, and another module related to the immune system, blood-oxygen transport, and endothelial cells. Both were highly expressed in regions involved in the preclinical stages of the disease. Finally, alterations in genes underlying these region-specific functions may contribute to the selective regional vulnerability in Parkinson's disease brains.

Targeted sequencing of the LRRTM gene family in suicide attempters with bipolar disorder

Glutamatergic signaling is the primary excitatory neurotransmission pathway in the brain, and its relationship to neuropsychiatric disorders is of considerable interest. Our previous attempted suicide genome-wide association study, and numerous studies investigating gene expression, genetic variation, and DNA methylation have implicated aberrant glutamatergic signaling in suicide risk. The glutamatergic pathway gene LRRTM4 was an associated gene identified in our attempted suicide genome-wide association study, with association support seen primarily in females. Recent evidence has also shown that glutamatergic signaling is partly regulated by sex-related hormones. The LRRTM gene family encodes neuronal leucine-rich transmembrane proteins that localize to and promote glutamatergic synapse development. In this study, we sequenced the coding and regulatory regions of all four LRRTM gene members plus a large intronic region of LRRTM4 in 476 bipolar disorder suicide attempters and 473 bipolar disorder nonattempters. We identified two male-specific variants, one female- and five male-specific haplotypes significantly associated with attempted suicide in LRRTM4. Furthermore, variants within significant haplotypes may be brain expression quantitative trait loci for LRRTM4 and some of these variants overlap with predicted hormone response elements. Overall, these results provide supporting evidence for a sex-specific association of genetic variation in LRRTM4 with attempted suicide.

Male-specific association of the 2p25 region with suicide attempt in bipolar disorder

We previously conducted a genome-wide association study (GWAS) of attempted suicide within bipolar disorder, which implicated common variation in the 2p25 region primarily in males. The top association signal from our GWAS occurred in an intergenic region of 2p25 (p = 5.07 × 10^-8) and was supported by two independent studies. In the current study, to better characterize the association of the 2p25 region with attempted suicide, we sequenced the entire 350kb 2p25 region in 476 bipolar suicide attempters and 473 bipolar non-attempters using targeted next-generation sequencing. This fine-mapping project identified 4,681 variants in the 2p25 region. We performed both gene-level and individual-variant tests on our sequencing results and identified 375 variants which were nominally significant (p < 0.05) and three common variants that were significantly associated with attempted suicide in males (corrected p = 0.035, odds ratio (OR) = 2.13). These three variants are in strong linkage disequilibrium with the top variant from our GWAS. Our top five variants are also predicted expression quantitative trait loci (eQTL) for three genes in the 2p25 region based on publicly available brain expression databases. Our sequencing and eQTL data implicate these three genes - SH3YL1, ACP1, and FAM150B - and three additional pathways - androgen receptor, Wnt signaling, and glutamatergic/GABAergic signaling - in the association of the 2p25 region with suicide. The current study provides additional support for an association of the 2p25 region with attempted suicide in males and identifies several candidate genes and pathways that warrant further investigation to understand their role in suicidal behavior.

Genome-wide human brain eQTLs: In-depth analysis and insights using the UKBEC dataset

Understanding the complexity of the human brain transcriptome architecture is one of the most important human genetics study areas. Previous studies have applied expression quantitative trait loci (eQTL) analysis at the genome-wide level of the brain to understand the underlying mechanisms relating to neurodegenerative diseases, primarily at the transcript level. To increase the resolution of our understanding, the current study investigates multi/single-region, transcript/exon-level and cis versus trans-acting eQTL, across 10 regions of the human brain. Some of the key findings of this study are: (i) only a relatively small proportion of eQTLs will be detected, where the sensitivity is under 5%; (ii) when an eQTL is acting in multiple regions (MR-eQTL), it tends to have very similar effects on gene expression in each of these regions, as well as being cis-acting; (iii) trans-acting eQTLs tend to have larger effects on expression compared to cis-acting eQTLs and tend to be specific to a single region (SR-eQTL) of the brain; (iv) the cerebellum has a very large number of eQTLs that function exclusively in this region, compared with other regions of the brain; (v) importantly, an interactive visualisation tool (Shiny app) was developed to visualise the MR/SR-eQTL at transcript and exon levels.

Association of Schizophrenia Risk With Disordered Niacin Metabolism in an Indian Genome-wide Association Study

Importance

Genome-wide association studies (GWASs) in European populations have identified more than 100 schizophrenia-associated loci. A schizophrenia GWAS in a unique Indian population offers novel findings.

Objective

To discover and functionally evaluate genetic loci for schizophrenia in a GWAS of a unique Indian population.

Design, Setting, and Participants

This GWAS included a sample of affected individuals, family members, and unrelated cases and controls. Three thousand ninety-two individuals were recruited and diagnostically ascertained via medical records, hospitals, clinics, and clinical networks in Chennai and surrounding regions. Affected participants fulfilled DSM-IV diagnostic criteria for schizophrenia. Unrelated control participants had no personal or family history of psychotic disorder. Recruitment, genotyping, and analysis occurred in consecutive phases beginning January 1, 2001. Recruitment was completed on February 28, 2018, and genotyping and analysis are ongoing.

Main Outcomes and Measures

Associations of single-nucleotide polymorphisms and gene expression with schizophrenia.

Results

The study population included 1321 participants with schizophrenia, 885 family controls, and 886 unrelated controls. Among participants with schizophrenia, mean (SD) age was 39.1 (11.4) years, and 52.7% were male. This sample demonstrated uniform ethnicity, a degree of inbreeding, and negligible rates of substance abuse. A novel genome-wide significant association was observed between schizophrenia and a chromosome 8q24.3 locus (rs10866912, allele A; odds ratio [OR], 1.27 [95% CI, 1.17-1.38]; P = 4.35 × 10-8) that attracted support in the schizophrenia Psychiatric Genomics Consortium 2 data (rs10866912, allele A; OR, 1.04 [95% CI, 1.02-1.06]; P = 7.56 × 10-4). This locus has undergone natural selection, with the risk allele A declining in frequency from India (approximately 72%) to Europe (approximately 43%). rs10866912 directly modifies the abundance of the nicotinate phosphoribosyltransferase gene (NAPRT1) transcript in brain cortex (normalized effect size, 0.79; 95% CI, 0.6-1.0; P = 5.8 × 10-13). NAPRT1 encodes a key enzyme for niacin metabolism. In Indian lymphoblastoid cell lines, (risk) allele A of rs10866912 was associated with NAPRT1 downregulation (AA: 0.74, n = 21; CC: 1.56, n = 17; P = .004). Preliminary zebrafish data further suggest that partial loss of function of NAPRT1 leads to abnormal brain development.

Conclusions and Relevance

Bioinformatic analyses and cellular and zebrafish gene expression studies implicate NAPRT1 as a novel susceptibility gene. Given this gene's role in niacin metabolism and the evidence for niacin deficiency provoking schizophrenialike symptoms in neuropsychiatric diseases such as pellagra and Hartnup disease, these results suggest that the rs10866912 genotype and niacin status may have implications for schizophrenia susceptibility and treatment.