NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

In search of epigenetic marks in testes and sperm cells of differentially fed boars

In search of transmittable epigenetic marks we investigated gene expression in testes and sperm cells of differentially fed F0 boars from a three generation pig feeding experiment that showed phenotypic differences in the F2 generation. RNA samples from 8 testes of boars that received either a diet enriched in methylating micronutrients or a control diet were analyzed by microarray analysis. We found moderate differential expression between testes of differentially fed boars with a high FDR of 0.82 indicating that most of the differentially expressed genes were false positives. Nevertheless, we performed a pathway analysis and found disparate pathway maps of development_A2B receptor: action via G-protein alpha s, cell adhesion_Tight junctions and cell adhesion_Endothelial cell contacts by junctional mechanisms which show inconclusive relation to epigenetic inheritance. Four RNA samples from sperm cells of these differentially fed boars were analyzed by RNA-Seq methodology. We found no differential gene expression in sperm cells of the two groups (adjusted P-value>0.05). Nevertheless, we also explored gene expression in sperm by a pathway analysis showing that genes were enriched for the pathway maps of bacterial infections in cystic fibrosis (CF) airways, glycolysis and gluconeogenesis p.3 and cell cycle_Initiation of mitosis. Again, these pathway maps are miscellaneous without an obvious relationship to epigenetic inheritance. It is concluded that the methylating micronutrients moderately if at all affects RNA expression in testes of differentially fed boars. Furthermore, gene expression in sperm cells is not significantly affected by extensive supplementation of methylating micronutrients and thus RNA molecules could not be established as the epigenetic mark in this feeding experiment.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

Metabolomic tissue signature in human non-alcoholic fatty liver disease identifies protective candidate metabolites

Background Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disorder in industrialized countries, yet its pathophysiology is incompletely understood. Small-molecule metabolite screens may offer new insights into disease mechanisms and reveal new treatment targets. Methods Discovery (N = 33) and replication (N = 66) of liver biopsies spanning the range from normal liver histology to non-alcoholic steatohepatitis (NASH) were ascertained ensuring rapid freezing under 30 s in patients. 252 metabolites were assessed using GC/MS. Replicated metabolites were evaluated in a murine high-fat diet model of NAFLD. Results In a two-stage metabolic screening, hydroquinone (HQ, pcombined = 3.0 × 10−4) and nicotinic acid (NA, pcombined = 3.9 × 10−9) were inversely correlated with histological NAFLD severity. A murine high-fat diet model of NAFLD demonstrated a protective effect of these two substances against NAFLD: Supplementation with 1% HQ reduced only liver steatosis, whereas 0.6% NA reduced both liver fat content and serum transaminase levels and induced a complex regulatory network of genes linked to NALFD pathogenesis in a global expression pathway analysis. Human nutritional intake of NA equivalent was also consistent with a protective effect of NA against NASH progression. Conclusion This first small-molecular screen of human liver tissue identified two replicated protective metabolites. Either the use of NA or targeting its regulatory pathways might be explored to treat or prevent human NAFLD.

Generation of a murine hepatic angiosarcoma cell line and reproducible mouse tumor model

Hepatic angiosarcoma (AS) is a rare and highly aggressive tumor of endothelial origin with dismal prognosis. Studies of the molecular biology of AS and treatment options are limited as animal models are rare. We have previously shown that inducible knockout of Notch1 in mice leads to spontaneous formation of hepatic AS. The aims of this study were to: (1) establish and characterize a cell line derived from this murine AS, (2) identify molecular pathways involved in the pathogenesis and potential therapeutic targets, and (3) generate a tumor transplantation model. AS cells retained specific endothelial properties such as tube formation activity, as well as expression of CD31 and Von Willebrand factor. However, electron microscopy analysis revealed signs of dedifferentiation with loss of fenestrae and loss of contact inhibition. Microarray and pathway analysis showed substantial changes in gene expression and revealed activation of the Myc pathway. Exposing the AS cells to sorafenib reduced migration, filopodia dynamics, and cell proliferation but did not induce apoptosis. In addition, sorafenib suppressed ERK phosphorylation and expression of cyclin D2. Injection of AS cells into NOD/SCID mice resulted in formation of undifferentiated tumors, confirming the tumorigenic potential of these cells. In summary, we established and characterized a murine model of spontaneous AS formation and hepatic AS cell lines as a useful in vitro tool. Our data demonstrate antitumor activity of sorafenib in AS cells with potent inhibition of migration, filopodia formation, and cell proliferation, supporting further evaluation of sorafenib as a novel treatment strategy. In addition, AS cell transplantation provides a subcutaneous tumor model useful for in vivo preclinical drug testing.Laboratory Investigation advance online publication, 24 November 2014; doi:10.1038/labinvest.2014.141.

