Annotating cancer variants and anti-cancer therapeutics in reactome

Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes

Germline and somatic mutations can give rise to proteins with altered activity, including both gain and loss-of-function. The effects of these variants can be captured in disease-specific reactions and pathways that highlight the resulting changes to normal biology. A disease reaction is defined as an aberrant reaction in which a variant protein participates. A disease pathway is defined as a pathway that contains a disease reaction. Annotation of disease variants as participants of disease reactions and disease pathways can provide a standardized overview of molecular phenotypes of pathogenic variants that is amenable to computational mining and mathematical modeling. Reactome (https://reactome.org/), an open source, manually curated, peer-reviewed database of human biological pathways, in addition to providing annotations for >11 000 unique human proteins in the context of ∼15 000 wild-type reactions within more than 2000 wild-type pathways, also provides annotations for >4000 disease variants of close to 400 genes as participants of ∼800 disease reactions in the context of ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, described in wild-type reactions and pathways, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Reactome's data model enables mapping of disease variant datasets to specific disease reactions within disease pathways, providing a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity. Database URL: https://reactome.org/.

Cardiopulmonary Phenotypes and Protein Signatures in Children With Down Syndrome

Pulmonary disease, lower respiratory tract infection, and pneumonia are the largest causes of morbidity and mortality in individuals with Down syndrome (DS), but whether pulmonary diagnoses in children with DS are common and occur independently of cardiac disease and pulmonary hypertension (PH) is unknown. Cardiopulmonary phenotypes were examined in a cohort of 1248 children with DS. Aptamer-based proteomic analysis of blood was performed in a subset (n = 120) of these children. By the age of 10 years, half of the patients in this cohort (n = 634, 50.8%) had co-occurring pulmonary diagnoses. That proteins and related pathways were distinct between children with pulmonary diagnoses and those with cardiac disease and/or PH may indicate that pulmonary diagnoses appear to occur independently of cardiac disease and PH. Heparin sulfate-glycosaminoglycandegradation, nicotinate metabolism, and elastic fiber formation were ranked highest in the group with pulmonary diagnoses.

Predisposing deleterious variants in the cancer-associated human kinases in the global populations

Human kinases play essential and diverse roles in the cellular activities of maintaining homeostasis and growth. Genetic mutations in the genes encoding the kinases (or phosphotransferases) have been linked with various types of cancers. In this study, we cataloged mutations in 500 kinases genes in >65,000 individuals of global populations from the Human Genetic Diversity Project (HGDP) and ExAC databases, and assessed their potentially deleterious impact by using the in silico tools SIFT, Polyphen2, and CADD. The analysis highlighted 35 deleterious non-synonymous SNVs in the ExAC and 5 SNVs in the HGDP project. Notably, a higher number of deleterious mutations was observed in the Non-Finnish Europeans (26 SNVs), followed by the Africans (14 SNVs), East Asians (13 SNVs), and South Asians (12 SNVs). The gene set enrichment analysis highlighted NTRK1 and FGFR3 being most significantly enriched among the kinases. The gene expression analysis revealed over-expression of NTRK1 in liver cancer, whereas, FGFR3 was found over-expressed in lung, breast, and liver cancers compared to their expression in the respective normal tissues. Also, 13 potential drugs were identified that target the NTRK1 protein, whereas 6 potential drugs for the FGFR3 target were identified. Taken together, the study provides a framework for exploring the predisposing germline mutations in kinases to suggest the underlying pathogenic mechanisms in cancers. The potential drugs are also suggested for personalized cancer management.

Single nucleotide variants in nuclear pore complex disassembly pathway associated with poor survival in osteosarcoma

Introduction

The bone tumor, osteosarcoma, remains challenging to treat in children and young adults, especially when patients present with metastatic disease. Developing new therapies based on genomic data from sequencing projects has proven difficult given the lack of recurrent genetic lesions across tumors. MYC overexpression has been associated with poor outcomes in osteosarcoma. However, other genomic markers of disease severity are lacking.

Materials and Methods

We utilized whole genome sequencing of 106 tumors and matched normal controls in order to define genomic characteristics that correlate with overall survival. Single nucleotide variants were overlaid onto annotated molecular pathways in order to define aberrant pathway signatures specific to aggressive osteosarcoma. Additionally, we calculated differential gene expression in a subsample of 71 tumors. Differentially expressed genes were then queried for known MYC-responsive genes.

Results

Molecular pathways specific to nuclear pore complex disassembly (NPCD) show significant correlation with poor overall survival in osteosarcoma when mutations were present. Genes involved in immune response and immune regulation are enriched in the differential expression analysis of samples with and without NPCD pathway aberrations. Furthermore, neither MYC nor MYC-responsive genes show differential expression between NPCD-aberrant and non-aberrant groups. The NPCD pathway mutations are dominated by regulatory region variants rather than protein-altering mutations, suggesting that dysregulation of genetic regulatory networks may be the underlying mechanism for their relation to osteosarcoma phenotype.

Discussion

Overall survival is significantly worse in patients whose tumors show aberrations in the NPCD pathway. Moreover, this difference in survival is not driven by MYC-overexpression, suggesting a novel mechanism for some aggressive osteosarcomas. These findings add light to the evolving understanding of the drivers of osteosarcoma and may aid in the search for new treatments based on patient-specific genetic data.

Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes

Disease variant annotation in the context of biological reactions and pathways can provide a standardized overview of molecular phenotypes of pathogenic mutations that is amenable to computational mining and mathematical modeling. Reactome, an open source, manually curated, peer-reviewed database of human biological pathways, provides annotations for over 4000 disease variants of close to 400 genes in the context of ∼800 disease reactions constituting ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics (ACMG). Reactome's pathway-based, reaction-specific disease variant dataset and data model provide a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity.

Cellular signaling modulated by miRNA-3652 in ovarian cancer: unveiling mechanistic pathways for future therapeutic strategies

MicroRNAs (miRNAs) are small non-coding RNA molecules that play pivotal roles in regulating gene expression and have been implicated in the pathogenesis of numerous cancers. miRNA-3652, though relatively less explored, has recently emerged as a potential key player in ovarian cancer's molecular landscape. This review aims to delineate the functional significance and tumor progression role of miRNA-3652 in ovarian cancer, shedding light on its potential as both a diagnostic biomarker and therapeutic target. A comprehensive literature search was carried out using established databases, the focus was on articles that reported the role of miRNA-3652 in ovarian cancer, encompassing mechanistic insights, functional studies, and its association with clinical outcomes. This updated review highlighted that miRNA-3652 is intricately involved in ovarian cancer cell proliferation, migration, and invasion, its dysregulation was linked to altered expression of critical genes involved in tumor growth and metastasis; furthermore, miRNA-3652 expression levels were found to correlate with clinical stages, prognosis, and response to therapy in ovarian cancer patients. miRNA-3652 holds significant promise as a vital molecular player in ovarian cancer's pathophysiology. Its functional role and impact on tumor progression make it a potential candidate for diagnostic and therapeutic applications in ovarian cancer. Given the pivotal role of miRNA-3652 in ovarian cancer, future studies should emphasize in-depth mechanistic explorations, utilizing advanced genomic and proteomic tools. Collaboration between basic scientists and clinicians will be vital to translating these findings into innovative diagnostic and therapeutic strategies, ultimately benefiting ovarian cancer patients. Video Abstract.

Integrative systems biology characterizes immune-mediated neurodevelopmental changes in murine Zika virus microcephaly

Characterizing perturbation of molecular pathways in congenital Zika virus (ZIKV) infection is critical for improved therapeutic approaches. Leveraging integrative systems biology, proteomics, and RNA-seq, we analyzed embryonic brain tissues from an immunocompetent, wild-type congenital ZIKV infection mouse model. ZIKV induced a robust immune response accompanied by the downregulation of critical neurodevelopmental gene programs. We identified a negative correlation between ZIKV polyprotein abundance and host cell cycle-inducing proteins. We further captured the downregulation of genes/proteins, many of which are known to be causative for human microcephaly, including Eomesodermin/T-box Brain Protein 2 (EOMES/TBR2) and Neuronal Differentiation 2 (NEUROD2). Disturbances of distinct molecular pathways in neural progenitors and post-mitotic neurons may contribute to complex brain phenotype of congenital ZIKV infection. Overall, this report on protein- and transcript-level dynamics enhances understanding of the ZIKV immunopathological landscape through characterization of fetal immune response in the developing brain.

Systems-level analyses of protein-protein interaction network dysfunctions via epichaperomics identify cancer-specific mechanisms of stress adaptation

Systems-level assessments of protein-protein interaction (PPI) network dysfunctions are currently out-of-reach because approaches enabling proteome-wide identification, analysis, and modulation of context-specific PPI changes in native (unengineered) cells and tissues are lacking. Herein, we take advantage of chemical binders of maladaptive scaffolding structures termed epichaperomes and develop an epichaperome-based 'omics platform, epichaperomics, to identify PPI alterations in disease. We provide multiple lines of evidence, at both biochemical and functional levels, demonstrating the importance of these probes to identify and study PPI network dysfunctions and provide mechanistically and therapeutically relevant proteome-wide insights. As proof-of-principle, we derive systems-level insight into PPI dysfunctions of cancer cells which enabled the discovery of a context-dependent mechanism by which cancer cells enhance the fitness of mitotic protein networks. Importantly, our systems levels analyses support the use of epichaperome chemical binders as therapeutic strategies aimed at normalizing PPI networks.

Genetic Networks of Alzheimer's Disease, Aging, and Longevity in Humans

Human genomic analysis and genome-wide association studies (GWAS) have identified genes that are risk factors for early and late-onset Alzheimer's disease (AD genes). Although the genetics of aging and longevity have been extensively studied, previous studies have focused on a specific set of genes that have been shown to contribute to or are a risk factor for AD. Thus, the connections among the genes involved in AD, aging, and longevity are not well understood. Here, we identified the genetic interaction networks (referred to as pathways) of aging and longevity within the context of AD by using a gene set enrichment analysis by Reactome that cross-references more than 100 bioinformatic databases to allow interpretation of the biological functions of gene sets through a wide variety of gene networks. We validated the pathways with a threshold of p-value < 1.00 × 10^-5 using the databases to extract lists of 356 AD genes, 307 aging-related (AR) genes, and 357 longevity genes. There was a broad range of biological pathways involved in AR and longevity genes shared with AD genes. AR genes identified 261 pathways within the threshold of p < 1.00 × 10^-5, of which 26 pathways (10% of AR gene pathways) were further identified by overlapping genes among AD and AR genes. The overlapped pathways included gene expression (p = 4.05 × 10^-11) including ApoE, SOD2, TP53, and TGFB1 (p = 2.84 × 10^-10); protein metabolism and SUMOylation, including E3 ligases and target proteins (p = 1.08 × 10^-7); ERBB4 signal transduction (p = 2.69 × 10^-6); the immune system, including IL-3 and IL-13 (p = 3.83 × 10^-6); programmed cell death (p = 4.36 × 10^-6); and platelet degranulation (p = 8.16 × 10^-6), among others. Longevity genes identified 49 pathways within the threshold, of which 12 pathways (24% of longevity gene pathways) were further identified by overlapping genes among AD and longevity genes. They include the immune system, including IL-3 and IL-13 (p = 7.64 × 10^-8), plasma lipoprotein assembly, remodeling and clearance (p < 4.02 × 10^-6), and the metabolism of fat-soluble vitamins (p = 1.96 × 10^-5). Thus, this study provides shared genetic hallmarks of aging, longevity, and AD backed up by statistical significance. We discuss the significant genes involved in these pathways, including TP53, FOXO, SUMOylation, IL4, IL6, APOE, and CEPT, and suggest that mapping the gene network pathways provide a useful basis for further medical research on AD and healthy aging.

Microglial TNFα orchestrates protein phosphorylation in the cortex during the sleep period and controls homeostatic sleep

Sleep intensity is adjusted by the length of previous awake time, and under tight homeostatic control by protein phosphorylation. Here, we establish microglia as a new cellular component of the sleep homeostasis circuit. Using quantitative phosphoproteomics of the mouse frontal cortex, we demonstrate that microglia-specific deletion of TNFα perturbs thousands of phosphorylation sites during the sleep period. Substrates of microglial TNFα comprise sleep-related kinases such as MAPKs and MARKs, and numerous synaptic proteins, including a subset whose phosphorylation status encodes sleep need and determines sleep duration. As a result, microglial TNFα loss attenuates the build-up of sleep need, as measured by electroencephalogram slow-wave activity and prevents immediate compensation for loss of sleep. Our data suggest that microglia control sleep homeostasis by releasing TNFα which acts on neuronal circuitry through dynamic control of phosphorylation.

Effects of ivacaftor on systemic inflammation and the plasma proteome in people with CF and G551D

BACKGROUND

Ivacaftor is a cystic fibrosis transmembrane conductance regulator (CFTR) potentiator for people with CF and the G551D mutation. We aimed to investigate the biology of CFTR modulation and systemic effects of CFTR restoration by examining changes in circulating measurements of inflammation and growth and novel proteins with ivacaftor treatment.

METHODS

Blood samples from 64 CF subjects with G551D-CFTR were analyzed for inflammatory and growth-related proteins at baseline, 1 and 6 months after ivacaftor initiation. In 30 subjects, plasma was assayed for 1,322 proteins using the SomaScan proteomic platform at baseline and 6 months post-ivacaftor. Correlations with clinical outcomes were assessed.

MEASUREMENTS AND MAIN RESULTS

Significant reductions in high mobility group box-1 protein (HMGB-1), calprotectin, serum amyloid A, and granulocyte colony-stimulating factor (G-CSF), and an increase in insulin-like growth factor (IGF-1) occurred 1 month after ivacaftor. This treatment effect was sustained at 6 months for HMGB-1 and calprotectin. Correcting for multiple comparisons in the proteomic analysis, 9 proteins (albumin, afamin, leptin, trypsin, pancreatic stone protein [PSP], pituitary adenylate cyclase-activating polypeptide-38, repulsive guidance molecule A [RGMA], calreticulin, GTPase KRas) changed significantly with ivacaftor. Proteins changing with treatment are involved in lipid digestion and transport and extracellular matrix organization biological processes. Reductions in calprotectin and G-CSF and increases in calreticulin, and RGMA correlated with improved lung function, while increasing IGF-1, leptin and afamin and decreasing PSP correlated with increased weight.

CONCLUSIONS

Ivacaftor led to changes in inflammatory, lipid digestion, and extracellular matrix proteins, lending insights into the extrapulmonary effects of CFTR modulation.

mSep: investigating physiological and immune-metabolic biomarkers in septic and healthy pregnant women to predict feto-maternal immune health - a prospective observational cohort study protocol

INTRODUCTION

Maternal sepsis remains a leading cause of death in pregnancy. Physiological adaptations to pregnancy obscure early signs of sepsis and can result in delays in recognition and treatment. Identifying biomarkers that can reliably diagnose sepsis will reduce morbidity and mortality and antibiotic overuse. We have previously identified an immune-metabolic biomarker network comprising three pathways with a >99% accuracy for detecting bacterial neonatal sepsis. In this prospective study, we will describe physiological parameters and novel biomarkers in two cohorts-healthy pregnant women and pregnant women with suspected sepsis-with the aim of mapping pathophysiological drivers and evaluating predictive biomarkers for diagnosing maternal sepsis.

METHODS AND ANALYSIS

Women aged over 18 with an ultrasound-confirmed pregnancy will be recruited to a pilot and two main study cohorts. The pilot will involve blood sample collection from 30 pregnant women undergoing an elective caesarean section. Cohort A will follow 100 healthy pregnant women throughout their pregnancy journey, with collection of blood samples from participants at routine time points in their pregnancy: week 12 'booking', week 28 and during labour. Cohort B will follow 100 pregnant women who present with suspected sepsis in pregnancy or labour and will have at least two blood samples taken during their care pathway. Study blood samples will be collected during routine clinical blood sampling. Detailed medical history and physiological parameters at the time of blood sampling will be recorded, along with the results of routine biochemical tests, including C reactive protein, lactate and white blood cell count. In addition, study blood samples will be processed and analysed for transcriptomic, lipidomic and metabolomic analyses and both qualitative and functional immunophenotyping.

ETHICS AND DISSEMINATION

Ethical approval has been obtained from the Wales Research Ethics Committee 2 (SPON1752-19, 30 October 2019).

TRIAL REGISTRATION NUMBER

NCT05023954.

Identification of the Potential Molecular Mechanisms Linking RUNX1 Activity with Nonalcoholic Fatty Liver Disease, by Means of Systems Biology

Nonalcoholic fatty liver disease (NAFLD) is the most prevalent chronic hepatic disease; nevertheless, no definitive diagnostic method exists yet, apart from invasive liver biopsy, and nor is there a specific approved treatment. Runt-related transcription factor 1 (RUNX1) plays a major role in angiogenesis and inflammation; however, its link with NAFLD is unclear as controversial results have been reported. Thus, the objective of this work was to determine the proteins involved in the molecular mechanisms between RUNX1 and NAFLD, by means of systems biology. First, a mathematical model that simulates NAFLD pathophysiology was generated by analyzing Anaxomics databases and reviewing available scientific literature. Artificial neural networks established NAFLD pathophysiological processes functionally related to RUNX1: hepatic insulin resistance, lipotoxicity, and hepatic injury-liver fibrosis. Our study indicated that RUNX1 might have a high relationship with hepatic injury-liver fibrosis, and a medium relationship with lipotoxicity and insulin resistance motives. Additionally, we found five RUNX1-regulated proteins with a direct involvement in NAFLD motives, which were NFκB1, NFκB2, TNF, ADIPOQ, and IL-6. In conclusion, we suggested a relationship between RUNX1 and NAFLD since RUNX1 seems to regulate NAFLD molecular pathways, posing it as a potential therapeutic target of NAFLD, although more studies in this field are needed.

Focal adhesion kinase inhibitor TAE226 combined with Sorafenib slows down hepatocellular carcinoma by multiple epigenetic effects

Background

Hepatocellular carcinoma (HCC) is one of the most common and lethal malignant tumours worldwide. Sorafenib (SOR) is one of the most effective single-drug systemic therapy against advanced HCC, but the identification of novel combination regimens for a continued improvement in overall survival is a big challenge. Recent studies highlighted the crucial role of focal adhesion kinase (FAK) in HCC growth. The aim of this study was to investigate the antitumor effects of three different FAK inhibitors (FAKi), alone or in combination with SOR, using in vitro and in vivo models of HCC.

Methods

The effect of PND1186, PF431396, TAE226 on cell viability was compared to SOR. Among them TAE226, emerging as the most effective FAKi, was tested alone or in combination with SOR using 2D/3D human HCC cell line cultures and HCC xenograft murine models. The mechanisms of action were assessed by gene/protein expression and imaging approaches, combined with high-throughput methods.

Results

TAE226 was the more effective FAKi to be combined with SOR against HCC. Combined TAE226 and SOR treatment reduced HCC growth both in vitro and in vivo by affecting tumour-promoting gene expression and inducing epigenetic changes via dysregulation of FAK nuclear interactome. We characterized a novel nuclear functional interaction between FAK and the NuRD complex. TAE226-mediated FAK depletion and SOR-promoted MAPK down-modulation caused a decrease in the nuclear amount of HDAC1/2 and a consequent increase of the histone H3 lysine 27 acetylation, thus counteracting histone H3 lysine 27 trimethylation.

Conclusions

Altogether, our findings provide the first evidence that TAE226 combined with SOR efficiently reduces HCC growth in vitro and in vivo. Also, our data highlight that deep analysis of FAK nuclear interactome may lead to the identification of new promising targets for HCC therapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02154-8.

Identifying digenic disease genes via machine learning in the Undiagnosed Diseases Network

Rare diseases affect millions of people worldwide, and discovering their genetic causes is challenging. More than half of the individuals analyzed by the Undiagnosed Diseases Network (UDN) remain undiagnosed. The central hypothesis of this work is that many of these rare genetic disorders are caused by multiple variants in more than one gene. However, given the large number of variants in each individual genome, experimentally evaluating combinations of variants for potential to cause disease is currently infeasible. To address this challenge, we developed the digenic predictor (DiGePred), a random forest classifier for identifying candidate digenic disease gene pairs by features derived from biological networks, genomics, evolutionary history, and functional annotations. We trained the DiGePred classifier by using DIDA, the largest available database of known digenic-disease-causing gene pairs, and several sets of non-digenic gene pairs, including variant pairs derived from unaffected relatives of UDN individuals. DiGePred achieved high precision and recall in cross-validation and on a held-out test set (PR area under the curve > 77%), and we further demonstrate its utility by using digenic pairs from the recent literature. In contrast to other approaches, DiGePred also appropriately controls the number of false positives when applied in realistic clinical settings. Finally, to enable the rapid screening of variant gene pairs for digenic disease potential, we freely provide the predictions of DiGePred on all human gene pairs. Our work enables the discovery of genetic causes for rare non-monogenic diseases by providing a means to rapidly evaluate variant gene pairs for the potential to cause digenic disease.

Interaction Network Provides Clues on the Role of BCAR1 in Cellular Response to Changes in Gravity

When culturing cells in space or under altered gravity conditions on Earth to investigate the impact of gravity, their adhesion and organoid formation capabilities change. In search of a target where the alteration of gravity force could have this impact, we investigated p130cas/BCAR1 and its interactions more thoroughly, particularly as its activity is sensitive to applied forces. This protein is well characterized regarding its role in growth stimulation and adhesion processes. To better understand BCAR1′s force-dependent scaffolding of other proteins, we studied its interactions with proteins we had detected by proteome analyses of MCF-7 breast cancer and FTC-133 thyroid cancer cells, which are both sensitive to exposure to microgravity and express BCAR1. Using linked open data resources and our experiments, we collected comprehensive information to establish a semantic knowledgebase and analyzed identified proteins belonging to signaling pathways and their networks. The results show that the force-dependent phosphorylation and scaffolding of BCAR1 influence the structure, function, and degradation of intracellular proteins as well as the growth, adhesion and apoptosis of cells similarly to exposure of whole cells to altered gravity. As BCAR1 evidently plays a significant role in cell responses to gravity changes, this study reveals a clear path to future research performing phosphorylation experiments on BCAR1.

miRNome and Functional Network Analysis of PGRMC1 Regulated miRNA Target Genes Identify Pathways and Biological Functions Associated With Triple Negative Breast Cancer

Background

Increased expression of the progesterone receptor membrane component 1, a heme and progesterone binding protein, is frequently found in triple negative breast cancer tissue. The basis for the expression of PGRMC1 and its regulation on cellular signaling mechanisms remain largely unknown. Therefore, we aim to study microRNAs that target selective genes and mechanisms that are regulated by PGRMC1 in TNBCs.

Methods

To identify altered miRNAs, whole human miRNome profiling was performed following AG-205 treatment and PGRMC1 silencing. Network analysis identified miRNA target genes while KEGG, REACTOME and Gene ontology were used to explore altered signaling pathways, biological processes, and molecular functions.

Results

KEGG term pathway analysis revealed that upregulated miRNAs target specific genes that are involved in signaling pathways that play a major role in carcinogenesis. While multiple downregulated miRNAs are known oncogenes and have been previously demonstrated to be overexpressed in a variety of cancers. Overlapping miRNA target genes associated with KEGG term pathways were identified and overexpression/amplification of these genes was observed in invasive breast carcinoma tissue from TCGA. Further, the top two genes (CCND1 and YWHAZ) which are highly genetically altered are also associated with poorer overall survival.

Conclusions

Thus, our data demonstrates that therapeutic targeting of PGRMC1 in aggressive breast cancers leads to the activation of miRNAs that target overexpressed genes and deactivation of miRNAs that have oncogenic potential.

Interferon regulatory factor 7 impairs cellular metabolism in aging adipose-derived stromal cells

Dysregulated immunity and widespread metabolic dysfunctions are the most relevant hallmarks of the passing of time over the course of adult life, and their combination at midlife is strongly related to increased vulnerability to diseases; however, the causal connection between them remains largely unclear. By combining multi-omics and functional analyses of adipose-derived stromal cells established from young (1 month) and midlife (12 months) mice, we show that an increase in expression of interferon regulatory factor 7 (IRF7) during adult life drives major metabolic changes, which include impaired mitochondrial function, altered amino acid biogenesis and reduced expression of genes involved in branched-chain amino acid (BCAA) degradation. Our results draw a new paradigm of aging as the 'sterile' activation of a cell-autonomous pathway of self-defense and identify a crucial mediator of this pathway, IRF7, as driver of metabolic dysfunction with age.

Out-of-Field Hippocampus from Partial-Body Irradiated Mice Displays Changes in Multi-Omics Profile and Defects in Neurogenesis

The brain undergoes ionizing radiation exposure in many clinical situations, particularly during radiotherapy for brain tumors. The critical role of the hippocampus in the pathogenesis of radiation-induced neurocognitive dysfunction is well recognized. The goal of this study is to test the potential contribution of non-targeted effects in the detrimental response of the hippocampus to irradiation and to elucidate the mechanisms involved. C57Bl/6 mice were whole body (WBI) or partial body (PBI) irradiated with 0.1 or 2.0 Gy of X-rays or sham irradiated. PBI consisted of the exposure of the lower third of the mouse body, whilst the upper two thirds were shielded. Hippocampi were collected 15 days or 6 months post-irradiation and a multi-omics approach was adopted to assess the molecular changes in non-coding RNAs, proteins and metabolic levels, as well as histological changes in the rate of hippocampal neurogenesis. Notably, at 2.0 Gy the pattern of early molecular and histopathological changes induced in the hippocampus at 15 days following PBI were similar in quality and quantity to the effects induced by WBI, thus providing a proof of principle of the existence of out-of-target radiation response in the hippocampus of conventional mice. We detected major alterations in DAG/IP3 and TGF-β signaling pathways as well as in the expression of proteins involved in the regulation of long-term neuronal synaptic plasticity and synapse organization, coupled with defects in neural stem cells self-renewal in the hippocampal dentate gyrus. However, compared to the persistence of the WBI effects, most of the PBI effects were only transient and tended to decrease at 6 months post-irradiation, indicating important mechanistic difference. On the contrary, at low dose we identified a progressive accumulation of molecular defects that tended to manifest at later post-irradiation times. These data, indicating that both targeted and non-targeted radiation effects might contribute to the pathogenesis of hippocampal radiation-damage, have general implications for human health.

A developmental lineage-based gene co-expression network for mouse pancreatic β-cells reveals a role for Zfp800 in pancreas development

To gain a deeper understanding of pancreatic β-cell development, we used iterative weighted gene correlation network analysis to calculate a gene co-expression network (GCN) from 11 temporally and genetically defined murine cell populations. The GCN, which contained 91 distinct modules, was then used to gain three new biological insights. First, we found that the clustered protocadherin genes are differentially expressed during pancreas development. Pcdhγ genes are preferentially expressed in pancreatic endoderm, Pcdhβ genes in nascent islets, and Pcdhα genes in mature β-cells. Second, after extracting sub-networks of transcriptional regulators for each developmental stage, we identified 81 zinc finger protein (ZFP) genes that are preferentially expressed during endocrine specification and β-cell maturation. Third, we used the GCN to select three ZFPs for further analysis by CRISPR mutagenesis of mice. Zfp800 null mice exhibited early postnatal lethality, and at E18.5 their pancreata exhibited a reduced number of pancreatic endocrine cells, alterations in exocrine cell morphology, and marked changes in expression of genes involved in protein translation, hormone secretion and developmental pathways in the pancreas. Together, our results suggest that developmentally oriented GCNs have utility for gaining new insights into gene regulation during organogenesis.

Differential Effects of Trp53 Alterations in Murine Colorectal Cancer

BACKGROUND

Colorectal cancer (CRC) development is a multi-step process resulting in the accumulation of genetic alterations. Despite its high incidence, there are currently no mouse models that accurately recapitulate this process and mimic sporadic CRC. We aimed to develop and characterize a genetically engineered mouse model (GEMM) of Apc/Kras/Trp53 mutant CRC, the most frequent genetic subtype of CRC.

METHODS

Tumors were induced in mice with conditional mutations or knockouts in Apc, Kras, and Trp53 by a segmental adeno-cre viral infection, monitored via colonoscopy and characterized on multiple levels via immunohistochemistry and next-generation sequencing.

RESULTS

The model accurately recapitulates human colorectal carcinogenesis clinically, histologically and genetically. The Trp53 R172H hotspot mutation leads to significantly increased metastatic capacity. The effects of Trp53 alterations, as well as the response to treatment of this model, are similar to human CRC. Exome sequencing revealed spontaneous protein-modifying alterations in multiple CRC-related genes and oncogenic pathways, resulting in a genetic landscape resembling human CRC.

CONCLUSIONS

This model realistically mimics human CRC in many aspects, allows new insights into the role of TP53 in CRC, enables highly predictive preclinical studies and demonstrates the value of GEMMs in current translational cancer research and drug development.

Developmental effects of maternal smoking during pregnancy on the human frontal cortex transcriptome

Cigarette smoking during pregnancy is a major public health concern. While there are well-described consequences in early child development, there is very little known about the effects of maternal smoking on human cortical biology during prenatal life. We therefore performed a genome-wide differential gene expression analysis using RNA sequencing (RNA-seq) on prenatal (N = 33; 16 smoking-exposed) as well as adult (N = 207; 57 active smokers) human postmortem prefrontal cortices. Smoking exposure during the prenatal period was directly associated with differential expression of 14 genes; in contrast, during adulthood, despite a much larger sample size, only two genes showed significant differential expression (FDR < 10%). Moreover, 1,315 genes showed significantly different exposure effects between maternal smoking during pregnancy and direct exposure in adulthood (FDR < 10%)-these differences were largely driven by prenatal differences that were enriched for pathways previously implicated in addiction and synaptic function. Furthermore, prenatal and age-dependent differentially expressed genes were enriched for genes implicated in non-syndromic autism spectrum disorder (ASD) and were differentially expressed as a set between patients with ASD and controls in postmortem cortical regions. These results underscore the enhanced sensitivity to the biological effect of smoking exposure in the developing brain and offer insight into how maternal smoking during pregnancy affects gene expression in the prenatal human cortex. They also begin to address the relationship between in utero exposure to smoking and the heightened risks for the subsequent development of neuropsychiatric disorders.

Plasma Biomarkers of Reticular Pseudodrusen and the Risk of Progression to Advanced Age-Related Macular Degeneration

Purpose

To determine, using an aptamer-based technology in patients with intermediate age-related macular degeneration (AMD), (1) if there is a difference in plasma levels of 4979 proteins in patients with and without reticular pseudodrusen (RPD), and (2) if plasma levels of proteins are related to time to conversion to advanced AMD.

Methods

Patients with intermediate AMD and RPD were identified from an AMD registry. Relative concentrations of each protein were log (base 2) transformed and compared between patients with and without RPD using linear regression. A Cox proportional hazards survival model was fit to each aptamer to quantify associations with time to conversion. A pathway analysis was conducted in converters versus non-converters using the Reactome database.

Results

Of the 109 intermediate AMD patients, 39 had bilateral RPD (36%). Two proteins, TCL1A and CNDP1, were lower in patients in the intermediate AMD group with RPD. Twenty-one patients converted to advanced AMD with a median time to conversion of 25.2 months (range, 2.3-48.5 months) and median follow-up time in non-converters of 26.4 months (range, 0.03-49.7 months). Several proteins (lysozyme C, TFF3, RNAS6, and SAP3) distinguished patients who converted from those who did not convert to advanced AMD. The top conversion pathways included tumor necrosis factors bind their physiological receptors, digestion and absorption, signaling by activin, and signaling by TGF-β family members.

Conclusions

We identified a protein signature related to RPD, as well as to conversion to advanced AMD. The pathway analysis suggests that dysfunction of critical systemic pathways may have links to conversion to advanced AMD.

Translational Relevance

Biomarkers identified in plasma likely reflect systemic alterations in protein expression in patients with intermediate AMD.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

The Molecular Processes in the Trabecular Meshwork After Exposure to Corticosteroids and in Corticosteroid-Induced Ocular Hypertension

Purpose

To identify processes that contribute to corticosteroid-induced ocular hypertension and candidate target genes for treatment.

Methods

A systematic search identified five human microarray datasets investigating the effect of dexamethasone versus a control medium on trabecular meshwork (TM) tissue. After thorough quality control, samples of low quality were removed, and the datasets were integrated. Additionally, a bovine RNA-sequencing dataset allowed to investigate differences in gene expression profiling between cows with and without corticosteroid-induced ocular hypertension (responders vs. nonresponders). The obtained datasets were used as input for parallel pathway analyses. Significantly changed pathways were clustered into functional categories and the results were further investigated. A network visualizing the differences between the responders and nonresponders was created.

Results

Seven functional pathway clusters were found to be significantly changed in TM cells exposed to dexamethasone versus a control medium and in TM cells of responders versus nonresponders: collagen, extracellular matrix, adhesion, WNT-signaling, inflammation, adipogenesis, and glucose metabolism. In addition, cell cycle and senescence were only significantly changed in responders versus nonresponders. The network of the differential gene expression between responders and nonresponders shows many connections between the identified processes via shared genes.

Conclusions

Nine functional pathway clusters synthesize the molecular response to dexamethasone exposure in TM cells and are likely to be involved in the pathogenesis of corticosteroid-induced ocular hypertension.

Laser capture microdissection coupled mass spectrometry (LCM-MS) for spatially resolved analysis of formalin-fixed and stained human lung tissues

BACKGROUND

Haematoxylin and eosin (H&E)-which respectively stain nuclei blue and other cellular and stromal material pink-are routinely used for clinical diagnosis based on the identification of morphological features. A richer characterization can be achieved by laser capture microdissection coupled to mass spectrometry (LCM-MS), giving an unbiased assay of the proteins that make up the tissue. However, the process of fixing and H&E staining of tissues provides challenges with standard sample preparation methods for mass spectrometry, resulting in low protein yield. Here we describe a microproteomics technique to analyse H&E-stained, formalin-fixed paraffin-embedded (FFPE) tissues.

METHODS

Herein, we utilize heat extraction, physical disruption, and in column digestion for the analysis of H&E stained FFPE tissues. Micro-dissected morphologically normal human lung alveoli (0.082 mm3) and human lung blood vessels (0.094 mm3) from FFPE-fixed H&E-stained sections from Idiopathic Pulmonary Fibrosis (IPF) specimens (n = 3 IPF specimens) were then subject to a qualitative and then quantitative proteomics approach using BayesENproteomics. In addition, we tested the sensitivity of this method by processing and analysing a range of micro-dissected human lung blood vessel tissue volumes.

RESULTS

This approach yields 1252 uniquely expressed proteins (at a protein identification threshold of 3 unique peptides) with 892 differentially expressed proteins between these regions. In accord with prior knowledge, our methodology approach confirms that human lung blood vessels are enriched with smoothelin, CNN1, ITGA7, MYH11, TAGLN, and PTGIS; whereas morphologically normal human lung alveoli are enriched with cytokeratin-7, -8, -18, -19, 14, and -17. In addition, we identify a total of 137 extracellular matrix (ECM) proteins and immunohistologically validate that laminin subunit beta-1 localizes to morphologically normal human lung alveoli and tenascin localizes to human lung blood vessels. Lastly, we show that this micro-proteomics technique can be applied to tissue volumes as low as 0.0125 mm3.

CONCLUSION

Herein we show that our multistep sample preparation methodology of LCM-MS can identify distinct, characteristic proteomic compositions of anatomical features within complex fixed and stained tissues.

BayesENproteomics: Bayesian Elastic Nets for Quantification of Peptidoforms in Complex Samples

Multivariate regression modelling provides a statistically powerful means of quantifying the effects of a given treatment while compensating for sources of variation and noise, such as variability between human donors and the behavior of different peptides during mass spectrometry. However, methods to quantify endogenous post-translational modifications (PTMs) are typically reliant on summary statistical methods that fail to consider sources of variability such as changes in the levels of the parent protein. Here, we compare three multivariate regression methods, including a novel Bayesian elastic net algorithm (BayesENproteomics) that enables assessment of relative protein abundances while also quantifying identified PTMs for each protein. We tested the ability of these methods to accurately quantify expression of proteins in a mixed-species benchmark experiment and to quantify synthetic PTMs induced by stable isotope labelling. Finally, we extended our regression pipeline to calculate fold changes at the pathway level, providing a complement to commonly used enrichment analysis. Our results show that BayesENproteomics can quantify changes to protein levels across a broad dynamic range while also accurately quantifying PTM and pathway-level fold changes.

The recent progress and therapy in endometriosis-associated ovarian cancer

Endometriosis-associated ovarian cancers (EAOCs) including endometrioid and clear cell ovarian carcinoma are subgroups of epithelial ovarian carcinomas (EOCs), which is generally acknowledged as the most lethal gynecological malignancy. Endometriosis (ES), a common clinical disease among women, presents with clinical symptoms of pelvic pain, infertility, or adnexal masses with the formation of endometrioma. It has long been considered to be a potential risk factor for developing EOCs, mainly of endometrioid and clear cell subtypes. Here, we compiled data from previous researches on deregulated molecular functions among ES and EOCs using gene set-based integrative analysis to decipher molecular and genetic relationships between ovarian ES and EOCs, especially EAOCs. We conclude that epidermal growth factor receptor (ERBB) and Phosphoinositide 3-kinases (PI3K)-related pathways are important in the carcinogenesis of type I EOCs, including clear cell, endometrioid, and mucinous ovarian carcinoma. Dysfunctional molecular pathways, such as deregulated oxidoreductase activity, metabolism, hormone activity, inflammatory response, innate immune response, and cell-cell signaling, played key roles in the malignant transformation of EAOCs. Nine genes related to inflammasome complex and inflammasome-related pathway were identified, indicating the importance of inflammation/immunity in EAOC transformation. We also collect progressive treatments of EAOC focused on targeted therapies and immunotherapy so far. This summarized information can contribute toward effective detection and treatment of EAOCs in the future.

The influence of rare variants in circulating metabolic biomarkers

Circulating metabolite levels are biomarkers for cardiovascular disease (CVD). Here we studied, association of rare variants and 226 serum lipoproteins, lipids and amino acids in 7,142 (discovery plus follow-up) healthy participants. We leveraged the information from multiple metabolite measurements on the same participants to improve discovery in rare variant association analyses for gene-based and gene-set tests by incorporating correlated metabolites as covariates in the validation stage. Gene-based analysis corrected for the effective number of tests performed, confirmed established associations at APOB, APOC3, PAH, HAL and PCSK (p<1.32x10-7) and identified novel gene-trait associations at a lower stringency threshold with ACSL1, MYCN, FBXO36 and B4GALNT3 (p<2.5x10-6). Regulation of the pyruvate dehydrogenase (PDH) complex was associated for the first time, in gene-set analyses also corrected for effective number of tests, with IDL and LDL parameters, as well as circulating cholesterol (pMETASKAT<2.41x10-6). In conclusion, using an approach that leverages metabolite measurements obtained in the same participants, we identified novel loci and pathways involved in the regulation of these important metabolic biomarkers. As large-scale biobanks continue to amass sequencing and phenotypic information, analytical approaches such as ours will be useful to fully exploit the copious amounts of biological data generated in these efforts.

Comparison and evaluation of integrative methods for the analysis of multilevel omics data: a study based on simulated and experimental cancer data

Integrative analysis aims to identify the driving factors of a biological process by the joint exploration of data from multiple cellular levels. The volume of omics data produced is constantly increasing, and so too does the collection of tools for its analysis. Comparative studies assessing performance and the biological value of results, however, are rare but in great demand. We present a comprehensive comparison of three integrative analysis approaches, sparse canonical correlation analysis (sCCA), non-negative matrix factorization (NMF) and logic data mining MicroArray Logic Analyzer (MALA), by applying them to simulated and experimental omics data. We find that sCCA and NMF are able to identify differential features in simulated data, while the Logic Data Mining method, MALA, falls short. Applied to experimental data, we show that MALA performs best in terms of sample classification accuracy, and in general, the classification power of prioritized feature sets is high (97.1-99.5% accuracy). The proportion of features identified by at least one of the other methods, however, is approximately 60% for sCCA and NMF and nearly 30% for MALA, and the proportion of features jointly identified by all methods is only around 16%. Similarly, the congruence on functional levels (Gene Ontology, Reactome) is low. Furthermore, the agreement of identified feature sets with curated gene signatures relevant to the investigated disease is modest. We discuss possible reasons for the moderate overlap of identified feature sets with each other and with curated cancer signatures. The R code to create simulated data, results and figures is provided at https://github.com/ThallingerLab/IamComparison.

Comprehensive bioinformatics analysis of trabecular meshwork gene expression data to unravel the molecular pathogenesis of primary open-angle glaucoma

PURPOSE

Performing bioinformatics analyses using trabecular meshwork (TM) gene expression data in order to further elucidate the molecular pathogenesis of primary open-angle glaucoma (POAG), and to identify candidate target genes.

METHODS

A systematic search in Gene Expression Omnibus and ArrayExpress was conducted, and quality control and preprocessing of the data was performed with ArrayAnalysis.org. Molecular pathway overrepresentation analysis was performed with PathVisio using pathway content from three pathway databases: WikiPathways, KEGG and Reactome. In addition, Gene Ontology (GO) analysis was performed on the gene expression data. The significantly changed pathways were clustered into functional categories which were combined into a network of connected genes.

RESULTS

Ninety-two significantly changed pathways were clustered into five functional categories: extracellular matrix (ECM), inflammation, complement activation, senescence and Rho GTPase signalling. ECM included pathways involved in collagen, actin and cell-matrix interactions. Inflammation included pathways entailing NF-κB and arachidonic acid. The network analysis showed that several genes overlap between the inflammation cluster on the one hand, and the ECM, complement activation and senescence clusters on the other hand. GO analysis, identified additional clusters, related to development and corticosteroids.

CONCLUSION

This study provides an overview of the processes involved in the molecular pathogenesis of POAG in the TM. The results show good face validity and confirm findings from histological, biochemical, genome-wide association and transcriptomics studies. The identification of known points of action for drugs, such as Rho GTPase, arachidonic acid, NF-κB, prostaglandins and corticosteroid clusters, supports the value of this approach to identify potential drug targets.

Early angiogenic proteins associated with high risk for bronchopulmonary dysplasia and pulmonary hypertension in preterm infants

Early pulmonary vascular disease in preterm infants is associated with the subsequent development of bronchopulmonary dysplasia (BPD) and pulmonary hypertension (PH); however, mechanisms that contribute to or identify infants with increased susceptibility for BPD and/or PH are incompletely understood. Therefore, we tested if changes in circulating angiogenic peptides during the first week of life are associated with the later development of BPD and/or PH. We further sought to determine alternate peptides and related signaling pathways with the risk for BPD or PH. We prospectively enrolled infants with gestational age <34 wk and collected blood samples during their first week of life. BPD and PH were assessed at 36 wk postmenstrual age. Samples were assayed for each of the 1,121 peptides included in the SOMAscan scan technology, with subsequent pathway analysis. Of 102 infants in the study, 82 had BPD, and 13 had PH. Multiple angiogenic proteins (PF-4, VEGF121, ANG-1, bone morphogenetic protein 10 [BMP10], hepatocyte growth factor (HGF), ANG-2) were associated with the subsequent diagnosis of BPD; and FGF-19, PF-4, connective tissue activating peptide (CTAP)-III, and PDGF-AA levels were associated with BPD severity. Early increases in BMP10 was strongly associated with the late risk for BPD and PH. We found that early alterations of circulating angiogenic peptides and others were associated with the subsequent development of BPD. We further identified peptides that were associated with BPD severity and BPD-associated PH, including BMP10. We speculate that proteomic biomarkers during the first week of life may identify infants at risk for BPD and/or PH to enhance care and research.

Impact of short- and long-term electrically induced muscle exercise on gene signaling pathways, gene expression, and PGC1a methylation in men with spinal cord injury

Exercise attenuates the development of chronic noncommunicable diseases (NCDs). Gene signaling pathway analysis offers an opportunity to discover if electrically induced muscle exercise regulates key pathways among people living with spinal cord injury (SCI). We examined short-term and long-term durations of electrically induced skeletal muscle exercise on complex gene signaling pathways, specific gene regulation, and epigenetic tagging of PGC1a, a major transcription factor in skeletal muscle of men with SCI. After short- or long-term electrically induced exercise training, participants underwent biopsies of the trained and untrained muscles. RNA was hybridized to an exon microarray and analyzed by a gene set enrichment analysis. We discovered that long-term exercise training regulated the Reactome gene sets for metabolism (38 gene sets), cell cycle (36 gene sets), disease (27 gene sets), gene expression and transcription (22 gene sets), organelle biogenesis (4 gene sets), cellular response to stimuli (8 gene sets), immune system (8 gene sets), vesicle-mediated transport (4 gene sets), and transport of small molecules (3 gene sets). Specific gene expression included: oxidative catabolism of glucose including PDHB (P < 0.001), PDHX (P < 0.001), MPC1 (P < 0.009), and MPC2 (P < 0.007); Oxidative phosphorylation genes including SDHA (P < 0.006), SDHB (P < 0.001), NDUFB1 (P < 0.002), NDUFA2 (P < 0.001); transcription genes including PGC1α (P < 0.030) and PRKAB2 (P < 0.011); hypertrophy gene MSTN (P < 0.001); and the myokine generating FNDC5 gene (P < 0.008). Long-term electrically induced exercise demethylated the major transcription factor PGC1a. Taken together, these findings support that long-term electrically induced muscle activity regulates key pathways associated with muscle health and systemic metabolism.

Pulmonary Aptamer Signatures in Children's Interstitial and Diffuse Lung Disease

Rationale: Biomarker signatures are needed in children with children's interstitial and diffuse lung disease (chILD) to improve diagnostic approaches, increase our understanding of disease pathogenesis, monitor disease progression, and develop new treatment strategies. Proteomic technology using SOMAmer (Slow Off-rate Modified Aptamer) nucleic acid-based protein-binding reagents allows for biomarker discovery.Objectives: We hypothesized that proteins and protein pathways in BAL fluid (BALF) would distinguish children with neuroendocrine cell hyperplasia of infancy (NEHI), surfactant dysfunction mutations, and other chILD diagnoses and control subjects.Methods: BALF was collected for clinical indications and banked in patients with chILD and disease control subjects using standardized protocols over 10 years. BALF supernatant was analyzed using an aptamer assay to measure 1,129 protein levels. Protein levels were compared between groups using an ANOVA and adjusted for multiple comparisons using false discovery rate. Proteins were classified into pathways. Hierarchical clustering was used to define endotypes in the group of children with NEHI.Measurements and Main Results: After correcting for multiple testing, children with NEHI (n = 22) had 202 aptamers that were significantly different (P < 0.05) in BALF compared with control subjects (n = 9). Children with surfactant mutation (n = 8) had 51 aptamers significantly different (P < 0.05) in BALF compared with control subjects (n = 9). Proteins associated with pulmonary fibrosis and inflammation were associated with the surfactant dysfunction group but not the NEHI group. Using hierarchical clustering analysis, two distinct NEHI endotypes were identified.Conclusions: Distinct proteins and protein pathways can be determined from BALF of children with chILD, and these hold promise to further our understanding of chILD.

Novel TMEM173 Mutation and the Role of Disease Modifying Alleles

Upon binding to pathogen or self-derived cytosolic nucleic acids cyclic GMP-AMP synthase (cGAS) triggers the production of cGAMP that further activates transmembrane protein STING. Upon activation STING translocates from ER via Golgi to vesicles. Monogenic STING gain-of-function mutations cause early-onset type I interferonopathy, with disease presentation ranging from fatal vasculopathy to mild chilblain lupus. Molecular mechanisms underlying the variable phenotype-genotype correlation are presently unclear. Here, we report a novel gain-of-function G207E STING mutation causing a distinct phenotype with alopecia, photosensitivity, thyroid dysfunction, and features of STING-associated vasculopathy with onset in infancy (SAVI), such as livedo reticularis, skin vasculitis, nasal septum perforation, facial erythema, and bacterial infections. Polymorphism in TMEM173 and IFIH1 showed variable penetrance in the affected family, implying contribution to varying phenotype spectrum. The G207E mutation constitutively activates inflammation-related pathways in vitro, and causes aberrant interferon signature and inflammasome activation in patient PBMCs. Treatment with Janus kinase 1 and 2 (JAK1/2) inhibitor baricitinib was beneficiary for a vasculitic ulcer, induced hair regrowth and improved overall well-being in one patient. Protein-protein interactions propose impaired cellular trafficking of G207E mutant. These findings reveal the molecular landscape of STING and propose common polymorphisms in TMEM173 and IFIH1 as likely modifiers of the phenotype.

Prader-Willi syndrome and Angelman syndrome: Visualisation of the molecular pathways for two chromosomal disorders

Objectives: Prader-Willi syndrome (PWS) and Angelman syndrome (AS) are two syndromes that are caused by the same chromosomal deletion on 15q11.2-q13. Due to methylation patterns, different genes are responsible for the two distinct phenotypes resulting in the disorders. Patients of both disorders exhibit hypotonia in neonatal stage, delay in development and hypopigmentation. Typical features for PWS include hyperphagia, which leads to obesity, the major cause of mortality, and hypogonadism. In AS, patients suffer from a more severe developmental delay, they have a distinctive behaviour that is often described as unnaturally happy, and a tendency for epileptic seizures. For both syndromes, we identified and visualised molecular downstream pathways of the deleted genes that could give insight on the development of the clinical features.Methods: This was done by consulting literature, genome browsers and pathway databases to identify molecular interactions and to construct downstream pathways.Results: A pathway visualisation was created and uploaded to the open pathway database WikiPathways covering all molecular pathways that were found.Conclusions: The visualisation of the downstream pathways of PWS- and AS-deleted genes shows that some of the typical symptoms are caused by multiple genes and reveals critical gaps in the current knowledge.

Ionizing radiation induces endothelial transdifferentiation of glioblastoma stem-like cells through the Tie2 signaling pathway

Glioblastomas (GBM) are brain tumors with a poor prognosis despite treatment that combines surgical resection and radio-chemotherapy. These tumors are characterized by abundant vascularization and significant cellular heterogeneity including GBM stem-like cells (GSC) which contribute to tumor aggressiveness, resistance, and recurrence. Recent data has demonstrated that GSC are directly involved in the formation of new vessels via their transdifferentiation into Tumor Derived Endothelial Cells (TDEC). We postulate that cellular stress such as ionizing radiation (IR) could enhance the transdifferentiation of GSC into TDEC. GSC neurospheres isolated from 3 different patients were irradiated or not and were then transdifferentiated into TDEC. In fact, TDEC obtained from irradiated GSC (TDEC IR+) migrate more towards VEGF, form more pseudotubes in Matrigel^TM in vitro and develop more functional blood vessels in Matrigel^TM plugs implanted in Nude mice than TDEC obtained from non-irradiated GSC. Transcriptomic analysis allows us to highlight an overexpression of Tie2 in TDEC IR+. All IR-induced effects on TDEC were abolished by using a Tie2 kinase inhibitor, which confirms the role of the Tie2 signaling pathway in this process. Finally, by analyzing Tie2 expression in patient GBMs by immunohistochemistry, we demonstrated that the number of Tie2+ vessels increases in recurrent GBM compared with matched untreated tumors. In conclusion, we demonstrate that IR potentiates proangiogenic features of TDEC through the Tie2 signaling pathway, which indicates a new pathway of treatment-induced tumor adaptation. New therapeutic strategies that associate standard treatment and a Tie2 signaling pathway inhibitor should be considered for future trials.

Transfusion Strategies are Associated with Epigenetic Changes Following Blunt Trauma

INTRODUCTION

Epigenetics has been identified in multiple diseases. The effect of transfusion strategy on epigenetics is unknown. We hypothesized that expression of epigenetic regulating genes would be associated with resuscitation strategy following blunt trauma.

METHODS

Retrospective study using the inflammation in host response to injury (glue grant) dataset. Volume transfused over 24 h of packed red blood cells (PRBC), fresh frozen plasma (FFP), platelets (PLT) as well as crystalloids was extracted along with leucocyte microarray data of genes with known epigenetic modulating activity from day 1 after injury.Principal component analysis (PCA) was used to identify principal components (PC) within the transcriptomic dataset. Multiple regression associated these with volume of blood products and crystalloids while controlling for confounders. Genes co-expressed with genes central genes in the identified PCs were furthermore subjected to pathway analysis using the reactome database.

RESULTS

PCA identified seven components. PRBC and crystalloid volumes were positively associated with PC1, implicating histone acetylation (HAT1), DNA and histone methylation (KDM6B, SET1DB) and histone phosphorylation (RPS6KA5). Conversely, PLT volume was negatively associated with PC1,5 and 6, implicating DNA methylation (DNMT1) as well as histone acetylation (HAT1) and phosphorylation (RPSKA3).FFP was negatively associated with PC3, implicating histone methylation (SETD1B) and phosphorylation (RPS6KA5).Co-expression network analysis identified downstream pathways relevant to inflammation and the innate immune response, including Toll-like receptor, interleukin and mitogen activated protein kinase activation.

CONCLUSIONS

Resuscitation strategy was associated with epigenetic transcriptomic alterations, especially for PRBC and PLT transfusions. Downstream effects may include pertubations of pathways involved in immune signaling and cellular survival.

Petri net-based model of the human DNA base excision repair pathway

Cellular DNA is daily exposed to several damaging agents causing a plethora of DNA lesions. As a first aid to restore DNA integrity, several enzymes got specialized in damage recognition and lesion removal during the process called base excision repair (BER). A large number of DNA damage types and several different readers of nucleic acids lesions during BER pathway as well as two sub-pathways were considered in the definition of a model using the Petri net framework. The intuitive graphical representation in combination with precise mathematical analysis methods are the strong advantages of the Petri net-based representation of biological processes and make Petri nets a promising approach for modeling and analysis of human BER. The reported results provide new information that will aid efforts to characterize in silico knockouts as well as help to predict the sensitivity of the cell with inactivated repair proteins to different types of DNA damage. The results can also help in identifying the by-passing pathways that may lead to lack of pronounced phenotypes associated with mutations in some of the proteins. This knowledge is very useful when DNA damage-inducing drugs are introduced for cancer therapy, and lack of DNA repair is desirable for tumor cell death.

Altered gating of Kv1.4 in the nucleus accumbens suppresses motivation for reward

Deficient motivation contributes to numerous psychiatric disorders, including withdrawal from drug use, depression, schizophrenia, and others. Nucleus accumbens (NAc) has been implicated in motivated behavior, but it remains unclear whether motivational drive is linked to discrete neurobiological mechanisms within the NAc. To examine this, we profiled cohorts of Sprague-Dawley rats in a test of motivation to consume sucrose. We found that substantial variability in willingness to exert effort for reward was not associated with operant responding under low-effort conditions or stress levels. Instead, effort-based motivation was mirrored by a divergent NAc shell transcriptome with differential regulation at potassium and dopamine signaling genes. Functionally, motivation was inversely related to excitability of NAc principal neurons. Furthermore, neuronal and behavioral outputs associated with low motivation were linked to faster inactivation of a voltage-gated potassium channel, Kv1.4. These results raise the prospect of targeting Kv1.4 gating in psychiatric conditions associated with motivational dysfunction.

A STAT3-based gene signature stratifies glioma patients for targeted therapy

Intratumoral heterogeneity is a hallmark of glioblastoma (GBM) tumors, thought to negatively influence therapeutic outcome. Previous studies showed that mesenchymal tumors have a worse outcome than the proneural subtype. Here we focus on STAT3 as its activation precedes the proneural-mesenchymal transition. We first establish a STAT3 gene signature that stratifies GBM patients into STAT3-high and -low cohorts. STAT3 inhibitor treatment selectively mitigates STAT3-high cell viability and tumorigenicity in orthotopic mouse xenograft models. We show the mechanism underlying resistance in STAT3-low cells by combining STAT3 signature analysis with kinome screen data on STAT3 inhibitor-treated cells. This allows us to draw connections between kinases affected by STAT3 inhibitors, their associated transcription factors and target genes. We demonstrate that dual inhibition of IGF-1R and STAT3 sensitizes STAT3-low cells and improves survival in mice. Our study underscores the importance of serially profiling tumors so as to accurately target individuals who may demonstrate molecular subtype switching.

The histone chaperone HIRA promotes the induction of host innate immune defences in response to HSV-1 infection

Host innate immune defences play a critical role in restricting the intracellular propagation and pathogenesis of invading viral pathogens. Here we show that the histone H3.3 chaperone HIRA (histone cell cycle regulator) associates with promyelocytic leukaemia nuclear bodies (PML-NBs) to stimulate the induction of innate immune defences against herpes simplex virus 1 (HSV-1) infection. Following the activation of innate immune signalling, HIRA localized at PML-NBs in a Janus-Associated Kinase (JAK), Cyclin Dependent Kinase (CDK), and Sp100-dependent manner. RNA-seq analysis revealed that HIRA promoted the transcriptional upregulation of a broad repertoire of host genes that regulate innate immunity to HSV-1 infection, including those involved in MHC-I antigen presentation, cytokine signalling, and interferon stimulated gene (ISG) expression. ChIP-seq analysis revealed that PML, the principle scaffolding protein of PML-NBs, was required for the enrichment of HIRA onto ISGs, identifying a role for PML in the HIRA-dependent regulation of innate immunity to virus infection. Our data identifies independent roles for HIRA in the intrinsic silencing of viral gene expression and the induction of innate immune defences to restrict the initiation and propagation of HSV-1 infection, respectively. These intracellular host defences are antagonized by the HSV-1 ubiquitin ligase ICP0, which disrupts the stable recruitment of HIRA to infecting viral genomes and PML-NBs at spatiotemporally distinct phases of infection. Our study highlights the importance of histone chaperones to regulate multiple phases of intracellular immunity to virus infection, findings that are likely to be highly pertinent in the cellular restriction of many clinically important viral pathogens.

Host innate immune defences play critical roles in the cellular restriction of invading viral pathogens and the stimulation of adaptive immune responses. A key component in the regulation of this arm of host immunity is the rapid induction of cytokine signalling and the expression of interferon stimulated gene products (ISGs), which confer a refractory antiviral state to limit virus propagation and pathogenesis. While the signal transduction cascades that activate innate immune defences are well established, little is known about the cellular host factors that expedite the expression of this broad repertoire of antiviral host genes in response to pathogen invasion. Here we show that HIRA, a histone H3.3 chaperone, associates with PML-NBs to stimulate the induction of innate immune defences in response to HSV-1 infection. Our study highlights the importance of histone chaperones in the coordinated regulation of multiple phases of host immunity in response to pathogen invasion and identifies a key role for HIRA in the induction of innate immunity to virus infection.

Gene expression profiling meta-analysis reveals novel gene signatures and pathways shared between tuberculosis and rheumatoid arthritis

Tuberculosis (TB) is among the leading causes of death by infectious diseases. An epidemiological association between Mycobacterium tuberculosis infection and autoimmune diseases like rheumatoid arthritis (RA) has been reported but it remains unclear if there is a causal relationship, and if so, which molecular pathways and regulatory mechanisms contribute to it. Here we used a computational biology approach by global gene expression meta-analysis to identify candidate genes and pathways that may link TB and RA. Data were collected from public expression databases such as NCBI GEO. Studies were selected that analyzed mRNA-expression in whole blood or blood cell populations in human case control studies at comparable conditions. Six TB and RA datasets (41 active TB patients, 33 RA patients, and 67 healthy controls) were included in the downstream analysis. This approach allowed the identification of deregulated genes that had not been identified in the single analysis of TB or RA patients and that were co-regulated in TB and RA patients compared to healthy subjects. The genes encoding TLR5, TNFSF10/TRAIL, PPP1R16B/TIMAP, SIAH1, PIK3IP1, and IL17RA were among the genes that were most significantly deregulated in TB and RA. Pathway enrichment analysis revealed 'T cell receptor signaling pathway', 'Toll-like receptor signaling pathway,' and 'virus defense related pathways' among the pathways most strongly associated with both diseases. The identification of a common gene signature and pathways substantiates the observation of an epidemiological association of TB and RA and provides clues on the mechanistic basis of this association. Newly identified genes may be a basis for future functional and epidemiological studies.

Proteomics identifies new therapeutic targets of early-stage hepatocellular carcinoma

Hepatocellular carcinoma is the third leading cause of deaths from cancer worldwide. Infection with the hepatitis B virus is one of the leading risk factors for developing hepatocellular carcinoma, particularly in East Asia¹. Although surgical treatment may be effective in the early stages, the five-year overall rate of survival after developing this cancer is only 50-70%². Here, using proteomic and phospho-proteomic profiling, we characterize 110 paired tumour and non-tumour tissues of clinical early-stage hepatocellular carcinoma related to hepatitis B virus infection. Our quantitative proteomic data highlight heterogeneity in early-stage hepatocellular carcinoma: we used this to stratify the cohort into the subtypes S-I, S-II and S-III, each of which has a different clinical outcome. S-III, which is characterized by disrupted cholesterol homeostasis, is associated with the lowest overall rate of survival and the greatest risk of a poor prognosis after first-line surgery. The knockdown of sterol O-acyltransferase 1 (SOAT1)-high expression of which is a signature specific to the S-III subtype-alters the distribution of cellular cholesterol, and effectively suppresses the proliferation and migration of hepatocellular carcinoma. Finally, on the basis of a patient-derived tumour xenograft mouse model of hepatocellular carcinoma, we found that treatment with avasimibe, an inhibitor of SOAT1, markedly reduced the size of tumours that had high levels of SOAT1 expression. The proteomic stratification of early-stage hepatocellular carcinoma presented in this study provides insight into the tumour biology of this cancer, and suggests opportunities for personalized therapies that target it.

Proteomic Profiles in Advanced Age-Related Macular Degeneration Using an Aptamer-Based Proteomic Technology

Purpose

To explore top-ranked plasma proteins related to neovascular age-related macular degeneration (AMD) and geographic atrophy (GA), and explore pathways related to neovascular AMD and GA.

Methods

We conducted a pilot study of patients with neovascular AMD (n = 10), GA (n = 10), and age-matched cataract controls (n = 10) who were recruited into an AMD registry. We measured 4001 proteins in ethylenediaminetetraacetic acid plasma samples using an aptamer-based proteomic technology. Relative concentrations of each of 4001 proteins were log (base 2) transformed and compared between cases of neovascular AMD and GA versus controls using linear regression. Pathway analysis was conducted using pathways downloaded from Reactome.

Results

In this pilot study, higher levels of vinculin and lower levels of CD177 were found in patients with neovascular AMD compared with controls. Neuregulin-4 was higher and soluble intercellular adhesion molecule-1 was lower in patients with GA compared with controls. For neovascular AMD, cargo trafficking to the periciliary membrane, fibroblast growth factor receptor 3b ligand binding and activation, and vascular endothelial growth factor–related pathways were in the top ranked pathways. The top-ranked pathways for GA included several related to ErbB4 signaling.

Conclusions

We found different proteins and different pathways associated with neovascular AMD and GA. Vinculin and some of the top-ranked pathways have been previously associated with AMD, whereas others have not been described.

Translational Relevance

Biomarkers identified in plasma likely reflect systemic alterations in protein expression and may improve our understanding of the mechanisms leading to AMD.

Novel Application of Aptamer Proteomic Analysis in Cystic Fibrosis Bronchoalveolar Lavage Fluid

PURPOSE

Biomarkers are needed in cystic fibrosis (CF) to understand disease progression, assess response to therapy, and enrich enrollment for clinical trials. Aptamer-based proteomics have proven useful in blood samples. The aim is to evaluate proteins in bronchoalveolar lavage fluid (BALF) in CF children compared to controls and identify endotypes during CF exacerbations.

EXPERIMENTAL DESIGN

BALF is collected clinically from 50 patients with CF and nine disease controls, processed, and stored per protocol. BALF supernatants are analyzed for 1129 proteins by aptamer approach (SOMAscan proteomics platform). Proteins are compared across groups and used for pathway analysis. Endotypes are identified within the CF group.

RESULTS

CF BALF has increased concentrations of neutrophil elastase, myeloperoxidase, and decreased concentration of protein folding and host defense proteins. Pathways that distinguished CF subjects included interferon gamma signaling, membrane trafficking, and phospholipid metabolism. In the CF group, unbiased analysis of proteins identified two distinct endotypes that differed based on BALF white blood cell and neutrophil counts and detection of CF pathogens.

CONCLUSIONS AND CLINICAL RELEVANCE

Proteomic analysis of the CF airway demonstrates a complex environment of proteins and pathways. This work provides evidence that aptamer-based proteomics can differentiate between groups and can determine endotypes within CF.

Obesity modifies the stoichiometry of mitochondrial proteins in a way that is distinct to the subcellular localization of the mitochondria in skeletal muscle

BACKGROUND

Skeletal muscle mitochondrial content and function appear to be altered in obesity. Mitochondria in muscle are found in well-defined regions within cells, and they are arranged in a way that form distinct subpopulations of subsarcolemmal (SS) and intermyofibrillar (IMF) mitochondria. We sought to investigate differences in the proteomes of SS and IMF mitochondria between lean subjects and subjects with obesity.

METHODS

We performed comparative proteomic analyses on SS and IMF mitochondria isolated from muscle samples obtained from lean subjects and subjects with obesity. Mitochondria were isolated using differential centrifugation, and proteins were subjected to label-free quantitative tandem mass spectrometry analyses. Collected data were evaluated for abundance of mitochondrial proteins using spectral counting. The Reactome pathway database was used to determine metabolic pathways that are altered in obesity.

RESULTS

Among proteins, 73 and 41 proteins showed different (mostly lower) expression in subjects with obesity in the SS and IMF mitochondria, respectively (false discovery rate-adjusted P ≤ 0.05). We specifically found an increase in proteins forming the tricarboxylic acid cycle and electron transport chain (ETC) complex II, but a decrease in proteins forming protein complexes I and III of the ETC and adenosine triphosphate (ATP) synthase in subjects with obesity in the IMF, but not SS, mitochondria. Obesity was associated with differential effects on metabolic pathways linked to protein translation in the SS mitochondria and ATP formation in the IMF mitochondria.

CONCLUSIONS

Obesity alters the expression of mitochondrial proteins regulating key metabolic processes in skeletal muscle, and these effects are distinct to mitochondrial subpopulations located in different regions of the muscle fibers.

TRIAL REGISTRATION

ClinicalTrials.gov (NCT01824173).

In Silico Profiling of Clinical Phenotypes for Human Targets Using Adverse Event Data

We present a novel approach for the molecular transformation and analysis of patient clinical phenotypes. Building on the fact that drugs perturb the function of targets/genes, we integrated data from 8.2 million clinical reports detailing drug-induced side effects with the molecular world of drug-target information. Using this dataset, we extracted 1.8 million associations of clinical phenotypes to 770 human drug-targets. This collection is perhaps the largest phenotypic profiling reference of human targets to-date, and unique in that it enables rapid development of testable molecular hypotheses directly from human-specific information. We also present validation results demonstrating analytical utilities of the approach, including drug safety prediction, and the design of novel combination therapies. Challenging the long-standing notion that molecular perturbation studies cannot be performed in humans, our data allows researchers to capitalize on the vast tomes of clinical information available throughout the healthcare system.