A Hidradenitis Suppurativa molecular disease signature derived from patient samples by high-throughput RNA sequencing and re-analysis of previously reported transcriptomic data sets.

Hidradenitis suppurativa (HS) is a common, debilitating inflammatory skin disease linked to immune dysregulation and abnormalities in follicular structure and function. Several studies have characterized the transcriptomic profile of affected and unaffected skin in small populations. In this study of 20 patients, RNA from lesional and matching non-lesional skin biopsies in 20 subjects were used to identify an expression-based HS disease signature. This was followed by differential expression and pathway enrichment analyses, as well as jointly reanalyzing our findings with previously published transcriptomic profiles. We establish an RNA-Seq based HS expression disease signature that is mostly consistent with previous reports. Bulk-RNA profiles from 104 subjects in 7 previously reported data sets identified a disease signature of 118 differentially regulated genes compared to three control data sets from non-lesional skin. We confirmed previously reported expression profiles and further characterized dysregulation in complement activation and host response to bacteria in disease pathogenesis. Changes in the transcriptome of lesional skin in this cohort of HS patients is consistent with smaller previously reported populations. The findings further support the significance of immune dysregulation, in particular with regard to bacterial response mechanisms. Joint analysis of this and previously reported cohorts indicate a remarkably consistent expression profile.

The essentiality of drug targets: an analysis of current literature and genomic databases.

Essential genes encode proteins thought to be crucial for the survival of an organism. However, the role of essential genes in human disease and their suitability as drug targets is less clear. Here, we use a recent catalog of nearly 9000 known essential and nonessential genes to evaluate their involvement as therapeutic targets in human diseases. We find that essential genes are more likely than nonessential genes to play a part in specific therapeutic areas such as cardiovascular diseases and neoplasms. We also find significant differences between essential and nonessential genes among protein classes relevant to drug discovery. Taken together, our analyses suggest that the essentiality status of a potential new target is an important consideration in drug discovery.

Integrating Biological Networks for Drug Target Prediction and Prioritization.

Computational prediction of the clinical success or failure of a potential drug target for therapeutic use is a challenging problem. Novel network propagation algorithms that integrate heterogeneous biological networks are proving useful for drug target identification and prioritization. These approaches typically utilize a network describing relationships between targets, a method to disseminate the relevant information through the network, and a method to elucidate new associations between targets and diseases. Here, we utilize one such network propagation-based approach, DTINet, which starts with diffusion component analysis of networks of both potential drug targets and diseases. Then an inductive matrix completion algorithm is applied to identify novel disease targets based on their network topological similarities with known disease targets with successfully launched drugs. DTINet performed well as assessed with area under the precision-recall curve (AUPR = 0.88 ± 0.007) and area under the receiver operating characteristic curve (AUROC = 0.86 ± 0.008). These metrics improved when we combined data from multiple networks in the target space but reduced significantly when we used a more conservative method to define negative controls (AUPR = 0.56 ± 0.007, AUROC = 0.57 ± 0.007). We are optimistic that integration of more relevant and cleaner datasets and networks, careful calibration of model parameters, as well as algorithmic improvements will improve prediction accuracy. However, we also recognize that predicting drug targets that are likely to be successful is an extremely challenging problem due to its complex nature and sparsity of known disease targets.

Correction to: HIF prolyl hydroxylase inhibition protects skeletal muscle from eccentric contraction induced injury.

Following publication of the original article [1], the authors flagged that there is a discrepancy with the Availability of data and materials statement on page 12 of the article.

HIF prolyl hydroxylase inhibition protects skeletal muscle from eccentric contraction-induced injury.

BACKGROUND: In muscular dystrophy and old age, skeletal muscle repair is compromised leading to fibrosis and fatty tissue accumulation. Therefore, therapies that protect skeletal muscle or enhance repair would be valuable medical treatments. Hypoxia-inducible factors (HIFs) regulate gene transcription under conditions of low oxygen, and HIF target genes EPO and VEGF have been associated with muscle protection and repair. We tested the importance of HIF activation following skeletal muscle injury, in both a murine model and human volunteers, using prolyl hydroxylase inhibitors that stabilize and activate HIF. METHODS: Using a mouse eccentric limb injury model, we characterized the protective effects of prolyl hydroxylase inhibitor, GSK1120360A. We then extended these studies to examine the impact of EPO modulation and infiltrating immune cell populations on muscle protection. Finally, we extended this study with an experimental medicine approach using eccentric arm exercise in untrained volunteers to measure the muscle-protective effects of a clinical prolyl hydroxylase inhibitor, daprodustat. RESULTS: GSK1120360A dramatically prevented functional deficits and histological damage, while accelerating recovery after eccentric limb injury in mice. Surprisingly, this effect was independent of EPO, but required myeloid HIF1α-mediated iNOS activity. Treatment of healthy human volunteers with high-dose daprodustat reduced accumulation of circulating damage markers following eccentric arm exercise, although we did not observe any diminution of functional deficits with compound treatment. CONCLUSION: The results of these experiments highlight a novel skeletal muscle protective effect of prolyl hydroxylase inhibition via HIF-mediated expression of iNOS in macrophages. Partial recapitulation of these findings in healthy volunteers suggests elements of consistent pharmacology compared to responses in mice although there are clear differences between these two systems.

