Genetic and multi-omic risk assessment of Alzheimer's disease implicates core associated biological domains.

INTRODUCTION: Alzheimer's disease (AD) is the predominant dementia globally, with heterogeneous presentation and penetrance of clinical symptoms, variable presence of mixed pathologies, potential disease subtypes, and numerous associated endophenotypes. Beyond the difficulty of designing treatments that address the core pathological characteristics of the disease, therapeutic development is challenged by the uncertainty of which endophenotypic areas and specific targets implicated by those endophenotypes to prioritize for further translational research. However, publicly funded consortia driving large-scale open science efforts have produced multiple omic analyses that address both disease risk relevance and biological process involvement of genes across the genome. METHODS: Here we report the development of an informatic pipeline that draws from genetic association studies, predicted variant impact, and linkage with dementia associated phenotypes to create a genetic risk score. This is paired with a multi-omic risk score utilizing extensive sets of both transcriptomic and proteomic studies to identify system-level changes in expression associated with AD. These two elements combined constitute our target risk score that ranks AD risk genome-wide. The ranked genes are organized into endophenotypic space through the development of 19 biological domains associated with AD in the described genetics and genomics studies and accompanying literature. The biological domains are constructed from exhaustive Gene Ontology (GO) term compilations, allowing automated assignment of genes into objectively defined disease-associated biology. This rank-and-organize approach, performed genome-wide, allows the characterization of aggregations of AD risk across biological domains. RESULTS: The top AD-risk-associated biological domains are Synapse, Immune Response, Lipid Metabolism, Mitochondrial Metabolism, Structural Stabilization, and Proteostasis, with slightly lower levels of risk enrichment present within the other 13 biological domains. DISCUSSION: This provides an objective methodology to localize risk within specific biological endophenotypes and drill down into the most significantly associated sets of GO terms and annotated genes for potential therapeutic targets.

STOP-AD portal: Selecting the optimal pharmaceutical for preclinical drug testing in Alzheimer's disease.

We propose an unbiased methodology to rank compounds for advancement into comprehensive preclinical testing for Alzheimer's disease (AD). Translation of compounds to the clinic in AD has been hampered by poor predictive validity of models, compounds with limited pharmaceutical properties, and studies that lack rigor. To overcome this, MODEL-AD's Preclinical Testing Core developed a standardized pipeline for assessing efficacy in AD mouse models. We hypothesize that rank-ordering compounds based upon pharmacokinetic, efficacy, and toxicity properties in preclinical models will enhance successful translation to the clinic. Previously compound selection was based solely on physiochemical properties, with arbitrary cutoff limits, making ranking challenging. Since no gold standard exists for systematic prioritization, validating a selection criteria has remained elusive. The STOP-AD framework evaluates the drug-like properties to rank compounds for in vivo studies, and uses an unbiased approach that overcomes the validation limitation by performing Monte-Carlo simulations. HIGHLIGHTS: Promising preclinical studies for AD drugs have not translated to clinical success. Systematic assessment of AD drug candidates may increase clinical translatability. We describe a well-defined framework for compound selection with clear selection metrics.

Open drug discovery in Alzheimer's disease.

Alzheimer's disease (AD) drug discovery has focused on a set of highly studied therapeutic hypotheses, with limited success. The heterogeneous nature of AD processes suggests that a more diverse, systems-integrated strategy may identify new therapeutic hypotheses. Although many target hypotheses have arisen from systems-level modeling of human disease, in practice and for many reasons, it has proven challenging to translate them into drug discovery pipelines. First, many hypotheses implicate protein targets and/or biological mechanisms that are under-studied, meaning there is a paucity of evidence to inform experimental strategies as well as high-quality reagents to perform them. Second, systems-level targets are predicted to act in concert, requiring adaptations in how we characterize new drug targets. Here we posit that the development and open distribution of high-quality experimental reagents and informatic outputs-termed target enabling packages (TEPs)-will catalyze rapid evaluation of emerging systems-integrated targets in AD by enabling parallel, independent, and unencumbered research.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent.

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

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.

