Novel TCF4:TCF12 heterodimer inhibits glioblastoma growth

TWIST1 (TW) is a pro-oncogenic basic helix-loop-helix (bHLH) transcription factor and promotes the hallmark features of malignancy (e.g., cell invasion, cancer cell stemness, and treatment resistance), which contribute to poor prognoses of glioblastoma (GBM). We previously reported that specific TW dimerization motifs regulate unique cellular phenotypes in GBM. For example, the TW:E12 heterodimer increases periostin (POSTN) expression and promotes cell invasion. TW dimer-specific transcriptional regulation requires binding to the regulatory E-box consensus sequences, but alternative bHLH dimers that balance TW dimer activity in regulating pro-oncogenic TW target genes are unknown. We leveraged the ENCODE DNase I hypersensitivity data to identify E-box sites and tethered TW:E12 and TW:TW proteins to validate dimer binding to E-boxes in vitro. Subsequently, TW knockdown revealed a novel TCF4:TCF12 bHLH dimer occupying the same TW E-box site that, when expressed as a tethered TCF4:TCF12 dimer, markedly repressed POSTN expression and extended animal survival. These observations support TCF4:TCF12 as a novel dimer with tumor-suppressor activity in GBM that functions in part through displacement of and/or competitive inhibition of pro-oncogenic TW dimers at E-box sites.

Aβ Amyloid Scaffolds the Accumulation of Matrisome and Additional Proteins in Alzheimer's Disease

We report a highly significant correlation in brain proteome changes between Alzheimers disease (AD) and CRND8 APP695NL/F transgenic mice. However, integrating protein changes observed in the CRND8 mice with co-expression networks derived from human AD, reveals both conserved and divergent module changes. For the most highly conserved module (M42, matrisome) we find many proteins accumulate in plaques, cerebrovascular amyloid (CAA), dystrophic processes, or a combination thereof. Overexpression of two M42 proteins, midkine (Mdk) and pleiotrophin (PTN), in CRND8 mice brains leads to increased accumulation of A β ; in plaques and in CAA; further, recombinant MDK and PTN enhance A β ; aggregation into amyloid. Multiple M42 proteins, annotated as heparan sulfate binding proteins, bind to fibrillar A β 42 and a non-human amyloid fibril in vitro. Supporting this binding data, MDK and PTN co-accumulate with transthyretin (TTR) amyloid in the heart and islet amyloid polypeptide (IAPP) amyloid in the pancreas. Our findings establish several critical insights. Proteomic changes in modules observed in human AD brains define an A β ; amyloid responsome that is well conserved from mouse model to human. Further, distinct amyloid structures may serve as scaffolds, facilitating the co-accumulation of proteins with signaling functions. We hypothesize that this co-accumulation may contribute to downstream pathological sequalae. Overall, this contextualized understanding of proteomic changes and their interplay with amyloid deposition provides valuable insights into the complexity of AD pathogenesis and potential biomarkers and therapeutic targets.

A metabolomic signature of the APOE2 allele

With the goal of identifying metabolites that significantly correlate with the protective e2 allele of the apolipoprotein E (APOE) gene, we established a consortium of five studies of healthy aging and extreme human longevity with 3545 participants. This consortium includes the New England Centenarian Study, the Baltimore Longitudinal Study of Aging, the Arivale study, the Longevity Genes Project/LonGenity studies, and the Long Life Family Study. We analyzed the association between APOE genotype groups E2 (e2e2 and e2e3 genotypes, N = 544), E3 (e3e3 genotypes, N = 2299), and E4 (e3e4 and e4e4 genotypes, N = 702) with metabolite profiles in the five studies and used fixed effect meta-analysis to aggregate the results. Our meta-analysis identified a signature of 19 metabolites that are significantly associated with the E2 genotype group at FDR < 10%. The group includes 10 glycerolipids and 4 glycerophospholipids that were all higher in E2 carriers compared to E3, with fold change ranging from 1.08 to 1.25. The organic acid 6-hydroxyindole sulfate, previously linked to changes in gut microbiome that were reflective of healthy aging and longevity, was also higher in E2 carriers compared to E3 carriers. Three sterol lipids and one sphingolipid species were significantly lower in carriers of the E2 genotype group. For some of these metabolites, the effect of the E2 genotype opposed the age effect. No metabolites reached a statistically significant association with the E4 group. This work confirms and expands previous results connecting the APOE gene to lipid regulation and suggests new links between the e2 allele, lipid metabolism, aging, and the gut-brain axis.

Risk factors for severe COVID-19 differ by age for hospitalized adults

Risk stratification for hospitalized adults with COVID-19 is essential to inform decisions about individual patients and allocation of resources. So far, risk models for severe COVID outcomes have included age but have not been optimized to best serve the needs of either older or younger adults. Models also need to be updated to reflect improvements in COVID-19 treatments. This retrospective study analyzed data from 6906 hospitalized adults with COVID-19 from a community health system across five states in the western United States. Risk models were developed to predict mechanical ventilation illness or death across one to 56 days of hospitalization, using clinical data available within the first hour after either admission with COVID-19 or a first positive SARS-CoV-2 test. For the seven-day interval, models for age ≥ 18 and < 50 years reached AUROC 0.81 (95% CI 0.71-0.91) and models for age ≥ 50 years reached AUROC 0.82 (95% CI 0.77-0.86). Models revealed differences in the statistical significance and relative predictive value of risk factors between older and younger patients including age, BMI, vital signs, and laboratory results. In addition, for hospitalized patients, sex and chronic comorbidities had lower predictive value than vital signs and laboratory results.

Manifestations of Alzheimer's disease genetic risk in the blood are evident in a multiomic analysis in healthy adults aged 18 to 90

