Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism.

Autism spectrum disorder (ASD) is a complex developmental syndrome of unknown etiology. Recent studies employing exome- and genome-wide sequencing have identified nine high-confidence ASD (hcASD) genes. Working from the hypothesis that ASD-associated mutations in these biologically pleiotropic genes will disrupt intersecting developmental processes to contribute to a common phenotype, we have attempted to identify time periods, brain regions, and cell types in which these genes converge. We have constructed coexpression networks based on the hcASD "seed" genes, leveraging a rich expression data set encompassing multiple human brain regions across human development and into adulthood. By assessing enrichment of an independent set of probable ASD (pASD) genes, derived from the same sequencing studies, we demonstrate a key point of convergence in midfetal layer 5/6 cortical projection neurons. This approach informs when, where, and in what cell types mutations in these specific genes may be productively studied to clarify ASD pathophysiology.

Massively parallel discovery of human-specific substitutions that alter enhancer activity.

Genetic changes that altered the function of gene regulatory elements have been implicated in the evolution of human traits such as the expansion of the cerebral cortex. However, identifying the particular changes that modified regulatory activity during human evolution remain challenging. Here we used massively parallel enhancer assays in neural stem cells to quantify the functional impact of >32,000 human-specific substitutions in >4,300 human accelerated regions (HARs) and human gain enhancers (HGEs), which include enhancers with novel activities in humans. We found that >30% of active HARs and HGEs exhibited differential activity between human and chimpanzee. We isolated the effects of human-specific substitutions from background genetic variation to identify the effects of genetic changes most relevant to human evolution. We found that substitutions interacted in both additive and nonadditive ways to modify enhancer function. Substitutions within HARs, which are highly constrained compared to HGEs, showed smaller effects on enhancer activity, suggesting that the impact of human-specific substitutions is buffered in enhancers with constrained ancestral functions. Our findings yield insight into how human-specific genetic changes altered enhancer function and provide a rich set of candidates for studies of regulatory evolution in humans.

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.

Bone marrow microenvironments that contribute to patient outcomes in newly diagnosed multiple myeloma: A cohort study of patients in the Total Therapy clinical trials.

BACKGROUND: The tumor microenvironment (TME) is increasingly appreciated as an important determinant of cancer outcome, including in multiple myeloma (MM). However, most myeloma microenvironment studies have been based on bone marrow (BM) aspirates, which often do not fully reflect the cellular content of BM tissue itself. To address this limitation in myeloma research, we systematically characterized the whole bone marrow (WBM) microenvironment during premalignant, baseline, on treatment, and post-treatment phases. METHODS AND FINDINGS: Between 2004 and 2019, 998 BM samples were taken from 436 patients with newly diagnosed MM (NDMM) at the University of Arkansas for Medical Sciences in Little Rock, Arkansas, United States of America. These patients were 61% male and 39% female, 89% White, 8% Black, and 3% other/refused, with a mean age of 58 years. Using WBM and matched cluster of differentiation (CD)138-selected tumor gene expression to control for tumor burden, we identified a subgroup of patients with an adverse TME associated with 17 fewer months of progression-free survival (PFS) (95% confidence interval [CI] 5-29, 49-69 versus 70-82 months, χ2 p = 0.001) and 15 fewer months of overall survival (OS; 95% CI -1 to 31, 92-120 versus 113-129 months, χ2 p = 0.036). Using immunohistochemistry-validated computational tools that identify distinct cell types from bulk gene expression, we showed that the adverse outcome was correlated with elevated CD8+ T cell and reduced granulocytic cell proportions. This microenvironment develops during the progression of premalignant to malignant disease and becomes less prevalent after therapy, in which it is associated with improved outcomes. In patients with quantified International Staging System (ISS) stage and 70-gene Prognostic Risk Score (GEP-70) scores, taking the microenvironment into consideration would have identified an additional 40 out of 290 patients (14%, premutation p = 0.001) with significantly worse outcomes (PFS, 95% CI 6-36, 49-73 versus 74-90 months) who were not identified by existing clinical (ISS stage III) and tumor (GEP-70) criteria as high risk. The main limitations of this study are that it relies on computationally identified cell types and that patients were treated with thalidomide rather than current therapies. CONCLUSIONS: In this study, we observe that granulocyte signatures in the MM TME contribute to a more accurate prognosis. This implies that future researchers and clinicians treating patients should quantify TME components, in particular monocytes and granulocytes, which are often ignored in microenvironment studies.

Origin and evolution of developmental enhancers in the mammalian neocortex.

Morphological innovations such as the mammalian neocortex may involve the evolution of novel regulatory sequences. However, de novo birth of regulatory elements active during morphogenesis has not been extensively studied in mammals. Here, we use H3K27ac-defined regulatory elements active during human and mouse corticogenesis to identify enhancers that were likely active in the ancient mammalian forebrain. We infer the phylogenetic origins of these enhancers and find that ∼20% arose in the mammalian stem lineage, coincident with the emergence of the neocortex. Implementing a permutation strategy that controls for the nonrandom variation in the ages of background genomic sequences, we find that mammal-specific enhancers are overrepresented near genes involved in cell migration, cell signaling, and axon guidance. Mammal-specific enhancers are also overrepresented in modules of coexpressed genes in the cortex that are associated with these pathways, notably ephrin and semaphorin signaling. Our results also provide insight into the mechanisms of regulatory innovation in mammals. We find that most neocortical enhancers did not originate by en bloc exaptation of transposons. Young neocortical enhancers exhibit smaller H3K27ac footprints and weaker evolutionary constraint in eutherian mammals than older neocortical enhancers. Based on these observations, we present a model of the enhancer life cycle in which neocortical enhancers initially emerge from genomic background as short, weakly constrained "proto-enhancers." Many proto-enhancers are likely lost, but some may serve as nucleation points for complex enhancers to evolve.