Author Correction: Molecular estimation of neurodegeneration pseudotime in older brains.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20261-6.

Molecular estimation of neurodegeneration pseudotime in older brains.

The temporal molecular changes that lead to disease onset and progression in Alzheimer's disease (AD) are still unknown. Here we develop a temporal model for these unobserved molecular changes with a manifold learning method applied to RNA-Seq data collected from human postmortem brain samples collected within the ROS/MAP and Mayo Clinic RNA-Seq studies. We define an ordering across samples based on their similarity in gene expression and use this ordering to estimate the molecular disease stage-or disease pseudotime-for each sample. Disease pseudotime is strongly concordant with the burden of tau (Braak score, P = 1.0 × 10-5), Aß (CERAD score, P = 1.8 × 10-5), and cognitive diagnosis (P = 3.5 × 10-7) of late-onset (LO) AD. Early stage disease pseudotime samples are enriched for controls and show changes in basic cellular functions. Late stage disease pseudotime samples are enriched for late stage AD cases and show changes in neuroinflammation and amyloid pathologic processes. We also identify a set of late stage pseudotime samples that are controls and show changes in genes enriched for protein trafficking, splicing, regulation of apoptosis, and prevention of amyloid cleavage pathways. In summary, we present a method for ordering patients along a trajectory of LOAD disease progression from brain transcriptomic data.

The AD Knowledge Portal: A Repository for Multi-Omic Data on Alzheimer's Disease and Aging.

The AD Knowledge Portal (adknowledgeportal.org) is a public data repository that shares data and other resources generated by multiple collaborative research programs focused on aging, dementia, and Alzheimer's disease (AD). In this article, we highlight how to use the Portal to discover and download genomic variant and transcriptomic data from the same individuals. First, we show how to use the web interface to browse and search for data of interest using relevant file annotations. We demonstrate how to learn more about the context surrounding the data, including diagnostic criteria and methodological details about sample preparation and data analysis. We present two primary ways to download data-using a web interface, and using a programmatic method that provides access using the command line. Finally, we show how to merge separate sources of metadata into a comprehensive file that contains factors and covariates necessary in downstream analyses. © 2020 The Authors. Basic Protocol 1: Find and download files associated with a selected study Basic Protocol 2: Download files in bulk using the command line client Basic Protocol 3: Working with file annotations and metadata.

Large eQTL meta-analysis reveals differing patterns between cerebral cortical and cerebellar brain regions.

The availability of high-quality RNA-sequencing and genotyping data of post-mortem brain collections from consortia such as CommonMind Consortium (CMC) and the Accelerating Medicines Partnership for Alzheimer's Disease (AMP-AD) Consortium enable the generation of a large-scale brain cis-eQTL meta-analysis. Here we generate cerebral cortical eQTL from 1433 samples available from four cohorts (identifying >4.1 million significant eQTL for >18,000 genes), as well as cerebellar eQTL from 261 samples (identifying 874,836 significant eQTL for >10,000 genes). We find substantially improved power in the meta-analysis over individual cohort analyses, particularly in comparison to the Genotype-Tissue Expression (GTEx) Project eQTL. Additionally, we observed differences in eQTL patterns between cerebral and cerebellar brain regions. We provide these brain eQTL as a resource for use by the research community. As a proof of principle for their utility, we apply a colocalization analysis to identify genes underlying the GWAS association peaks for schizophrenia and identify a potentially novel gene colocalization with lncRNA RP11-677M14.2 (posterior probability of colocalization 0.975).

Dopaminergic inhibition of gonadotropin-releasing hormone neurons in the cichlid fish Astatotilapia burtoni.

