Search

The cancer syndrome polymerase proofreading-associated polyposis results from germline mutations in the POLE and POLD1 genes. Mutations in the exonuclease domain of these genes are associated with hyper- and ultra-mutated tumors with a predominance of base substitutions resulting from faulty proofreading during DNA replication. When a new variant is identified by gene testing of POLE and POLD1, it is important to verify whether the variant is associated with PPAP or not, to guide genetic counseling of mutation carriers. In 2015, we reported the likely pathogenic (class 4) germline POLE c.1373A > T p.(Tyr458Phe) variant and we have now characterized this variant to verify that it is a class 5 pathogenic variant. For this purpose, we investigated (1) mutator phenotype in tumors from two carriers, (2) mutation frequency in cell-based mutagenesis assays, and (3) structural consequences based on protein modeling. Whole-exome sequencing of two tumors identified an ultra-mutator phenotype with a predominance of base substitutions, the majority of which are C > T. A SupF mutagenesis assay revealed increased mutation frequency in cells overexpressing the variant of interest as well as in isogenic cells encoding the variant. Moreover, exonuclease repair yeast-based assay supported defect in proofreading activity. Lastly, we present a homology model of human POLE to demonstrate structural consequences leading to pathogenic impact of the p.(Tyr458Phe) mutation. The three lines of evidence, taken together with updated co-segregation and previously published data, allow the germline variant POLE c.1373A > T p.(Tyr458Phe) to be reclassified as a class 5 variant. That means the variant is associated with PPAP.

We present an update of EpiFactors, a manually curated database providing information about epigenetic regulators, their complexes, targets, and products which is openly accessible at http://epifactors.autosome.org. An updated version of the EpiFactors contains information on 902 proteins, including 101 histones and protamines, and, as a main update, a newly curated collection of 124 lncRNAs involved in epigenetic regulation. The amount of publications concerning the role of lncRNA in epigenetics is rapidly growing. Yet, the resource that compiles, integrates, organizes, and presents curated information on lncRNAs in epigenetics is missing. EpiFactors fills this gap and provides data on epigenetic regulators in an accessible and user-friendly form. For 820 of the genes in EpiFactors, we include expression estimates across multiple cell types assessed by CAGE-Seq in the FANTOM5 project. In addition, the updated EpiFactors contains information on 73 protein complexes involved in epigenetic regulation. Our resource is practical for a wide range of users, including biologists, bioinformaticians and molecular/systems biologists.

The protein Vitellogenin (Vg) plays a central role in lipid transportation in most egg-laying animals. High Vg levels correlate with stress resistance and lifespan potential in honey bees (Apis mellifera). Vg is the primary circulating zinc-carrying protein in honey bees. Zinc is an essential metal ion in numerous biological processes, including the function and structure of many proteins. Measurements of Zn2+ suggest a variable number of ions per Vg molecule in different animal species, but the molecular implications of zinc-binding by this protein are not well understood. We used inductively coupled plasma mass spectrometry (ICP-MS) to determine that, on average, each honey bee Vg molecule binds 3 Zn2+-ions. Our full-length protein structure and sequence analysis revealed seven potential zinc-binding sites. These are located in the β-barrel and α-helical subdomains of the N-terminal domain, the lipid binding site, and the cysteine-rich C-terminal region of unknown function. Interestingly, two potential zinc-binding sites in the β-barrel can support a proposed role for this structure in DNA-binding. Overall, our findings illustrate the capacity of honey bee Vg to bind zinc at several functional regions, indicating that Zn2+-ions are important for many of the activities of this protein. In addition to being potentially relevant for other egg-laying species, these insights provide a platform for studies of metal ions in bee health, which is of global interest due to recent declines in pollinator numbers.