Genetics play an important role in late-onset Alzheimer's Disease (AD) etiology and dozens of genetic variants have been implicated in AD risk through large-scale GWAS meta-analyses. However, the precise mechanistic effects of most of these variants have yet to be determined. Deeply phenotyped cohort data can reveal physiological changes associated with genetic risk for AD across an age spectrum that may provide clues to the biology of the disease. We utilized over 2000 high-quality quantitative measurements obtained from blood of 2831 cognitively normal adult clients of a consumer-based scientific wellness company, each with CLIA-certified whole-genome sequencing data. Measurements included: clinical laboratory blood tests, targeted chip-based proteomics, and metabolomics. We performed a phenome-wide association study utilizing this diverse blood marker data and 25 known AD genetic variants and an AD-specific polygenic risk score (PGRS), adjusting for sex, age, vendor (for clinical labs), and the first four genetic principal components; sex-SNP interactions were also assessed. We observed statistically significant SNP-analyte associations for five genetic variants after correction for multiple testing (for SNPs in or near NYAP1, ABCA7, INPP5D, and APOE), with effects detectable from early adulthood. The ABCA7 SNP and the APOE2 and APOE4 encoding alleles were associated with lipid variability, as seen in previous studies; in addition, six novel proteins were associated with the e2 allele. The most statistically significant finding was between the NYAP1 variant and PILRA and PILRB protein levels, supporting previous functional genomic studies in the identification of a putative causal variant within the PILRA gene. We did not observe associations between the PGRS and any analyte. Sex modified the effects of four genetic variants, with multiple interrelated immune-modulating effects associated with the PICALM variant. In post-hoc analysis, sex-stratified GWAS results from an independent AD case-control meta-analysis supported sex-specific disease effects of the PICALM variant, highlighting the importance of sex as a biological variable. Known AD genetic variation influenced lipid metabolism and immune response systems in a population of non-AD individuals, with associations observed from early adulthood onward. Further research is needed to determine whether and how these effects are implicated in early-stage biological pathways to AD. These analyses aim to complement ongoing work on the functional interpretation of AD-associated genetic variants.

The viral hypothesis: how herpesviruses may contribute to Alzheimer's disease

The hypothesis that infectious agents, particularly herpesviruses, contribute to Alzheimer's disease (AD) pathogenesis has been investigated for decades but has long engendered controversy. In the past 3 years, several studies in mouse models, human tissue models, and population cohorts have reignited interest in this hypothesis. Collectively, these studies suggest that many of the hallmarks of AD, like amyloid beta production and neuroinflammation, can arise as a protective response to acute infection that becomes maladaptive in the case of chronic infection. We place this work in its historical context and explore its etiological implications.

Systems modeling of metabolic dysregulation in neurodegenerative diseases

Neurodegenerative diseases (NDDs) encompass a wide range of conditions that arise owing to progressive degeneration and the ultimate loss of nerve cells in the brain and peripheral nervous system. NDDs such as Alzheimer's, Parkinson's, and Huntington's diseases negatively impact both length and quality of life, due to lack of effective disease-modifying treatments. Herein, we review the use of genome-scale metabolic models, network-based approaches, and integration with multiomics data to identify key biological processes that characterize NDDs. We describe powerful systems biology approaches for modeling NDD pathophysiology by leveraging in silico models that are informed by patient-derived multiomics data. These approaches can enable mechanistic insights into NDD-specific metabolic dysregulations that can be leveraged to identify potential metabolic markers of disease and predisease states.

Case Study: A Precision Medicine Approach to Multifactorial Dementia and Alzheimer's Disease

We report a case of a patient with mixed dementia successfully treated with a personalized multimodal therapy. Monotherapeutics are inadequate for the treatment of Alzheimer's disease (AD) and mixed dementia; therefore, we approach treatment through an adaptive personalized multimodal program. Many multimodal programs are pre-determined, and thus may not address the underlying contributors to cognitive decline in each particular individual. The combination of a targeted, personalized, precision medicine approach using a multimodal program promises advantages over monotherapies and untargeted multimodal therapies for multifactorial dementia. In this case study, we describe successful treatment for a patient diagnosed with AD, using a multimodal, programmatic, precision medicine intervention encompassing therapies targeting multiple dementia diastheses. We describe specific interventions used in this case that are derived from a comprehensive protocol for AD precision medicine. After treatment, our patient demonstrated improvements in quantitative neuropsychological testing, volumetric neuroimaging, PET scans, and serum chemistries, accompanied by symptomatic improvement over a 3.5-year period. This case outcome supports the need for rigorous trials of comprehensive, targeted combination therapies to stabilize, restore, and prevent cognitive decline in individuals with potentially many underlying causes of such decline and dementia. Our multimodal therapy included personalized treatments to address each potential perturbation to neuroplasticity. In particular, neuroinflammation and metabolic subsystems influence cognitive function and hippocampal volume. In this patient with a primary biliary cholangitis (PBC) multimorbidity component, we introduced a personalized diet that helped reduce liver inflammation. Together, all these components of multimodal therapy showed a sustained functional and cognitive benefit. Multimodal therapies may have systemwide benefits on all dementias, particularly in the context of multimorbidity. Furthermore, these therapies provide generalized health benefits, as many of the factors - such as inflammation - that impact cognitive function also impact other systems.

A systems-biology clinical trial of a personalized multimodal lifestyle intervention for early Alzheimer's disease

INTRODUCTION

There is an urgent need to develop effective interventional treatments for people with Alzheimer's disease (AD). AD results from a complex multi-decade interplay of multiple interacting dysfunctional biological systems that have not yet been fully elucidated. Epidemiological studies have linked several modifiable lifestyle factors with increased incidence for AD. Because monotherapies have failed to prevent or ameliorate AD, interventional studies should deploy multiple, targeted interventions that address the dysfunctional systems that give rise to AD.

METHODS

