The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

Massive mining of publicly available RNA-seq data from human and mouse

RNA sequencing (RNA-seq) is the leading technology for genome-wide transcript quantification. However, publicly available RNA-seq data is currently provided mostly in raw form, a significant barrier for global and integrative retrospective analyses. ARCHS4 is a web resource that makes the majority of published RNA-seq data from human and mouse available at the gene and transcript levels. For developing ARCHS4, available FASTQ files from RNA-seq experiments from the Gene Expression Omnibus (GEO) were aligned using a cloud-based infrastructure. In total 187,946 samples are accessible through ARCHS4 with 103,083 mouse and 84,863 human. Additionally, the ARCHS4 web interface provides intuitive exploration of the processed data through querying tools, interactive visualization, and gene pages that provide average expression across cell lines and tissues, top co-expressed genes for each gene, and predicted biological functions and protein-protein interactions for each gene based on prior knowledge combined with co-expression.

lncRNA-based study of epigenetic regulations in diabetic peripheral neuropathy

Diabetes remains one of the most prevalent non-communicable diseases in the world, affecting over 400 million of people worldwide, causing serious complications leading to amputations and even death. Over the years, researchers have found that, in addition to genomic mutations, epigenetic mechanisms also play a role in the development of diabetes-specifically type-2 diabetes. Long noncoding RNAs (lncRNAs) have been linked to mediate epigenetic mechanisms, including those in late-stage diabetes. This study attempts to assess the unexplored topic of how lncRNAs could be used to assess the epigenetic mechanisms present in diabetic peripheral neuropathy (DPN); a serious complication of the disease often leading to amputation. Differential lncRNA expression analysis was done with a dataset containing DPN and healthy patients. Standard and corrected t test, and also LIMMA was applied. Results of this study indicates the usefulness of lncRNAs as an exploratory tool to elucidate the complexity of the epigenetic mechanisms of human DPN.

Multi-omics analysis reveals neoantigen-independent immune cell infiltration in copy-number driven cancers

To realize the full potential of immunotherapy, it is critical to understand the drivers of tumor infiltration by immune cells. Previous studies have linked immune infiltration with tumor neoantigen levels, but the broad applicability of this concept remains unknown. Here, we find that while this observation is true across cancers characterized by recurrent mutations, it does not hold for cancers driven by recurrent copy number alterations, such as breast and pancreatic tumors. To understand immune invasion in these cancers, we developed an integrative multi-omics framework, identifying the DNA damage response protein ATM as a driver of cytokine production leading to increased immune infiltration. This prediction was validated in numerous orthogonal datasets, as well as experimentally in vitro and in vivo by cytokine release and immune cell migration. These findings demonstrate diverse drivers of immune cell infiltration across cancer lineages and may facilitate the clinical adaption of immunotherapies across diverse malignancies.

Ikaros family zinc finger 1 regulates dendritic cell development and function in humans

Ikaros family zinc finger 1 (IKZF1) is a haematopoietic transcription factor required for mammalian B-cell development. IKZF1 deficiency also reduces plasmacytoid dendritic cell (pDC) numbers in mice, but its effects on human DC development are unknown. Here we show that heterozygous mutation of IKZF1 in human decreases pDC numbers and expands conventional DC1 (cDC1). Lenalidomide, a drug that induces proteosomal degradation of IKZF1, also decreases pDC numbers in vivo, and reduces the ratio of pDC/cDC1 differentiated from progenitor cells in vitro in a dose-dependent manner. In addition, non-classical monocytes are reduced by IKZF1 deficiency in vivo. DC and monocytes from patients with IKZF1 deficiency or lenalidomide-treated cultures secrete less IFN-α, TNF and IL-12. These results indicate that human DC development and function are regulated by IKZF1, providing further insights into the consequences of IKZF1 mutation on immune function and the mechanism of immunomodulation by lenalidomide.

IKZF1 is a transcription factor known to regulate mammalian B-cell development. Here the authors show that IKZF1 is required for human pDC development and regulation of DC cytokine production in patients with IKZF1 haploinsufficiency, findings which are recapitulated in lenalidomide-induced IKZF1 deficiency.

Hydrophilic Auristatin Glycoside Payload Enables Improved Antibody-Drug Conjugate Efficacy and Biocompatibility

Antibody-drug conjugates (ADCs) offer a combination of antibody therapy and specific delivery of potent small-molecule payloads to target cells. The properties of the ADC molecule are determined by the balance of its components. The efficacy of the payload component increases with higher drug-to-antibody ratio (DAR), while homogeneous DAR = 8 ADCs are easily prepared by conjugation to the four accessible antibody hinge cystines. However, use of hydrophobic payloads has permitted only DAR = 2-4, due to poor pharmacokinetics and aggregation problems. Here, we describe generation and characterization of homogeneous DAR = 8 ADCs carrying a novel auristatin β-D-glucuronide, MMAU. The glycoside payload contributed to overall hydrophilicity of the ADC reducing aggregation. Compared to standard DAR = 2-4 ADCs, cytotoxicity of the homogeneous DAR = 8 ADCs was improved to low-picomolar IC₅₀ values against cancer cells in vitro. Bystander efficacy was restored after ADC internalization and subsequent cleavage of the glycoside, although unconjugated MMAU was relatively non-toxic to cells. DAR = 8 MMAU ADCs were effective against target antigen-expressing xenograft tumors. The ADCs were also studied in 3D in vitro patient-derived xenograft (PDX) assays where they outperformed clinically used ADC. In conclusion, increased hydrophilicity of the payload contributed to the ADC's hydrophilicity, stability and safety to non-target cells, while significantly improving cytotoxicity and enabling bystander efficacy.

Mutations in THAP1/DYT6 reveal that diverse dystonia genes disrupt similar neuronal pathways and functions

