Proteomic Network Analysis of Alzheimer's Disease Cerebrospinal Fluid Reveals Alterations Associated with APOE ε4 Genotype and Atomoxetine Treatment.

Alzheimer's disease (AD) is currently defined at the research level by the aggregation of amyloid-ß (Aß) and tau proteins in brain. While biofluid biomarkers are available to measure Aß and tau pathology, few biomarkers are available to measure the complex pathophysiology that is associated with these two cardinal neuropathologies. Here we describe the proteomic landscape of cerebrospinal fluid (CSF) changes associated with Aß and tau pathology in 300 individuals as assessed by two different proteomic technologies-tandem mass tag (TMT) mass spectrometry and SomaScan. Harmonization and integration of both data types allowed for generation of a robust protein co-expression network consisting of 34 modules derived from 5242 protein measurements, including disease-relevant modules associated with autophagy, ubiquitination, endocytosis, and glycolysis. Three modules strongly associated with the apolipoprotein E ε4 (APOE ε4) AD risk genotype mapped to oxidant detoxification, mitogen associated protein kinase (MAPK) signaling, neddylation, and mitochondrial biology, and overlapped with a previously described lipoprotein module in serum. Neddylation and oxidant detoxification/MAPK signaling modules had a negative association with APOE ε4 whereas the mitochondrion module had a positive association with APOE ε4. The directions of association were consistent between CSF and blood in two independent longitudinal cohorts, and altered levels of all three modules in blood were associated with dementia over 20 years prior to diagnosis. Dual-proteomic platform analysis of CSF samples from an AD phase 2 clinical trial of atomoxetine (ATX) demonstrated that abnormal elevations in the glycolysis CSF module-the network module most strongly correlated to cognitive function-were reduced by ATX treatment. Individuals who had more severe glycolytic changes at baseline responded better to ATX. Clustering of individuals based on their CSF proteomic network profiles revealed ten groups that did not cleanly stratify by Aß and tau status, underscoring the heterogeneity of pathological changes not fully reflected by Aß and tau. AD biofluid proteomics holds promise for the development of biomarkers that reflect diverse pathologies for use in clinical trials and precision medicine.

Sex differences in brain protein expression and disease.

Most complex human traits differ by sex, but we have limited insight into the underlying mechanisms. Here, we investigated the influence of biological sex on protein expression and its genetic regulation in 1,277 human brain proteomes. We found that 13.2% (1,354) of brain proteins had sex-differentiated abundance and 1.5% (150) of proteins had sex-biased protein quantitative trait loci (sb-pQTLs). Among genes with sex-biased expression, we found 67% concordance between sex-differentiated protein and transcript levels; however, sex effects on the genetic regulation of expression were more evident at the protein level. Considering 24 psychiatric, neurologic and brain morphologic traits, we found that an average of 25% of their putatively causal genes had sex-differentiated protein abundance and 12 putatively causal proteins had sb-pQTLs. Furthermore, integrating sex-specific pQTLs with sex-stratified genome-wide association studies of six psychiatric and neurologic conditions, we uncovered another 23 proteins contributing to these traits in one sex but not the other. Together, these findings begin to provide insights into mechanisms underlying sex differences in brain protein expression and disease.

Shared mechanisms across the major psychiatric and neurodegenerative diseases.

Several common psychiatric and neurodegenerative diseases share epidemiologic risk; however, whether they share pathophysiology is unclear and is the focus of our investigation. Using 25 GWAS results and LD score regression, we find eight significant genetic correlations between psychiatric and neurodegenerative diseases. We integrate the GWAS results with human brain transcriptomes (n = 888) and proteomes (n = 722) to identify cis- and trans- transcripts and proteins that are consistent with a pleiotropic or causal role in each disease, referred to as causal proteins for brevity. Within each disease group, we find many distinct and shared causal proteins. Remarkably, 30% (13 of 42) of the neurodegenerative disease causal proteins are shared with psychiatric disorders. Furthermore, we find 2.6-fold more protein-protein interactions among the psychiatric and neurodegenerative causal proteins than expected by chance. Together, our findings suggest these psychiatric and neurodegenerative diseases have shared genetic and molecular pathophysiology, which has important ramifications for early treatment and therapeutic development.