Uncovering new disease indications for G-protein coupled receptors and their endogenous ligands.

BACKGROUND: The Open Targets Platform integrates different data sources in order to facilitate identification of potential therapeutic drug targets to treat human diseases. It currently provides evidence for nearly 2.6 million potential target-disease pairs. G-protein coupled receptors are a drug target class of high interest because of the number of successful drugs being developed against them over many years. Here we describe a systematic approach utilizing the Open Targets Platform data to uncover and prioritize potential new disease indications for the G-protein coupled receptors and their ligands. RESULTS: Utilizing the data available in the Open Targets platform, potential G-protein coupled receptor and endogenous ligand disease association pairs were systematically identified. Intriguing examples such as GPR35 for inflammatory bowel disease and CXCR4 for viral infection are used as illustrations of how a systematic approach can aid in the prioritization of interesting drug discovery hypotheses. Combining evidences for G-protein coupled receptors and their corresponding endogenous peptidergic ligands increases confidence and provides supportive evidence for potential new target-disease hypotheses. Comparing such hypotheses to the global pharma drug discovery pipeline to validate the approach showed that more than 93% of G-protein coupled receptor-disease pairs with a high overall Open Targets score involved receptors with an existing drug discovery program. CONCLUSIONS: The Open Targets gene-disease score can be used to prioritize potential G-protein coupled receptors-indication hypotheses. In addition, availability of multiple different evidence types markedly increases confidence as does combining evidence from known receptor-ligand pairs. Comparing the top-ranked hypotheses to the current global pharma pipeline serves validation of our approach and identifies and prioritizes new therapeutic opportunities.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes.

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

The epigenetic signature of systemic insulin resistance in obese women.

AIMS/HYPOTHESIS: Insulin resistance (IR) links obesity to type 2 diabetes. The aim of this study was to explore whether white adipose tissue (WAT) epigenetic dysregulation is associated with systemic IR by genome-wide CG dinucleotide (CpG) methylation and gene expression profiling in WAT from insulin-resistant and insulin-sensitive women. A secondary aim was to determine whether the DNA methylation signature in peripheral blood mononuclear cells (PBMCs) reflects WAT methylation and, if so, can be used as a marker for systemic IR. METHODS: From 220 obese women, we selected a total of 80 individuals from either of the extreme ends of the distribution curve of HOMA-IR, an indirect measure of systemic insulin sensitivity. Genome-wide transcriptome and DNA CpG methylation profiling by array was performed on subcutaneous (SAT) and visceral (omental) adipose tissue (VAT). CpG methylation in PBMCs was assayed in the same cohort. RESULTS: There were 647 differentially expressed genes (false discovery rate [FDR] 10%) in SAT, all of which displayed directionally consistent associations in VAT. This suggests that IR is associated with dysregulated expression of a common set of genes in SAT and VAT. The average degree of DNA methylation did not differ between the insulin-resistant and insulin-sensitive group in any of the analysed tissues/cells. There were 223 IR-associated genes in SAT containing a total of 336 nominally significant differentially methylated sites (DMS). The 223 IR-associated genes were over-represented in pathways related to integrin cell surface interactions and insulin signalling and included COL5A1, GAB1, IRS2, PFKFB3 and PTPRJ. In VAT there were a total of 51 differentially expressed genes (FDR 10%); 18 IR-associated genes contained a total of 29 DMS. CONCLUSIONS/INTERPRETATION: In individuals discordant for insulin sensitivity, the average DNA CpG methylation in SAT and VAT is similar, although specific genes, particularly in SAT, display significantly altered expression and DMS in IR, possibly indicating that epigenetic regulation of these genes influences metabolism.