Insulin, CCAAT/enhancer-binding proteins and lactate regulate the human 11β-hydroxysteroid dehydrogenase type 2 gene expression in colon cancer cell lines.

11β-Hydroxysteroid dehydrogenases (11beta-HSD) modulate mineralocorticoid receptor transactivation by glucocorticoids and regulate access to the glucocorticoid receptor. The isozyme 11beta-HSD2 is selectively expressed in mineralocorticoid target tissues and its activity is reduced in various disease states with abnormal sodium retention and hypertension, including the apparent mineralocorticoid excess. As 50% of patients with essential hypertension are insulin resistant and hyperinsulinemic, we hypothesized that insulin downregulates the 11beta-HSD2 activity. In the present study we show that insulin reduced the 11beta-HSD2 activity in cancer colon cell lines (HCT116, SW620 and HT-29) at the transcriptional level, in a time and dose dependent manner. The downregulation was reversible and required new protein synthesis. Pathway analysis using mRNA profiling revealed that insulin treatment modified the expression of the transcription factor family C/EBPs (CCAAT/enhancer-binding proteins) but also of glycolysis related enzymes. Western blot and real time PCR confirmed an upregulation of C/EBP beta isoforms (LAP and LIP) with a more pronounced increase in the inhibitory isoform LIP. EMSA and reporter gene assays demonstrated the role of C/EBP beta isoforms in HSD11B2 gene expression regulation. In addition, secretion of lactate, a byproduct of glycolysis, was shown to mediate insulin-dependent HSD11B2 downregulation. In summary, we demonstrate that insulin downregulates HSD11B2 through increased LIP expression and augmented lactate secretion. Such mechanisms are of interest and potential significance for sodium reabsorption in the colon.

miRNA-26b Overexpression in Ulcerative Colitis-associated Carcinogenesis

BACKGROUND Longstanding ulcerative colitis (UC) bears a high risk for development of UC-associated colorectal carcinoma (UCC). The inflammatory microenvironment influences microRNA expression, which in turn deregulates target gene expression. microRNA-26b (miR-26b) was shown to be instrumental in normal tissue growth and differentiation. Thus, we aimed to investigate the impact of miR-26b in inflammation-associated colorectal carcinogenesis. METHODS Two different cohorts of patients were investigated. In the retrospective group, a tissue microarray with 38 samples from 17 UC/UCC patients was used for miR-26b in situ hybridization and quantitative reverse transcription polymerase chain reaction analyses. In the prospective group, we investigated miR-26b expression in 25 fresh-frozen colon biopsies and corresponding serum samples of 6 UC and 15 non-UC patients, respectively. In silico analysis, Ago2-RNA immunoprecipitation, luciferase reporter assay, quantitative reverse transcription polymerase chain reaction examination, and miR-26b mimic overexpression were employed for target validation. RESULTS miR-26b expression was shown to be upregulated with disease progression in tissues and serum of UC and UCC patients. Using miR-26b and Ki-67 expression levels, an UCC was predicted with high accuracy. We identified 4 novel miR-26b targets (DIP1, MDM2, CREBBP, BRCA1). Among them, the downregulation of the E3 ubiquitin ligase DIP1 was closely related to death-associated protein kinase stabilization along the normal mucosa-UC-UCC sequence. In silico functional pathway analysis revealed that the common cellular pathways affected by miR-26b are highly related to cancerogenesis and the development of gastrointestinal diseases. CONCLUSIONS We suggest that miR-26b could serve as a biomarker for inflammation-associated processes in the gastrointestinal system. Because miR-26b expression is downregulated in sporadic colon cancer, it could discriminate between UCC and the sporadic cancer type.

Common target genes of palatal and gingival fibroblasts for EMD: the microarray approach.