Alzheimer's disease genetic burden is associated with mid-life depression among persons with normal cognition.

INTRODUCTION: Despite an established link between depression and higher Alzheimer's disease (AD) risk, it is unclear whether the conditions share pathophysiology. Here, we investigated whether depression manifesting after age 50 is associated with a genetic predisposition to AD. METHODS: From the population-based Health and Retirement Study cohort with biennial assessments of depressive symptoms and cognitive performance, we studied 6656 individuals of European ancestry with whole-genome genotyping. Polygenic risk scores (PRSs) for AD were estimated and examined for an association with depression in cognitively normal participants using regression modeling. RESULTS: Among cognitively normal participants, those with a higher AD PRS were more likely to experience depression after age 50 after accounting for the effects of genetic predisposition to depression, sex, age, and education. DISCUSSION: Genetic predisposition to AD may be one of the factors contributing to the pathogenesis of mid-life depression. Whether there is a shared genetic basis between mid-life depression and AD merits further study.

Integrating human brain proteomes with genome-wide association data implicates novel proteins in post-traumatic stress disorder.

Genome-wide association studies (GWAS) have identified several risk loci for post-traumatic stress disorder (PTSD); however, how they confer PTSD risk remains unclear. We aimed to identify genes that confer PTSD risk through their effects on brain protein abundance to provide new insights into PTSD pathogenesis. To that end, we integrated human brain proteomes with PTSD GWAS results to perform a proteome-wide association study (PWAS) of PTSD, followed by Mendelian randomization, using a discovery and confirmatory study design. Brain proteomes (N = 525) were profiled from the dorsolateral prefrontal cortex using mass spectrometry. The Million Veteran Program (MVP) PTSD GWAS (n = 186,689) was used for the discovery PWAS, and the Psychiatric Genomics Consortium PTSD GWAS (n = 174,659) was used for the confirmatory PWAS. To understand whether genes identified at the protein-level were also evident at the transcript-level, we performed a transcriptome-wide association study (TWAS) using human brain transcriptomes (N = 888) and the MVP PTSD GWAS results. We identified 11 genes that contribute to PTSD pathogenesis via their respective cis-regulated brain protein abundance. Seven of 11 genes (64%) replicated in the confirmatory PWAS and 4 of 11 also had their cis-regulated brain mRNA levels associated with PTSD. High confidence level was assigned to 9 of 11 genes after considering evidence from the confirmatory PWAS and TWAS. Most of the identified genes are expressed in other PTSD-relevant brain regions and several are preferentially expressed in excitatory neurons, astrocytes, and oligodendrocyte precursor cells. These genes are novel, promising targets for mechanistic and therapeutic studies to find new treatments for PTSD.

Genetic Evidence Supporting a Causal Role of Depression in Alzheimer's Disease.

BACKGROUND: Depression has been associated with a higher risk of Alzheimer's disease (AD) in several prospective studies; however, mechanisms underlying this association remain unclear. METHODS: We examined genetic correlation between depression and AD using linkage disequilibrium score regression. We then tested for evidence of causality between depression and AD using Mendelian randomization and genome-wide association study results. Subsequently, cis and trans quantitative trait locus analyses for the depression genome-wide association study signals were performed to resolve the genetic signals to specific DNA methylation sites, brain transcripts, and proteins. These transcripts and proteins were then examined for associations with AD and its endophenotypes. Finally, the associations between depression polygenic risk score and AD endophenotypes were examined. RESULTS: We detected a significant genetic correlation between depression and AD, suggesting that they have a shared genetic basis. Furthermore, we found that depression had a causal role in AD through Mendelian randomization but did not find evidence for a causal role of AD on depression. Moreover, we identified 75 brain transcripts and 28 brain proteins regulated by the depression genome-wide association study signals through quantitative trait locus analyses. Of these, 46 transcripts and seven proteins were associated with rates of cognitive decline over time, AD pathologies, and AD diagnosis in two separate cohorts, thus implicating them in AD. In addition, we found that a higher depression polygenic risk score was associated with a faster decline of episodic memory over time. CONCLUSIONS: Depression appears to have a causal role in AD, and this causal relationship is likely driven, in part, by the 53 brain transcripts and proteins identified in this study.