Significant obesity-associated gene expression changes occur in the stomach but not intestines in obese mice.

The gastrointestinal (GI) tract can have significant impact on the regulation of the whole-body metabolism and may contribute to the development of obesity and diabetes. To systemically elucidate the role of the GI tract in obesity, we performed a transcriptomic analysis in different parts of the GI tract of two obese mouse models: ob/ob and high-fat diet (HFD) fed mice. Compared to their lean controls, significant changes in the gene expression were observed in both obese mouse groups in the stomach (ob/ob: 959; HFD: 542). In addition, these changes were quantitatively much higher than in the intestine. Despite the difference in genetic background, the two mouse models shared 296 similar gene expression changes in the stomach. Among those genes, some had known associations to obesity, diabetes, and insulin resistance. In addition, the gene expression profiles strongly suggested an increased gastric acid secretion in both obese mouse models, probably through an activation of the gastrin pathway. In conclusion, our data reveal a previously unknown dominant connection between the stomach and obesity in murine models extensively used in research.

PHARMACOLOGICAL SIRT1 ACTIVATION IMPROVES MORTALITY AND MARKEDLY ALTERS TRANSCRIPTIONAL PROFILES THAT ACCOMPANY EXPERIMENTAL SEPSIS.

The sirtuin family consists of seven NAD+-dependent enzymes affecting a broad array of regulatory protein networks by primarily catalyzing the deacetylation of key lysine residues in regulatory proteins. The enzymatic activity of SIRT1 can be enhanced by small molecule activators known as SIRT1 activator compounds (STACs). We tested the therapeutic potential of the STAC SRT3025 in two preclinical models of severe infection, the murine cecal ligation and puncture (CLP) model to induce peritonitis and intratracheal installation of Streptococcus pneumoniae to induce severe bacterial pneumonia. SRT3025 provided significant survival benefits over vehicle control in both the peritonitis and pneumococcal pneumonia models when administered with appropriate antimicrobial agents. The survival benefit of SRT3025 in the CLP model was absent in SIRT1 knockout showing the SIRT1 dependency of SRT3025's effects. SRT3025 administration promoted bacterial clearance and significantly reduced inflammatory cytokines from the lungs of animals challenged with S. pneumoniae. SRT3025 treatment was also accompanied by striking changes in the transcription profiles in multiple inflammatory and metabolic pathways in liver, spleen, small bowel, and lung tissue. Remarkably, these organ-specific changes in the transcriptome analyses were similar following CLP or pneumococcal challenge despite different sets of pathogens at disparate sites of infection. Pharmacologic activation of SIRT1 modulates the innate host response and could represent a novel treatment strategy for severe infection.

Network-Based Identification and Prioritization of Key Regulators of Coronary Artery Disease Loci.

OBJECTIVE: Recent genome-wide association studies of coronary artery disease (CAD) have revealed 58 genome-wide significant and 148 suggestive genetic loci. However, the molecular mechanisms through which they contribute to CAD and the clinical implications of these findings remain largely unknown. We aim to retrieve gene subnetworks of the 206 CAD loci and identify and prioritize candidate regulators to better understand the biological mechanisms underlying the genetic associations. APPROACH AND RESULTS: We devised a new integrative genomics approach that incorporated (1) candidate genes from the top CAD loci, (2) the complete genetic association results from the 1000 genomes-based CAD genome-wide association studies from the Coronary Artery Disease Genome Wide Replication and Meta-Analysis Plus the Coronary Artery Disease consortium, (3) tissue-specific gene regulatory networks that depict the potential relationship and interactions between genes, and (4) tissue-specific gene expression patterns between CAD patients and controls. The networks and top-ranked regulators according to these data-driven criteria were further queried against literature, experimental evidence, and drug information to evaluate their disease relevance and potential as drug targets. Our analysis uncovered several potential novel regulators of CAD such as LUM and STAT3, which possess properties suitable as drug targets. We also revealed molecular relations and potential mechanisms through which the top CAD loci operate. Furthermore, we found that multiple CAD-relevant biological processes such as extracellular matrix, inflammatory and immune pathways, complement and coagulation cascades, and lipid metabolism interact in the CAD networks. CONCLUSIONS: Our data-driven integrative genomics framework unraveled tissue-specific relations among the candidate genes of the CAD genome-wide association studies loci and prioritized novel network regulatory genes orchestrating biological processes relevant to CAD.

Assessing Colonic Exposure, Safety, and Clinical Activity of SRT2104, a Novel Oral SIRT1 Activator, in Patients with Mild to Moderate Ulcerative Colitis.

