Control of nutrient stress-induced metabolic reprogramming by PKCζ in tumorigenesis.

Tumor cells have high-energetic and anabolic needs and are known to adapt their metabolism to be able to survive and keep proliferating under conditions of nutrient stress. We show that PKCζ deficiency promotes the plasticity necessary for cancer cells to reprogram their metabolism to utilize glutamine through the serine biosynthetic pathway in the absence of glucose. PKCζ represses the expression of two key enzymes of the pathway, PHGDH and PSAT1, and phosphorylates PHGDH at key residues to inhibit its enzymatic activity. Interestingly, the loss of PKCζ in mice results in enhanced intestinal tumorigenesis and increased levels of these two metabolic enzymes, whereas patients with low levels of PKCζ have a poor prognosis. Furthermore, PKCζ and caspase-3 activities are correlated with PHGDH levels in human intestinal tumors. Taken together, this demonstrates that PKCζ is a critical metabolic tumor suppressor in mouse and human cancer.

Crosstalk of MAP3K1 and EGFR signaling mediates gene-environment interactions that block developmental tissue closure.

Aberrant regulation of signal transduction pathways can adversely derail biological processes for tissue development. One such process is the embryonic eyelid closure that is dependent on the mitogen-activated protein kinase kinase kinase 1 (MAP3K1). Map3k1 KO in mice results in defective eyelid closure and an autosomal recessive eye-open at birth phenotype. We have shown that in utero exposure to dioxin, a persistent environmental toxicant, induces the same eye defect in Map3k1+/- heterozygous but not WT pups. Here, we explore the mechanisms of the Map3k1 (gene) and dioxin (environment) interactions (GxE) underlying defective eyelid closure. We show that, acting through the aryl hydrocarbon receptor, dioxin activates epidermal growth factor receptor signaling, which in turn depresses MAP3K1-dependent Jun N-terminal kinase (JNK) activity. The dioxin-mediated JNK repression is moderate but is exacerbated by Map3k1 heterozygosity. Therefore, dioxin exposed Map3k1+/- embryonic eyelids have a marked reduction of JNK activity, accelerated differentiation and impeded polarization in the epithelial cells. Knocking out Ahr or Egfr in eyelid epithelium attenuates the open-eye defects in dioxin-treated Map3k1+/- pups, whereas knockout of Jnk1 and S1pr that encodes the sphigosin-1-phosphate (S1P) receptors upstream of the MAP3K1-JNK pathway potentiates the dioxin toxicity. Our novel findings show that the crosstalk of aryl hydrocarbon receptor, epidermal growth factor receptor, and S1P-MAP3K1-JNK pathways determines the outcome of dioxin exposure. Thus, gene mutations targeting these pathways are potential risk factors for the toxicity of environmental chemicals.

Microtubule-associated protein MAP1LC3C regulates lysosomal exocytosis and induces zinc reprogramming in renal cancer cells.

Microtubule-associated protein 1 light chain 3 gamma (MAP1LC3C or LC3C) is a member of the microtubule-associated family of proteins that are essential in the formation of autophagosomes and lysosomal degradation of cargo. LC3C has tumor-suppressing activity, and its expression is dependent on kidney cancer tumor suppressors, such as von Hippel-Lindau protein and folliculin. Recently, we demonstrated that LC3C autophagy is regulated by noncanonical upstream regulatory complexes and targets for degradation postdivision midbody rings associated with cancer cell stemness. Here, we show that loss of LC3C leads to peripheral positioning of the lysosomes and lysosomal exocytosis (LE). This process is independent of the autophagic activity of LC3C. Analysis of isogenic cells with low and high LE shows substantial transcriptomic reprogramming with altered expression of zinc (Zn)-related genes and activity of polycomb repressor complex 2, accompanied by a robust decrease in intracellular Zn. In addition, metabolomic analysis revealed alterations in amino acid steady-state levels. Cells with augmented LE show increased tumor initiation properties and form aggressive tumors in xenograft models. Immunocytochemistry identified high levels of lysosomal-associated membrane protein 1 on the plasma membrane of cancer cells in human clear cell renal cell carcinoma and reduced levels of Zn, suggesting that LE occurs in clear cell renal cell carcinoma, potentially contributing to the loss of Zn. These data indicate that the reprogramming of lysosomal localization and Zn metabolism with implication for epigenetic remodeling in a subpopulation of tumor-propagating cancer cells is an important aspect of tumor-suppressing activity of LC3C.