This randomized controlled trial (RCT) will examine the efficacy of a 12-month personalized, multimodal, lifestyle intervention in 60 mild cognitive impairment (MCI) and early stage AD patients (aged 50+, amyloid positivity). Both groups receive data-driven, lifestyle recommendations designed to target multiple systemic pathways implicated in AD. One group receives these personalized recommendations without coaching. The other group receives personalized recommendations with health coaching, dietary counseling, exercise training, cognitive stimulation, and nutritional supplements. We collect clinical, proteomic, metabolomic, neuroimaging, and genetic data to fuel systems-biology analyses. We will examine effects on cognition and hippocampal volume. The overarching goal of the study is to longitudinally track biological systems implicated in AD to reveal the dynamics between these systems during the intervention to understand differences in treatment response.

RESULTS

We have developed and implemented a protocol for a personalized, multimodal intervention program for early AD patients. We began enrollment in September 2019; we have enrolled a third of our target (20 of 60) with a 95% retention and 86% compliance rate.

DISCUSSION

This study presents a paradigm shift in designing multimodal, lifestyle interventions to reduce cognitive decline, and how to elucidate the biological systems being targeted. Analytical efforts to explain mechanistic or causal underpinnings of individual trajectories and the interplay between multi-omic variables will inform the design of future hypotheses and development of effective precision medicine trials.

Microglia show differential transcriptomic response to Aβ peptide aggregates ex vivo and in vivo

Aggregation and accumulation of amyloid-β (Aβ) is a defining feature of Alzheimer's disease pathology. To study microglial responses to Aβ, we applied exogenous Aβ peptide, in either oligomeric or fibrillar conformation, to primary mouse microglial cultures and evaluated system-level transcriptional changes and then compared these with transcriptomic changes in the brains of CRND8 APP mice. We find that primary microglial cultures have rapid and massive transcriptional change in response to Aβ. Transcriptomic responses to oligomeric or fibrillar Aβ in primary microglia, although partially overlapping, are distinct and are not recapitulated in vivo where Aβ progressively accumulates. Furthermore, although classic immune mediators show massive transcriptional changes in the primary microglial cultures, these changes are not observed in the mouse model. Together, these data extend previous studies which demonstrate that microglia responses ex vivo are poor proxies for in vivo responses. Finally, these data demonstrate the potential utility of using microglia as biosensors of different aggregate conformation, as the transcriptional responses to oligomeric and fibrillar Aβ can be distinguished.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.

Partial inhibition of mitochondrial complex I ameliorates Alzheimer's disease pathology and cognition in APP/PS1 female mice

Alzheimer's Disease (AD) is a devastating neurodegenerative disorder without a cure. Here we show that mitochondrial respiratory chain complex I is an important small molecule druggable target in AD. Partial inhibition of complex I triggers the AMP-activated protein kinase-dependent signaling network leading to neuroprotection in symptomatic APP/PS1 female mice, a translational model of AD. Treatment of symptomatic APP/PS1 mice with complex I inhibitor improved energy homeostasis, synaptic activity, long-term potentiation, dendritic spine maturation, cognitive function and proteostasis, and reduced oxidative stress and inflammation in brain and periphery, ultimately blocking the ongoing neurodegeneration. Therapeutic efficacy in vivo was monitored using translational biomarkers FDG-PET, 31P NMR, and metabolomics. Cross-validation of the mouse and the human transcriptomic data from the NIH Accelerating Medicines Partnership-AD database demonstrated that pathways improved by the treatment in APP/PS1 mice, including the immune system response and neurotransmission, represent mechanisms essential for therapeutic efficacy in AD patients.

Metabolic Network Analysis Reveals Altered Bile Acid Synthesis and Metabolism in Alzheimer's Disease

Increasing evidence suggests Alzheimer's disease (AD) pathophysiology is influenced by primary and secondary bile acids, the end product of cholesterol metabolism. We analyze 2,114 post-mortem brain transcriptomes and identify genes in the alternative bile acid synthesis pathway to be expressed in the brain. A targeted metabolomic analysis of primary and secondary bile acids measured from post-mortem brain samples of 111 individuals supports these results. Our metabolic network analysis suggests that taurine transport, bile acid synthesis, and cholesterol metabolism differ in AD and cognitively normal individuals. We also identify putative transcription factors regulating metabolic genes and influencing altered metabolism in AD. Intriguingly, some bile acids measured in brain tissue cannot be explained by the presence of enzymes responsible for their synthesis, suggesting that they may originate from the gut microbiome and are transported to the brain. These findings motivate further research into bile acid metabolism in AD to elucidate their possible connection to cognitive decline.

Meta-Analysis of the Alzheimer's Disease Human Brain Transcriptome and Functional Dissection in Mouse Models

We present a consensus atlas of the human brain transcriptome in Alzheimer's disease (AD), based on meta-analysis of differential gene expression in 2,114 postmortem samples. We discover 30 brain coexpression modules from seven regions as the major source of AD transcriptional perturbations. We next examine overlap with 251 brain differentially expressed gene sets from mouse models of AD and other neurodegenerative disorders. Human-mouse overlaps highlight responses to amyloid versus tau pathology and reveal age- and sex-dependent expression signatures for disease progression. Human coexpression modules enriched for neuronal and/or microglial genes broadly overlap with mouse models of AD, Huntington's disease, amyotrophic lateral sclerosis, and aging. Other human coexpression modules, including those implicated in proteostasis, are not activated in AD models but rather following other, unexpected genetic manipulations. Our results comprise a cross-species resource, highlighting transcriptional networks altered by human brain pathophysiology and identifying correspondences with mouse models for AD preclinical studies.

iREAD: a tool for intron retention detection from RNA-seq data

BACKGROUND

Intron retention (IR) has been traditionally overlooked as 'noise' and received negligible attention in the field of gene expression analysis. In recent years, IR has become an emerging field for interrogating transcriptomes because it has been recognized to carry out important biological functions such as gene expression regulation and it has been found to be associated with complex diseases such as cancers. However, methods for detecting IR today are limited. Thus, there is a need to develop novel methods to improve IR detection.