Dopamine regulates reproduction in part by modulating neuronal activity within the hypothalamic-pituitary-gonadal (HPG) axis. Previous studies suggested numerous mechanisms by which dopamine exerts inhibitory control over the HPG axis, ultimately changing the levels of sex steroids that regulate reproductive behaviors. However, it is not known whether these mechanisms are conserved across vertebrate species. In particular, it is unknown whether mechanisms underlying dopaminergic control of reproduction are shared between mammals and teleost fish. In mammals, dopamine directly inhibits gonadotropin-releasing hormone (GnRH1) hypothalamic neurons, the gatekeepers for activation of the HPG axis. Here, we demonstrate, for the first time in teleost fish, dopaminergic control of GnRH1 neurons via direct dopamine type-2-like receptor (D2R)-mediated inhibition within the hypothalamus. These results suggest that direct dopaminergic control of GnRH1 neurons via interactions in the hypothalamus is not exclusive to tetrapod reproductive control, but is likely conserved across vertebrate species.

Evolution of Schooling Behavior in Threespine Sticklebacks Is Shaped by the Eda Gene.

Despite longstanding interest in the genetic mechanisms that underlie behavioral evolution, very few genes that underlie naturally occurring variation in behavior between individuals or species are known, particularly in vertebrates. Here, we build on our previous forward genetic mapping experiments and use transgenic approaches to identify Ectodysplasin as a gene that causes differences in schooling behavior between wild populations of threespine stickleback (Gasterosteus aculeatus) fish. This work provides rare insight into the proximate mechanisms that have shaped the evolution of vertebrate behavior.

Social Regulation of Gene Expression in Threespine Sticklebacks.

Identifying genes that are differentially expressed in response to social interactions is informative for understanding the molecular basis of social behavior. To address this question, we described changes in gene expression as a result of differences in the extent of social interactions. We housed threespine stickleback (Gasterosteus aculeatus) females in either group conditions or individually for one week, then measured levels of gene expression in three brain regions using RNA-sequencing. We found that numerous genes in the hindbrain/cerebellum had altered expression in response to group or individual housing. However, relatively few genes were differentially expressed in either the diencephalon or telencephalon. The list of genes upregulated in fish from social groups included many genes related to neural development and cell adhesion as well as genes with functions in sensory signaling, stress, and social and reproductive behavior. The list of genes expressed at higher levels in individually-housed fish included several genes previously identified as regulated by social interactions in other animals. The identified genes are interesting targets for future research on the molecular mechanisms of normal social interactions.

Genetic mapping of natural variation in schooling tendency in the threespine stickleback.

Although there is a heritable basis for many animal behaviors, the genetic architecture of behavioral variation in natural populations remains mostly unknown, particularly in vertebrates. We sought to identify the genetic basis for social affiliation in two populations of threespine sticklebacks (Gasterosteus aculeatus) that differ in their propensity to school. Marine sticklebacks from Japan school strongly whereas benthic sticklebacks from a lake in Canada are more solitary. Here, we expanded on our previous efforts to identify quantitative trait loci (QTL) for differences in schooling tendency. We tested fish multiple times in two assays that test different aspects of schooling tendency: 1) the model school assay, which presents fish with a school of eight model sticklebacks; and 2) the choice assay, in which fish are given a choice between the model school and a stationary artificial plant. We found low-to-moderate levels of repeatability, ranging from 0.1 to 0.5, in schooling phenotypes. To identify the genomic regions that contribute to differences in schooling tendency, we used QTL mapping in two types of crosses: benthic × marine backcrosses and an F2 intercross. We found two QTL for time spent with the school in the model school assay, and one QTL for number of approaches to the school in the choice assay. These QTL were on three different linkage groups, not previously linked to behavioral differences in sticklebacks. Our results highlight the importance of using multiple crosses and robust behavioral assays to uncover the genetic basis of behavioral variation in natural populations.

Pleiotropic effects of a single gene on skeletal development and sensory system patterning in sticklebacks.

