Characterization and predictive functional profiles on metagenomic 16S rRNA data of liver transplant recipients: A longitudinal study.

Long-term survival after Liver Transplantation (LT) is often compromised by infectious and metabolic complications. We aimed to delineate alterations in intestinal microbiome (IM) over time that could contribute to medical complications compromising long-term survival following LT. Fecal samples from LT recipients were collected at 3 months (3 M) and 6 months (6 M) post-LT. The bacterial DNA was extracted using E.Z.N.A. Stool DNA Kit and 16S rRNA gene sequencing at V4 hypervariable region was performed. DADA2 and Phyloseq was implemented to analyze the taxonomic composition. Differentially abundant taxa were identified by metagenomeSeq and LEfSe. Piphillin, an Inferred functional metagenomic analysis tool was used to study the bacterial functional content. For comparison, healthy samples were extracted from NCBI and analyzed similarly. The taxonomic & functional profiles of LT recipients were validated with metagenomic sequencing data from animals exposed to immunosuppressants using Venny. Our findings provide a new perspective on longitudinal increase in specific IM communities post-LT along with an increase in bacterial genes associated with metabolic and infectious disease.

Estrogen Receptor 1 Inhibition of Wnt/ß-Catenin Signaling Contributes to Sex Differences in Hepatocarcinogenesis.

BACKGROUND: Hepatocellular Carcinoma (HCC) is a sexually dimorphic cancer, with female sex being independently protective against HCC incidence and progression. The aim of our study was to understand the mechanism of estrogen receptor signaling in driving sex differences in hepatocarcinogenesis. METHODS: We integrated 1,268 HCC patient sample profiles from publicly available gene expression data to identify the most differentially expressed genes (DEGs). We mapped DEGs into a physical protein interaction network and performed network topology analysis to identify the most important proteins. Experimental validation was performed in vitro on HCC cell lines, in and in vivo, using HCC mouse model. RESULTS: We showed that the most central protein, ESR1, is HCC prognostic, as increased ESR1 expression was protective for overall survival, with HR=0.45 (95%CI 0.32-0.64, p=4.4E-06), and was more pronounced in women. Transfection of HCC cell lines with ESR1 and exposure to estradiol affected expression of genes involved in the Wnt/ß-catenin signaling pathway. ER-α (protein product of ESR1) agonist treatment in a mouse model of HCC resulted in significantly longer survival and decreased tumor burden (p<0.0001), with inhibition of Wnt/ß-Catenin signaling. In vitro experiments confirmed colocalization of ß-catenin with ER-α, leading to inhibition of ß-catenin-mediated transcription of target genes c-Myc and Cyclin D1. CONCLUSION: Combined, the centrality of ESR1 and its inhibition of the Wnt/ß-catenin signaling axis provide a biological rationale for protection against HCC incidence and progression in women.

Combined proteomic/transcriptomic signature of recurrence post-liver transplantation for hepatocellular carcinoma beyond Milan.

BACKGROUND AND AIMS: Liver transplantation (LT) can be offered to patients with Hepatocellular carcinoma (HCC) beyond Milan criteria. However, there are currently limited molecular markers on HCC explant histology to predict recurrence, which arises in up to 20% of LT recipients. The goal of our study was to derive a combined proteomic/transcriptomic signature on HCC explant predictive of recurrence post-transplant using unbiased, high-throughput approaches. METHODS: Patients who received a LT for HCC beyond Milan criteria in the context of hepatitis B cirrhosis were identified. Tumor explants from patients with post-transplant HCC recurrence (N = 7) versus those without recurrence (N = 4) were analyzed by mass spectrometry and gene expression array. Univariate analysis was used to generate a combined proteomic/transcriptomic signature linked to recurrence. Significantly predictive genes and proteins were verified and internally validated by immunoblotting and immunohistochemistry. RESULTS: Seventy-nine proteins and 636 genes were significantly differentially expressed in HCC tumors with subsequent recurrence (p < 0.05). Univariate survival analysis identified Aldehyde Dehydrogenase 1 Family Member A1 (ALDH1A1) gene (HR = 0.084, 95%CI 0.01-0.68, p = 0.0152), ALDH1A1 protein (HR = 0.039, 95%CI 0.16-0.91, p = 0.03), Galectin 3 Binding Protein (LGALS3BP) gene (HR = 7.14, 95%CI 1.20-432.96, p = 0.03), LGALS3BP protein (HR = 2.6, 95%CI 1.1-6.1, p = 0.036), Galectin 3 (LGALS3) gene (HR = 2.89, 95%CI 1.01-8.3, p = 0.049) and LGALS3 protein (HR = 2.6, 95%CI 1.2-5.5, p = 0.015) as key dysregulated analytes in recurrent HCC. In concordance with our proteome findings, HCC recurrence was linked to decreased ALDH1A1 and increased LGALS3 protein expression by Western Blot. LGALS3BP protein expression was validated in 29 independent HCC samples. CONCLUSIONS: Significantly increased LGALS3 and LGALS3BP gene and protein expression on explant were associated with post-transplant recurrence, whereas increased ALDH1A1 was associated with absence of recurrence in patients transplanted for HCC beyond Milan criteria. This combined proteomic/transcriptomic signature could help in predicting HCC recurrence risk and guide post-transplant surveillance.