RESULTS

Here we present iREAD (intron REtention Analysis and Detector), a tool to detect IR events genome-wide from high-throughput RNA-seq data. The command line interface for iREAD is implemented in Python. iREAD takes as input a BAM file, representing the transcriptome, and a text file containing the intron coordinates of a genome. It then 1) counts all reads that overlap intron regions, 2) detects IR events by analyzing the features of reads such as depth and distribution patterns, and 3) outputs a list of retained introns into a tab-delimited text file. iREAD provides significant added value in detecting IR compared with output from IRFinder with a higher AUC on all datasets tested. Both methods showed low false positive rates and high false negative rates in different regimes, indicating that use together is generally beneficial. The output from iREAD can be directly used for further exploratory analysis such as differential intron expression and functional enrichment. The software is freely available at https://github.com/genemine/iread.

CONCLUSION

Being complementary to existing tools, iREAD provides a new and generic tool to interrogate poly-A enriched transcriptomic data of intron regions. Intron retention analysis provides a complementary approach for understanding transcriptome.

TLR5 decoy receptor as a novel anti-amyloid therapeutic for Alzheimer's disease

Chakrabarty et al. show that human TLR5 ectodomain reduces amyloid β (Aβ) plaques by direct interaction with Aβ, demonstrating the feasibility of such immune decoy receptor strategies as potential biotherapies in Alzheimer’s disease.

There is considerable interest in harnessing innate immunity to treat Alzheimer’s disease (AD). Here, we explore whether a decoy receptor strategy using the ectodomain of select TLRs has therapeutic potential in AD. AAV-mediated expression of human TLR5 ectodomain (sTLR5) alone or fused to human IgG4 Fc (sTLR5Fc) results in robust attenuation of amyloid β (Aβ) accumulation in a mouse model of Alzheimer-type Aβ pathology. sTLR5Fc binds to oligomeric and fibrillar Aβ with high affinity, forms complexes with Aβ, and blocks Aβ toxicity. Oligomeric and fibrillar Aβ modulates flagellin-mediated activation of human TLR5 but does not, by itself, activate TLR5 signaling. Genetic analysis shows that rare protein coding variants in human TLR5 may be associated with a reduced risk of AD. Further, transcriptome analysis shows altered TLR gene expression in human AD. Collectively, our data suggest that TLR5 decoy receptor–based biologics represent a novel and safe Aβ-selective class of biotherapy in AD.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

Targeting TWIST1 through loss of function inhibits tumorigenicity of human glioblastoma

TWIST1 (TW) is a bHLH transcription factor (TF) and master regulator of the epithelial-to-mesenchymal transition (EMT). In vitro, TW promotes mesenchymal change, invasion, and self-renewal in glioblastoma (GBM) cells. However, the potential therapeutic relevance of TW has not been established through loss-of-function studies in human GBM cell xenograft models. The effects of TW loss of function (gene editing and knockdown) on inhibition of tumorigenicity of U87MG and GBM4 glioma stem cells were tested in orthotopic xenograft models and conditional knockdown in established flank xenograft tumors. RNAseq and the analysis of tumors investigated putative TW-associated mechanisms. Multiple bioinformatic tools revealed significant alteration of ECM, membrane receptors, signaling transduction kinases, and cytoskeleton dynamics leading to identification of PI3K/AKT signaling. We experimentally show alteration of AKT activity and periostin (POSTN) expression in vivo and/or in vitro. For the first time, we show that effect of TW knockout inhibits AKT activity in U87MG cells in vivo independent of PTEN mutation. The clinical relevance of TW and candidate mechanisms was established by analysis of the TCGA and ENCODE databases. TW expression was associated with decreased patient survival and LASSO regression analysis identified POSTN as one of top targets of TW in human GBM. While we previously demonstrated the role of TW in promoting EMT and invasion of glioma cells, these studies provide direct experimental evidence supporting protumorigenic role of TW independent of invasion in vivo and the therapeutic relevance of targeting TW in human GBM. Further, the role of TW driving POSTN expression and AKT signaling suggests actionable targets, which could be leveraged to mitigate the oncogenic effects of TW in GBM.

Rare coding variants in PLCG2, ABI3, and TREM2 implicate microglial-mediated innate immunity in Alzheimer's disease

We identified rare coding variants associated with Alzheimer's disease in a three-stage case-control study of 85,133 subjects. In stage 1, we genotyped 34,174 samples using a whole-exome microarray. In stage 2, we tested associated variants (P < 1 × 10-4) in 35,962 independent samples using de novo genotyping and imputed genotypes. In stage 3, we used an additional 14,997 samples to test the most significant stage 2 associations (P < 5 × 10-8) using imputed genotypes. We observed three new genome-wide significant nonsynonymous variants associated with Alzheimer's disease: a protective variant in PLCG2 (rs72824905: p.Pro522Arg, P = 5.38 × 10-10, odds ratio (OR) = 0.68, minor allele frequency (MAF)cases = 0.0059, MAFcontrols = 0.0093), a risk variant in ABI3 (rs616338: p.Ser209Phe, P = 4.56 × 10-10, OR = 1.43, MAFcases = 0.011, MAFcontrols = 0.008), and a new genome-wide significant variant in TREM2 (rs143332484: p.Arg62His, P = 1.55 × 10-14, OR = 1.67, MAFcases = 0.0143, MAFcontrols = 0.0089), a known susceptibility gene for Alzheimer's disease. These protein-altering changes are in genes highly expressed in microglia and highlight an immune-related protein-protein interaction network enriched for previously identified risk genes in Alzheimer's disease. These genetic findings provide additional evidence that the microglia-mediated innate immune response contributes directly to the development of Alzheimer's disease.