Dystonia is characterized by involuntary muscle contractions. Its many forms are genetically, phenotypically and etiologically diverse and it is unknown whether their pathogenesis converges on shared pathways. Mutations in THAP1 [THAP (Thanatos-associated protein) domain containing, apoptosis associated protein 1], a ubiquitously expressed transcription factor with DNA binding and protein-interaction domains, cause dystonia, DYT6. There is a unique, neuronal 50-kDa Thap1-like immunoreactive species, and Thap1 levels are auto-regulated on the mRNA level. However, THAP1 downstream targets in neurons, and the mechanism via which it causes dystonia are largely unknown. We used RNA-Seq to assay the in vivo effect of a heterozygote Thap1 C54Y or ΔExon2 allele on the gene transcription signatures in neonatal mouse striatum and cerebellum. Enriched pathways and gene ontology terms include eIF2α Signaling, Mitochondrial Dysfunction, Neuron Projection Development, Axonal Guidance Signaling, and Synaptic LongTerm Depression, which are dysregulated in a genotype and tissue-dependent manner. Electrophysiological and neurite outgrowth assays were consistent with those enrichments, and the plasticity defects were partially corrected by salubrinal. Notably, several of these pathways were recently implicated in other forms of inherited dystonia, including DYT1. We conclude that dysfunction of these pathways may represent a point of convergence in the pathophysiology of several forms of inherited dystonia.

Dystonia is a brain disorder that causes disabling involuntary muscle contractions and abnormal postures. Mutations in THAP1, a zinc-finger transcription factor, cause DYT6, but its neuronal targets and functions are unknown. In this study, we sought to determine the effects of Thap1C54Y and ΔExon2 alleles on the gene transcription signatures at postnatal day 1 (P1) in the mouse striatum and cerebellum in order to correlate function with specific genes or pathways. Our unbiased transcriptomics approach showed that Thap1 mutants revealed multiple signaling pathways involved in neuronal plasticity, axonal guidance, and oxidative stress response, which are also present in other forms of dystonia, particularly DYT1. We conclude that dysfunction of these pathways may represent a point of convergence on the pathogenesis of unrelated forms of inherited dystonia.

The Library of Integrated Network-Based Cellular Signatures NIH Program: System-Level Cataloging of Human Cells Response to Perturbations

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program that catalogs how human cells globally respond to chemical, genetic, and disease perturbations. Resources generated by LINCS include experimental and computational methods, visualization tools, molecular and imaging data, and signatures. By assembling an integrated picture of the range of responses of human cells exposed to many perturbations, the LINCS program aims to better understand human disease and to advance the development of new therapies. Perturbations under study include drugs, genetic perturbations, tissue micro-environments, antibodies, and disease-causing mutations. Responses to perturbations are measured by transcript profiling, mass spectrometry, cell imaging, and biochemical methods, among other assays. The LINCS program focuses on cellular physiology shared among tissues and cell types relevant to an array of diseases, including cancer, heart disease, and neurodegenerative disorders. This Perspective describes LINCS technologies, datasets, tools, and approaches to data accessibility and reusability.

Profiling and Co-expression Network Analysis of Learned Helplessness Regulated mRNAs and lncRNAs in the Mouse Hippocampus

Although studies provide insights into the neurobiology of stress and depression, the exact molecular mechanisms underlying their pathologies remain largely unknown. Long non-coding RNA (lncRNA) has been implicated in brain functions and behavior. A potential link between lncRNA and psychiatric disorders has been proposed. However, it remains undetermined whether IncRNA regulation, in the brain, contributes to stress or depression pathologies. In this study, we used a valid animal model of depression-like symptoms; namely learned helplessness, RNA-seq, Gene Ontology and co-expression network analyses to profile the expression pattern of lncRNA and mRNA in the hippocampus of mice. We identified 6346 differentially expressed transcripts. Among them, 340 lncRNAs and 3559 protein coding mRNAs were differentially expressed in helpless mice in comparison with control and/or non-helpless mice (inescapable stress resilient mice). Gene Ontology and pathway enrichment analyses indicated that induction of helplessness altered expression of mRNAs enriched in fundamental biological functions implicated in stress/depression neurobiology such as synaptic, metabolic, cell survival and proliferation, developmental and chromatin modification functions. To explore the possible regulatory roles of the altered lncRNAs, we constructed co-expression networks composed of the lncRNAs and mRNAs. Among our differentially expressed lncRNAs, 17% showed significant correlation with genes. Functional co-expression analysis linked the identified lncRNAs to several cellular mechanisms implicated in stress/depression neurobiology. Importantly, 57% of the identified regulatory lncRNAs significantly correlated with 18 different synapse-related functions. Thus, the current study identifies for the first time distinct groups of lncRNAs regulated by induction of learned helplessness in the mouse brain. Our results suggest that lncRNA-directed regulatory mechanisms might contribute to stress-induced pathologies; in particular, to inescapable stress-induced synaptic modifications.

(Pro)renin Receptor Inhibition Reprograms Hepatic Lipid Metabolism and Protects Mice From Diet-Induced Obesity and Hepatosteatosis

RATIONALE

An elevated level of plasma LDL (low-density lipoprotein) is an established risk factor for cardiovascular disease. Recently, we reported that the (pro)renin receptor ([P]RR) regulates LDL metabolism in vitro via the LDLR (LDL receptor) and SORT1 (sortilin-1), independently of the renin-angiotensin system.

OBJECTIVES

To investigate the physiological role of (P)RR in lipid metabolism in vivo.

METHODS AND RESULTS

We used N-acetylgalactosamine modified antisense oligonucleotides to specifically inhibit hepatic (P)RR expression in C57BL/6 mice and studied the consequences this has on lipid metabolism. In line with our earlier report, hepatic (P)RR silencing increased plasma LDL-C (LDL cholesterol). Unexpectedly, this also resulted in markedly reduced plasma triglycerides in a SORT1-independent manner in C57BL/6 mice fed a normal- or high-fat diet. In LDLR-deficient mice, hepatic (P)RR inhibition reduced both plasma cholesterol and triglycerides, in a diet-independent manner. Mechanistically, we found that (P)RR inhibition decreased protein abundance of ACC (acetyl-CoA carboxylase) and PDH (pyruvate dehydrogenase). This alteration reprograms hepatic metabolism, leading to reduced lipid synthesis and increased fatty acid oxidation. As a result, hepatic (P)RR inhibition attenuated diet-induced obesity and hepatosteatosis.

CONCLUSIONS

Collectively, our study suggests that (P)RR plays a key role in energy homeostasis and regulation of plasma lipids by integrating hepatic glucose and lipid metabolism.

BioDataome: a collection of uniformly preprocessed and automatically annotated datasets for data-driven biology

