Rare variant associations with plasma protein levels in the UK Biobank.

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

Functional evaluation of epilepsy-associated KCNT1 variants in multiple cellular systems reveals a predominant gain of function impact on channel properties.

OBJECTIVE: Gain of function variants in the sodium-activated potassium channel KCNT1 have been associated with pediatric epilepsy disorders. Here, we systematically examine a spectrum of KCNT1 variants and establish their impact on channel function in multiple cellular systems. METHODS: KCNT1 variants identified from published reports and genetic screening of pediatric epilepsy patients were expressed in Xenopus oocytes and HEK cell lines. Variant impact on current magnitude, current-voltage relationships, and sodium ion modulation were examined. RESULTS: We determined basic properties of KCNT1 in Xenopus oocyte and HEK systems, including the role of extra- and intracellular sodium in regulating KCNT1 activity. The most common six KCNT1 variants demonstrated strong gain of function (GOF) effects on one or more channel properties. Analysis of 36 total variants identified phenotypic heterogeneity but a strong tendency for pathogenic variants to exert GOF effects on channel properties. By controlling intracellular sodium, we demonstrate that multiple pathogenic KCNT1 variants modulate channel voltage dependence by altering the sensitivity to sodium ions. SIGNIFICANCE: This study represents the largest systematic functional examination of KCNT1 variants to date. We both confirm previously reported GOF channel phenotypes and expand the number of variants with in vitro GOF effects. Our data provide further evidence that novel KCNT1 variants identified in epilepsy patients lead to disease through generalizable GOF mechanisms including increases in current magnitude and/or current-voltage relationships.

FreeSurfer-based segmentation of hippocampal subfields: A review of methods and applications, with a novel quality control procedure for ENIGMA studies and other collaborative efforts.

Structural hippocampal abnormalities are common in many neurological and psychiatric disorders, and variation in hippocampal measures is related to cognitive performance and other complex phenotypes such as stress sensitivity. Hippocampal subregions are increasingly studied, as automated algorithms have become available for mapping and volume quantification. In the context of the Enhancing Neuro Imaging Genetics through Meta Analysis Consortium, several Disease Working Groups are using the FreeSurfer software to analyze hippocampal subregion (subfield) volumes in patients with neurological and psychiatric conditions along with data from matched controls. In this overview, we explain the algorithm's principles, summarize measurement reliability studies, and demonstrate two additional aspects (subfield autocorrelation and volume/reliability correlation) with illustrative data. We then explain the rationale for a standardized hippocampal subfield segmentation quality control (QC) procedure for improved pipeline harmonization. To guide researchers to make optimal use of the algorithm, we discuss how global size and age effects can be modeled, how QC steps can be incorporated and how subfields may be aggregated into composite volumes. This discussion is based on a synopsis of 162 published neuroimaging studies (01/2013-12/2019) that applied the FreeSurfer hippocampal subfield segmentation in a broad range of domains including cognition and healthy aging, brain development and neurodegeneration, affective disorders, psychosis, stress regulation, neurotoxicity, epilepsy, inflammatory disease, childhood adversity and posttraumatic stress disorder, and candidate and whole genome (epi-)genetics. Finally, we highlight points where FreeSurfer-based hippocampal subfield studies may be optimized.

SCN1A overexpression, associated with a genomic region marked by a risk variant for a common epilepsy, raises seizure susceptibility.