BACKGROUND AND OBJECTIVE Connective tissue grafts are frequently applied, together with Emdogain(®) , for root coverage. However, it is unknown whether fibroblasts from the gingiva and from the palate respond similarly to Emdogain. The aim of this study was therefore to evaluate the effect of Emdogain(®) on fibroblasts from palatal and gingival connective tissue using a genome-wide microarray approach. MATERIAL AND METHODS Human palatal and gingival fibroblasts were exposed to Emdogain(®) and RNA was subjected to microarray analysis followed by gene ontology screening with Database for Annotation, Visualization and Integrated Discovery functional annotation clustering, Kyoto Encyclopedia of Genes and Genomes pathway analysis and the Search Tool for the Retrieval of Interacting Genes/Proteins functional protein association network. Microarray results were confirmed by quantitative RT-PCR analysis. RESULTS The transcription levels of 106 genes were up-/down-regulated by at least five-fold in both gingival and palatal fibroblasts upon exposure to Emdogain(®) . Gene ontology screening assigned the respective genes into 118 biological processes, six cellular components, eight molecular functions and five pathways. Among the striking patterns observed were the changing expression of ligands targeting the transforming growth factor-beta and gp130 receptor family as well as the transition of mesenchymal epithelial cells. Moreover, Emdogain(®) caused changes in expression of receptors for chemokines, lipids and hormones, and for transcription factors such as SMAD3, peroxisome proliferator-activated receptor gamma and those of the ETS family. CONCLUSION The present data suggest that Emdogain(®) causes substantial alterations in gene expression, with similar patterns observed in palatal and gingival fibroblasts.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.

Potential REST Targets in Neural Stem Cells and Glioblastoma-Derived Stem Cells

Glioblastoma multiforme (GBM) tumors are the most common malignant primary brain tumors in adults. The current theory is that these tumors are caused by self-renewing glioblastoma-derived stem cells (GSCs). At the current time, the mechanisms that regulate self-renewal and other oncogenic properties of GSCs remain unknown. Recently, we found transcriptional repressor REST maintains self-renewal in neural stem cells (NSCs) and in GSCs. REST also regulates other oncogenic properties, such as apoptosis, invasion and proliferation. However, the mechanisms by which REST regulates these oncogenic properties are unknown. In an attempt to determine these mechanisms, we performed loss and gain-of-function experiments and genome-wide mRNA expression analysis in GSCs, and we were able to identify REST-regulated genes in GSCs. This was accomplished, after screening concordantly regulated genes in NSCs and GSCs, utilizing two RE1 databases, and setting two-fold expression as filters on the resulting genes. These results received further validation by qRT-PCR. Ingenuity Pathway Analysis (IPA) analysis further revealed the top REST target genes in GSCs were downstream targets of REST and/or involved in other cancers in other cell lines. IPA also revealed that many of the differentially-regulated genes identified in this study are involved in oncogenic properties seen in GBM, and which we believe are related to REST expression.

Experiment on marine fungus Microascus brevicaulis strain

The marine fungus Microascus brevicaulis strain LF580 is a non-model secondary metabolite producer with high yields of the two secondary metabolites scopularide A and B, which exhibit distinct activities against tumour cell lines. A mutant strain was obtained using UV mutagenesis, showing besides higher production levels faster growth and differences in pellet formation. Comparative proteomics were applied to gain deeper understanding of the regulation of production and of the physiology of this fungus and its mutant. For this purpose, an optimised protein extraction protocol was established. Here, we show the first proteome study of a marine fungus. In total, 4759 proteins were identified. The central metabolic pathway of LF580 could be mapped by using KEGG pathway analysis and GO annotation. Using iTRAQ labelling, 318 proteins were shown to be significantly regulated in the mutant strain: 189 were down- and 129 upregulated. Proteomics are a powerful tool for the understanding of regulatory aspects: The differences on proteome level could be attributed to a limited nutrient availability in wild type strain due to a strong pellet formation. This information can be applied to optimisation on strain and process level. The linkage between nutrient limitation and pellet formation in the non-model fungus M. brevicaulis is in consensus with the knowledge on model organisms like Aspergillus niger and Penicillium chrysogenum.

Enhanced Glycosylation and Sulfation of Secretory Proteoglycans Is Coupled to the Expression of a Basic Secretory Protein