BackgroundSecretion of the metabolites citrate and spermine into prostate lumen is a unique hallmark for normal prostate epithelial cells. However, the identity of the genes controlling citrate and spermine secretion remains mostly unknown despite their obvious relevance for progression to aggressive prostate cancer.Materials & MethodsIn this study, we have correlated simultaneous measurement of citrate/spermine and transcriptomics data. We have refined these gene correlations in 12 prostate cancer cohorts containing 2915 tissue samples to create a novel gene signature of 150 genes connected with citrate and spermine secretion. We further explored the signature in public data, interrogating over 18 000 samples from various tissues and model systems, including 3826 samples from prostate and prostate cancer.ResultsIn prostate cancer, the expression of this gene signature is gradually lost in tissue from normal epithelial cells through PIN, low grade (Gleason = 8) and metastatic lesions. The accuracy of the signature is validated by its unique enrichment in prostate compared to other tissues, and its strong enrichment in epithelial tissue compartments compared to stroma. Several zinc-binding proteins that are not previously investigated in the prostate are present in the gene signature, suggesting new mechanisms for controlling zinc homeostasis in citrate/spermine secretion. However, the absence of the gene signature in all common prostate normal and cancer cell-lines, as well as prostate organoids, underlines the challenge to study the role of these genes during prostate cancer progression in model systems.ConclusionsA large collection of transcriptomics data integrated with metabolomics identifies the genes related to citrate and spermine secretion in the prostate, and show that the expression of these genes gradually decreases on the path towards aggressive prostate cancer. In addition, the study questions the relevance of currently available model systems to study metabolism in prostate cancer development.

Background Detection of copy number variation (CNV) in genes associated with disease is important in genetic diagnostics, and next generation sequencing (NGS) technology provides data that can be used for CNV detection. However, CNV detection based on NGS data is in general not often used in diagnostic labs as the data analysis is challenging, especially with data from targeted gene panels. Wet lab methods like MLPA (MRC Holland) are widely used, but are expensive, time consuming and have gene-specific limitations. Our aim has been to develop a bioinformatic tool for CNV detection from NGS data in medical genetic diagnostic samples. Results Our computational pipeline for detection of CNVs in NGS data from targeted gene panels utilizes coverage depth of the captured regions and calculates a copy number ratio score for each region. This is computed by comparing the mean coverage of the sample with the mean coverage of the same region in other samples, defined as a pool. The pipeline selects pools for comparison dynamically from previously sequenced samples, using the pool with an average coverage depth that is nearest to the one of the samples. A sliding window-based approach is used to analyze each region, where length of sliding window and sliding distance can be chosen dynamically to increase or decrease the resolution. This helps in detecting CNVs in small or partial exons. With this pipeline we have correctly identified the CNVs in 36 positive control samples, with sensitivity of 100% and specificity of 91%. We have detected whole gene level deletion/duplication, single/multi exonic level deletion/duplication, partial exonic deletion and mosaic deletion. Since its implementation in mid-2018 it has proven its diagnostic value with more than 45 CNV findings in routine tests. Conclusions With this pipeline as part of our diagnostic practices it is now possible to detect partial, single or multi-exonic, and intragenic CNVs in all genes in our target panel. This has helped our diagnostic lab to expand the portfolio of genes where we offer CNV detection, which previously was limited by the availability of MLPA kits.