SigCom LINCS: data and metadata search engine for a million gene expression signatures.

Millions of transcriptome samples were generated by the Library of Integrated Network-based Cellular Signatures (LINCS) program. When these data are processed into searchable signatures along with signatures extracted from Genotype-Tissue Expression (GTEx) and Gene Expression Omnibus (GEO), connections between drugs, genes, pathways and diseases can be illuminated. SigCom LINCS is a webserver that serves over a million gene expression signatures processed, analyzed, and visualized from LINCS, GTEx, and GEO. SigCom LINCS is built with Signature Commons, a cloud-agnostic skeleton Data Commons with a focus on serving searchable signatures. SigCom LINCS provides a rapid signature similarity search for mimickers and reversers given sets of up and down genes, a gene set, a single gene, or any search term. Additionally, users of SigCom LINCS can perform a metadata search to find and analyze subsets of signatures and find information about genes and drugs. SigCom LINCS is findable, accessible, interoperable, and reusable (FAIR) with metadata linked to standard ontologies and vocabularies. In addition, all the data and signatures within SigCom LINCS are available via a well-documented API. In summary, SigCom LINCS, available at https://maayanlab.cloud/sigcom-lincs, is a rich webserver resource for accelerating drug and target discovery in systems pharmacology.

The aryl hydrocarbon receptor directs the differentiation of murine progenitor blastomeres.

Key regulatory decisions during cleavage divisions in mammalian embryogenesis determine the fate of preimplantation embryonic cells. Single-cell RNA sequencing of early-stage-2-cell, 4-cell, and 8-cell-blastomeres show that the aryl hydrocarbon receptor (AHR), traditionally considered as an environmental sensor, directs blastomere differentiation. Disruption of AHR functions in Ahr knockout embryos or in embryos from dams exposed to dioxin, the prototypic xenobiotic AHR agonist, significantly impairs blastocyst formation, causing repression and loss of transcriptional heterogeneity of OCT4 and CDX2 and incidence of nonspecific downregulation of pluripotency. Trajectory-the path of differentiation-and gene variability analyses further confirm that deregulation of OCT4 functions and changes of transcriptional heterogeneity resulting from disruption of AHR functions restrict the emergence of differentiating blastomeres in 4-cell embryos. It appears that AHR directs the differentiation of progenitor blastomeres and that disruption of preimplantation AHR functions may significantly perturb embryogenesis leading to long-lasting conditions at the heart of disease in offspring's adulthood.

piNET: a versatile web platform for downstream analysis and visualization of proteomics data.

Rapid progress in proteomics and large-scale profiling of biological systems at the protein level necessitates the continued development of efficient computational tools for the analysis and interpretation of proteomics data. Here, we present the piNET server that facilitates integrated annotation, analysis and visualization of quantitative proteomics data, with emphasis on PTM networks and integration with the LINCS library of chemical and genetic perturbation signatures in order to provide further mechanistic and functional insights. The primary input for the server consists of a set of peptides or proteins, optionally with PTM sites, and their corresponding abundance values. Several interconnected workflows can be used to generate: (i) interactive graphs and tables providing comprehensive annotation and mapping between peptides and proteins with PTM sites; (ii) high resolution and interactive visualization for enzyme-substrate networks, including kinases and their phospho-peptide targets; (iii) mapping and visualization of LINCS signature connectivity for chemical inhibitors or genetic knockdown of enzymes upstream of their target PTM sites. piNET has been built using a modular Spring-Boot JAVA platform as a fast, versatile and easy to use tool. The Apache Lucene indexing is used for fast mapping of peptides into UniProt entries for the human, mouse and other commonly used model organism proteomes. PTM-centric network analyses combine PhosphoSitePlus, iPTMnet and SIGNOR databases of validated enzyme-substrate relationships, for kinase networks augmented by DeepPhos predictions and sequence-based mapping of PhosphoSitePlus consensus motifs. Concordant LINCS signatures are mapped using iLINCS. For each workflow, a RESTful API counterpart can be used to generate the results programmatically in the json format. The server is available at http://pinet-server.org, and it is free and open to all users without login requirement.