BACKGROUND: Sirtuins are a class of proteins with important physiologic roles in metabolism and inflammation. Sirtuin (silent mating type information regulation 2 homolog) 1, or SIRT1, activation is an unexplored therapeutic approach for the treatment of ulcerative colitis (UC). METHODS: Patients with mild to moderately active UC were blindly randomized to 50 mg or 500 mg daily of SRT2104, a selective activator of SIRT1, for 8 weeks. Colonic exposure and safety were assessed, as well as blinded endoscopic scoring and disease activity by Mayo score, Simple Clinical Colitis Activity Index and fecal calprotectin. RESULTS: Across both SRT2104 groups, only 3 of 26 evaluable subjects achieved remission on blinded endoscopic assessment. Clinical remission (Mayo score ≤2, no subscore >1) was achieved in 4 patients (2 of 13 evaluable patients in each dose group). Fecal calprotectin levels declined with treatment in both groups, but after 56 days of treatment subjects were still found to have levels approximately 4-fold elevated above normal. One subject experienced an SAE requiring study withdrawal and another was withdrawn for a severe UC flare; 19 subjects (61%) across both treatment groups experienced at least 1 treatment emergent adverse event. Average drug exposure increased in a dose-dependent manner for escalating doses of SRT2104, and colonic exposure was 140 to 160 times higher than plasma exposures. CONCLUSIONS: SRT2104 did not demonstrate significant clinical activity in mild to moderately active UC. This suggests that further evaluation of SRT2104 as a therapeutic strategy for the treatment of UC is not warranted.

Effect of Roux-en-Y gastric bypass surgery on bile acid metabolism in normal and obese diabetic rats.

In addition to classic functions of facilitating hepatobiliary secretion and intestinal absorption of lipophilic nutrients, bile acids (BA) are also endocrine factors and regulate glucose and lipid metabolism. Recent data indicate that antiobesity bariatric procedures e.g. Roux-en-Y gastric bypass surgery (RYGB), which also remit diabetes, increase plasma BAs in humans, leading to the hypothesis that BAs may play a role in diabetes resolution following surgery. To investigate the effect of RYGB on BA physiology and its relationship with glucose homeostasis, we undertook RYGB and SHAM surgery in Zucker diabetic fatty (ZDF) and normoglycemic Sprague Dawley (SD) rats and measured plasma and fecal BA levels, as well as plasma glucose, insulin, Glucagon like peptide 1 (GLP-1) and Peptide YY (PYY), 2 days before and 3, 7, 14 and 28 days after surgery. RYGB decreased body weight and increased plasma GLP-1 in both SD and ZDF rats while decreasing plasma insulin and glucose in ZDF rats starting from the first week. Compared to SHAM groups, both SD-RYGB and ZDF-RYGB groups started to have increases in plasma total BAs in the second week, which might not contribute to early post-surgery metabolic changes. While there was no significant difference in fecal BA excretion between SD-RYGB and SD-SHAM groups, the ZDF-RYGB group had a transient 4.2-fold increase (P<0.001) in 24-hour fecal BA excretion on post-operative day 3 compared to ZDF-SHAM, which paralleled a significant increase in plasma PYY. Ratios of plasma and fecal cholic acid/chenodeoxycholic acid derived BAs were decreased in RYGB groups. In addition, tissue mRNA expression analysis suggested early intestinal BA reabsorption and potentially reduced hepatic cholic acid production in RYGB groups. In summary, we present novel data on RYGB-mediated changes in BA metabolism to further understand the role of BAs in RYGB-induced metabolic effects in humans.

Histone-fold domain protein NF-Y promotes chromatin accessibility for cell type-specific master transcription factors.

Cell type-specific master transcription factors (TFs) play vital roles in defining cell identity and function. However, the roles ubiquitous factors play in the specification of cell identity remain underappreciated. Here we show that the ubiquitous CCAAT-binding NF-Y complex is required for the maintenance of embryonic stem cell (ESC) identity and is an essential component of the core pluripotency network. Genome-wide studies in ESCs and neurons reveal that NF-Y regulates not only genes with housekeeping functions through cell type-invariant promoter-proximal binding, but also genes required for cell identity by binding to cell type-specific enhancers with master TFs. Mechanistically, NF-Y's distinct DNA-binding mode promotes master/pioneer TF binding at enhancers by facilitating a permissive chromatin conformation. Our studies unearth a conceptually unique function for histone-fold domain (HFD) protein NF-Y in promoting chromatin accessibility and suggest that other HFD proteins with analogous structural and DNA-binding properties may function in similar ways.