Objective Properties of gene products can be described or annotated with Gene Ontology (GO) terms. But for many genes we have limited information about their products, for example with respect to function. This is particularly true for long non-coding RNAs (lncRNAs), where the function in most cases is unknown. However, it has been shown that annotation as described by GO terms to some extent can be predicted by enrichment analysis on properties of co-expressed genes. Results GAPGOM integrates two relevant algorithms, lncRNA2GOA and TopoICSim, into a user-friendly R package. Here lncRNA2GOA does annotation prediction by co-expression, whereas TopoICSim estimates similarity between GO graphs, which can be used for benchmarking of prediction performance, but also for comparison of GO graphs in general. The package provides an improved implementation of the original tools, with substantial improvements in performance and documentation, unified interfaces, and additional features. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Background: Many types of data from genomic analyses can be represented as genomic tracks, i.e. features linked to the genomic coordinates of a reference genome. Examples of such data are epigenetic DNA methylation data, ChIP-seq peaks, germline or somatic DNA variants, as well as RNA-seq expression levels. Researchers often face difficulties in locating, accessing and combining relevant tracks from external sources, as well as locating the raw data, reducing the value of the generated information. Description of work: We propose to advance the application of FAIR data principles (Findable, Accessible, Interoperable, and Reusable) to produce searchable metadata for genomic tracks. Findability and Accessibility of metadata can then be ensured by a track search service that integrates globally identifiable metadata from various track hubs in the Track Hub Registry and other relevant repositories. Interoperability and Reusability need to be ensured by the specification and implementation of a basic set of recommendations for metadata. We have tested this concept by developing such a specification in a JSON Schema, called FAIRtracks, and have integrated it into a novel track search service, called TrackFind. We demonstrate practical usage by importing datasets through TrackFind into existing examples of relevant analytical tools for genomic tracks: EPICO and the GSuite HyperBrowser. Conclusion: We here provide a first iteration of a draft standard for genomic track metadata, as well as the accompanying software ecosystem. It can easily be adapted or extended to future needs of the research community regarding data, methods and tools, balancing the requirements of both data submitters and analytical end-users.

The k-Nearest Neighbor (kNN) classifier represents a simple and very general approach to classification. Still, the performance of kNN classifiers can often compete with more complex machine-learning algorithms. The core of kNN depends on a “guilt by association” principle where classification is performed by measuring the similarity between a query and a set of training patterns, often computed as distances. The relative performance of kNN classifiers is closely linked to the choice of distance or similarity measure, and it is therefore relevant to investigate the effect of using different distance measures when comparing biomedical data. In this study on classification of cancer data sets, we have used both common and novel distance measures, including the novel distance measures Sobolev and Fisher, and we have evaluated the performance of kNN with these distances on 4 cancer data sets of different type. We find that the performance when using the novel distance measures is comparable to the performance with more well-established measures, in particular for the Sobolev distance. We define a robust ranking of all the distance measures according to overall performance. Several distance measures show robust performance in kNN over several data sets, in particular the Hassanat, Sobolev, and Manhattan measures. Some of the other measures show good performance on selected data sets but seem to be more sensitive to the nature of the classification data. It is therefore important to benchmark distance measures on similar data prior to classification to identify the most suitable measure in each case. https://creativecommons.org/licenses/by/4.0/This article is distributed under the terms of the Creative Commons Attribution 4.0 License (https://creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access page (https://us.sagepub.com/en-us/nam/open-access-at-sage).

Background Diseases like cancer will lead to changes in gene expression, and it is relevant to identify key regulatory genes that can be linked directly to these changes. This can be done by computing a Regulatory Impact Factor (RIF) score for relevant regulators. However, this computation is based on estimating correlated patterns of gene expression, often Pearson correlation, and an assumption about a set of specific regulators, normally transcription factors. This study explores alternative measures of correlation, using the Fisher and Sobolev metrics, and an extended set of regulators, including epigenetic regulators and long non-coding RNAs (lncRNAs). Data on prostate cancer have been used to explore the effect of these modifications. Results A tool for computation of RIF scores with alternative correlation measures and extended sets of regulators was developed and tested on gene expression data for prostate cancer. The study showed that the Fisher and Sobolev metrics lead to improved identification of well-documented regulators of gene expression in prostate cancer, and the sets of identified key regulators showed improved overlap with previously defined gene sets of relevance to cancer. The extended set of regulators lead to identification of several interesting candidates for further studies, including lncRNAs. Several key processes were identified as important, including spindle assembly and the epithelial-mesenchymal transition (EMT). Conclusions The study has shown that using alternative metrics of correlation can improve the performance of tools based on correlation of gene expression in genomic data. The Fisher and Sobolev metrics should be considered also in other correlation-based applications.