LINCS Data Portal 2.0: next generation access point for perturbation-response signatures.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program with the goal of generating a large-scale and comprehensive catalogue of perturbation-response signatures by utilizing a diverse collection of perturbations across many model systems and assay types. The LINCS Data Portal (LDP) has been the primary access point for the compendium of LINCS data and has been widely utilized. Here, we report the first major update of LDP (http://lincsportal.ccs.miami.edu/signatures) with substantial changes in the data architecture and APIs, a completely redesigned user interface, and enhanced curated metadata annotations to support more advanced, intuitive and deeper querying, exploration and analysis capabilities. The cornerstone of this update has been the decision to reprocess all high-level LINCS datasets and make them accessible at the data point level enabling users to directly access and download any subset of signatures across the entire library independent from the originating source, project or assay. Access to the individual signatures also enables the newly implemented signature search functionality, which utilizes the iLINCS platform to identify conditions that mimic or reverse gene set queries. A newly designed query interface enables global metadata search with autosuggest across all annotations associated with perturbations, model systems, and signatures.

Predicting mechanism of action of cellular perturbations with pathway activity signatures.

MOTIVATION: Misregulation of signaling pathway activity is etiologic for many human diseases, and modulating activity of signaling pathways is often the preferred therapeutic strategy. Understanding the mechanism of action (MOA) of bioactive chemicals in terms of targeted signaling pathways is the essential first step in evaluating their therapeutic potential. Changes in signaling pathway activity are often not reflected in changes in expression of pathway genes which makes MOA inferences from transcriptional signatures (TSeses) a difficult problem. RESULTS: We developed a new computational method for implicating pathway targets of bioactive chemicals and other cellular perturbations by integrated analysis of pathway network topology, the Library of Integrated Network-based Cellular Signature TSes of genetic perturbations of pathway genes and the TS of the perturbation. Our methodology accurately predicts signaling pathways targeted by the perturbation when current pathway analysis approaches utilizing only the TS of the perturbation fail. AVAILABILITY AND IMPLEMENTATION: Open source R package paslincs is available at https://github.com/uc-bd2k/paslincs. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Monocyte and bone marrow macrophage transcriptional phenotypes in systemic juvenile idiopathic arthritis reveal TRIM8 as a mediator of IFN-γ hyper-responsiveness and risk for macrophage activation syndrome.

OBJECTIVES: Systemic juvenile idiopathic arthritis (SJIA) confers high risk for macrophage activation syndrome (MAS), a life-threatening cytokine storm driven by interferon (IFN)-γ. SJIA monocytes display IFN-γ hyper-responsiveness, but the molecular basis of this remains unclear. The objective of this study is to identify circulating monocyte and bone marrow macrophage (BMM) polarisation phenotypes in SJIA including molecular features contributing to IFN response. METHODS: Bulk RNA-seq was performed on peripheral blood monocytes (n=26 SJIA patients) and single cell (sc) RNA-seq was performed on BMM (n=1). Cultured macrophages were used to define consequences of tripartite motif containing 8 (TRIM8) knockdown on IFN-γ signalling. RESULTS: Bulk RNA-seq of SJIA monocytes revealed marked transcriptional changes in patients with elevated ferritin levels. We identified substantial overlap with multiple polarisation states but little evidence of IFN-induced signature. Interestingly, among the most highly upregulated genes was TRIM8, a positive regulator of IFN-γ signalling. In contrast to PBMC from SJIA patients without MAS, scRNA-seq of BMM from a patient with SJIA and MAS identified distinct subpopulations of BMM with altered transcriptomes, including upregulated IFN-γ response pathways. These BMM also showed significantly increased expression of TRIM8. In vitro knockdown of TRIM8 in macrophages significantly reduced IFN-γ responsiveness. CONCLUSIONS: Macrophages with an 'IFN-γ response' phenotype and TRIM8 overexpression were expanded in the bone marrow from an MAS patient. TRIM8 is also upregulated in SJIA monocytes, and augments macrophage IFN-γ response in vitro, providing both a candidate molecular mechanism and potential therapeutic target for monocyte hyper-responsiveness to IFNγ in cytokine storms including MAS.

ChIAPoP: a new tool for ChIA-PET data analysis.

Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) is a popular assay method for studying genome-wide chromatin interactions mediated by a protein of interest. The main goal of ChIA-PET data analysis is to detect interactions between DNA regions. Here, we propose a new method and the associated data analysis pipeline, ChIAPoP, to detect chromatin interactions from ChIA-PET data. We compared ChIAPoP with other popular methods, including a hypergeometric model (used in ChIA-PET tool), MICC (used in ChIA-PET2), ChiaSig and mango. The results showed that ChIA-PoP performed better than or at least as well as these top existing methods in detecting true chromatin interactions. ChIAPoP is freely available to the public at https://github.com/wh90999/ChIAPoP.