Biotechnology revolution generates a plethora of omics data with an exponential growth pace. Therefore, biological data mining demands automatic, 'high quality' curation efforts to organize biomedical knowledge into online databases. BioDataome is a database of uniformly preprocessed and disease-annotated omics data with the aim to promote and accelerate the reuse of public data. We followed the same preprocessing pipeline for each biological mart (microarray gene expression, RNA-Seq gene expression and DNA methylation) to produce ready for downstream analysis datasets and automatically annotated them with disease-ontology terms. We also designate datasets that share common samples and automatically discover control samples in case-control studies. Currently, BioDataome includes ∼5600 datasets, ∼260 000 samples spanning ∼500 diseases and can be easily used in large-scale massive experiments and meta-analysis. All datasets are publicly available for querying and downloading via BioDataome web application. We demonstrate BioDataome's utility by presenting exploratory data analysis examples. We have also developed BioDataome R package found in: https://github.com/mensxmachina/BioDataome/.Database URL: http://dataome.mensxmachina.org/.

Suppression of Adaptive Responses to Targeted Cancer Therapy by Transcriptional Repression

Acquired drug resistance is a major factor limiting the effectiveness of targeted cancer therapies. Targeting tumors with kinase inhibitors induces complex adaptive programs that promote the persistence of a fraction of the original cell population, facilitating the eventual outgrowth of inhibitor-resistant tumor clones. We show that the addition of a newly identified CDK7/12 inhibitor, THZ1, to targeted therapy enhances cell killing and impedes the emergence of drug-resistant cell populations in diverse cellular and in vivo cancer models. We propose that targeted therapy induces a state of transcriptional dependency in a subpopulation of cells poised to become drug tolerant, which THZ1 can exploit by blocking dynamic transcriptional responses, promoting remodeling of enhancers and key signaling outputs required for tumor cell survival in the setting of targeted therapy. These findings suggest that the addition of THZ1 to targeted therapies is a promising broad-based strategy to hinder the emergence of drug-resistant cancer cell populations.Significance: CDK7/12 inhibition prevents active enhancer formation at genes, promoting resistance emergence in response to targeted therapy, and impedes the engagement of transcriptional programs required for tumor cell survival. CDK7/12 inhibition in combination with targeted cancer therapies may serve as a therapeutic paradigm for enhancing the effectiveness of targeted therapies. Cancer Discov; 8(1); 59-73. ©2017 AACR.See related commentary by Carugo and Draetta, p. 17This article is highlighted in the In This Issue feature, p. 1.

Reactive oxygen species as signaling molecules in the development of lung fibrosis

Pulmonary fibrosis is a relatively rare but devastating disease characterized by the excessive deposition of extracellular matrix. The increased matrix results in reduced lung compliance and increased work of breathing, while the obliteration of alveolar-capillary structures can result in hypoxemia and pulmonary hypertension, which manifests clinically as worsening shortness of breath, respiratory failure, and death. Unbiased genome-wide association studies combined with animal models suggest that damage to the alveolar epithelium is the initiating factor in pulmonary fibrosis. This epithelial injury leads to the activation and proliferation of myofibroblasts that secrete extracellular matrix proteins characteristic of fibrosis. The best described molecular link between alveolar epithelial dysfunction and myofibroblast activation and proliferation is the profibrotic cytokine transforming growth factor-β (TGF-β). We and others have found that mitochondrial and NAD(P)H oxidase-generated reactive oxygen species (ROS) play a signaling role to enhance TGF-β signaling and promote fibrosis. The purpose of this article is to review how ROS signaling leads to the activation of TGF-β. We suggest that an improved understanding of these pathways might explain the failure of nonselective antioxidants to improve outcomes in patients with pulmonary fibrosis and might identify novel targets for therapy.

CHARGE syndrome modeling using patient-iPSCs reveals defective migration of neural crest cells harboring CHD7 mutations

CHARGE syndrome is caused by heterozygous mutations in the chromatin remodeler, CHD7, and is characterized by a set of malformations that, on clinical grounds, were historically postulated to arise from defects in neural crest formation during embryogenesis. To better delineate neural crest defects in CHARGE syndrome, we generated induced pluripotent stem cells (iPSCs) from two patients with typical syndrome manifestations, and characterized neural crest cells differentiated in vitro from these iPSCs (iPSC-NCCs). We found that expression of genes associated with cell migration was altered in CHARGE iPSC-NCCs compared to control iPSC-NCCs. Consistently, CHARGE iPSC-NCCs showed defective delamination, migration and motility in vitro, and their transplantation in ovo revealed overall defective migratory activity in the chick embryo. These results support the historical inference that CHARGE syndrome patients exhibit defects in neural crest migration, and provide the first successful application of patient-derived iPSCs in modeling craniofacial disorders.

CHARGE syndrome is a disease in which organs including the heart, eyes and ears may not develop properly. The cells that form the tissues affected by CHARGE syndrome develop in embryos from precursor cells called neural crest cells. Individuals with CHARGE syndrome also have mutations in a gene called CHD7. However, it is difficult to examine how CHD7 mutations affect neural crest cells in embryos.

In recent years, cell reprogramming techniques have made it possible to create induced pluripotent stem cells (iPSCs) from the specialized somatic cells found in the human body. These iPSCs can be developed into many different cell types, including neural crest cells.

Okuno et al. created iPSCs from the skin cells of people with CHARGE syndrome, developed these cells into neural crest cells, and compared them with neural crest cells that were developed from the skin cells of people without CHARGE syndrome. The neural crest cells developed from people with CHARGE syndrome showed multiple abnormalities. For example, they were not able to move around correctly. This is an important observation because neural crest cells must move through tissues to form the various organs affected by CHARGE syndrome.

Okuno et al. also observed changes in the activity of many genes other than CHD7 in the neural crest cells developed from CHARGE patients. Further research is now needed to find out which genes are the most important for restoring the normal activity of neural crest cells.

Pathways of cellular proteostasis in aging and disease

Ensuring cellular protein homeostasis, or proteostasis, requires precise control of protein synthesis, folding, conformational maintenance, and degradation. A complex and adaptive proteostasis network coordinates these processes with molecular chaperones of different classes and their regulators functioning as major players. This network serves to ensure that cells have the proteins they need while minimizing misfolding or aggregation events that are hallmarks of age-associated proteinopathies, including neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. It is now clear that the capacity of cells to maintain proteostasis undergoes a decline during aging, rendering the organism susceptible to these pathologies. Here we discuss the major proteostasis pathways in light of recent research suggesting that their age-dependent failure can both contribute to and result from disease. We consider different strategies to modulate proteostasis capacity, which may help develop urgently needed therapies for neurodegeneration and other age-dependent pathologies.