Mesial temporal lobe epilepsy with hippocampal sclerosis and a history of febrile seizures is associated with common variation at rs7587026, located in the promoter region of SCN1A. We sought to explore possible underlying mechanisms. SCN1A expression was analysed in hippocampal biopsy specimens of individuals with mesial temporal lobe epilepsy with hippocampal sclerosis who underwent surgical treatment, and hippocampal neuronal cell loss was quantitatively assessed using immunohistochemistry. In healthy individuals, hippocampal volume was measured using MRI. Analyses were performed stratified by rs7587026 type. To study the functional consequences of increased SCN1A expression, we generated, using transposon-mediated bacterial artificial chromosome transgenesis, a zebrafish line expressing exogenous scn1a, and performed EEG analysis on larval optic tecta at 4 day post-fertilization. Finally, we used an in vitro promoter analysis to study whether the genetic motif containing rs7587026 influences promoter activity. Hippocampal SCN1A expression differed by rs7587026 genotype (Kruskal-Wallis test P = 0.004). Individuals homozygous for the minor allele showed significantly increased expression compared to those homozygous for the major allele (Dunn's test P = 0.003), and to heterozygotes (Dunn's test P = 0.035). No statistically significant differences in hippocampal neuronal cell loss were observed between the three genotypes. Among 597 healthy participants, individuals homozygous for the minor allele at rs7587026 displayed significantly reduced mean hippocampal volume compared to major allele homozygotes (Cohen's D = - 0.28, P = 0.02), and to heterozygotes (Cohen's D = - 0.36, P = 0.009). Compared to wild type, scn1lab-overexpressing zebrafish larvae exhibited more frequent spontaneous seizures [one-way ANOVA F(4,54) = 6.95 (P < 0.001)]. The number of EEG discharges correlated with the level of scn1lab overexpression [one-way ANOVA F(4,15) = 10.75 (P < 0.001]. Finally, we showed that a 50 bp promoter motif containing rs7587026 exerts a strong regulatory role on SCN1A expression, though we could not directly link this to rs7587026 itself. Our results develop the mechanistic link between rs7587026 and mesial temporal lobe epilepsy with hippocampal sclerosis and a history of febrile seizures. Furthermore, we propose that quantitative precision may be important when increasing SCN1A expression in current strategies aiming to treat seizures in conditions involving SCN1A haploinsufficiency, such as Dravet syndrome.

The ENIGMA-Epilepsy working group: Mapping disease from large data sets.

Epilepsy is a common and serious neurological disorder, with many different constituent conditions characterized by their electro clinical, imaging, and genetic features. MRI has been fundamental in advancing our understanding of brain processes in the epilepsies. Smaller-scale studies have identified many interesting imaging phenomena, with implications both for understanding pathophysiology and improving clinical care. Through the infrastructure and concepts now well-established by the ENIGMA Consortium, ENIGMA-Epilepsy was established to strengthen epilepsy neuroscience by greatly increasing sample sizes, leveraging ideas and methods established in other ENIGMA projects, and generating a body of collaborating scientists and clinicians to drive forward robust research. Here we review published, current, and future projects, that include structural MRI, diffusion tensor imaging (DTI), and resting state functional MRI (rsfMRI), and that employ advanced methods including structural covariance, and event-based modeling analysis. We explore age of onset- and duration-related features, as well as phenomena-specific work focusing on particular epilepsy syndromes or phenotypes, multimodal analyses focused on understanding the biology of disease progression, and deep learning approaches. We encourage groups who may be interested in participating to make contact to further grow and develop ENIGMA-Epilepsy.

Structural brain abnormalities in the common epilepsies assessed in a worldwide ENIGMA study.

Progressive functional decline in the epilepsies is largely unexplained. We formed the ENIGMA-Epilepsy consortium to understand factors that influence brain measures in epilepsy, pooling data from 24 research centres in 14 countries across Europe, North and South America, Asia, and Australia. Structural brain measures were extracted from MRI brain scans across 2149 individuals with epilepsy, divided into four epilepsy subgroups including idiopathic generalized epilepsies (n =367), mesial temporal lobe epilepsies with hippocampal sclerosis (MTLE; left, n = 415; right, n = 339), and all other epilepsies in aggregate (n = 1026), and compared to 1727 matched healthy controls. We ranked brain structures in order of greatest differences between patients and controls, by meta-analysing effect sizes across 16 subcortical and 68 cortical brain regions. We also tested effects of duration of disease, age at onset, and age-by-diagnosis interactions on structural measures. We observed widespread patterns of altered subcortical volume and reduced cortical grey matter thickness. Compared to controls, all epilepsy groups showed lower volume in the right thalamus (Cohen's d = -0.24 to -0.73; P < 1.49 × 10-4), and lower thickness in the precentral gyri bilaterally (d = -0.34 to -0.52; P < 4.31 × 10-6). Both MTLE subgroups showed profound volume reduction in the ipsilateral hippocampus (d = -1.73 to -1.91, P < 1.4 × 10-19), and lower thickness in extrahippocampal cortical regions, including the precentral and paracentral gyri, compared to controls (d = -0.36 to -0.52; P < 1.49 × 10-4). Thickness differences of the ipsilateral temporopolar, parahippocampal, entorhinal, and fusiform gyri, contralateral pars triangularis, and bilateral precuneus, superior frontal and caudal middle frontal gyri were observed in left, but not right, MTLE (d = -0.29 to -0.54; P < 1.49 × 10-4). Contrastingly, thickness differences of the ipsilateral pars opercularis, and contralateral transverse temporal gyrus, were observed in right, but not left, MTLE (d = -0.27 to -0.51; P < 1.49 × 10-4). Lower subcortical volume and cortical thickness associated with a longer duration of epilepsy in the all-epilepsies, all-other-epilepsies, and right MTLE groups (beta, b < -0.0018; P < 1.49 × 10-4). In the largest neuroimaging study of epilepsy to date, we provide information on the common epilepsies that could not be realistically acquired in any other way. Our study provides a robust ranking of brain measures that can be further targeted for study in genetic and neuropathological studies. This worldwide initiative identifies patterns of shared grey matter reduction across epilepsy syndromes, and distinctive abnormalities between epilepsy syndromes, which inform our understanding of epilepsy as a network disorder, and indicate that certain epilepsy syndromes involve more widespread structural compromise than previously assumed.

Genome-wide association analysis links multiple psychiatric liability genes to oscillatory brain activity.

Oscillatory activity is crucial for information processing in the brain, and has a long history as a biomarker for psychopathology. Variation in oscillatory activity is highly heritable, but current understanding of specific genetic influences remains limited. We performed the largest genome-wide association study to date of oscillatory power during eyes-closed resting electroencephalogram (EEG) across a range of frequencies (delta 1-3.75 Hz, theta 4-7.75 Hz, alpha 8-12.75 Hz, and beta 13-30 Hz) in 8,425 subjects. Additionally, we performed KGG positional gene-based analysis and brain-expression analyses. GABRA2-a known genetic marker for alcohol use disorder and epilepsy-significantly affected beta power, consistent with the known relation between GABAA interneuron activity and beta oscillations. Tissue-specific SNP-based imputation of gene-expression levels based on the GTEx database revealed that hippocampal GABRA2 expression may mediate this effect. Twenty-four genes at 3p21.1 were significant for alpha power (FDR q < .05). SNPs in this region were linked to expression of GLYCTK in hippocampal tissue, and GNL3 and ITIH4 in the frontal cortex-genes that were previously implicated in schizophrenia and bipolar disorder. In sum, we identified several novel genetic variants associated with oscillatory brain activity; furthermore, we replicated and advanced understanding of previously known genes associated with psychopathology (i.e., schizophrenia and alcohol use disorders). Importantly, these psychopathological liability genes affect brain functioning, linking the genes' expression to specific cortical/subcortical brain regions.

Smaller hippocampal CA1 subfield volume in posttraumatic stress disorder.

BACKGROUND: Smaller hippocampal volume in patients with posttraumatic stress disorder (PTSD) represents the most consistently reported structural alteration in the brain. Subfields of the hippocampus play distinct roles in encoding and processing of memories, which are disrupted in PTSD. We examined PTSD-associated alterations in 12 hippocampal subfields in relation to global hippocampal shape, and clinical features. METHODS: Case-control cross-sectional studies of U.S. military veterans (n = 282) from the Iraq and Afghanistan era were grouped into PTSD (n = 142) and trauma-exposed controls (n = 140). Participants underwent clinical evaluation for PTSD and associated clinical parameters followed by MRI at 3 T. Segmentation with FreeSurfer v6.0 produced hippocampal subfield volumes for the left and right CA1, CA3, CA4, DG, fimbria, fissure, hippocampus-amygdala transition area, molecular layer, parasubiculum, presubiculum, subiculum, and tail, as well as hippocampal meshes. Covariates included age, gender, trauma exposure, alcohol use, depressive symptoms, antidepressant medication use, total hippocampal volume, and MRI scanner model. RESULTS: Significantly lower subfield volumes were associated with PTSD in left CA1 (P = 0.01; d = 0.21; uncorrected), CA3 (P = 0.04; d = 0.08; uncorrected), and right CA3 (P = 0.02; d = 0.07; uncorrected) only if ipsilateral whole hippocampal volume was included as a covariate. A trend level association of L-CA1 with PTSD (F4, 221  = 3.32, P = 0.07) is present and the other subfield findings are nonsignificant if ipsilateral whole hippocampal volume is not included as a covariate. PTSD-associated differences in global hippocampal shape were nonsignificant. CONCLUSIONS: The present finding of smaller hippocampal CA1 in PTSD is consistent with model systems in rodents that exhibit increased anxiety-like behavior from repeated exposure to acute stress. Behavioral correlations with hippocampal subfield volume differences in PTSD will elucidate their relevance to PTSD, particularly behaviors of associative fear learning, extinction training, and formation of false memories.

Temporal Cortex Morphology in Mesial Temporal Lobe Epilepsy Patients and Their Asymptomatic Siblings.

Temporal cortex abnormalities are common in patients with mesial temporal lobe epilepsy due to hippocampal sclerosis (MTLE+HS) and believed to be relevant to the underlying mechanisms. In the present study, we set out to determine the familiarity of temporal cortex morphologic alterations in a cohort of MTLE+HS patients and their asymptomatic siblings. A surface-based morphometry (SBM) method was applied to process MRI data acquired from 140 individuals (50 patients with unilateral MTLE+HS, 50 asymptomatic siblings of patients, and 40 healthy controls). Using a region-of-interest approach, alterations in temporal cortex morphology were determined in patients and their asymptomatic siblings by comparing with the controls. Alterations in temporal cortex morphology were identified in MTLE+HS patients ipsilaterally within the anterio-medial regions, including the entorhinal cortex, parahippocampal gyrus, and temporal pole. Subtle but similar pattern of morphology changes with a medium effect size were also noted in the asymptomatic siblings. These localized alterations were related to volume loss that appeared driven by shared contractions in cerebral cortex surface area. These findings indicate that temporal cortex morphologic alterations are common to patients and their asymptomatic siblings and suggest that such localized traits are possibly heritable.

Heritability and reliability of automatically segmented human hippocampal formation subregions.

The human hippocampal formation can be divided into a set of cytoarchitecturally and functionally distinct subregions, involved in different aspects of memory formation. Neuroanatomical disruptions within these subregions are associated with several debilitating brain disorders including Alzheimer's disease, major depression, schizophrenia, and bipolar disorder. Multi-center brain imaging consortia, such as the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) consortium, are interested in studying disease effects on these subregions, and in the genetic factors that affect them. For large-scale studies, automated extraction and subsequent genomic association studies of these hippocampal subregion measures may provide additional insight. Here, we evaluated the test-retest reliability and transplatform reliability (1.5T versus 3T) of the subregion segmentation module in the FreeSurfer software package using three independent cohorts of healthy adults, one young (Queensland Twins Imaging Study, N=39), another elderly (Alzheimer's Disease Neuroimaging Initiative, ADNI-2, N=163) and another mixed cohort of healthy and depressed participants (Max Planck Institute, MPIP, N=598). We also investigated agreement between the most recent version of this algorithm (v6.0) and an older version (v5.3), again using the ADNI-2 and MPIP cohorts in addition to a sample from the Netherlands Study for Depression and Anxiety (NESDA) (N=221). Finally, we estimated the heritability (h(2)) of the segmented subregion volumes using the full sample of young, healthy QTIM twins (N=728). Test-retest reliability was high for all twelve subregions in the 3T ADNI-2 sample (intraclass correlation coefficient (ICC)=0.70-0.97) and moderate-to-high in the 4T QTIM sample (ICC=0.5-0.89). Transplatform reliability was strong for eleven of the twelve subregions (ICC=0.66-0.96); however, the hippocampal fissure was not consistently reconstructed across 1.5T and 3T field strengths (ICC=0.47-0.57). Between-version agreement was moderate for the hippocampal tail, subiculum and presubiculum (ICC=0.78-0.84; Dice Similarity Coefficient (DSC)=0.55-0.70), and poor for all other subregions (ICC=0.34-0.81; DSC=0.28-0.51). All hippocampal subregion volumes were highly heritable (h(2)=0.67-0.91). Our findings indicate that eleven of the twelve human hippocampal subregions segmented using FreeSurfer version 6.0 may serve as reliable and informative quantitative phenotypes for future multi-site imaging genetics initiatives such as those of the ENIGMA consortium.

White matter alterations in patients with MRI-negative temporal lobe epilepsy and their asymptomatic siblings.

OBJECTIVE: The identification of "endophenotypes"-measurable variations along the pathways between genes and distal disease state-may help deconstruct focal epilepsies into more sensitive phenomena and improve future efforts to map the genetic underpinnings of the disorder. In this study, we set out to determine if diffusion tensor imaging (DTI)-inferred white matter (WM) alterations represent a suitable structural endophenotype for focal epilepsy. METHODS: We recruited 25 patients with sporadic mesial temporal lobe epilepsy (MTLE) with normal magnetic resonance imaging (MRI) findings, 25 of their gender-matched, asymptomatic siblings, and 60 control subjects. Whole-brain, voxelwise statistics were conducted to identify regions of microstructural degeneration in patients with MTLE and/or their asymptomatic siblings. WM tracts exhibiting evidence of microstructural disruption were then reconstructed using deterministic tractography. Diffusion metrics including fractional anisotropy (FA) and mean diffusivity (MD) were compared across groups using a series of one-way multivariate analyses of covariance (MANCOVAs). RESULTS: Voxelwise statistics revealed significant FA reductions in the corpus callosum (CC), bilateral superior longitudinal fasciculi (SLF), bilateral inferior longitudinal fasciculi (ILF), and left corticospinal tract (CST) in MTLE patients only. MD increases were observed in MTLE patients and their asymptomatic siblings in the left SLF and left CST. Deterministic tractography supported the voxelwise results, revealing significant FA alterations in the left SLF and CST in patients only and significant MD alterations in MTLE patients and their unaffected siblings. The diffusion scalars of MTLE patients and their asymptomatic siblings were highly correlated in the SLF and CST ipsilateral to patients' sides of seizure onset. SIGNIFICANCE: These findings confirm the presence of microstructural WM alterations in patients with MRI-negative MTLE and provide preliminary support for a diffusion-based endophenotype in the disorder. Further studies of narrow-sense heritability in larger cohorts of first-degree relatives of MTLE patients are required to confirm these results.

Cerebellum and cognition: evidence for the encoding of higher order rules.

Converging anatomical and functional evidence suggests that the cerebellum processes both motor and nonmotor information originating from the primary motor cortex and prefrontal cortex, respectively. However, it has not been established whether the cerebellum only processes prefrontal information where rules specify actions or whether the cerebellum processes any form of prefrontal information no matter how abstract. Using functional magnetic resonance imaging, we distinguish between two competing hypotheses: (1) activity within prefrontal-projecting cerebellar lobules (Crus I and II) will only be evoked by rules that specify action (i.e. first-order rules; arbitrary S-R mappings) and (2) activity will be evoked in these lobules by both first-order rules and second-order rules that govern the application of lower order rules. The results showed that prefrontal-projecting cerebellar lobules Crus I and II were commonly activated by processing both first- and second-order rules. We demonstrate for the first time that cerebellar circuits engage both first- and second-order rules and in doing so show that the cerebellum can contribute to cognitive control independent of motor control.

Genome-Wide Association Studies of ARIA From the Aducanumab Phase 3 ENGAGE and EMERGE Studies.

BACKGROUND AND OBJECTIVES: Amyloid-related imaging abnormalities (ARIA) were the most common adverse events reported in the phase 3 ENGAGE and EMERGE trials of aducanumab, an anti-amyloid monoclonal antibody. APOE ε4 carrier status has been shown to increase risk of ARIA in prior trials of aducanumab and other anti-amyloid therapies; however, the remainder of the human genome has not been evaluated for ARIA risk factors. Therefore, we sought to determine in a hypothesis-free manner whether genetic variants beyond APOE influence risk of ARIA in aducanumab-treated patients. METHODS: We performed genome-wide association studies (GWAS) of ARIA in participants in the ENGAGE and EMERGE trials. Participants had mild cognitive impairment due to Alzheimer disease or mild Alzheimer disease dementia and were amyloid-positive on PET scans. All participants underwent regular MRI monitoring to detect and diagnose ARIA. RESULTS: Of the 3,285 participants in the intent-to-treat population, this analysis included 1,691 with genotyping array data who received at least one dose of aducanumab with at least one post-baseline MRI. All participants in the study cohort were of European ancestry; 51% were female. The mean age was 70.3 years. 31% had ARIA-E, 19% had ARIA-H microhemorrhage, and 14% had ARIA-H superficial siderosis. We identified one genome-wide significant (p < 5.0 × 10-8) association at the chromosome 19 locus encompassing APOE. The APOE association with ARIA was stronger in ε4/ε4 homozygotes (OR = 4.28, 4.58, 7.84; p < 2.9 × 10-14 for ARIA-E, ARIA-H microhemorrhage, and ARIA-H superficial siderosis, respectively) than in ε3/ε4 heterozygotes (OR = 1.74, 1.46, 3.14; p ≤ 0.03). We found greater odds of radiographically severe ARIA (OR = 7.04-24.64, p ≤ 2.72 × 10-5) than radiographically mild ARIA (OR = 3.19-5.00, p ≤ 1.37 × 10-5) among ε4/ε4 homozygotes. APOE ε4 was also significantly associated with both symptomatic (ε4/ε4 OR = 3.64-9.52; p < 0.004) and asymptomatic (ε4/ε4 OR = 4.20-7.94, p < 1.7 × 10-11) cases, although among ARIA cases, APOE did not appear to modulate symptomatic status. No other genome-wide significant associations were found. DISCUSSION: We identified a strong, genome-wide significant association between APOE and risk of ARIA. Future, larger studies may be better powered to detect associations beyond APOE. These findings indicate that APOE is the strongest genetic risk factor of ARIA incidence, with implications for patient management and risk-benefit treatment decisions. TRIAL REGISTRATION INFORMATION: Both trials (ENGAGE [221AD301]: NCT02477800 and EMERGE [221AD302]: NCT02484547) were registered in June 2015 at clinicaltrials.gov and enrolled patients from August 2015 to July 2018.

Regional increase of cerebral cortex thickness in juvenile myoclonic epilepsy.

The goal of this study was to characterize cerebral cortex thickness patterns in juvenile myoclonic epilepsy (JME). Surface-based morphometry (SBM) was applied to process brain magnetic resonance images acquired from 24 patients with JME and 40 healthy controls and quantify cerebral cortex thickness. Differences in cortical thickness between patients and controls were determined using generalized linear model (covariates: age and gender). In patients with JME, thickness increase was detected bilaterally within localized regions in the orbitofrontal and mesial frontal cortices. Such thickness patterns coexisted with significant bilateral reduction in thalamic volume. These findings confirm that the underlying mechanisms in JME are related to aberrant corticothalamic structure and indicate that frontal cortex abnormalities are possibly linked to regional increase in cerebral cortical thickness.

The genetic regulation of protein expression in cerebrospinal fluid.

Studies of the genetic regulation of cerebrospinal fluid (CSF) proteins may reveal pathways for treatment of neurological diseases. 398 proteins in CSF were measured in 1,591 participants from the BioFINDER study. Protein quantitative trait loci (pQTL) were identified as associations between genetic variants and proteins, with 176 pQTLs for 145 CSF proteins (P < 1.25 × 10-10 , 117 cis-pQTLs and 59 trans-pQTLs). Ventricular volume (measured with brain magnetic resonance imaging) was a confounder for several pQTLs. pQTLs for CSF and plasma proteins were overall correlated, but CSF-specific pQTLs were also observed. Mendelian randomization analyses suggested causal roles for several proteins, for example, ApoE, CD33, and GRN in Alzheimer's disease, MMP-10 in preclinical Alzheimer's disease, SIGLEC9 in amyotrophic lateral sclerosis, and CD38, GPNMB, and ADAM15 in Parkinson's disease. CSF levels of GRN, MMP-10, and GPNMB were altered in Alzheimer's disease, preclinical Alzheimer's disease, and Parkinson's disease, respectively. These findings point to pathways to be explored for novel therapies. The novel finding that ventricular volume confounded pQTLs has implications for design of future studies of the genetic regulation of the CSF proteome.

The burden of rare protein-truncating genetic variants on human lifespan.

Genetic predisposition has been shown to contribute substantially to the age at which we die. Genome-wide association studies (GWASs) have linked more than 20 loci to phenotypes related to human lifespan1. However, little is known about how lifespan is impacted by gene loss of function. Through whole-exome sequencing of 352,338 UK Biobank participants of European ancestry, we assessed the relevance of protein-truncating variant (PTV) gene burden on individual and parental survival. We identified four exome-wide significant (P < 4.2 × 10-7) human lifespan genes, BRCA1, BRCA2, ATM and TET2. Gene and gene-set, PTV-burden, phenome-wide association studies support known roles of these genes in cancer to impact lifespan at the population level. The TET2 PTV burden was associated with a lifespan through somatic mutation events presumably due to clonal hematopoiesis. The overlap between PTV burden and common variant-based lifespan GWASs was modest, underscoring the value of exome sequencing in well-powered biobank cohorts to complement GWASs for identifying genes underlying complex traits.

Microglial activation protects against accumulation of tau aggregates in nondemented individuals with underlying Alzheimer's disease pathology.

The role of microglia in tau accumulation is currently unclear but could provide an important insight into the mechanisms underlying Alzheimer's disease (AD)1. Here, we measured the microglial marker soluble TREM2 and the disease-associated microglial activation stage 2 markers AXL, MERTK, GAS6, LPL, CST7, SPP1 and CSF1 in nondemented individuals from the Swedish BioFINDER-2 cohort who underwent longitudinal tau-positron emission tomography (PET), amyloid-PET and global cognitive assessment. To assess whether baseline microglial markers had an effect on AD-related changes, we studied three sub-groups of individuals: 121 with evidence of amyloid-PET pathology (A+), 64 with additional evidence of tau-PET pathology (A+T+) and 159 without amyloid- or tau-PET pathology (A-T-). Our results showed that increased levels of TREM2 were associated with slower amyloid accumulation in A+ individuals in addition to slower tau deposition and cognitive decline in A+T+ subjects. Similarly, higher levels of AXL, MERTK, GAS6, LPL, CST7 and CSF1 predicted slower tau accumulation and/or cognitive decline in the A+T+ group. These findings have important implications for future therapeutic strategies aiming to boost microglial protective functions in AD.

Genetic map of regional sulcal morphology in the human brain from UK biobank data.

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Genetic variants associated with longitudinal changes in brain structure across the lifespan.

Human brain structure changes throughout the lifespan. Altered brain growth or rates of decline are implicated in a vast range of psychiatric, developmental and neurodegenerative diseases. In this study, we identified common genetic variants that affect rates of brain growth or atrophy in what is, to our knowledge, the first genome-wide association meta-analysis of changes in brain morphology across the lifespan. Longitudinal magnetic resonance imaging data from 15,640 individuals were used to compute rates of change for 15 brain structures. The most robustly identified genes GPR139, DACH1 and APOE are associated with metabolic processes. We demonstrate global genetic overlap with depression, schizophrenia, cognitive functioning, insomnia, height, body mass index and smoking. Gene set findings implicate both early brain development and neurodegenerative processes in the rates of brain changes. Identifying variants involved in structural brain changes may help to determine biological pathways underlying optimal and dysfunctional brain development and aging.

BioInfograph: An Online Tool to Design and Display Multi-Panel Scientific Figure Interactively.

Many fit-for-purpose bioinformatics tools generate plots to interpret complex biological data and illustrate findings. However, assembling individual plots in different formats from various sources into one high-resolution figure in the desired layout requires mastery of commercial tools or even programming skills. In addition, it is a time-consuming and sometimes frustrating process even for a computationally savvy scientist who frequently takes a trial-and-error iterative approach to get satisfactory results. To address the challenge, we developed bioInfograph, a web-based tool that allows users to interactively arrange high-resolution images in diversified formats, mainly Scalable Vector Graphics (SVG), to produce one multi-panel publication-quality composite figure in both PDF and HTML formats in a user-friendly manner, requiring no programming skills. It solves stylesheet conflicts of coexisting SVG plots, integrates a rich-text editor, and allows creative design by providing advanced functionalities like image transparency, controlled vertical stacking of plots, versatile image formats, and layout templates. To highlight, the sharable interactive HTML output with zoom-in function is a unique feature not seen in any other similar tools. In the end, we make the online tool publicly available at https://baohongz.github.io/bioInfograph while releasing the source code at https://github.com/baohongz/bioInfograph under MIT open-source license.