Connectivity Map Analysis of a Single-Cell RNA-Sequencing -Derived Transcriptional Signature of mTOR Signaling.

In the connectivity map (CMap) approach to drug repositioning and development, transcriptional signature of disease is constructed by differential gene expression analysis between the diseased tissue or cells and the control. The negative correlation between the transcriptional disease signature and the transcriptional signature of the drug, or a bioactive compound, is assumed to indicate its ability to "reverse" the disease process. A major limitation of traditional CMaP analysis is the use of signatures derived from bulk disease tissues. Since the key driver pathways are most likely dysregulated in only a subset of cells, the "averaged" transcriptional signatures resulting from bulk analysis lack the resolution to effectively identify effective therapeutic agents. The use of single-cell RNA-seq (scRNA-seq) transcriptomic assay facilitates construction of disease signatures that are specific to individual cell types, but methods for using scRNA-seq data in the context of CMaP analysis are lacking. Lymphangioleiomyomatosis (LAM) mutations in TSC1 or TSC2 genes result in the activation of the mTOR complex 1 (mTORC1). The mTORC1 inhibitor Sirolimus is the only FDA-approved drug to treat LAM. Novel therapies for LAM are urgently needed as the disease recurs with discontinuation of the treatment and some patients are insensitive to the drug. We developed methods for constructing disease transcriptional signatures and CMaP analysis using scRNA-seq profiling and applied them in the analysis of scRNA-seq data of lung tissue from naïve and sirolimus-treated LAM patients. New methods successfully implicated mTORC1 inhibitors, including Sirolimus, as capable of reverting the LAM transcriptional signatures. The CMaP analysis mimicking standard bulk-tissue approach failed to detect any connection between the LAM signature and mTORC1 signaling. This indicates that the precise signature derived from scRNA-seq data using our methods is the crucial difference between the success and the failure to identify effective therapeutic treatments in CMaP analysis.

The N6-Methyladenosine mRNA Methylase METTL3 Controls Cardiac Homeostasis and Hypertrophy.

BACKGROUND: N6-Methyladenosine (m6A) methylation is the most prevalent internal posttranscriptional modification on mammalian mRNA. The role of m6A mRNA methylation in the heart is not known. METHODS: To determine the role of m6A methylation in the heart, we isolated primary cardiomyocytes and performed m6A immunoprecipitation followed by RNA sequencing. We then generated genetic tools to modulate m6A levels in cardiomyocytes by manipulating the levels of the m6A RNA methylase methyltransferase-like 3 (METTL3) both in culture and in vivo. We generated cardiac-restricted gain- and loss-of-function mouse models to allow assessment of the METTL3-m6A pathway in cardiac homeostasis and function. RESULTS: We measured the level of m6A methylation on cardiomyocyte mRNA, and found a significant increase in response to hypertrophic stimulation, suggesting a potential role for m6A methylation in the development of cardiomyocyte hypertrophy. Analysis of m6A methylation showed significant enrichment in genes that regulate kinases and intracellular signaling pathways. Inhibition of METTL3 completely abrogated the ability of cardiomyocytes to undergo hypertrophy when stimulated to grow, whereas increased expression of the m6A RNA methylase METTL3 was sufficient to promote cardiomyocyte hypertrophy both in vitro and in vivo. Finally, cardiac-specific METTL3 knockout mice exhibit morphological and functional signs of heart failure with aging and stress, showing the necessity of RNA methylation for the maintenance of cardiac homeostasis. CONCLUSIONS: Our study identified METTL3-mediated methylation of mRNA on N6-adenosines as a dynamic modification that is enhanced in response to hypertrophic stimuli and is necessary for a normal hypertrophic response in cardiomyocytes. Enhanced m6A RNA methylation results in compensated cardiac hypertrophy, whereas diminished m6A drives eccentric cardiomyocyte remodeling and dysfunction, highlighting the critical importance of this novel stress-response mechanism in the heart for maintaining normal cardiac function.