A Cell-Surface Membrane Protein Signature for Glioblastoma

We present a systems strategy that facilitated the development of a molecular signature for glioblastoma (GBM), composed of 33 cell-surface transmembrane proteins. This molecular signature, GBMSig, was developed through the integration of cell-surface proteomics and transcriptomics from patient tumors in the REMBRANDT (n = 228) and TCGA datasets (n = 547) and can separate GBM patients from control individuals with a Matthew's correlation coefficient value of 0.87 in a lock-down test. Functionally, 17/33 GBMSig proteins are associated with transforming growth factor β signaling pathways, including CD47, SLC16A1, HMOX1, and MRC2. Knockdown of these genes impaired GBM invasion, reflecting their role in disease-perturbed changes in GBM. ELISA assays for a subset of GBMSig (CD44, VCAM1, HMOX1, and BIGH3) on 84 plasma specimens from multiple clinical sites revealed a high degree of separation of GBM patients from healthy control individuals (area under the curve is 0.98 in receiver operating characteristic). In addition, a classifier based on these four proteins differentiated the blood of pre- and post-tumor resections, demonstrating potential clinical value as biomarkers.

Ion channel expression patterns in glioblastoma stem cells with functional and therapeutic implications for malignancy

Ion channels and transporters have increasingly recognized roles in cancer progression through the regulation of cell proliferation, migration, and death. Glioblastoma stem-like cells (GSCs) are a source of tumor formation and recurrence in glioblastoma multiforme, a highly aggressive brain cancer, suggesting that ion channel expression may be perturbed in this population. However, little is known about the expression and functional relevance of ion channels that may contribute to GSC malignancy. Using RNA sequencing, we assessed the enrichment of ion channels in GSC isolates and non-tumor neural cell types. We identified a unique set of GSC-enriched ion channels using differential expression analysis that is also associated with distinct gene mutation signatures. In support of potential clinical relevance, expression of selected GSC-enriched ion channels evaluated in human glioblastoma databases of The Cancer Genome Atlas and Ivy Glioblastoma Atlas Project correlated with patient survival times. Finally, genetic knockdown as well as pharmacological inhibition of individual or classes of GSC-enriched ion channels constrained growth of GSCs compared to normal neural stem cells. This first-in-kind global examination characterizes ion channels enriched in GSCs and explores their potential clinical relevance to glioblastoma molecular subtypes, gene mutations, survival outcomes, regional tumor expression, and experimental responses to loss-of-function. Together, the data support the potential biological and therapeutic impact of ion channels on GSC malignancy and provide strong rationale for further examination of their mechanistic and therapeutic importance.

Motivational, proteostatic and transcriptional deficits precede synapse loss, gliosis and neurodegeneration in the B6.HttQ111/+ model of Huntington's disease

We investigated the appearance and progression of disease-relevant signs in the B6.Htt^Q111/+ mouse, a genetically precise model of the mutation that causes Huntington’s disease (HD). We find that B6.Htt^Q111/+ mice are healthy, show no overt signs of central or peripheral inflammation, and no gross motor impairment as late as 12 months of age. Behaviorally, we find that 4–9 month old B6.Htt^Q111/+ mice have normal activity levels and show no clear signs of anxiety or depression, but do show clear signs of reduced motivation. The neuronal density, neuronal size, synaptic density and number of glia is normal in B6.Htt^Q111/+ striatum, the most vulnerable brain region in HD, up to 12 months of age. Despite this preservation of the synaptic and cellular composition of the striatum, we observe clear progressive, striatal-specific transcriptional dysregulation and accumulation of neuronal intranuclear inclusions (NIIs). Simulation studies suggest these molecular endpoints are sufficiently robust for future preclinical studies, and that B6.Htt^Q111/+ mice are a useful tool for modeling disease-modifying or neuroprotective strategies for disease processes before the onset of overt phenotypes.

Isoforms of ADAR1 independently control MDA5-driven autoimmunity and multi-organ development

Mutations in the ADAR gene that encodes the ADAR1 RNA editing enzyme cause Aicardi-Goutières Syndrome (AGS), a severe autoimmune disease associated with an aberrant type I interferon response. How ADAR1 prevents autoimmunity remains incompletely defined. Here, we demonstrate that ADAR1 is a specific and essential negative regulator of the MDA5-MAVS RNA sensing pathway. Moreover, we uncovered a MDA5-MAVS-independent function for ADAR1 in the development of multiple organs. We showed that the p150 isoform of ADAR1 uniquely regulated the MDA5 pathway, whereas both the p150 and p110 isoforms contributed to development. Abrupt deletion of ADAR1 in adult mice revealed that both of these functions were required throughout life. Our findings delineate genetically distinct roles for both ADAR1 isoforms in vivo, with implications for the human diseases caused by ADAR mutations.

Isoforms of RNA-Editing Enzyme ADAR1 Independently Control Nucleic Acid Sensor MDA5-Driven Autoimmunity and Multi-organ Development

Mutations in ADAR, which encodes the ADAR1 RNA-editing enzyme, cause Aicardi-Goutières syndrome (AGS), a severe autoimmune disease associated with an aberrant type I interferon response. How ADAR1 prevents autoimmunity remains incompletely defined. Here, we demonstrate that ADAR1 is a specific and essential negative regulator of the MDA5-MAVS RNA sensing pathway. Moreover, we uncovered a MDA5-MAVS-independent function for ADAR1 in the development of multiple organs. We showed that the p150 isoform of ADAR1 uniquely regulated the MDA5 pathway, whereas both the p150 and p110 isoforms contributed to development. Abrupt deletion of ADAR1 in adult mice revealed that both of these functions were required throughout life. Our findings delineate genetically separable roles for both ADAR1 isoforms in vivo, with implications for the human diseases caused by ADAR mutations.

SNAPR: a bioinformatics pipeline for efficient and accurate RNA-seq alignment and analysis