Linking site-specific loss of histone acetylation to repression of gene expression by the mycotoxin ochratoxin A

Ochratoxin A (OTA) is a potent renal carcinogen but its mechanism has not been fully resolved. In vitro and in vivo gene expression studies consistently revealed down-regulation of gene expression as the predominant transcriptional response to OTA. Based on the importance of specific histone acetylation marks in regulating gene transcription and our recent finding that OTA inhibits histone acetyltransferases (HATs), leading to loss of acetylation of histones and non-histone proteins, we hypothesized that OTA-mediated repression of gene expression may be causally linked to HAT inhibition and loss of histone acetylation. In this study, we used a novel mass spectrometry approach employing chemical ¹³C-acetylation of unmodified lysine residues for quantification of post-translational acetylation sites to identify site-specific alterations in histone acetylation in human kidney epithelial cells (HK-2) exposed to OTA. These results showed OTA-mediated hypoacetylation at almost all lysine residues of core histones, including loss of acetylation at H3K9 and H3K14, which are hallmarks of gene activation. ChIP-qPCR used to establish a possible link between H3K9 or H3K14 hypoacetylation and OTA-mediated down-regulation of selected genes (AMIGO2, CLASP2, CTNND1) confirmed OTA-mediated H3K9 hypoacetylation at promoter regions of these genes. Integrated analysis of OTA-mediated genome-wide changes in H3K9 acetylation identified by ChIP-Seq with published gene expression data further demonstrated that among OTA-responsive genes almost 80% of hypoacetylated genes were down-regulated, thus confirming an association between H3K9 acetylation status and gene expression of these genes. However, only 7% of OTA repressed genes showed loss of H3K9 acetylation within promoter regions. Interestingly, however, GO analysis and functional enrichment of down-regulated genes showing loss of H3K9 acetylation at their respective promoter regions revealed enrichment of genes involved in the regulation of transcription, including a number of transcription factors that are predicted to directly or indirectly regulate the expression of 98% of OTA repressed genes. Thus, it is possible that histone acetylation changes in a fairly small set of genes but with key function in transcriptional regulation may trigger a cascade of events that may lead to overall repression of gene expression. Taken together, our data provide evidence for a mechanistic link between loss of H3K9 acetylation as a consequence of OTA-mediated inhibition of HATs and repression of gene expression by OTA, thereby affecting cellular processes critical to tumorigenesis.

Controversial Effects of D-Amino Acid Oxidase Activator (DAOA)/G72 on D-Amino Acid Oxidase (DAO) Activity in Human Neuronal, Astrocyte and Kidney Cell Lines: The N-methyl D-aspartate (NMDA) Receptor Hypofunction Point of View

Dysfunction of D-amino acid oxidase (DAO) and DAO activator (DAOA)/G72 genes have been linked to neuropsychiatric disorders. The glutamate hypothesis of schizophrenia has proposed that increased DAO activity leads to decreased D-serine, which subsequently may lead to N-methyl-D-aspartate (NMDA) receptor hypofunction. It has been shown that DAOA binds to DAO and increases its activity. However, there are also studies showing DAOA decreases DAO activity. Thus, the effect of DAOA on DAO is controversial. We aimed to understand the effect of DAOA on DAO activity in neuron-like (SH-SY5Y), astrocyte-like (1321N1) and kidney-like (HEK293) human cell lines. DAO activity was measured based on the release of hydrogen peroxide and its interaction with Amplex Red reagent. We found that DAOA increases DAO activity only in HEK293 cells, but has no effect on DAO activity in SH-SY5Y and 1321N1 cells. This might be because of different signaling pathways, or due to lower DAO and DAOA expression in SH-SY5Y and 1321N1 cells compared to HEK293 cells, but also due to different compartmentalization of the proteins. The lower DAO and DAOA expression in neuron-like SH-SY5Y and astrocyte-like 1321N1 cells might be due to tightly regulated expression, as previously reported in the human post-mortem brain. Our simulation experiments to demonstrate the interaction between DAOA and human DAO (hDAO) showed that hDAO holoenzyme [hDAO with flavine adenine dinucleotide (FAD)] becomes more flexible and misfolded in the presence of DAOA, whereas DAOA had no effect on hDAO apoprotein (hDAO without FAD), which indicate that DAOA inactivates hDAO holoenzyme. Furthermore, patch-clamp analysis demonstrated no effect of DAOA on NMDA receptor activity in NR1/NR2A HEK293 cells. In summary, the interaction between DAO and DAOA seems to be cell type and its biochemical characteristics dependent which still needs to be elucidated.

Early Transcriptional Divergence Marks Virus-Specific Primary Human CD8+ T Cells in Chronic versus Acute Infection

Distinct molecular pathways govern the differentiation of CD8+ effector T cells into memory or exhausted T cells during acute and chronic viral infection, but these are not well studied in humans. Here, we employed an integrative systems immunology approach to identify transcriptional commonalities and differences between virus-specific CD8+ T cells from patients with persistent and spontaneously resolving hepatitis C virus (HCV) infection during the acute phase. We observed dysregulation of metabolic processes during early persistent infection that was linked to changes in expression of genes related to nucleosomal regulation of transcription, T cell differentiation, and the inflammatory response and correlated with subject age, sex, and the presence of HCV-specific CD4+ T cell populations. These early changes in HCV-specific CD8+ T cell transcription preceded the overt establishment of T cell exhaustion, making this signature a prime target in the search for the regulatory origins of T cell dysfunction in chronic viral infection.

Structural basis for regulation of the nucleo-cytoplasmic distribution of Bag6 by TRC35