Background Prostate cancer (PCa) has the highest incidence rates of cancers in men in western countries. Unlike several other types of cancer, PCa has few genetic drivers, which has led researchers to look for additional epigenetic and transcriptomic contributors to PCa development and progression. Especially datasets on DNA methylation, the most commonly studied epigenetic marker, have recently been measured and analysed in several PCa patient cohorts. DNA methylation is most commonly associated with downregulation of gene expression. However, positive associations of DNA methylation to gene expression have also been reported, suggesting a more diverse mechanism of epigenetic regulation. Such additional complexity could have important implications for understanding prostate cancer development but has not been studied at a genome-wide scale. Results In this study, we have compared three sets of genome-wide single-site DNA methylation data from 870 PCa and normal tissue samples with multi-cohort gene expression data from 1117 samples, including 532 samples where DNA methylation and gene expression have been measured on the exact same samples. Genes were classified according to their corresponding methylation and expression profiles. A large group of hypermethylated genes was robustly associated with increased gene expression (UPUP group) in all three methylation datasets. These genes demonstrated distinct patterns of correlation between DNA methylation and gene expression compared to the genes showing the canonical negative association between methylation and expression (UPDOWN group). This indicates a more diversified role of DNA methylation in regulating gene expression than previously appreciated. Moreover, UPUP and UPDOWN genes were associated with different compartments — UPUP genes were related to the structures in nucleus, while UPDOWN genes were linked to extracellular features. Conclusion We identified a robust association between hypermethylation and upregulation of gene expression when comparing samples from prostate cancer and normal tissue. These results challenge the classical view where DNA methylation is always associated with suppression of gene expression, which underlines the importance of considering corresponding expression data when assessing the downstream regulatory effect of DNA methylation. © The Author(s). 2020 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Germline variants inactivating the mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2 cause Lynch syndrome that implies an increased cancer risk, where colon and endometrial cancer are the most frequent. Identification of these pathogenic variants is important to identify endometrial cancer patients with inherited increased risk of new cancers, in order to offer them lifesaving surveillance. However, several other genes are also part of the MMR pathway. It is therefore relevant to search for variants in additional genes that may be associated with cancer risk by including all known genes involved in the MMR pathway. Next-generation sequencing was used to screen 22 genes involved in the MMR pathway in constitutional DNA extracted from full blood from 199 unselected endometrial cancer patients. Bioinformatic pipelines were developed for identification and functional annotation of variants, using several different software tools and custom programs. This facilitated identification of 22 exonic, 4 UTR and 9 intronic variants that could be classified according to pathogenicity. This study has identified several germline variants in genes of the MMR pathway that potentially may be associated with an increased risk for cancer, in particular endometrial cancer, and therefore are relevant for further investigation. We have also developed bioinformatics strategies to analyse targeted sequencing data, including low quality data and genomic regions outside of the protein coding exons of the relevant genes. Copyright: © 2020 Singh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Most genes coding for long non-coding RNAs (lncRNAs) have unknown function. There is therefore a strong interest in approaches for function prediction of such genes. One possible approach is to hypothesize that genes with similar patterns of gene expression across multiple conditions may be involved in similar processes, which can be identified by enrichment analysis on gene annotations. This paper describes an implementation of this approach, tested on lncRNAs associated with cancer. The paper describes a submission for The LncRNA Function Prediction Challenge at the BIBM 2019 workshop on lncRNAs.

Proliferating cell nuclear antigen (PCNA), a member of the highly conserved DNA sliding clamp family, is an essential protein for cellular processes including DNA replication and repair. A large number of proteins from higher eukaryotes contain one of two PCNA‐interacting motifs: PCNA‐interacting protein box (PIP box) and AlkB homologue 2 PCNA‐interacting motif (APIM). APIM has been shown to be especially important during cellular stress. PIP box is known to be functionally conserved in yeast, and here, we show that this is also the case for APIM. Several of the 84 APIM‐containing yeast proteins are associated with cellular signaling as hub proteins, which are able to interact with a large number of other proteins. Cellular signaling is highly conserved throughout evolution, and we recently suggested a novel role for PCNA as a scaffold protein in cellular signaling in human cells. A cell‐penetrating peptide containing the APIM sequence increases the sensitivity toward the chemotherapeutic agent cisplatin in both yeast and human cells, and both yeast and human cells become hypersensitive when the Hog1/p38 MAPK pathway is blocked. These results suggest that the interactions between APIM‐containing signaling proteins and PCNA during the DNA damage response is evolutionary conserved between yeast and mammals and that PCNA has a role in cellular signaling also in yeast. © 2018 The Authors. Published by FEBS Press and John Wiley & Sons Ltd. This is an open access article under the terms of the Creative Commons Attribution License