The process of converting raw RNA sequencing data to interpretable results can be circuitous and time consuming, requiring multiple steps. We present an RNA-seq mapping algorithm that streamlines this process. Our algorithm utilizes a hash table approach to leverage the availability and power of high memory machines. SNAPR, which can be run on a single library or thousands of libraries, can take compressed or uncompressed FASTQ and BAM files as inputs, and can output a sorted BAM file, individual read counts, gene fusions and identify exogenous RNA species in a single step. SNAPR also does native Phred score filtering of reads. SNAPR is also well suited for future sequencing platforms that generate longer reads. Using SNAPR, we show how we can analyze data from hundreds of TCGA samples in a matter of hours, while identifying gene fusions and viral events at the same time. With the references genome and transcriptome undergoing periodic updates, and the need for uniform parameters when integrating multiple data sets, there is great need for a streamlined process for RNA-seq analysis. We demonstrate how SNAPR does this efficiently and accurately, with the high-throughput capacity needed to do high-volume analyses.

Fostering synergy between cell biology and systems biology

In the shared pursuit of elucidating detailed mechanisms of cell function, systems biology presents a natural complement to ongoing efforts in cell biology. Systems biology aims to characterize biological systems through integrated and quantitative modeling of cellular information. The process of model building and analysis provides value through synthesizing and cataloging information about cells and molecules, predicting mechanisms and identifying generalizable themes, generating hypotheses and guiding experimental design, and highlighting knowledge gaps and refining understanding. In turn, incorporating domain expertise and experimental data is crucial for building towards whole cell models. An iterative cycle of interaction between cell and systems biologists advances the goals of both fields and establishes a framework for mechanistic understanding of the genome-to-phenome relationship.

IL-10 alters immunoproteostasis in APP mice, increasing plaque burden and worsening cognitive behavior

Anti-inflammatory strategies are proposed to have beneficial effects in Alzheimer's disease. To explore how anti-inflammatory cytokine signaling affects Aβ pathology, we investigated the effects of adeno-associated virus (AAV2/1)-mediated expression of Interleukin (IL)-10 in the brains of APP transgenic mouse models. IL-10 expression resulted in increased Aβ accumulation and impaired memory in APP mice. A focused transcriptome analysis revealed changes consistent with enhanced IL-10 signaling and increased ApoE expression in IL-10-expressing APP mice. ApoE protein was selectively increased in the plaque-associated insoluble cellular fraction, likely because of direct interaction with aggregated Aβ in the IL-10-expressing APP mice. Ex vivo studies also show that IL-10 and ApoE can individually impair glial Aβ phagocytosis. Our observations that IL-10 has an unexpected negative effect on Aβ proteostasis and cognition in APP mouse models demonstrate the complex interplay between innate immunity and proteostasis in neurodegenerative diseases, an interaction we call immunoproteostasis.

Quantitative proteomic analysis reveals effects of epidermal growth factor receptor (EGFR) on invasion-promoting proteins secreted by glioblastoma cells

Glioblastoma multiforme is a highly invasive and aggressive brain tumor with an invariably poor prognosis. The overexpression of epidermal growth factor receptor (EGFR) is a primary influencer of invasion and proliferation in tumor cells and the constitutively active EGFRvIII mutant, found in 30-65% of Glioblastoma multiforme, confers more aggressive invasion. To better understand how EGFR contributes to tumor aggressiveness, we investigated the effect of EGFR on the secreted levels of 65 rationally selected proteins involved in invasion. We employed selected reaction monitoring targeted mass spectrometry using stable isotope labeled internal peptide standards to quantity proteins in the secretome from five GBM (U87) isogenic cell lines in which EGFR, EGFRvIII, and/or PTEN were expressed. Our results show that cell lines with EGFR overexpression and constitutive EGFRvIII expression differ remarkably in the expression profiles for both secreted and intracellular signaling proteins, and alterations in EGFR signaling result in reproducible changes in concentrations of secreted proteins. Furthermore, the EGFRvIII-expressing mutant cell line secretes the majority of the selected invasion-promoting proteins at higher levels than other cell lines tested. Additionally, the intracellular and extracellular protein measurements indicate elevated oxidative stress in the EGFRvIII-expressing cell line. In conclusion, the results of our study demonstrate that EGFR signaling has a significant effect on the levels of secreted invasion-promoting proteins, likely contributing to the aggressiveness of Glioblastoma multiforme. Further characterization of these proteins may provide candidates for new therapeutic strategies and targets as well as biomarkers for this aggressive disease.

Detection of 1α,25-dihydroxyvitamin D-regulated miRNAs in zebrafish by whole transcriptome sequencing

The sterol hormone, 1α,25-dihydroxyvitamin D₃ (1α,25(OH)₂D₃), regulates gene expression and messenger RNA (mRNA) concentrations in zebrafish in vivo. Since mRNA concentrations and translation are influenced by micro-RNAs (miRNAs), we examined the influence of 1α,25(OH)₂D₃ on miRNA expression in zebrafish in vivo with whole transcriptome RNA sequencing, searched for miRNA binding sites in 1α,25(OH)₂D₃-sensitive genes, and performed correlation analyses between 1α,25(OH)₂D₃-sensitive miRNAs and mRNAs. In vehicle- and 1α,25(OH)₂D₃-treated, 7-day postfertilization larvae, between 282 and 295 known precursor miRNAs were expressed, and in vehicle- and 1α,25(OH)₂D₃-treated fish, between 83 and 122 novel miRNAs were detected. Following 1α,25(OH)₂D₃ treatment, 31 precursor miRNAs were differentially expressed (p<0.05). The differentially expressed miRNAs are predicted to potentially alter mRNAs for metabolic enzymes, transcription factors, growth factors, and Jak-STAT signaling. We verified the role of a 1α,25(OH)₂D₃-sensitive miRNA, miR125b, by demonstrating alterations in the concentrations of the mRNA of a 1α,25(OH)₂D₃-regulated gene, Cyp24a1, following transfection of renal cells with a miR125b miRNA mimic. Changes in the Cyp24a1 mRNA concentration by the miR125b miRNA mimic were associated with changes in the protein for Cyp24a1. Our data show that 1α,25(OH)₂D₃ regulates miRNA in zebrafish larvae in vivo and could thereby influence vitamin D-sensitive mRNA concentrations.