Integrative analysis of layers of data in hepatocellular carcinoma reveals pathway dependencies.

BACKGROUND: The broader use of high-throughput technologies has led to improved molecular characterization of hepatocellular carcinoma (HCC). AIM: To comprehensively analyze and characterize all publicly available genomic, gene expression, methylation, miRNA and proteomic data in HCC, covering 85 studies and 3355 patient sample profiles, to identify the key dysregulated genes and pathways they affect. METHODS: We collected and curated all well-annotated and publicly available high-throughput datasets from PubMed and Gene Expression Omnibus derived from human HCC tissue. Comprehensive pathway enrichment analysis was performed using pathDIP for each data type (genomic, gene expression, methylation, miRNA and proteomic), and the overlap of pathways was assessed to elucidate pathway dependencies in HCC. RESULTS: We identified a total of 8733 abstracts retrieved by the search on PubMed on HCC for the different layers of data on human HCC samples, published until December 2016. The common key dysregulated pathways in HCC tissue across different layers of data included epidermal growth factor (EGFR) and ß1-integrin pathways. Genes along these pathways were significantly and consistently dysregulated across the different types of high-throughput data and had prognostic value with respect to overall survival. Using CTD database, estradiol would best modulate and revert these genes appropriately. CONCLUSION: By analyzing and integrating all available high-throughput genomic, transcriptomic, miRNA, methylation and proteomic data from human HCC tissue, we identified EGFR, ß1-integrin and axon guidance as pathway dependencies in HCC. These are master regulators of key pathways in HCC, such as the mTOR, Ras/Raf/MAPK and p53 pathways. The genes implicated in these pathways had prognostic value in HCC, with Netrin and Slit3 being novel proteins of prognostic importance to HCC. Based on this integrative analysis, EGFR, and ß1-integrin are master regulators that could serve as potential therapeutic targets in HCC.

Diabetogenic Effects of Immunosuppression: An Integrative Analysis.

BACKGROUND: Posttransplant diabetes mellitus (PTDM) affects up to 50% of solid organ transplant recipients and compromises long-term outcomes. The goal of this study was to investigate how immunosuppressants affect gene expression in a manner that increases diabetes risk, by performing integrative analysis on publicly available, high-throughput gene expression data. METHODS: All high-throughput gene expression datasets of solid organ transplant recipients were retrieved from the Gene Expression Omnibus. Significantly dysregulated genes and pathways were determined, and those in common with type 2 diabetes were identified. THP-1 and HepG2 cells were exposed in vitro to tacrolimus, and validation of genes involved in insulin signaling and glucose metabolism was performed using specific arrays. These cells were then treated with the hypoglycemic agents, metformin, and insulin to assess for appropriate reversion of specific diabetogenic genes. RESULTS: Insulin signaling and secretion were the most commonly dysregulated pathways that overlapped with diabetes in transplant recipients. KRAS, GRB2, PCK2, BCL2L1, INSL3, DOK3, and PTPN1 were among the most significantly upregulated genes in both immunosuppression and diabetes subsets and were appropriately reverted by metformin as confirmed in vitro. CONCLUSIONS: We discovered that the significantly dysregulated genes in the context of immunosuppression are implicated in insulin signaling and insulin secretion, as a manifestation of pancreatic ß-cell function. In vitro validation confirmed key diabetes-related genes in the context of immunosuppression. Further analysis and in vitro validation revealed that metformin optimally reverts diabetogenic genes dysregulated in the context of immunosuppression. The optimal therapeutic management of posttransplant diabetes mellitus needs to be further investigated, taking into account the mechanistic impact of immunosuppressants.