BACKGROUND: Adaptation to a new environment can be facilitated by co-inheritance of a suite of phenotypes that are all advantageous in the new habitat. Although experimental evidence demonstrates that multiple phenotypes often map to the same genomic regions, it is challenging to determine whether phenotypes are associated due to pleiotropic effects of a single gene or to multiple tightly linked genes. In the threespine stickleback fish (Gasterosteus aculeatus), multiple phenotypes are associated with a genomic region around the gene Ectodysplasin (Eda), but only the presence of bony lateral plates has been conclusively shown to be caused by Eda. RESULTS: Here, we ask whether pleiotropy or linkage is responsible for the association between lateral plates and the distribution of the neuromasts of the lateral line. We first identify a strong correlation between plate appearance and changes in the spatial distribution of neuromasts through development. We then use an Eda transgene to induce the formation of ectopic plates in low plated fish, which also results in alterations to neuromast distribution. Our results also show that other loci may modify the effects of Eda on plate formation and neuromast distribution. CONCLUSIONS: Together, these results demonstrate that Eda has pleiotropic effects on at least two phenotypes, highlighting the role of pleiotropy in the genetic basis of adaptation.

Genetic and neural modularity underlie the evolution of schooling behavior in threespine sticklebacks.

Although descriptions of striking diversity in animal behavior are plentiful, little is known about the mechanisms by which behaviors change and evolve between groups. To fully understand behavioral evolution, it will be necessary to identify the genetic mechanisms that mediate behavioral change in a natural context. Genetic analysis of behavior can also reveal associations between behavior and morphological or neural phenotypes, providing insight into the proximate mechanisms that control behavior. Relatively few studies to date have successfully identified genes or genomic regions that contribute to behavioral variation among natural populations or species, particularly in vertebrates. Here, we apply genetic approaches to dissect a complex social behavior that has long fascinated biologists, schooling behavior. We performed quantitative trait locus (QTL) analysis of schooling in an F2 intercross between strongly schooling marine and weakly schooling benthic sticklebacks (Gasterosteus aculeatus) and found that distinct genetic modules control different aspects of schooling behavior. Two key components of the behavior, tendency to school and body position when schooling, are uncorrelated in hybrids and map to different genomic regions. Our results further point to a genetic link between one behavioral component, schooling position, and variation in the neurosensory lateral line.

Molecular and developmental contributions to divergent pigment patterns in marine and freshwater sticklebacks.

Pigment pattern variation across species or populations offers a tractable framework in which to investigate the evolution of development. Juvenile threespine sticklebacks (Gasterosteus aculeatus) from marine and freshwater environments exhibit divergent pigment patterns that are associated with ecological differences. Juvenile marine sticklebacks have a silvery appearance, whereas sticklebacks from freshwater environments exhibit a pattern of vertical bars. We investigated both the developmental and molecular basis of this population-level variation in pigment pattern. Time course imaging during the transition from larval to juvenile stages revealed differences between marine and freshwater fish in spatial patterns of chromatophore differentiation as well as in pigment amount and dispersal. In freshwater fish, melanophores appear primarily within dark bars whereas iridophores appear within light bars. By contrast, in marine fish, these chromatophores are interspersed across the flank. In addition to spatially segregated chromatophore differentiation, pigment amount and dispersal within melanophores varies spatially across the flank of freshwater, but not marine fish. To gain insight into the molecular pathways that underlie the differences in pigment pattern development, we evaluated differential gene expression in the flanks of developing fish using high-throughput cDNA sequencing (RNA-seq) and quantitative PCR. We identified several genes that were differentially expressed across dark and light bars of freshwater fish, and between freshwater and marine fish. Together, these experiments begin to shed light on the process of pigment pattern evolution in sticklebacks.

Heritable differences in schooling behavior among threespine stickleback populations revealed by a novel assay.