Transcriptional analysis of aggressiveness and heterogeneity across grades of astrocytomas

Astrocytoma is the most common glioma, accounting for half of all primary brain and spinal cord tumors. Late detection and the aggressive nature of high-grade astrocytomas contribute to high mortality rates. Though many studies identify candidate biomarkers using high-throughput transcriptomic profiling to stratify grades and subtypes, few have resulted in clinically actionable results. This shortcoming can be attributed, in part, to pronounced lab effects that reduce signature robustness and varied individual gene expression among patients with the same tumor. We addressed these issues by uniformly preprocessing publicly available transcriptomic data, comprising 306 tumor samples from three astrocytoma grades (Grade 2, 3, and 4) and 30 non-tumor samples (normal brain as control tissues). Utilizing Differential Rank Conservation (DIRAC), a network-based classification approach, we examined the global and individual patterns of network regulation across tumor grades. Additionally, we applied gene-based approaches to identify genes whose expression changed consistently with increasing tumor grade and evaluated their robustness across multiple studies using statistical sampling. Applying DIRAC, we observed a global trend of greater network dysregulation with increasing tumor aggressiveness. Individual networks displaying greater differences in regulation between adjacent grades play well-known roles in calcium/PKC, EGF, and transcription signaling. Interestingly, many of the 90 individual genes found to monotonically increase or decrease with astrocytoma grade are implicated in cancer-affected processes such as calcium signaling, mitochondrial metabolism, and apoptosis. The fact that specific genes monotonically increase or decrease with increasing astrocytoma grade may reflect shared oncogenic mechanisms among phenotypically similar tumors. This work presents statistically significant results that enable better characterization of different human astrocytoma grades and hopefully can contribute towards improvements in diagnosis and therapy choices. Our results also identify a number of testable hypotheses relating to astrocytoma etiology that may prove helpful in developing much-needed biomarkers for earlier disease detection.

Multi-study integration of brain cancer transcriptomes reveals organ-level molecular signatures

We utilized abundant transcriptomic data for the primary classes of brain cancers to study the feasibility of separating all of these diseases simultaneously based on molecular data alone. These signatures were based on a new method reported herein – Identification of Structured Signatures and Classifiers (ISSAC) – that resulted in a brain cancer marker panel of 44 unique genes. Many of these genes have established relevance to the brain cancers examined herein, with others having known roles in cancer biology. Analyses on large-scale data from multiple sources must deal with significant challenges associated with heterogeneity between different published studies, for it was observed that the variation among individual studies often had a larger effect on the transcriptome than did phenotype differences, as is typical. For this reason, we restricted ourselves to studying only cases where we had at least two independent studies performed for each phenotype, and also reprocessed all the raw data from the studies using a unified pre-processing pipeline. We found that learning signatures across multiple datasets greatly enhanced reproducibility and accuracy in predictive performance on truly independent validation sets, even when keeping the size of the training set the same. This was most likely due to the meta-signature encompassing more of the heterogeneity across different sources and conditions, while amplifying signal from the repeated global characteristics of the phenotype. When molecular signatures of brain cancers were constructed from all currently available microarray data, 90% phenotype prediction accuracy, or the accuracy of identifying a particular brain cancer from the background of all phenotypes, was found. Looking forward, we discuss our approach in the context of the eventual development of organ-specific molecular signatures from peripheral fluids such as the blood.

From a multi-study, integrated transcriptomic dataset, we identified a marker panel for differentiating major human brain cancers at the gene-expression level. The ISSAC molecular signatures for brain cancers, composed of 44 unique genes, are based on comparing expression levels of pairs of genes, and phenotype prediction follows a diagnostic hierarchy. We found that sufficient dataset integration across multiple studies greatly enhanced diagnostic performance on truly independent validation sets, whereas signatures learned from only one dataset typically led to high error rate. Molecular signatures of brain cancers, when obtained using all currently available gene-expression data, achieved 90% phenotype prediction accuracy. Thus, our integrative approach holds significant promise for developing organ-level, comprehensive, molecular signatures of disease.

AUREA: an open-source software system for accurate and user-friendly identification of relative expression molecular signatures

Background

Public databases such as the NCBI Gene Expression Omnibus contain extensive and exponentially increasing amounts of high-throughput data that can be applied to molecular phenotype characterization. Collectively, these data can be analyzed for such purposes as disease diagnosis or phenotype classification. One family of algorithms that has proven useful for disease classification is based on relative expression analysis and includes the Top-Scoring Pair (TSP), k-Top-Scoring Pairs (k-TSP), Top-Scoring Triplet (TST) and Differential Rank Conservation (DIRAC) algorithms. These relative expression analysis algorithms hold significant advantages for identifying interpretable molecular signatures for disease classification, and have been implemented previously on a variety of computational platforms with varying degrees of usability. To increase the user-base and maximize the utility of these methods, we developed the program AUREA (Adaptive Unified Relative Expression Analyzer)—a cross-platform tool that has a consistent application programming interface (API), an easy-to-use graphical user interface (GUI), fast running times and automated parameter discovery.

Results

Herein, we describe AUREA, an efficient, cohesive, and user-friendly open-source software system that comprises a suite of methods for relative expression analysis. AUREA incorporates existing methods, while extending their capabilities and bringing uniformity to their interfaces. We demonstrate that combining these algorithms and adaptively tuning parameters on the training sets makes these algorithms more consistent in their performance and demonstrate the effectiveness of our adaptive parameter tuner by comparing accuracy across diverse datasets.