Mining emerging biomedical literature for understanding disease associations in drug discovery.

Systematically evaluating the exponentially growing body of scientific literature has become a critical task that every drug discovery organization must engage in in order to understand emerging trends for scientific investment and strategy development. Developing trends analysis uses the number of publications within a 3-year window to determine concepts derived from well-established disease and gene ontologies to aid in recognizing and predicting emerging areas of scientific discoveries relevant to that space. In this chapter, we describe such a method and use obesity and psoriasis as use-case examples by analyzing the frequency of disease-related MeSH terms in PubMed abstracts over time. We share how our system can be used to predict emerging trends at a relatively early stage and we analyze the literature-identified genes for genetic associations, druggability, and biological pathways to explore any potential biological connections between the two diseases that could be utilized for drug discovery.

Integrative framework for identification of key cell identity genes uncovers determinants of ES cell identity and homeostasis.

Identification of genes associated with specific biological phenotypes is a fundamental step toward understanding the molecular basis underlying development and pathogenesis. Although RNAi-based high-throughput screens are routinely used for this task, false discovery and sensitivity remain a challenge. Here we describe a computational framework for systematic integration of published gene expression data to identify genes defining a phenotype of interest. We applied our approach to rank-order all genes based on their likelihood of determining ES cell (ESC) identity. RNAi-mediated loss-of-function experiments on top-ranked genes unearthed many novel determinants of ESC identity, thus validating the derived gene ranks to serve as a rich and valuable resource for those working to uncover novel ESC regulators. Underscoring the value of our gene ranks, functional studies of our top-hit Nucleolin (Ncl), abundant in stem and cancer cells, revealed Ncl's essential role in the maintenance of ESC homeostasis by shielding against differentiation-inducing redox imbalance-induced oxidative stress. Notably, we report a conceptually novel mechanism involving a Nucleolin-dependent Nanog-p53 bistable switch regulating the homeostatic balance between self-renewal and differentiation in ESCs. Our findings connect the dots on a previously unknown regulatory circuitry involving genes associated with traits in both ESCs and cancer and might have profound implications for understanding cell fate decisions in cancer stem cells. The proposed computational framework, by helping to prioritize and preselect candidate genes for tests using complex and expensive genetic screens, provides a powerful yet inexpensive means for identification of key cell identity genes.

Human and mouse skeletal muscle stem cells: convergent and divergent mechanisms of myogenesis.

Satellite cells are the chief contributor to skeletal muscle growth and regeneration. The study of mouse satellite cells has accelerated in recent years due to technical advancements in the isolation of these cells. The study of human satellite cells has lagged and thus little is known about how the biology of mouse and human satellite cells compare. We developed a flow cytometry-based method to prospectively isolate human skeletal muscle progenitors from the satellite cell pool using positive and negative selection markers. Results show that this pool is enriched in PAX7 expressing cells that possess robust myogenic potential including the ability to give rise to de novo muscle in vivo. We compared mouse and human satellite cells in culture and identify differences in the elaboration of the myogenic genetic program and in the sensitivity of the cells to cytokine stimulation. These results indicate that not all mechanisms regulating mouse satellite cell activation are conserved in human satellite cells and that such differences may impact the clinical translation of therapeutics validated in mouse models. Thus, the findings of this study are relevant to developing therapies to combat muscle disease.

Fip1 regulates mRNA alternative polyadenylation to promote stem cell self-renewal.

mRNA alternative polyadenylation (APA) plays a critical role in post-transcriptional gene control and is highly regulated during development and disease. However, the regulatory mechanisms and functional consequences of APA remain poorly understood. Here, we show that an mRNA 3' processing factor, Fip1, is essential for embryonic stem cell (ESC) self-renewal and somatic cell reprogramming. Fip1 promotes stem cell maintenance, in part, by activating the ESC-specific APA profiles to ensure the optimal expression of a specific set of genes, including critical self-renewal factors. Fip1 expression and the Fip1-dependent APA program change during ESC differentiation and are restored to an ESC-like state during somatic reprogramming. Mechanistically, we provide evidence that the specificity of Fip1-mediated APA regulation depends on multiple factors, including Fip1-RNA interactions and the distance between APA sites. Together, our data highlight the role for post-transcriptional control in stem cell self-renewal, provide mechanistic insight on APA regulation in development, and establish an important function for APA in cell fate specification.