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.

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.

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.

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.

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.

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.

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.

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/.

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.

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.

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.

Multi-Omic analyses characterize the ceramide/sphingomyelin pathway as a therapeutic target in Alzheimer's disease.

Dysregulation of sphingomyelin and ceramide metabolism have been implicated in Alzheimer's disease. Genome-wide and transcriptome-wide association studies have identified various genes and genetic variants in lipid metabolism that are associated with Alzheimer's disease. However, the molecular mechanisms of sphingomyelin and ceramide disruption remain to be determined. We focus on the sphingolipid pathway and carry out multi-omics analyses to identify central and peripheral metabolic changes in Alzheimer's patients, correlating them to imaging features. Our multi-omics approach is based on (a) 2114 human post-mortem brain transcriptomics to identify differentially expressed genes; (b) in silico metabolic flux analysis on context-specific metabolic networks identified differential reaction fluxes; (c) multimodal neuroimaging analysis on 1576 participants to associate genetic variants in sphingomyelin pathway with Alzheimer's disease pathogenesis; (d) plasma metabolomic and lipidomic analysis to identify associations of lipid species with dysregulation in Alzheimer's; and (e) metabolite genome-wide association studies to define receptors within the pathway as a potential drug target. We validate our hypothesis in amyloidogenic APP/PS1 mice and show prolonged exposure to fingolimod alleviated synaptic plasticity and cognitive impairment in mice. Our integrative multi-omics approach identifies potential targets in the sphingomyelin pathway and suggests modulators of S1P metabolism as possible candidates for Alzheimer's disease treatment.

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.

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.

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.

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.

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.