The metazoan protein BCL2-associated athanogene cochaperone 6 (Bag6) forms a hetero-trimeric complex with ubiquitin-like 4A and transmembrane domain recognition complex 35 (TRC35). This Bag6 complex is involved in tail-anchored protein targeting and various protein quality-control pathways in the cytosol as well as regulating transcription and histone methylation in the nucleus. Here we present a crystal structure of Bag6 and its cytoplasmic retention factor TRC35, revealing that TRC35 is remarkably conserved throughout the opisthokont lineage except at the C-terminal Bag6-binding groove, which evolved to accommodate Bag6, a unique metazoan factor. While TRC35 and its fungal homolog, guided entry of tail-anchored protein 4 (Get4), utilize a conserved hydrophobic patch to bind their respective partners, Bag6 wraps around TRC35 on the opposite face relative to the Get4-5 interface. We further demonstrate that TRC35 binding is critical not only for occluding the Bag6 nuclear localization sequence from karyopherin α to retain Bag6 in the cytosol but also for preventing TRC35 from succumbing to RNF126-mediated ubiquitylation and degradation. The results provide a mechanism for regulation of Bag6 nuclear localization and the functional integrity of the Bag6 complex in the cytosol.

Clustergrammer, a web-based heatmap visualization and analysis tool for high-dimensional biological data

Most tools developed to visualize hierarchically clustered heatmaps generate static images. Clustergrammer is a web-based visualization tool with interactive features such as: zooming, panning, filtering, reordering, sharing, performing enrichment analysis, and providing dynamic gene annotations. Clustergrammer can be used to generate shareable interactive visualizations by uploading a data table to a web-site, or by embedding Clustergrammer in Jupyter Notebooks. The Clustergrammer core libraries can also be used as a toolkit by developers to generate visualizations within their own applications. Clustergrammer is demonstrated using gene expression data from the cancer cell line encyclopedia (CCLE), original post-translational modification data collected from lung cancer cells lines by a mass spectrometry approach, and original cytometry by time of flight (CyTOF) single-cell proteomics data from blood. Clustergrammer enables producing interactive web based visualizations for the analysis of diverse biological data.

Genetic, Transcriptome, Proteomic, and Epidemiological Evidence for Blood-Brain Barrier Disruption and Polymicrobial Brain Invasion as Determinant Factors in Alzheimer's Disease

Diverse pathogens are detected in Alzheimer's disease (AD) brains. A bioinformatics survey showed that AD genome-wide association study (GWAS) genes (localized in bone marrow, immune locations and microglia) relate to multiple host/pathogen interactomes (Candida albicans, Cryptococcus neoformans, Bornavirus, Borrelia burgdorferri, cytomegalovirus, Ebola virus, HSV-1, HERV-W, HIV-1, Epstein-Barr, hepatitis C, influenza, Chlamydia pneumoniae, Porphyrymonas gingivalis, Helicobacter pylori, Toxoplasma gondii, Trypanosoma cruzi). These interactomes also relate to the AD hippocampal transcriptome and to plaque or tangle proteins. Upregulated AD hippocampal genes match those upregulated by multiple bacteria, viruses, fungi, or protozoa in immunocompetent cells. AD genes are enriched in GWAS datasets reflecting pathogen diversity, suggesting selection for pathogen resistance, as supported by the old age of AD patients, implying resistance to earlier infections. APOE4 is concentrated in regions of high parasitic burden and protects against childhood tropical infections and hepatitis C. Immune/inflammatory gain of function applies to APOE4, CR1, and TREM2 variants. AD genes are also expressed in the blood-brain barrier (BBB), which is disrupted by AD risk factors (age, alcohol, aluminum, concussion, cerebral hypoperfusion, diabetes, homocysteine, hypercholesterolemia, hypertension, obesity, pesticides, pollution, physical inactivity, sleep disruption, smoking) and by pathogens, directly or via olfactory routes to basal-forebrain BBB control centers. The BBB benefits from statins, NSAIDs, estrogen, melatonin, memantine, and the Mediterranean diet. Polymicrobial involvement is supported by upregulation of bacterial, viral, and fungal sensors/defenders in the AD brain, blood, or cerebrospinal fluid. AD serum amyloid-β autoantibodies may attenuate its antimicrobial effects favoring microbial survival and cerebral invasion leading to activation of neurodestructive immune/inflammatory processes, which may also be augmented by age-related immunosenescence. AD may thus respond to antibiotic, antifungal, or antiviral therapy.

Rapid Discovery of De Novo Deleterious Mutations in Cattle Enhances the Value of Livestock as Model Species

In humans, the clinical and molecular characterization of sporadic syndromes is often hindered by the small number of patients and the difficulty in developing animal models for severe dominant conditions. Here we show that the availability of large data sets of whole-genome sequences, high-density SNP chip genotypes and extensive recording of phenotype offers an unprecedented opportunity to quickly dissect the genetic architecture of severe dominant conditions in livestock. We report on the identification of seven dominant de novo mutations in CHD7, COL1A1, COL2A1, COPA, and MITF and exploit the structure of cattle populations to describe their clinical consequences and map modifier loci. Moreover, we demonstrate that the emergence of recessive genetic defects can be monitored by detecting de novo deleterious mutations in the genome of bulls used for artificial insemination. These results demonstrate the attractiveness of cattle as a model species in the post genomic era, particularly to confirm the genetic aetiology of isolated clinical case reports in humans.

Identification of high-confidence RNA regulatory elements by combinatorial classification of RNA-protein binding sites

Crosslinking immunoprecipitation sequencing (CLIP-seq) technologies have enabled researchers to characterize transcriptome-wide binding sites of RNA-binding protein (RBP) with high resolution. We apply a soft-clustering method, RBPgroup, to various CLIP-seq datasets to group together RBPs that specifically bind the same RNA sites. Such combinatorial clustering of RBPs helps interpret CLIP-seq data and suggests functional RNA regulatory elements. Furthermore, we validate two RBP–RBP interactions in cell lines. Our approach links proteins and RNA motifs known to possess similar biochemical and cellular properties and can, when used in conjunction with additional experimental data, identify high-confidence RBP groups and their associated RNA regulatory elements.

Cellular responses to human cytomegalovirus infection: Induction of a mesenchymal-to-epithelial transition (MET) phenotype

Human cytomegalovirus (HCMV) is the prototypical human β-herpes virus. Here we perform a systems analysis of the HCMV host-cell transcriptome, using gene set enrichment analysis (GSEA) as an engine to globally map the host-pathogen interaction across two cell types. Our analysis identified several previously unknown signatures of infection, such as induction of potassium channels and amino acid transporters, derepression of genes marked with histone H3 lysine 27 trimethylation (H3K27me3), and inhibition of genes related to epithelial-to-mesenchymal transition (EMT). The repression of EMT genes was dependent on early viral gene expression and correlated with induction E-cadherin (CDH1) and mesenchymal-to-epithelial transition (MET) genes. Infection of transformed breast carcinoma and glioma stem cells similarly inhibited EMT and induced MET, arguing that HCMV induces an epithelium-like cellular environment during infection.