The basis of liver regeneration: A systems biology approach.

INTRODUCTION: Liver regeneration is a normal response to liver injury. The aim of this study was to determine the molecular basis of liver regeneration, through an integrative analysis of high-throughput gene expression datasets. METHODS: We identified and curated datasets pertaining to liver regeneration from the Gene Expression Omnibus, where regenerating liver tissue was compared to healthy liver samples. The key dysregulated genes and pathways were identified using Ingenuity Pathway Analysis software. There were three eligible datasets in total. RESULTS: In the early phase after hepatectomy, inflammatory pathways such as Nrf2 oxidative stress-mediated response and cytokine signaling were significantly upregulated. At peak regeneration, we discovered that cell cycle genes were predominantly expressed to promote cell proliferation. Using the Betweenness centrality algorithm, we discovered that Jun is the key central gene in liver regeneration. Calcineurin inhibitors may inhibit liver regeneration, based on predictive modeling. CONCLUSION: There is a paucity of human literature in defining the molecular mechanisms of liver regeneration along a time continuum. Nonetheless, using an integrative computational analysis approach to the available high-throughput data, we determine that the oxidative stress response and cytokine signaling are key early after hepatectomy, whereas cell cycle control is important at peak regeneration. The transcription factor Jun is central to liver regeneration and a potential therapeutic target. Future studies of regeneration in humans along a time continuum are needed to better define the underlying mechanisms, and ultimately enhance care of patients with acute and chronic liver failure while awaiting transplant.

Epigenetic basis of hepatocellular carcinoma: A network-based integrative meta-analysis.

AIM: To identify the key epigenetically modulated genes and pathways in HCC by performing an integrative meta-analysis of all major, well-annotated and publicly available methylation datasets using tools of network analysis. METHODS: PubMed and Gene Expression Omnibus were searched for genome-wide DNA methylation datasets. Patient clinical and demographic characteristics were obtained. DNA methylation data were integrated using the Ingenuity Pathway Analysis, a software package for visualizing and analyzing biological networks. Pathway enrichment analysis was performed using IPA, which also provides literature-driven and computationally-predicted annotations for significant association of genes to curated molecular pathways. RESULTS: From an initial 928 potential abstracts, we identified and analyzed 11 eligible high-throughput methylation datasets representing 354 patients. A significant proportion of studies did not provide concomitant clinical data. In the promoter region, HIST1H2AJ and SPDYA were the most commonly methylated, whereas HRNBP3 gene was the most commonly hypomethylated. ESR1 and ERK were central genes in the principal networks. The pathways most associated with the frequently methylated genes were G-protein coupled receptor and cAMP-mediated signalling. CONCLUSION: Using an integrative network-based analysis approach of genome-wide DNA methylation data of both the promoter and body of genes, we identified G-protein coupled receptor signalling as the most highly associated with HCC. This encompasses a diverse range of cancer pathways, such as the PI3K/Akt/mTOR and Ras/Raf/MAPK pathways, and is therefore supportive of previous literature on gene expression in HCC. However, there are novel targetable genes such as HIST1H2AJ that are epigenetically modified, suggesting their potential as biomarkers and for therapeutic targeting of the HCC epigenome.

Systematic integrative analysis of gene expression identifies HNF4A as the central gene in pathogenesis of non-alcoholic steatohepatitis.

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in the Western world, and encompasses a spectrum from simple steatosis to steatohepatitis (NASH). There is currently no approved pharmacologic therapy against NASH, partly due to an incomplete understanding of its molecular basis. The goal of this study was to determine the key differentially expressed genes (DEGs), as well as those genes and pathways central to its pathogenesis. We performed an integrative computational analysis of publicly available gene expression data in NASH from GEO (GSE17470, GSE24807, GSE37031, GSE89632). The DEGs were identified using GEOquery, and only the genes present in at least three of the studies, to a total of 190 DEGs, were considered for further analyses. The pathways, networks, molecular interactions, functional analyses were generated through the use of Ingenuity Pathway Analysis (IPA). For selected networks, we computed the centrality using igraph package in R. Among the statistically significant predicted networks (p-val < 0.05), three were of most biological interest: the first is involved in antimicrobial response, inflammatory response and immunological disease, the second in cancer, organismal injury and development and the third in metabolic diseases. We discovered that HNF4A is the central gene in the network of NASH connected to metabolic diseases and that it regulates HNF1A, an additional transcription regulator also involved in lipid metabolism. Therefore, we show, for the first time to our knowledge, that HNF4A is central to the pathogenesis of NASH. This adds to previous literature demonstrating that HNF4A regulates the transcription of genes involved in the progression of NAFLD, and that HNF4A genetic variants play a potential role in NASH progression.