Identifying the proximate and ultimate mechanisms of social behavior remains a major goal of behavioral biology. In particular, the complex social interactions mediating schooling behavior have long fascinated biologists, leading to theoretical and empirical investigations that have focused on schooling as a group-level phenomenon. However, methods to examine the behavior of individual fish within a school are needed in order to investigate the mechanisms that underlie both the performance and the evolution of schooling behavior. We have developed a technique to quantify the schooling behavior of an individual in standardized but easily manipulated social circumstances. Using our model school assay, we show that threespine sticklebacks (Gasterosteus aculeatus) from alternative habitats differ in behavior when tested in identical social circumstances. Not only do marine sticklebacks show increased association with the model school relative to freshwater benthic sticklebacks, they also display a greater degree of parallel swimming with the models. Taken together, these data indicate that marine sticklebacks exhibit a stronger tendency to school than benthic sticklebacks. We demonstrate that these population-level differences in schooling tendency are heritable and are shared by individuals within a population even when they have experienced mixed-population housing conditions. Finally, we begin to explore the stimuli that elicit schooling behavior in these populations. Our data suggest that the difference in schooling tendency between marine and benthic sticklebacks is accompanied by differential preferences for social vs. non-social and moving vs. stationary shelter options. Our study thus provides novel insights into the evolution of schooling behavior, as well as a new experimental approach to investigate the genetic and neural mechanisms that underlie this complex social behavior.

Sensory evolution: picking up good vibrations.

A new study has identified the sensory basis for behavioral adaptation in cavefish, providing insight into the process of nervous system evolution.

Expression of arginine vasotocin in distinct preoptic regions is associated with dominant and subordinate behaviour in an African cichlid fish.

Neuropeptides have widespread modulatory effects on behaviour and physiology and are associated with phenotypic transitions in a variety of animals. Arginine vasotocin (AVT) is implicated in mediating alternative male phenotypes in teleost fish, but the direction of the association differs among species, with either higher or lower AVT related to more territorial behaviour in different fishes. To clarify the complex relationship between AVT and alternative phenotype, we evaluated AVT expression in an African cichlid in which social status is associated with divergent behaviour and physiology. We compared AVT mRNA expression between territorial and non-territorial (NT) males in both whole brains and microdissected anterior preoptic areas using transcription profiling, and in individual preoptic nuclei using in situ hybridization. These complementary methods revealed that in the posterior preoptic area (gigantocellular nucleus), territorial males exhibit higher levels of AVT expression than NT males. Conversely, in the anterior preoptic area (parvocellular nucleus), AVT expression is lower in territorial males than NT males. We further correlated AVT expression with behavioural and physiological characteristics of social status to gain insight into the divergent functions of individual AVT nuclei. Overall, our findings highlight a complex association between AVT and social behaviour.

Differential social regulation of two pituitary gonadotropin-releasing hormone receptors.

In many vertebrates, social interactions regulate reproductive capacity by altering the activity of the hypothalamic-pituitary-gonadal (HPG) axis. To better understand the mechanisms underlying social regulation of reproduction, we investigated the relationship between social status and one main component of the HPG axis: expression levels of gonadotropin-releasing hormone receptor (GnRH-R). Social interactions dictate reproductive capacity in the cichlid fish Astatotilapia burtoni. Reproductively active territory holders suppress the HPG axis of non-territorial males through repeated aggressive encounters. To determine whether the expression of GnRH-R is socially regulated, we quantified mRNA levels of two GnRH-R variants in the pituitaries and brains of territorial (T) and non-territorial (NT) A. burtoni males. We found that T males had significantly higher levels of pituitary GnRH-R1 mRNA than NT males. In contrast, GnRH-R2 mRNA levels in the pituitary did not vary with social status. Pituitaries from both T and NT males expressed significantly higher mRNA levels of GnRH-R1 than GnRH-R2. GnRH mRNA levels in the brain correlated positively with GnRH-R1 mRNA levels in the pituitary but did not correlate with pituitary GnRH-R2. Measurements of GnRH-R1 and GnRH-R2 mRNA levels across the whole brain revealed no social status differences. These results show that, in addition to the known effects of social status on other levels of the HPG axis, GnRH receptor in the pituitary is also a target of social regulation.