Size-based detection of sarcoma circulating tumor cells and cell clusters

Metastatic disease is the most important factor in determining the survival of sarcoma patients. Since sarcoma metastasis is predominantly hematogenous, we hypothesized that detection and quantification of circulating tumor cells (CTCs) could reflect response to therapy and risk of metastatic relapse. We evaluated the presence of CTCs using a novel animal model and in the blood of patients with high grade sarcomas utilizing the CellSieve™ size-based low pressure microfiltration system. Sarcoma CTCs were identified based on antibody staining patterns and nuclear morphology. Additionally, RNA was extracted from the CTCs for molecular analysis including demonstration of an EWS-FLI1 translocation, identification of a previously unrecognized p53 mutation in a patient with Ewing sarcoma, and single cell RNA sequencing of CTC from a child with alveolar rhabdomyosarcoma. In mouse xenograft models, the presence of CTC correlates with disease burden and with clinically silent metastases. In human patients, CTCs were readily detected at diagnosis, decreased with successful treatment, and were detectable in the blood of patients with no radiographic evidence of disease prior to the development of overt metastasis. Although evaluation of CTC is established in the care of patients with carcinomas, this technology has yet to be effectively applied to the evaluation and treatment of sarcoma patients. Our work demonstrates that the CellSieve™ microfiltration system can be used to study the biology of CTC in both mouse models and human sarcoma patients, with the potential for application to the monitoring of disease response and prediction of metastatic relapse.

The Tyr113His T/C rs1051740 and 'very slow' phenotype of the EPHX1 gene alters miR-26b-5p and miR-1207-5p expression in pregnancy

BACKGROUND

Environmental insults and microsomal epoxide hydrolase 1 (EPHX1) single nucleotide polymorphisms (SNPs), Tyr113His T/C rs1051740 and His139Arg A/G rs2234922, aberrantly alters microRNA (miR) expression and are linked to low birthweights (LBW).

OBJECTIVES

To investigate the interplay between pollution, EPHX1 SNPs and miRs during pregnancy and associated LBW outcomes.

METHODS

South African pregnant women (n=241) were recruited in the MACE birth cohort study in Durban, a city with high levels of industry and traffic related pollutants. EPHX1 SNPs were genotyped using PCR-RFLP and grouped into their respective phenotypes, i.e. normal (N), slow (S), very slow (VS) and fast (F). EPHX1, miR-26b-5p, miR-193b-3p and miR-1207-5p expression were determined using quantitative PCR.

RESULTS

Mothers with the Tyr113His SNP had low iron levels [TT vs. TC+CC: mean difference (MD)=0.67g/dl; p=0.0167], LBW [TT vs. TC+CC: MD=189.30g; p=0.0067], and low EPHX1 expression; p<0.0001. miR-26b-5p and miR-1207-5p expression were significantly higher in the CC genotypes compared to TT+TC groups; p<0.0001. The opposite trend occurred for miR-193b-3p; p=0.0045. Mothers with the VS phenotype had low iron levels [N vs. VS and VS vs. F: MD=2.03 and -1.96g/dl; p=0.0021, respectively], decreased gestational age [VS vs. F: MD=-2.14weeks; p=0.0051, respectively], and LBW [N vs. VS, VS vs. F and S vs. VS: MD=1000, -940.30 and 968.80g; p<0.0001, respectively]; F phenotype had the highest EPHX1 expression [N vs. F, VS vs. F and S vs. F: MD=-1.067, -1.854 and -1.379; p=0.0002, respectively]; and N phenotype had low miR-26b-5p [N vs. VS: MD=-0.6100; p=0.0159] and miR-1207-5p [N vs. VS and VS vs. F: MD=-0.834 and 1.103; p=0.0007, respectively] expression. miR-193b-3p expression between phenotypes remained unchanged.

CONCLUSION

The Tyr113His T/C variant of rs1051740 and VS phenotype alters EPHX1, miR-26b-5p and miR-1207-5p expression, and contributes towards low blood iron levels and LBW.

Necroptosis activation in Alzheimer's disease

Alzheimer's disease (AD) is characterized by severe neuronal loss; however, the mechanisms by which neurons die remain elusive. Necroptosis, a programmed form of necrosis, is executed by the mixed lineage kinase domain-like (MLKL) protein, which is triggered by receptor-interactive protein kinases (RIPK) 1 and 3. We found that necroptosis was activated in postmortem human AD brains, positively correlated with Braak stage, and inversely correlated with brain weight and cognitive scores. In addition, we found that the set of genes regulated by RIPK1 overlapped significantly with multiple independent AD transcriptomic signatures, indicating that RIPK1 activity could explain a substantial portion of transcriptomic changes in AD. Furthermore, we observed that lowering necroptosis activation reduced cell loss in a mouse model of AD. We anticipate that our findings will spur a new area of research in the AD field focused on developing new therapeutic strategies aimed at blocking its activation.

Reconstructing context-specific gene regulatory network and identifying modules and network rewiring through data integration

Reconstructing context-specific transcriptional regulatory network is crucial for deciphering principles of regulatory mechanisms underlying various conditions. Recently studies that reconstructed transcriptional networks have focused on individual organisms or cell types and relied on data repositories of context-free regulatory relationships. Here we present a comprehensive framework to systematically derive putative regulator-target pairs in any given context by integrating context-specific transcriptional profiling and public data repositories of gene regulatory networks. Moreover, our framework can identify core regulatory modules and signature genes underlying global regulatory circuitry, and detect network rewiring and core rewired modules in different contexts by considering gene modules and edge (gene interaction) modules collaboratively. We applied our methods to analyzing Autism RNA-seq experiment data and produced biologically meaningful results. In particular, all 11 hub genes in a predicted rewired autistic regulatory subnetwork have been linked to autism based on literature review. The predicted rewired autistic regulatory network may shed some new insight into disease mechanism.

Analysis of differential gene expression mediated by clozapine in human postmortem brains