We have used coexpression of a salivary basic proline-rich protein (PRP) along with a proline-rich proteoglycan (PRPg) in pituitary AtT-20 cells to examine the regulation of glycosaminoglycan (GAG) biosynthesis and the storage of these secretory products for regulated secretion. The basic PRP caused a dose-dependent increase in sulfation of PRPg and also increased the extent to which PRPg polypeptide backbones are modified by a GAG chain. The sulfation of an endogenous proteoglycan was similarly increased in the presence of basic PRP; however, other sulfated secretory products of AtT-20 cells were unaffected. These results imply that enzymes functioning in elongation and sulfation of proteoglycans are coordinately regulated and that their activities respond to a change in the milieu of the intracellular transport pathway. Analysis of the regulated secretion of both the basic PRP and PRPg has indicated that while the presence of the GAG chain improves the storage of PRPg, the presence of PRPg does not increase the storage of basic PRP. Therefore, sulfation of GAGs does not appear to be a primary factor in regulated secretory sorting.

The Arabidopsis dwarf1 Mutant Is Defective in the Conversion of 24-Methylenecholesterol to Campesterol in Brassinosteroid Biosynthesis1

Since the isolation and characterization of dwarf1-1 (dwf1-1) from a T-DNA insertion mutant population, phenotypically similar mutants, including deetiolated2 (det2), constitutive photomorphogenesis and dwarfism (cpd), brassinosteroid insensitive1 (bri1), and dwf4, have been reported to be defective in either the biosynthesis or the perception of brassinosteroids. We present further characterization of dwf1-1 and additional dwf1 alleles. Feeding tests with brassinosteroid-biosynthetic intermediates revealed that dwf1 can be rescued by 22α-hydroxycampesterol and downstream intermediates in the brassinosteroid pathway. Analysis of the endogenous levels of brassinosteroid intermediates showed that 24-methylenecholesterol in dwf1 accumulates to 12 times the level of the wild type, whereas the level of campesterol is greatly diminished, indicating that the defective step is in C-24 reduction. Furthermore, the deduced amino acid sequence of DWF1 shows significant similarity to a flavin adenine dinucleotide-binding domain conserved in various oxidoreductases, suggesting an enzymatic role for DWF1. In support of this, 7 of 10 dwf1 mutations directly affected the flavin adenine dinucleotide-binding domain. Our molecular characterization of dwf1 alleles, together with our biochemical data, suggest that the biosynthetic defect in dwf1 results in reduced synthesis of bioactive brassinosteroids, causing dwarfism.

Fungal metabolic model for human type I hereditary tyrosinaemia.

Type I hereditary tyrosinaemia (HT1) is a severe human inborn disease resulting from loss of fumaryl-acetoacetate hydrolase (Fah). Homozygous disruption of the gene encoding Fah in mice causes neonatal lethality, seriously limiting use of this animal as a model. We report here that fahA, the gene encoding Fah in the fungus Aspergillus nidulans, encodes a polypeptide showing 47.1% identity to its human homologue, fahA disruption results in secretion of succinylacetone (a diagnostic compound for human type I tyrosinaemia) and phenylalanine toxicity. We have isolated spontaneous suppressor mutations preventing this toxicity, presumably representing loss-of-function mutations in genes acting upstream of fahA in the phenylalanine catabolic pathway. Analysis of a class of these mutations demonstrates that loss of homogentisate dioxygenase (leading to alkaptonuria in humans) prevents the effects of a Fah deficiency. Our results strongly suggest human homogentisate dioxygenase as a target for HT1 therapy and illustrate the usefulness of this fungus as an alternative to animal models for certain aspects of human metabolic diseases.

Modeling hybridoma cell metabolism using a generic genome-scale metabolic model of Mus musculus

The reconstructed cellular metabolic network of Mus musculus, based on annotated genomic data, pathway databases, and currently available biochemical and physiological information, is presented. Although incomplete, it represents the first attempt to collect and characterize the metabolic network of a mammalian cell on the basis of genomic data. The reaction network is generic in nature and attempts to capture the carbon, energy, and nitrogen metabolism of the cell. The metabolic reactions were compartmentalized between the cytosol and the mitochondria, including transport reactions between the compartments and the extracellular medium. The reaction list consists of 872 internal metabolites involved in a total of 1220 reactions, whereof 473 relate to known open reading frames. Initial in silico analysis of the reconstructed model is presented.