Conclusions

We have integrated several relative expression analysis algorithms and provided a unified interface for their implementation while making data acquisition, parameter fixing, data merging, and results analysis ‘point-and-click’ simple. The unified interface and the adaptive parameter tuning of AUREA provide an effective framework in which to investigate the massive amounts of publically available data by both ‘in silico’ and ‘bench’ scientists. AUREA can be found at http://price.systemsbiology.net/AUREA/.

The estrogen-regulated transcription factor PITX1 coordinates gene-specific regulation by estrogen receptor-alpha in breast cancer cells

The estrogen receptor α (ERα) is a master regulator of gene expression and works along with cooperating transcription factors in mediating the actions of the hormone estradiol (E2) in ER-positive tissues and breast tumors. Here, we report that expression of paired-like homeodomain transcription factor (PITX1), a tumor suppressor and member of the homeobox family of transcription factors, is robustly up-regulated by E2 in several ERα-positive breast cancer cell lines via ERα-dependent interaction between the proximal promoter and an enhancer region 5' upstream of the PITX1 gene. Overexpression of PITX1 selectively inhibited the transcriptional activity of ERα and ERβ, while enhancing the activities of the glucocorticoid receptor and progesterone receptor. Reduction of PITX1 by small interfering RNA enhanced ERα-dependent transcriptional regulation of a subset of ERα target genes. The consensus PITX1 binding motif was found to be present in 28% of genome-wide ERα binding sites and was in close proximity to estrogen response elements in a subset of ERα binding sites, and E2 treatment enhanced PITX1 as well as ERα recruitment to these binding sites. These studies identify PITX1 as a new ERα transcriptional target that acts as a repressor to coordinate and fine tune target-specific, ERα-mediated transcriptional activity in human breast cancer cells.

Systems approaches to molecular cancer diagnostics

The search for improved molecular cancer diagnostics is a challenge for which systems approaches show great promise. As is becoming increasingly clear, cancer is a perpetually-evolving, highly multi-factorial disease. With next generation sequencing providing an ever-increasing amount of high-throughput data, the need for analytical tools that can provide meaningful context is critical. Systems approaches have demonstrated an ability to separate meaningful signal from noise that arises from population heterogeneity, heterogeneity within and across tumors, and multiple sources of technical variation when sufficient sample sizes are obtained and standardized measurement technologies are used. The ability to develop clinically useful molecular cancer diagnostics will be predicated on advancements on two major fronts: 1) more comprehensive and accurate measurements of multiple endpoints, and 2) more sophisticated analytical tools that synthesize high-throughput data into meaningful reflections of cellular states. To this end, systems approaches that have integrated transcriptomic data onto biomolecular networks have shown promise in their ability to classify tumor subtypes, predict clinical progression, and inform treatment options. Ultimately, the success of systems approaches will be measured by their ability to develop molecular cancer diagnostics through distilling complex, systems-wide information into actionable information in the clinic.

Estrogen Dendrimer Conjugates that Preferentially Activate Extranuclear, Nongenomic Versus Genomic Pathways of Estrogen Action

Estrogenic hormones are classically thought to exert their effects by binding to nuclear estrogen receptors and altering target gene transcription, but estrogens can also have nongenomic effects through rapid activation of membrane-initiated kinase cascades. The development of ligands that selectively activate only the nongenomic pathways would provide useful tools to investigate the significance of these pathways. We have prepared large, abiotic, nondegradable poly(amido)amine dendrimer macromolecules that are conjugated to multiple estrogen molecules through chemically robust linkages. Because of their charge and size, these estrogen-dendrimer conjugates (EDCs) remain outside the nucleus. They stimulate ERK, Shc, and Src phosphorylation in MCF-7 breast cancer cells at low concentrations, yet they are very ineffective in stimulating transcription of endogenous estrogen target genes, being approximately 10,000-fold less potent than estradiol in genomic actions. In contrast to estradiol, EDC was not effective in stimulating breast cancer cell proliferation. Because these EDC ligands activate nongenomic activity at concentrations at which they do not alter the transcription of estrogen target genes, they should be useful in studying extranuclear initiated pathways of estrogen action in a variety of target cells.

Estrogen down-regulation of the corepressor N-CoR: mechanism and implications for estrogen derepression of N-CoR-regulated genes

The nuclear receptor corepressor N-CoR plays a crucial role in the repressive activity of diverse transcription factors, yet little is known about what regulates its cellular level. We have found that estrogen markedly down-regulates N-CoR protein levels in estrogen receptor (ER)-positive breast cancer cells without affecting N-CoR mRNA levels, whereas levels of the related corepressor SMRT are unaffected. This effect is attributable to estrogen up-regulation of the ubiquitin ligase Siah2, which is a rapid and primary transcriptional response mediated by the ER, and precedes the loss of N-CoR. Treatment with proteasomal inhibitor or with small interfering RNA against Siah2 prevented the down-regulation of N-CoR by estrogen. Furthermore, the expression of 24-hydroxylase, a gene repressed by unliganded vitamin D receptor through its interaction with N-CoR, was up-regulated by estrogen and required Siah2. Our results illustrate a mechanism by which the estrogen-ER complex markedly reduces the level of N-CoR through a process involving the up-regulation of Siah2 and the subsequent targeting of N-CoR for proteasomal degradation. These findings reveal that, although estrogen directly regulates the transcription of many genes, by regulating a gene such as Siah2 it can exert profound "secondary" effects on cellular activity through mechanisms such as targeting regulatory proteins for degradation. This estrogen-evoked down-regulation of N-CoR could have a global derepressive effect on genes whose repression depends on N-CoR and thereby have broad impact on the activity of transcription factors and nuclear receptors whose actions involve N-CoR.