Large-scale deep multi-layer analysis of Alzheimer's disease brain reveals strong proteomic disease-related changes not observed at the RNA level.

The biological processes that are disrupted in the Alzheimer's disease (AD) brain remain incompletely understood. In this study, we analyzed the proteomes of more than 1,000 brain tissues to reveal new AD-related protein co-expression modules that were highly preserved across cohorts and brain regions. Nearly half of the protein co-expression modules, including modules significantly altered in AD, were not observed in RNA networks from the same cohorts and brain regions, highlighting the proteopathic nature of AD. Two such AD-associated modules unique to the proteomic network included a module related to MAPK signaling and metabolism and a module related to the matrisome. The matrisome module was influenced by the APOE ε4 allele but was not related to the rate of cognitive decline after adjustment for neuropathology. By contrast, the MAPK/metabolism module was strongly associated with the rate of cognitive decline. Disease-associated modules unique to the proteome are sources of promising therapeutic targets and biomarkers for AD.

A machine learning approach to brain epigenetic analysis reveals kinases associated with Alzheimer's disease.

Alzheimer's disease (AD) is influenced by both genetic and environmental factors; thus, brain epigenomic alterations may provide insights into AD pathogenesis. Multiple array-based Epigenome-Wide Association Studies (EWASs) have identified robust brain methylation changes in AD; however, array-based assays only test about 2% of all CpG sites in the genome. Here, we develop EWASplus, a computational method that uses a supervised machine learning strategy to extend EWAS coverage to the entire genome. Application to six AD-related traits predicts hundreds of new significant brain CpGs associated with AD, some of which are further validated experimentally. EWASplus also performs well on data collected from independent cohorts and different brain regions. Genes found near top EWASplus loci are enriched for kinases and for genes with evidence for physical interactions with known AD genes. In this work, we show that EWASplus implicates additional epigenetic loci for AD that are not found using array-based AD EWASs.

A Genetic Study of Cerebral Atherosclerosis Reveals Novel Associations with NTNG1 and CNOT3.

Cerebral atherosclerosis is a leading cause of stroke and an important contributor to dementia. Yet little is known about its genetic basis. To examine the association of common single nucleotide polymorphisms with cerebral atherosclerosis severity, we conducted a genomewide association study (GWAS) using data collected as part of two community-based cohort studies in the United States, the Religious Orders Study (ROS) and Rush Memory and Aging Project (MAP). Both studies enroll older individuals and exclude participants with signs of dementia at baseline. From our analysis of 1325 participants of European ancestry who had genotype and neuropathologically assessed cerebral atherosclerosis measures available, we found a novel locus for cerebral atherosclerosis in NTNG1. The locus comprises eight SNPs, including two independent significant SNPs: rs6664221 (ß = -0.27, 95% CI = (-0.35, -0.19), p = 1.29 × 10-10) and rs10881463 (ß = -0.20, 95% CI = (-0.27, -0.13), p = 3.40 × 10-8). We further found that the SNPs may influence cerebral atherosclerosis by regulating brain protein expression of CNOT3. CNOT3 is a subunit of CCR4-NOT, which has been shown to be a master regulator of mRNA stability and translation and an important complex for cholesterol homeostasis. In summary, we identify a novel genetic locus for cerebral atherosclerosis and a potential mechanism linking this variation to cerebral atherosclerosis progression. These findings offer insights into the genetic effects on cerebral atherosclerosis.

Transcriptome-wide association study of post-trauma symptom trajectories identified GRIN3B as a potential biomarker for PTSD development.