Retraction: Circulating plant miRNAs can regulate human gene expression in vitro.

This corrects the article DOI: 10.1038/srep32773.

Circulating plant miRNAs can regulate human gene expression in vitro.

While Brassica oleracea vegetables have been linked to cancer prevention, the exact mechanism remains unknown. Regulation of gene expression by cross-species microRNAs has been previously reported; however, its link to cancer suppression remains unexplored. In this study we address both issues. We confirm plant microRNAs in human blood in a large nutrigenomics study cohort and in a randomized dose-controlled trial, finding a significant positive correlation between the daily amount of broccoli consumed and the amount of microRNA in the blood. We also demonstrate that Brassica microRNAs regulate expression of human genes and proteins in vitro, and that microRNAs cooperate with other Brassica-specific compounds in a possible cancer-preventive mechanism. Combined, we provide strong evidence and a possible multimodal mechanism for broccoli in cancer prevention.

Robust quantitative scratch assay.

UNLABELLED: The wound healing assay (or scratch assay) is a technique frequently used to quantify the dependence of cell motility-a central process in tissue repair and evolution of disease-subject to various treatments conditions. However processing the resulting data is a laborious task due its high throughput and variability across images. This Robust Quantitative Scratch Assay algorithm introduced statistical outputs where migration rates are estimated, cellular behaviour is distinguished and outliers are identified among groups of unique experimental conditions. Furthermore, the RQSA decreased measurement errors and increased accuracy in the wound boundary at comparable processing times compared to previously developed method (TScratch). AVAILABILITY AND IMPLEMENTATION: The RQSA is freely available at: http://ophid.utoronto.ca/RQSA/RQSA_Scripts.zip The image sets used for training and validation and results are available at: (http://ophid.utoronto.ca/RQSA/trainingSet.zip, http://ophid.utoronto.ca/RQSA/validationSet.zip, http://ophid.utoronto.ca/RQSA/ValidationSetResults.zip, http://ophid.utoronto.ca/RQSA/ValidationSet_H1975.zip, http://ophid.utoronto.ca/RQSA/ValidationSet_H1975Results.zip, http://ophid.utoronto.ca/RQSA/RobustnessSet.zip, http://ophid.utoronto.ca/RQSA/RobustnessSet.zip). Supplementary Material is provided for detailed description of the development of the RQSA. CONTACT: juris@ai.utoronto.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prioritizing therapeutics for lung cancer: an integrative meta-analysis of cancer gene signatures and chemogenomic data.

Repurposing FDA-approved drugs with the aid of gene signatures of disease can accelerate the development of new therapeutics. A major challenge to developing reliable drug predictions is heterogeneity. Different gene signatures of the same disease or drug treatment often show poor overlap across studies, as a consequence of both biological and technical variability, and this can affect the quality and reproducibility of computational drug predictions. Existing algorithms for signature-based drug repurposing use only individual signatures as input. But for many diseases, there are dozens of signatures in the public domain. Methods that exploit all available transcriptional knowledge on a disease should produce improved drug predictions. Here, we adapt an established meta-analysis framework to address the problem of drug repurposing using an ensemble of disease signatures. Our computational pipeline takes as input a collection of disease signatures, and outputs a list of drugs predicted to consistently reverse pathological gene changes. We apply our method to conduct the largest and most systematic repurposing study on lung cancer transcriptomes, using 21 signatures. We show that scaling up transcriptional knowledge significantly increases the reproducibility of top drug hits, from 44% to 78%. We extensively characterize drug hits in silico, demonstrating that they slow growth significantly in nine lung cancer cell lines from the NCI-60 collection, and identify CALM1 and PLA2G4A as promising drug targets for lung cancer. Our meta-analysis pipeline is general, and applicable to any disease context; it can be applied to improve the results of signature-based drug repurposing by leveraging the large number of disease signatures in the public domain.