// Elise Sandsmark 1 , Ailin Falkmo Hansen 1 , Kirsten M. Selnaes 1 , Helena Bertilsson 2, 3 , Anna M. Bofin 4 , Alan J. Wright 5 , Trond Viset 6 , Elin Richardsen 7, 8 , Finn Drablos 3 , Tone F. Bathen 1 , May-Britt Tessem 1, * , Morten B. Rye 2, 3, 9, * 1 Department of Circulation and Medical Imaging, Faculty of Medicine, NTNU - Norwegian University of Science and Technology, Trondheim, Norway 2 Department of Urology, St. Olavs Hospital, Trondheim University Hospital, Norway 3 Department of Cancer Research and Molecular Medicine, Faculty of Medicine, NTNU - Norwegian University of Science and Technology, Trondheim, Norway 4 Department of Laboratory Medicine, Children’s and Women’s Health, Faculty of Medicine, NTNU - Norwegian University of Science and Technology, Trondheim, Norway 5 Cancer Research UK Cambridge Institute, University of Cambridge, United Kingdom 6 Department of Pathology and Medical Genetics, St. Olavs Hospital, Trondheim University Hospital, Norway 7 Department of Medical Biology, UiT - The Arctic University of Norway, Tromso, Norway 8 Department of Clinical Pathology, University Hospital of North Norway, Tromso, Norway 9 St. Olavs Hospital, Trondheim University Hospital, Norway * These authors have contributed equally to this work Correspondence to: Elise Sandsmark, email: elise.sandsmark@gmail.com Morten B. Rye, email: morten.rye@ntnu.no Keywords: EMT, gene expression signature, biochemical recurrence, spectroscopy, MRSI Received: August 26, 2016 Accepted: November 23, 2016 Published: December 24, 2016 ABSTRACT Activation of the Canonical Wnt pathway (CWP) has been linked to advanced and metastatic prostate cancer, whereas the Wnt5a-induced non-canonical Wnt pathway (NCWP) has been associated with both good and poor prognosis. A newly discovered NCWP, Wnt5/Fzd2, has been shown to induce epithelial-to-mesenchymal transition (EMT) in cancers, but has not been investigated in prostate cancer. The aim of this study was to investigate if the CWP and NCWP, in combination with EMT, are associated with metabolic alterations, aggressive disease and biochemical recurrence in prostate cancer. An initial analysis was performed using integrated transcriptomics, ex vivo and in vivo metabolomics, and histopathology of prostatectomy samples (n=129), combined with at least five-year follow-up. This analysis detected increased activation of NCWP through Wnt5a/ Fzd2 as the most common mode of Wnt activation in prostate cancer. This activation was associated with increased expression of EMT markers and higher Gleason score. The transcriptional association between NCWP and EMT was confirmed in five other publicly available patient cohorts (1519 samples in total). A novel gene expression signature of concordant activation of NCWP and EMT (NCWP-EMT) was developed, and this signature was significantly associated with metastasis and shown to be a significant predictor of biochemical recurrence. The NCWP-EMT signature was also associated with decreased concentrations of the metabolites citrate and spermine, which have previously been linked to aggressive prostate cancer. Our results demonstrate the importance of NCWP and EMT in prostate cancer aggressiveness, suggest a novel gene signature for improved risk stratification, and give new molecular insight.