Biomarkers that predict symptom trajectories after trauma can facilitate early detection or intervention for posttraumatic stress disorder (PTSD) and may also advance our understanding of its biology. Here, we aimed to identify trajectory-based biomarkers using blood transcriptomes collected in the immediate aftermath of trauma exposure. Participants were recruited from an Emergency Department in the immediate aftermath of trauma exposure and assessed for PTSD symptoms at baseline, 1, 3, 6, and 12 months. Three empirical symptom trajectories (chronic-PTSD, remitting, and resilient) were identified in 377 individuals based on longitudinal symptoms across four data points (1, 3, 6, and 12 months), using latent growth mixture modeling. Blood transcriptomes were examined for association with longitudinal symptom trajectories, followed by expression quantitative trait locus analysis. GRIN3B and AMOTL1 blood mRNA levels were associated with chronic vs. resilient post-trauma symptom trajectories at a transcriptome-wide significant level (N = 153, FDR-corrected p value = 0.0063 and 0.0253, respectively). We identified four genetic variants that regulate mRNA blood expression levels of GRIN3B. Among these, GRIN3B rs10401454 was associated with PTSD in an independent dataset (N = 3521, p = 0.04). Examination of the BrainCloud and GTEx databases revealed that rs10401454 was associated with brain mRNA expression levels of GRIN3B. While further replication and validation studies are needed, our data suggest that GRIN3B, a glutamate ionotropic receptor NMDA type subunit-3B, may be involved in the manifestation of PTSD. In addition, the blood mRNA level of GRIN3B may be a promising early biomarker for the PTSD manifestation and development.

Multi-tissue neocortical transcriptome-wide association study implicates 8 genes across 6 genomic loci in Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) is an incurable neurodegenerative disease currently affecting 1.75% of the US population, with projected growth to 3.46% by 2050. Identifying common genetic variants driving differences in transcript expression that confer AD risk is necessary to elucidate AD mechanism and develop therapeutic interventions. We modify the FUSION transcriptome-wide association study (TWAS) pipeline to ingest gene expression values from multiple neocortical regions. METHODS: A combined dataset of 2003 genotypes clustered to 1000 Genomes individuals from Utah with Northern and Western European ancestry (CEU) was used to construct a training set of 790 genotypes paired to 888 RNASeq profiles from temporal cortex (TCX = 248), prefrontal cortex (FP = 50), inferior frontal gyrus (IFG = 41), superior temporal gyrus (STG = 34), parahippocampal cortex (PHG = 34), and dorsolateral prefrontal cortex (DLPFC = 461). Following within-tissue normalization and covariate adjustment, predictive weights to impute expression components based on a gene's surrounding cis-variants were trained. The FUSION pipeline was modified to support input of pre-scaled expression values and support cross validation with a repeated measure design arising from the presence of multiple transcriptome samples from the same individual across different tissues. RESULTS: Cis-variant architecture alone was informative to train weights and impute expression for 6780 (49.67%) autosomal genes, the majority of which significantly correlated with gene expression; FDR < 5%: N = 6775 (99.92%), Bonferroni: N = 6716 (99.06%). Validation of weights in 515 matched genotype to RNASeq profiles from the CommonMind Consortium (CMC) was (72.14%) in DLPFC profiles. Association of imputed expression components from all 2003 genotype profiles yielded 8 genes significantly associated with AD (FDR < 0.05): APOC1, EED, CD2AP, CEACAM19, CLPTM1, MTCH2, TREM2, and KNOP1. CONCLUSIONS: We provide evidence of cis-genetic variation conferring AD risk through 8 genes across six distinct genomic loci. Moreover, we provide expression weights for 6780 genes as a valuable resource to the community, which can be abstracted across the neocortex and a wide range of neuronal phenotypes.

Brain proteome-wide association study implicates novel proteins in depression pathogenesis.

Depression is a common condition, but current treatments are only effective in a subset of individuals. To identify new treatment targets, we integrated depression genome-wide association study (GWAS) results (N = 500,199) with human brain proteomes (N = 376) to perform a proteome-wide association study of depression followed by Mendelian randomization. We identified 19 genes that were consistent with being causal in depression, acting via their respective cis-regulated brain protein abundance. We replicated nine of these genes using an independent depression GWAS (N = 307,353) and another human brain proteomic dataset (N = 152). Eleven of the 19 genes also had cis-regulated mRNA levels that were associated with depression, based on integration of the depression GWAS with human brain transcriptomes (N = 888). Meta-analysis of the discovery and replication proteome-wide association study analyses identified 25 brain proteins consistent with being causal in depression, 20 of which were not previously implicated in depression by GWAS. Together, these findings provide promising brain protein targets for further mechanistic and therapeutic studies.

Genetic control of the human brain proteome.