Massively parallel disruption of enhancers active in human neural stem cells.

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

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.

Characterization of covalent inhibitors that disrupt the interaction between the tandem SH2 domains of SYK and FCER1G phospho-ITAM.

RNA sequencing and genetic data support spleen tyrosine kinase (SYK) and high affinity immunoglobulin epsilon receptor subunit gamma (FCER1G) as putative targets to be modulated for Alzheimer's disease (AD) therapy. FCER1G is a component of Fc receptor complexes that contain an immunoreceptor tyrosine-based activation motif (ITAM). SYK interacts with the Fc receptor by binding to doubly phosphorylated ITAM (p-ITAM) via its two tandem SH2 domains (SYK-tSH2). Interaction of the FCER1G p-ITAM with SYK-tSH2 enables SYK activation via phosphorylation. Since SYK activation is reported to exacerbate AD pathology, we hypothesized that disruption of this interaction would be beneficial for AD patients. Herein, we developed biochemical and biophysical assays to enable the discovery of small molecules that perturb the interaction between the FCER1G p-ITAM and SYK-tSH2. We identified two distinct chemotypes using a high-throughput screen (HTS) and orthogonally assessed their binding. Both chemotypes covalently modify SYK-tSH2 and inhibit its interaction with FCER1G p-ITAM, however, these compounds lack selectivity and this limits their utility as chemical tools.

Characterization of covalent inhibitors that disrupt the interaction between the tandem SH2 domains of SYK and FCER1G phospho-ITAM.

RNA sequencing and genetic data support spleen tyrosine kinase (SYK) and high affinity immunoglobulin epsilon receptor subunit gamma (FCER1G) as putative targets to be modulated for Alzheimer's disease (AD) therapy. FCER1G is a component of Fc receptor complexes that contain an immunoreceptor tyrosine-based activation motif (ITAM). SYK interacts with the Fc receptor by binding to doubly phosphorylated ITAM (p-ITAM) via its two tandem SH2 domains (SYK-tSH2). Interaction of the FCER1G p-ITAM with SYK-tSH2 enables SYK activation via phosphorylation. Since SYK activation is reported to exacerbate AD pathology, we hypothesized that disruption of this interaction would be beneficial for AD patients. Herein, we developed biochemical and biophysical assays to enable the discovery of small molecules that perturb the interaction between the FCER1G p-ITAM and SYK-tSH2. We identified two distinct chemotypes using a high-throughput screen (HTS) and orthogonally assessed their binding. Both chemotypes covalently modify SYK-tSH2 and inhibit its interaction with FCER1G p-ITAM.

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.

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.

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.

The female protective effect in autism spectrum disorder is not mediated by a single genetic locus.

BACKGROUND: A 4:1 male to female sex bias has consistently been observed in autism spectrum disorder (ASD). Epidemiological and genetic studies suggest a female protective effect (FPE) may account for part of this bias; however, the mechanism of such protection is unknown. Quantitative assessment of ASD symptoms using the Social Responsiveness Scale (SRS) shows a bimodal distribution unique to females in multiplex families. This leads to the hypothesis that a single, common genetic locus on chromosome X might mediate the FPE and produce the ASD sex bias. Such a locus would represent a major therapeutic target and is likely to have been missed by conventional genome-wide association study (GWAS) analysis. METHODS: To explore this possibility, we performed an association study in affected versus unaffected females, considering three tiers of single nucleotide polymorphisms (SNPs) as follows: 1) regions of chromosome X that escape X-inactivation, 2) all of chromosome X, and 3) genome-wide. RESULTS: No evidence of a SNP meeting the criteria for a single FPE locus was observed, despite the analysis being well powered to detect this effect. CONCLUSIONS: The results do not support the hypothesis that the FPE is mediated by a single genetic locus; however, this does not exclude the possibility of multiple genetic loci playing a role in the FPE.

De novo insertions and deletions of predominantly paternal origin are associated with autism spectrum disorder.

Whole-exome sequencing (WES) studies have demonstrated the contribution of de novo loss-of-function single-nucleotide variants (SNVs) to autism spectrum disorder (ASD). However, challenges in the reliable detection of de novo insertions and deletions (indels) have limited inclusion of these variants in prior analyses. By applying a robust indel detection method to WES data from 787 ASD families (2,963 individuals), we demonstrate that de novo frameshift indels contribute to ASD risk (OR = 1.6; 95% CI = 1.0-2.7; p = 0.03), are more common in female probands (p = 0.02), are enriched among genes encoding FMRP targets (p = 6 × 10(-9)), and arise predominantly on the paternal chromosome (p < 0.001). On the basis of mutation rates in probands versus unaffected siblings, we conclude that de novo frameshift indels contribute to risk in approximately 3% of individuals with ASD. Finally, by observing clustering of mutations in unrelated probands, we uncover two ASD-associated genes: KMT2E (MLL5), a chromatin regulator, and RIMS1, a regulator of synaptic vesicle release.