Clozapine is the only medication indicated for treating refractory schizophrenia, due to its superior efficacy among all antipsychotic agents, but its mechanism of action is poorly understood. To date, no studies of human postmortem brain have characterized the gene expression response to clozapine. Therefore, we addressed this question by analyzing expression data extracted from published microarray studies involving brains of patients on antipsychotic therapy. We first performed a systematic review and identified four microarray studies of postmortem brains from antipsychotic-treated patients, then extracted the expression data. We then performed generalized linear model analysis on each study separately, and identified the genes differentially expressed in response to clozapine compared to other atypical antipsychotic medications, as well as their associated canonical pathways. We also found a number of genes common to all four studies that we analyzed: GCLM, ZNF652, and GYPC. In addition, pathway analysis highlighted the following processes in all four studies: clathrin-mediated endocytosis, SAPK/JNK signaling, 3-phosphoinositide synthesis, and paxillin signaling. Our analysis yielded the first comprehensive compendium of genes and pathways differentially expressed upon clozapine treatment in the human brain, which may provide insight into the mechanism and unique efficacy of clozapine, as well as the pathophysiology of schizophrenia.

Developing a framework for digital objects in the Big Data to Knowledge (BD2K) commons: Report from the Commons Framework Pilots workshop

The volume and diversity of data in biomedical research have been rapidly increasing in recent years. While such data hold significant promise for accelerating discovery, their use entails many challenges including: the need for adequate computational infrastructure, secure processes for data sharing and access, tools that allow researchers to find and integrate diverse datasets, and standardized methods of analysis. These are just some elements of a complex ecosystem that needs to be built to support the rapid accumulation of these data. The NIH Big Data to Knowledge (BD2K) initiative aims to facilitate digitally enabled biomedical research. Within the BD2K framework, the Commons initiative is intended to establish a virtual environment that will facilitate the use, interoperability, and discoverability of shared digital objects used for research. The BD2K Commons Framework Pilots Working Group (CFPWG) was established to clarify goals and work on pilot projects that address existing gaps toward realizing the vision of the BD2K Commons. This report reviews highlights from a two-day meeting involving the BD2K CFPWG to provide insights on trends and considerations in advancing Big Data science for biomedical research in the United States.

Envisioning the future of 'big data' biomedicine

Through the increasing availability of more efficient data collection procedures, biomedical scientists are now confronting ever larger sets of data, often finding themselves struggling to process and interpret what they have gathered. This, while still more data continues to accumulate. This torrent of biomedical information necessitates creative thinking about how the data are being generated, how they might be best managed, analyzed, and eventually how they can be transformed into further scientific understanding for improving patient care. Recognizing this as a major challenge, the National Institutes of Health (NIH) has spearheaded the "Big Data to Knowledge" (BD2K) program - the agency's most ambitious biomedical informatics effort ever undertaken to date. In this commentary, we describe how the NIH has taken on "big data" science head-on, how a consortium of leading research centers are developing the means for handling large-scale data, and how such activities are being marshalled for the training of a new generation of biomedical data scientists. All in all, the NIH BD2K program seeks to position data science at the heart of 21st Century biomedical research.

Altered expression of histamine signaling genes in autism spectrum disorder

The histaminergic system (HS) has a critical role in cognition, sleep and other behaviors. Although not well studied in autism spectrum disorder (ASD), the HS is implicated in many neurological disorders, some of which share comorbidity with ASD, including Tourette syndrome (TS). Preliminary studies suggest that antagonism of histamine receptors 1-3 reduces symptoms and specific behaviors in ASD patients and relevant animal models. In addition, the HS mediates neuroinflammation, which may be heightened in ASD. Together, this suggests that the HS may also be altered in ASD. Using RNA sequencing (RNA-seq), we investigated genome-wide expression, as well as a focused gene set analysis of key HS genes (HDC, HNMT, HRH1, HRH2, HRH3 and HRH4) in postmortem dorsolateral prefrontal cortex (DLPFC) initially in 13 subjects with ASD and 39 matched controls. At the genome level, eight transcripts were differentially expressed (false discovery rate <0.05), six of which were small nucleolar RNAs (snoRNAs). There was no significant diagnosis effect on any of the individual HS genes but expression of the gene set of HNMT, HRH1, HRH2 and HRH3 was significantly altered. Curated HS gene sets were also significantly differentially expressed. Differential expression analysis of these gene sets in an independent RNA-seq ASD data set from DLPFC of 47 additional subjects confirmed these findings. Understanding the physiological relevance of an altered HS may suggest new therapeutic options for the treatment of ASD.

Genomic Analysis of Thymic Epithelial Tumors Identifies Novel Subtypes Associated with Distinct Clinical Features

Purpose: To reconcile the heterogeneity of thymic epithelial tumors (TET) and gain deeper understanding of the molecular determinants of TETs, we set out to establish a clinically relevant molecular classification system for these tumors.Experimental Design: Molecular subgrouping of TETs was performed in 120 patients from The Cancer Genome Atlas using a multidimensional approach incorporating analyses of DNA mutations, mRNA expression, and somatic copy number alterations (SCNA), and validated in two independent cohorts.Results: Four distinct molecular subtypes of TETs were identified. The most commonly identified gene mutation was a missense mutation in General Transcription Factor II-I (GTF2I group), which was present in 38% of patients. The next group was identified by unsupervised mRNA clustering of GTF2I wild-type tumors and represented TETs enriched in expression of genes associated with T-cell signaling (TS group; 33%). The remaining two groups were distinguished by their degree of chromosomal stability (CS group; 8%) or instability (CIN group; 21%) based upon SCNA analyses. Disease-free survival and overall survival were favorable in the GTF2I group and unfavorable in the CIN group. These molecular subgroups were associated with TET histology and clinical features including disease-free survival. Finally, we demonstrate high expression of PD1 mRNA and correlation of PD1 and CD8A in the TS subgroup.Conclusions: Molecular subtyping of TETs is associated with disease-free and overall survival. Classification of TETs by a molecular framework could aid in the refinement of staging and in the discovery and development of rational treatment options for patients with TETs. Clin Cancer Res; 23(16); 4855-64. ©2017 AACR.

Construction of a Comprehensive Protein-Protein Interaction Map for Vitiligo Disease to Identify Key Regulatory Elements: A Systemic Approach