We generated an online brain pQTL resource for 7,376 proteins through the analysis of genetic and proteomic data derived from post-mortem samples of the dorsolateral prefrontal cortex of 330 older adults. The identified pQTLs tend to be non-synonymous variation, are over-represented among variants associated with brain diseases, and replicate well (77%) in an independent brain dataset. Comparison to a large study of brain eQTLs revealed that about 75% of pQTLs are also eQTLs. In contrast, about 40% of eQTLs were identified as pQTLs. These results are consistent with lower pQTL mapping power and greater evolutionary constraint on protein abundance. The latter is additionally supported by observations of pQTLs with large effects' tending to be rare, deleterious, and associated with proteins that have evidence for fewer protein-protein interactions. Mediation analyses using matched transcriptomic and proteomic data provided additional evidence that pQTL effects are often, but not always, mediated by mRNA. Specifically, we identified roughly 1.6 times more mRNA-mediated pQTLs than mRNA-independent pQTLs (550 versus 341). Our pQTL resource provides insight into the functional consequences of genetic variation in the human brain and a basis for novel investigations of genetics and disease.

Integrating human brain proteomes with genome-wide association data implicates new proteins in Alzheimer's disease pathogenesis.

Genome-wide association studies (GWAS) have identified many risk loci for Alzheimer's disease (AD)1,2, but how these loci confer AD risk is unclear. Here, we aimed to identify loci that confer AD risk through their effects on brain protein abundance to provide new insights into AD pathogenesis. To that end, we integrated AD GWAS results with human brain proteomes to perform a proteome-wide association study (PWAS) of AD, followed by Mendelian randomization and colocalization analysis. We identified 11 genes that are consistent with being causal in AD, acting via their cis-regulated brain protein abundance. Nine replicated in a confirmation PWAS and eight represent new AD risk genes not identified before by AD GWAS. Furthermore, we demonstrated that our results were independent of APOE e4. Together, our findings provide new insights into AD pathogenesis and promising targets for further mechanistic and therapeutic studies.

Publisher Correction: Shared proteomic effects of cerebral atherosclerosis and Alzheimer's disease on the human brain.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Shared proteomic effects of cerebral atherosclerosis and Alzheimer's disease on the human brain.

Cerebral atherosclerosis contributes to dementia via unclear processes. We performed proteomic sequencing of dorsolateral prefrontal cortex in 438 older individuals and found associations between cerebral atherosclerosis and reduced synaptic signaling and between RNA splicing and increased oligodendrocyte development and myelination. Consistently, single-cell RNA sequencing showed cerebral atherosclerosis associated with higher oligodendrocyte abundance. A subset of proteins and modules associated with cerebral atherosclerosis was also associated with Alzheimer's disease, suggesting shared mechanisms.

Benchmarking of computational error-correction methods for next-generation sequencing data.

BACKGROUND: Recent advancements in next-generation sequencing have rapidly improved our ability to study genomic material at an unprecedented scale. Despite substantial improvements in sequencing technologies, errors present in the data still risk confounding downstream analysis and limiting the applicability of sequencing technologies in clinical tools. Computational error correction promises to eliminate sequencing errors, but the relative accuracy of error correction algorithms remains unknown. RESULTS: In this paper, we evaluate the ability of error correction algorithms to fix errors across different types of datasets that contain various levels of heterogeneity. We highlight the advantages and limitations of computational error correction techniques across different domains of biology, including immunogenomics and virology. To demonstrate the efficacy of our technique, we apply the UMI-based high-fidelity sequencing protocol to eliminate sequencing errors from both simulated data and the raw reads. We then perform a realistic evaluation of error-correction methods. CONCLUSIONS: In terms of accuracy, we find that method performance varies substantially across different types of datasets with no single method performing best on all types of examined data. Finally, we also identify the techniques that offer a good balance between precision and sensitivity.

Brain microRNAs associated with late-life depressive symptoms are also associated with cognitive trajectory and dementia.