Sequencing technologies have changed not only our approaches to classical genetics, but also the field of epigenetics. Specific methods allow scientists to identify novel genome-wide epigenetic patterns of DNA methylation down to single-nucleotide resolution. DNA methylation is the most researched epigenetic mark involved in various processes in the human cell, including gene regulation and development of diseases, such as cancer. Increasing numbers of DNA methylation sequencing datasets from human genome are produced using various platforms—from methylated DNA precipitation to the whole genome bisulfite sequencing. Many of those datasets are fully accessible for repeated analyses. Sequencing experiments have become routine in laboratories around the world, while analysis of outcoming data is still a challenge among the majority of scientists, since in many cases it requires advanced computational skills. Even though various tools are being created and published, guidelines for their selection are often not clear, especially to non-bioinformaticians with limited experience in computational analyses. Separate tools are often used for individual steps in the analysis, and these can be challenging to manage and integrate. However, in some instances, tools are combined into pipelines that are capable to complete all the essential steps to achieve the result. In the case of DNA methylation sequencing analysis, the goal of such pipeline is to map sequencing reads, calculate methylation levels, and distinguish differentially methylated positions and/or regions. The objective of this review is to describe basic principles and steps in the analysis of DNA methylation sequencing data that in particular have been used for mammalian genomes, and more importantly to present and discuss the most pronounced computational pipelines that can be used to analyze such data. We aim to provide a good starting point for scientists with limited experience in computational analyses of DNA methylation and hydroxymethylation data, and recommend a few tools that are powerful, but still easy enough to use for their own data analysis. © The Author(s). 2019 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Upon the first publication of the fifth iteration of the Functional Annotation of Mammalian Genomes collaborative project, FANTOM5, we gathered a series of primary data and database systems into the FANTOM web resource (http://fantom.gsc.riken.jp) to facilitate researchers to explore transcriptional regulation and cellular states. In the course of the collaboration, primary data and analysis results have been expanded, and functionalities of the database systems enhanced. We believe that our data and web systems are invaluable resources, and we think the scientific community will benefit for this recent update to deepen their understanding of mammalian cellular organization. We introduce the contents of FANTOM5 here, report recent updates in the web resource and provide future perspectives. © The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/)

Immediate-early genes (IEGs) can be activated and transcribed within minutes after stimulation, without the need for de novo protein synthesis, and they are stimulated in response to both cell-extrinsic and cell-intrinsic signals. Extracellular signals are transduced from the cell surface, through receptors activating a chain of proteins in the cell, in particular extracellular-signal-regulated kinases (ERKs), mitogen-activated protein kinases (MAPKs) and members of the RhoA-actin pathway. These communicate through a signaling cascade by adding phosphate groups to neighboring proteins, and this will eventually activate and translocate TFs to the nucleus and thereby induce gene expression. The gene activation also involves proximal and distal enhancers that interact with promoters to simulate gene expression. The immediate-early genes have essential biological roles, in particular in stress response, like the immune system, and in differentiation. Therefore they also have important roles in various diseases, including cancer development. In this paper we summarize some recent advances on key aspects of the activation and regulation of immediate-early genes.

With the aim of identifying novel thermostable esterases, comprehensive sequence databases and cloned fosmid libraries of metagenomes derived from an offshore oil reservoir on the Norwegian Continental Shelf were screened for enzyme candidates using both sequence-and function-based screening. From several candidates identified in both approaches, one enzyme discovered by the functional approach was verified as a novel esterase and subjected to a deeper characterization. The enzyme was successfully over-produced in Escherichia coli and was shown to be thermostable up to 90°C, with the highest esterase activity on short-chain ester substrates and with tolerance to solvents and metal ions. The fact that the thermostable enzyme was solely found by functional screening of the oil reservoir metagenomes illustrates the importance of this approach as a complement to purely sequence-based screening, in which the enzyme candidate was not detected. In addition, this example indicates the large potential of deep-sub-surface oil reservoir metagenomes as a source of novel, thermostable enzymes of potential relevance for industrial applications.

Cytoscape is often used for visualization and analysis of metabolic pathways. For example, based on KEGG data, a reader for KEGG Markup Language (KGML) is used to load files into Cytoscape. However, although multiple genes can be responsible for the same reaction, the KGML-reader KEGGScape only presents the first listed gene in a network node for a given reaction. This can lead to incorrect interpretations of the pathways. Our new method, FunHoP, shows all possible genes in each node, making the pathways more complete. FunHoP collapses all genes in a node into one measurement using read counts from RNA-seq. Assuming that activity for an enzymatic reaction mainly depends upon the gene with the highest number of reads, and weighting the reads on gene length and ratio, a new expression value is calculated for the node as a whole. Differential expression at node level is then applied to the networks. Using prostate cancer as model, we integrate RNA-seq data from two patient cohorts with metabolism data from literature. Here we show that FunHoP gives more consistent pathways that are easier to interpret biologically. Code and documentation for running FunHoP can be found at https://github.com/kjerstirise/FunHoP.

The Norwegian e-Infrastructure for Life Sciences (NeLS) has been developed by ELIXIR Norway to provide its users with a system enabling data storage, sharing, and analysis in a project-oriented fashion. The system is available through easy-to-use web interfaces, including the Galaxy workbench for data analysis and workflow execution. Users confident with a command-line interface and programming may also access it through Secure Shell (SSH) and application programming interfaces (APIs). NeLS has been in production since 2015, with training and support provided by the help desk of ELIXIR Norway. Through collaboration with NorSeq, the national consortium for high-throughput sequencing, an integrated service is offered so that sequencing data generated in a research project is provided to the involved researchers through NeLS. Sensitive data, such as individual genomic sequencing data, are handled using the TSD (Services for Sensitive Data) platform provided by Sigma2 and the University of Oslo. NeLS integrates national e-infrastructure storage and computing resources, and is also integrated with the SEEK platform in order to store large data files produced by experiments described in SEEK. In this article, we outline the architecture of NeLS and discuss possible directions for further development Copyright: © 2018 Tekle KM et al. This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background The relationship between cholesterol and prostate cancer has been extensively studied for decades, where high levels of cellular cholesterol are generally associated with cancer progression and less favorable outcomes. However, the role of in vivo cellular cholesterol synthesis in this process is unclear, and data on the transcriptional activity of cholesterol synthesis pathway genes in tissue from prostate cancer patients are inconsistent. Methods A common problem with cancer tissue data from patient cohorts is the presence of heterogeneous tissue which confounds molecular analysis of the samples. In this study we present a general method to minimize systematic confounding from stroma tissue in any prostate cancer cohort comparing prostate cancer and normal samples. In particular we use samples assessed by histopathology to identify genes enriched and depleted in prostate stroma. These genes are then used to assess stroma content in tissue samples from other prostate cancer cohorts where no histopathology is available. Differential expression analysis is performed by comparing cancer and normal samples where the average stroma content has been balanced between the sample groups. In total we analyzed seven patient cohorts with prostate cancer consisting of 1713 prostate cancer and 230 normal tissue samples. Results When stroma confounding was minimized, differential gene expression analysis over all cohorts showed robust and consistent downregulation of nearly all genes in the cholesterol synthesis pathway. Additional Gene Ontology analysis also identified cholesterol synthesis as the most significantly altered metabolic pathway in prostate cancer at the transcriptional level. Conclusion The surprising observation that cholesterol synthesis genes are downregulated in prostate cancer is important for our understanding of how prostate cancer cells regulate cholesterol levels in vivo. Moreover, we show that tissue heterogeneity explains the lack of consistency in previous expression analysis of cholesterol synthesis genes in prostate cancer. © The Author(s). 2018 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( http://creativecommons.org/licenses/by/4.0/ ), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( http://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

We present model-based analysis for ChIA-PET (MACPET), which analyzes paired-end read sequences provided by ChIA-PET for finding binding sites of a protein of interest. MACPET uses information from both tags of each PET and searches for binding sites in a two-dimensional space, while taking into account different noise levels in different genomic regions. MACPET shows favorable results compared with MACS in terms of motif occurrence and spatial resolution. Furthermore, significant binding sites discovered by MACPET are involved in a higher number of significant three-dimensional interactions than those discovered by MACS. MACPET is freely available on Bioconductor. ChIA-PET; MACPET; Model-based clustering; Paired-end tags; Peak-calling algorithm. © The Author 2019. Published by Oxford University Press. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. Supplementary data

© 2018 Elsevier Ltd Oil biodegradation as a weathering process has been extensively investigated over the years, especially after the Deepwater Horizon blowout. In this study, we performed microcosm experiments at 5 °C with chemically dispersed oil in non-amended seawater. We link biodegradation processes with microbial community and metagenome dynamics and explain the succession based on substrate specialization. Reconstructed genomes and 16S rRNA gene analysis revealed that Bermanella and Zhongshania were the main contributors to initial n-alkane breakdown, while subsequent abundances of Colwellia and microorganisms closely related to Porticoccaceae were involved in secondary n‑alkane breakdown and beta‑oxidation. Cycloclasticus, Porticoccaceae and Spongiiabcteraceae were associated with degradation of mono- and poly-cyclic aromatics. Successional pattern of genes coding for hydrocarbon degrading enzymes at metagenome level, and reconstructed genomic content, revealed a high differentiation of bacteria involved in hydrocarbon biodegradation. A cooperation among oil degrading microorganisms is thus needed for the complete substrate transformation.

The relationship between cholesterol and prostate cancer has been extensively studied for decades, where high levels of cellular cholesterol are generally associated with cancer progression and less favorable outcomes. However, the role of in vivo cellular cholesterol synthesis in this process is unclear, and data on the transcriptional activity of cholesterol synthesis pathway genes in tissue from prostate cancer patients are inconsistent. A common problem with cancer tissue data from patient cohorts is the presence of heterogeneous tissue which confounds molecular analysis of the samples. In this study we present a method to minimize systematic confounding from stroma tissue in seven patient cohorts consisting of 1713 prostate cancer and 230 normal tissue samples. When confounding was minimized, differential gene expression analysis over all cohorts showed robust and consistent downregulation of nearly all genes in the cholesterol synthesis pathway. Additional analysis also identified cholesterol synthesis as the most significantly altered metabolic pathway in prostate cancer. This surprising observation is important for our understanding of how prostate cancer cells regulate cholesterol levels in vivo. Moreover, we show that tissue heterogeneity explains the lack of consistency in previous expression analysis of cholesterol synthesis genes in prostate cancer.

Maintenance of a genome requires DNA repair integrated with chromatin remodeling. We have analyzed six transcriptome data sets and one data set on translational regulation of known DNA repair and remodeling genes in synchronized human cells. These data are available through our new database: www.dnarepairgenes.com. Genes that have similar transcription profiles in at least two of our data sets generally agree well with known protein profiles. In brief, long patch base excision repair (BER) is enriched for S phase genes, whereas short patch BER uses genes essentially equally expressed in all cell cycle phases. Furthermore, most genes related to DNA mismatch repair, Fanconi anemia and homologous recombination have their highest expression in the S phase. In contrast, genes specific for direct repair, nucleotide excision repair, as well as non-homologous end joining do not show cell cycle-related expression. Cell cycle regulated chromatin remodeling genes were most frequently confined to G1/S and S. These include e.g. genes for chromatin assembly factor 1 (CAF-1) major subunits CHAF1A and CHAF1B; the putative helicases HELLS and ATAD2 that both co-activate E2F transcription factors central in G1/S-transition and recruit DNA repair and chromatin-modifying proteins and DNA double strand break repair proteins; and RAD54L and RAD54B involved in double strand break repair. TOP2A was consistently most highly expressed in G2, but also expressed in late S phase, supporting a role in regulating entry into mitosis. Translational regulation complements transcriptional regulation and appears to be a relatively common cell cycle regulatory mechanism for DNA repair genes. Our results identify cell cycle phases in which different pathways have highest activity, and demonstrate that periodically expressed genes in a pathway are frequently co-expressed. Furthermore, the data suggest that S phase expression and over-expression of some multifunctional chromatin remodeling proteins may set up feedback loops driving cancer cell proliferation.