Vitiligo is an idiopathic disorder characterized by depigmented patches on the skin due to progressive loss of melanocytes. Several genetic, immunological, and pathophysiological investigations have established vitiligo as a polygenetic disorder with multifactorial etiology. However, no definite model explaining the interplay between these causative factors has been established hitherto. Therefore, we studied the disorder at the system level to identify the key proteins involved by exploring their molecular connectivity in terms of topological parameters. The existing research data helped us in collating 215 proteins involved in vitiligo onset or progression. Interaction study of these proteins leads to a comprehensive vitiligo map with 4845 protein nodes linked with 107,416 edges. Based on centrality measures, a backbone network with 500 nodes has been derived. This has presented a clear overview of the proteins and processes involved and the crosstalk between them. Clustering backbone proteins revealed densely connected regions inferring major molecular interaction modules essential for vitiligo. Finally, a list of top order proteins that play a key role in the disease pathomechanism has been formulated. This includes SUMO2, ESR1, COPS5, MYC, SMAD3, and Cullin proteins. While this list is in fair agreement with the available literature, it also introduces new candidate proteins that can be further explored. A subnetwork of 64 vitiligo core proteins was built by analyzing the backbone and seed protein networks. Our finding suggests that the topology, along with functional clustering, provides a deep insight into the behavior of proteins. This in turn aids in the illustration of disease condition and discovery of significant proteins involved in vitiligo.

MAP7 Regulates Axon Collateral Branch Development in Dorsal Root Ganglion Neurons

Collateral branches from axons are key components of functional neural circuits that allow neurons to connect with multiple synaptic targets. Like axon growth and guidance, formation of collateral branches depends on the regulation of microtubules, but how such regulation is coordinated to ensure proper circuit development is not known. Based on microarray analysis, we have identified a role for microtubule-associated protein 7 (MAP7) during collateral branch development of dorsal root ganglion (DRG) sensory neurons. We show that MAP7 is expressed at the onset of collateral branch formation. Perturbation of its expression by overexpression or shRNA knockdown alters axon branching in cultured DRG neurons. Localization and time-lapse imaging analysis reveals that MAP7 is enriched at branch points and colocalizes with stable microtubules, but enters the new branch with a delay, suggesting a role in branch maturation. We have also investigated a spontaneous mutant mouse that expresses a truncated MAP7 and found a gain-of-function phenotype both in vitro and in vivo Further domain analysis suggests that the amino half of MAP7 is responsible for branch formation, suggesting a mechanism that is independent of its known interaction with kinesin. Moreover, this mouse exhibits increased pain sensitivity, a phenotype that is consistent with increased collateral branch formation. Therefore, our study not only uncovers the first neuronal function of MAP7, but also demonstrates the importance of proper microtubule regulation in neural circuit development. Furthermore, our data provide new insights into microtubule regulation during axonal morphogenesis and may shed light on MAP7 function in neurological disorders.SIGNIFICANCE STATEMENT Neurons communicate with multiple targets by forming axonal branches. In search of intrinsic factors that control collateral branch development, we identified a role for microtubule-associated protein 7 (MAP7) in dorsal root ganglion sensory neurons. We show that MAP7 expression is developmentally regulated and perturbation of this expression alters branch formation. Cell biological analysis indicates that MAP7 promotes branch maturation. Analysis of a spontaneous mouse mutant suggests a molecular mechanism for branch regulation and the potential influence of collateral branches on pain sensitivity. Our studies thus establish the first neuronal function of MAP7 and demonstrate its role in branch morphogenesis and neural circuit function. These findings may help in our understanding of the contribution of MAP7 to neurological disorders and nerve regeneration.

Governing roles for Trib3 pseudokinase during stress erythropoiesis

In response to anemia, the heightened production of erythropoietin (EPO) can sharply promote erythroid progenitor cell (EPC) formation. Specific mediators of such EPO- accelerated erythropoiesis, however, are not well understood. Presently, we first report that the expression of Trib3 in adult bone marrow EPCs in vivo is nominal at steady state, but strongly activated on EPO challenge. In a knockout mouse model, Trib3 disruption modestly increased steady-state erythrocyte numbers and decreased mean corpuscular volume. Following 5-fluorouracil myeloablation, however, rebound red blood cell production and hemoglobin levels were substantially (and selectively) compromised in Trib3^-/- mice versus Trib3^+/+ congenic controls. Erythrocytes from 5-fluorouracil-treated Trib3^-/- mice additionally were more prone to lysis and exhibited elevated peroxide-induced reactive oxygen species. Ex vivo, the development of CD71^posTer119^pos erythroblasts from Trib3^-/- bone marrow progenitors was attenuated, and this was associated with heightened EPO-dependent Erk1/2 activation and moderately increased Akt activation. For developmentally staged EPCs, gene profiling provided further initial insight into candidate mediators of EPO-induced Trib3 gene expression, including Cebp-beta, Atf4, Egr-1, and Nab1. Overall, Trib3 is indicated to act as a novel EPC-intrinsic governor of stress erythropoiesis.

Pharos: Collating protein information to shed light on the druggable genome

The 'druggable genome' encompasses several protein families, but only a subset of targets within them have attracted significant research attention and thus have information about them publicly available. The Illuminating the Druggable Genome (IDG) program was initiated in 2014, has the goal of developing experimental techniques and a Knowledge Management Center (KMC) that would collect and organize information about protein targets from four families, representing the most common druggable targets with an emphasis on understudied proteins. Here, we describe two resources developed by the KMC: the Target Central Resource Database (TCRD) which collates many heterogeneous gene/protein datasets and Pharos (https://pharos.nih.gov), a multimodal web interface that presents the data from TCRD. We briefly describe the types and sources of data considered by the KMC and then highlight features of the Pharos interface designed to enable intuitive access to the IDG knowledgebase. The aim of Pharos is to encourage 'serendipitous browsing', whereby related, relevant information is made easily discoverable. We conclude by describing two use cases that highlight the utility of Pharos and TCRD.

Open Targets: a platform for therapeutic target identification and validation

We have designed and developed a data integration and visualization platform that provides evidence about the association of known and potential drug targets with diseases. The platform is designed to support identification and prioritization of biological targets for follow-up. Each drug target is linked to a disease using integrated genome-wide data from a broad range of data sources. The platform provides either a target-centric workflow to identify diseases that may be associated with a specific target, or a disease-centric workflow to identify targets that may be associated with a specific disease. Users can easily transition between these target- and disease-centric workflows. The Open Targets Validation Platform is accessible at https://www.targetvalidation.org.