Late-life depression is associated with an increased risk for dementia but we have limited knowledge of the molecular mechanisms underlying this association. Here we investigated whether brain microRNAs, important posttranscriptional regulators of gene expression, contribute to this association. Late-life depressive symptoms were assessed annually in 300 participants of the Religious Orders Study and Rush Memory and Aging Project for a mean of 7 years. Participants underwent annual cognitive testing, clinical assessment of cognitive status, and uniform neuropathologic examination after death. microRNAs were profiled from the prefrontal cortex using NanoString platform in the discovery cohort and small RNA sequencing in the replication cohort. A global microRNA association study of late-life depressive symptoms was performed using linear mixed model adjusting for the potential confounding factors. Four brain microRNAs were associated with late-life depressive symptoms at adjusted p < 0.05: miR-484, miR-26b-5p, miR-30d-5p, and miR-197-3p. Lower expression levels of these miRNAs were associated having greater depressive symptoms. Furthermore, lower levels of miR-484 and miR-197-3p were associated with faster decline of cognition over time. Moreover, lower miR-484 level was associated with higher probability of having Alzheimer's dementia. Importantly, the associations between miR-484 and depressive symptoms and Alzheimer's dementia, respectively, were replicated in an independent cohort. Lastly, the predicted targets of miR-484 were enriched in a brain protein co-expression module involving synaptic transmission and regulation of synaptic plasticity. This study identified four brain microRNAs associated with late-life depressive symptoms assessed longitudinally. In addition, we found a molecular connection between late-life depression and dementia through miR-484.

Identification of Surface Epitopes Associated with Protection against Highly Immune-Evasive VlsE-Expressing Lyme Disease Spirochetes.

The tick-borne pathogen Borrelia burgdorferi is responsible for approximately 300,000 Lyme disease (LD) cases per year in the United States. Recent increases in the number of LD cases, in addition to the spread of the tick vector and a lack of a vaccine, highlight an urgent need for designing and developing an efficacious LD vaccine. Identification of protective epitopes that could be used to develop a second-generation (subunit) vaccine is therefore imperative. Despite the antigenicity of several lipoproteins and integral outer membrane proteins (OMPs) on the B. burgdorferi surface, the spirochetes successfully evade antibodies primarily due to the VlsE-mediated antigenic variation. VlsE is thought to sterically block antibody access to protective epitopes of B. burgdorferi However, it is highly unlikely that VlsE shields the entire surface epitome. Thus, identification of subdominant epitope targets that induce protection when they are made dominant is necessary to generate an efficacious vaccine. Toward the identification, we repeatedly immunized immunocompetent mice with live-attenuated VlsE-deleted B. burgdorferi and then challenged the animals with the VlsE-expressing (host-adapted) wild type. Passive immunization and Western blotting data suggested that the protection of 50% of repeatedly immunized animals against the highly immune-evasive B. burgdorferi was antibody mediated. Comparison of serum antibody repertoires identified in protected and nonprotected animals permitted the identification of several putative epitopes significantly associated with the protection. Most linear putative epitopes were conserved between the main pathogenic Borrelia genospecies and found within known subdominant regions of OMPs. Currently, we are performing immunization studies to test whether the identified protection-associated epitopes are protective for mice.

Identification of cancer-specific motifs in mimotope profiles of serum antibody repertoire.

BACKGROUND: For fighting cancer, earlier detection is crucial. Circulating auto-antibodies produced by the patient's own immune system after exposure to cancer proteins are promising bio-markers for the early detection of cancer. Since an antibody recognizes not the whole antigen but 4-7 critical amino acids within the antigenic determinant (epitope), the whole proteome can be represented by a random peptide phage display library. This opens the possibility to develop an early cancer detection test based on a set of peptide sequences identified by comparing cancer patients' and healthy donors' global peptide profiles of antibody specificities. RESULTS: Due to the enormously large number of peptide sequences contained in global peptide profiles generated by next generation sequencing, the large number of cancer and control sera is required to identify cancer-specific peptides with high degree of statistical significance. To decrease the number of peptides in profiles generated by nextgen sequencing without losing cancer-specific sequences we used for generation of profiles the phage library enriched by panning on the pool of cancer sera. To further decrease the complexity of profiles we used computational methods for transforming a list of peptides constituting the mimotope profiles to the list motifs formed by similar peptide sequences. CONCLUSION: We have shown that the amino-acid order is meaningful in mimotope motifs since they contain significantly more peptides than motifs among peptides where amino-acids are randomly permuted. Also the single sample motifs significantly differ from motifs in peptides drawn from multiple samples. Finally, multiple cancer-specific motifs have been identified.