Data Portal for the Library of Integrated Network-based Cellular Signatures (LINCS) program: integrated access to diverse large-scale cellular perturbation response data.

The Library of Integrated Network-based Cellular Signatures (LINCS) program is a national consortium funded by the NIH to generate a diverse and extensive reference library of cell-based perturbation-response signatures, along with novel data analytics tools to improve our understanding of human diseases at the systems level. In contrast to other large-scale data generation efforts, LINCS Data and Signature Generation Centers (DSGCs) employ a wide range of assay technologies cataloging diverse cellular responses. Integration of, and unified access to LINCS data has therefore been particularly challenging. The Big Data to Knowledge (BD2K) LINCS Data Coordination and Integration Center (DCIC) has developed data standards specifications, data processing pipelines, and a suite of end-user software tools to integrate and annotate LINCS-generated data, to make LINCS signatures searchable and usable for different types of users. Here, we describe the LINCS Data Portal (LDP) (http://lincsportal.ccs.miami.edu/), a unified web interface to access datasets generated by the LINCS DSGCs, and its underlying database, LINCS Data Registry (LDR). LINCS data served on the LDP contains extensive metadata and curated annotations. We highlight the features of the LDP user interface that is designed to enable search, browsing, exploration, download and analysis of LINCS data and related curated content.

Balancing yield, purity and practicality: a modified differential ultracentrifugation protocol for efficient isolation of small extracellular vesicles from human serum.

Ultracentrifugation remains the gold standard for isolation of small extracellular vesicles (sEV), particularly for cancer applications. The objective of this study was to determine if a widely used ultracentrifugation protocol for isolation of serum sEV could be modified to reduce the number of ultracentrifugation cycles and increase efficiency, while maintaining equal or better sample purity and yield. Serum was obtained from two healthy subjects. sEVs were isolated from 1 mL aliquots using three different ultracentrifugation protocols. Co-isolation of RNA carrier protein was assessed by performing Western blots for ApoA-I, ApoB, and Ago2. Small RNA-sequencing was performed on the sEV isolates, and differential detection of small ncRNA was compared across isolation protocols. Reduction from three- to two-ultracentrifuge cycles with no sucrose cushion resulted in a much higher sEV yield but also had the highest levels of lipoprotein and Ago2 contamination. However, the two-ultracentrifugation cycle protocol that incorporated a 30% sucrose cushion into the first cycle resulted in slightly higher sEV yields with lower levels of protein contamination compared to the lengthier three-ultracentrifugation cycle approach, therefore presenting a more efficient alternative approach for isolation of serum sEVs. It was also notable that there were some differences in sEV ncRNA cargo according to protocol, although it was less than expected given the differences in co-isolated RNA carrier proteins. Our results suggest that use of the modified serum sEV isolation protocol with two ultracentrifugation cycles and incorporating a 30% sucrose cushion offers a more efficient approach in terms of efficiency and purity.

The genetic fingerprint of susceptibility for transplant-associated thrombotic microangiopathy.

Transplant-associated thrombotic microangiopathy (TA-TMA) occurs frequently after hematopoietic stem cell transplantation (HSCT) and can lead to significant morbidity and mortality. There are no data addressing individual susceptibility to TA-TMA. We performed a hypothesis-driven analysis of 17 candidate genes known to play a role in complement activation as part of a prospective study of TMA in HSCT recipients. We examined the functional significance of gene variants by using gene expression profiling. Among 77 patients undergoing genetic testing, 34 had TMA. Sixty-five percent of patients with TMA had genetic variants in at least one gene compared with 9% of patients without TMA (P < .0001). Gene variants were increased in patients of all races with TMA, but nonwhites had more variants than whites (2.5 [range, 0-7] vs 0 [range, 0-2]; P < .0001). Variants in ≥3 genes were identified only in nonwhites with TMA and were associated with high mortality (71%). RNA sequencing analysis of pretransplantation samples showed upregulation of multiple complement pathways in patients with TMA who had gene variants, including variants predicted as possibly benign by computer algorithm, compared with those without TMA and without gene variants. Our data reveal important differences in genetic susceptibility to HSCT-associated TMA based on recipient genotype. These data will allow prospective risk assessment and intervention to prevent TMA in highly susceptible transplant recipients. Our findings may explain, at least in part, racial disparities previously reported in transplant recipients and may guide treatment strategies to improve outcomes.

CD44 variant isoform 9 emerges in response to injury and contributes to the regeneration of the gastric epithelium.

The CD44 gene encodes several protein isoforms due to alternative splicing and post translational modifications. Given that CD44 variant isoform 9 (CD44v9) is expressed within Spasmolytic Polypeptide/TFF2-Expressing Metaplasia (SPEM) glands during repair, CD44v9 may be play a funcitonal role during the process of regeneration of the gastric epithelium. Here we hypothesize that CD44v9 marks a regenerative cell lineage responsive to infiltrating macrophages during regeneration of the gastric epithelium. Ulcers were induced in CD44-deficient (CD44KO) and C57BL/6 (BL6) mice by a localized application of acetic acid to the serosal surface of the stomach. Gastric organoids expressing CD44v9 were derived from mouse stomachs and transplanted at the ulcer site of CD44KO mice. Ulcers, CD44v9 expression, proliferation and histology were measured 1, 3, 5 and 7-days post-injury. Human-derived gastric organoids were generated from stomach tissue collected from elderly (>55 years) or young (14-20 years) patients. Organoids were transplanted into the stomachs of NOD scid gamma (NSG) mice at the site of injury. Gastric injury was induced in NRG-SGM3 (NRGS) mice harboring human-derived immune cells (hnNRGS) and the immune profile anlayzed by CyTOF. CD44v9 expression emerged within regenerating glands the ulcer margin in response to injury. While ulcers in BL6 mice healed within 7-days post-injury, CD44KO mice exhibited loss of repair and epithelial regeneration. Ulcer healing was promoted in CD44KO mice by transplanted CD55v9-expressing gastric organoids. NSG mice exhibited loss of CD44v9 expression and gastric repair. Transplantation of human-derived gastric organoids from young, but not aged stomachs promoted repair in NSG mouse stomachs in response to injury. Finally, compared to NRGS mice, huNRGS animals exhibited reduced ulcer sizes, an infiltration of human CD162+ macrophages and an emergence of CD44v9 expression in SPEM. Thus, during repair of the gastic epithelium CD44v9 emerges within a regenerative cell lineage that coincides with macrophage inflitration within the injured mucosa. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Ultra-Deep Genomic Sequencing of HCV NS5A Resistance-Associated Substitutions in HCV/HIV Coinfected Patients.

BACKGROUND AND AIMS: The prevalence of naturally occurring HCV-NS5A resistance-associated substitutions (RAS) to DAA drugs might affect the response to treatment in HCV/HIV coinfected subjects. There are limited data on the frequency of HCV-NS5A naturally occurring drug-RAS at baseline in HCV/HIV coinfected patients when ultra-deep sequencing methodologies are applied. METHODS: HCV-NS5A-RAS were evaluated among 25 subjects in each group. Patients were matched by age, gender, and hepatic fibrosis stage category to control for selection bias. RESULTS: Within subtype 1a, RAS were observed in 28% of HCV monoinfected and 48% of HCV/HIV coinfected subjects. More patients in the HCV/HIV coinfected group had clinically relevant mutations to DAA directed at NS5A. CONCLUSION: While the clinical significance of this observation may be limited in highly drug adherent populations, some HCV/HIV coinfected persons may be at greater risk of viral resistance if suboptimal dosing occurs.

Transcriptional profiles of JIA patient blood with subsequent poor response to methotrexate.

Objective: The mechanisms that determine the efficacy or inefficacy of MTX in JIA are ill-defined. The objective of this study was to identify a gene expression transcriptional signature associated with poor response to MTX in patients with JIA. Methods: RNA sequencing was used to measure gene expression in peripheral blood mononuclear cells collected from 47 patients with JIA prior to MTX treatment and 14 age-matched controls. Differentially expressed baseline genes between responders and non-responders were evaluated. Biological differences between all JIA patients and controls were explored by constructing a signature of differentially expressed genes. Unsupervised clustering and pathway analysis was performed. Results: A signature of 99 differentially expressed genes (Bonferroni-corrected P < 0.05) capturing the biological differences between all JIA patients and controls was identified. Unsupervised clustering of samples based on this list of 99 genes produced subgroups enriched for MTX response status. Comparing this gene signature with reference signatures from sorted cell populations revealed high concordance between the expression signatures of monocytes and of MTX non-responders. CXCL8 (IL-8) was the most significantly differentially expressed gene transcript comparing all JIA patients with controls (Bonferroni-corrected P = 4.12 × 10-10). Conclusion: Variability in clinical response to MTX in JIA patients is associated with differences in gene transcripts modulated in monocytes. These gene expression profiles may provide a basis for biomarkers predictive of treatment response.

Genetic susceptibility to toxicologic lung responses among inbred mouse strains following exposure to carbon nanotubes and profiling of underlying gene networks.

The risk of human exposure to fiber nanoparticles has risen in recent years due to increases in the manufacture and utilization of carbon nanotubes (CNTs). CNTs are present as airborne particulates in occupational settings and their hazard potential has been demonstrated in experimental lung exposure studies using inbred mouse strains. However, it is not known whether different inbred strains differ in lung responses to CNTs by virtue of their genetics. In this work, common inbred strains (BALB/c, C57Bl/6, DBA/2, and C3H/He) were exposed to CNTs via oropharyngeal aspiration and lung histology and bronchoalveolar lavage (BAL) samples were evaluated over 28days with the objective of evaluating sensitivity/resistance among strains. C57Bl/6 mice developed significantly more extensive type II pneumocyte (T2P) hyperplasia and alveolar infiltrate compared to DBA/2 mice, which were resistant. Surprisingly, DBA/2 but not C57Bl/6 mice were extremely sensitive to increases in leukocytes recovered in BAL fluid. Underlying global gene expression patterns in the two strains were compared using mRNA sequencing to investigate regulatory networks associated with the different effects. The impact of exposure on gene networks regulating various aspects of immune response and cell survival was limited in DBA/2 mice compared to C57Bl/6. Investigation of B6D2F1 (C57Bl/6×DBA/2 hybrid) mice demonstrated inheritance of sensitivity to CNT exposures in regard to toxicologic lung pathology and BAL leukocyte accumulations. These findings demonstrate a genetic basis of susceptibility to CNT particle exposures and both inform the use of inbred mouse models and suggest the likelihood of differences in genetic susceptibility among humans.

GRcalculator: an online tool for calculating and mining dose-response data.

BACKGROUND: Quantifying the response of cell lines to drugs or other perturbagens is the cornerstone of pre-clinical drug development and pharmacogenomics as well as a means to study factors that contribute to sensitivity and resistance. In dividing cells, traditional metrics derived from dose-response curves such as IC 50 , AUC, and E max , are confounded by the number of cell divisions taking place during the assay, which varies widely for biological and experimental reasons. Hafner et al. (Nat Meth 13:521-627, 2016) recently proposed an alternative way to quantify drug response, normalized growth rate (GR) inhibition, that is robust to such confounders. Adoption of the GR method is expected to improve the reproducibility of dose-response assays and the reliability of pharmacogenomic associations (Hafner et al. 500-502, 2017). RESULTS: We describe here an interactive website ( www.grcalculator.org ) for calculation, analysis, and visualization of dose-response data using the GR approach and for comparison of GR and traditional metrics. Data can be user-supplied or derived from published datasets. The web tools are implemented in the form of three integrated Shiny applications (grcalculator, grbrowser, and grtutorial) deployed through a Shiny server. Intuitive graphical user interfaces (GUIs) allow for interactive analysis and visualization of data. The Shiny applications make use of two R packages (shinyLi and GRmetrics) specifically developed for this purpose. The GRmetrics R package is also available via Bioconductor and can be used for offline data analysis and visualization. Source code for the Shiny applications and associated packages (shinyLi and GRmetrics) can be accessed at www.github.com/uc-bd2k/grcalculator and www.github.com/datarail/gr_metrics . CONCLUSIONS: GRcalculator is a powerful, user-friendly, and free tool to facilitate analysis of dose-response data. It generates publication-ready figures and provides a unified platform for investigators to analyze dose-response data across diverse cell types and perturbagens (including drugs, biological ligands, RNAi, etc.). GRcalculator also provides access to data collected by the NIH LINCS Program ( http://www.lincsproject.org /) and other public domain datasets. The GRmetrics Bioconductor package provides computationally trained users with a platform for offline analysis of dose-response data and facilitates inclusion of GR metrics calculations within existing R analysis pipelines. These tools are therefore well suited to users in academia as well as industry.