<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Towards a Semantic Clinical Data Warehouse: A Case Study of Discovering Similar Genes</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Benedikt Kampgen</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Horst Werner</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Radwan Deeb</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Christof Bornhovd</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>FZI Research Center for Information Technology</institution>
          ,
          <addr-line>Karlsruhe</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>SAP AG</institution>
          ,
          <addr-line>Karlsruhe</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>Physicians nowadays have to consider a diverse range of data sources when treating a patient. Semantic clinical data warehouses allow to easily add new data and to pro-actively help the physician making sense of the data. In this work-in-progress paper we investigate an approach of using Linked Data as the access mechanism and a graph database for storage and query processing. We describe lessons learned from a case study of discovering similar genes where we use an existing similarity metric to derive new information, the Gene Ontology as a data source, and SAP HANA as an e cient graph database.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>In this work-in-progress paper we argue that overcoming such challenges
is possible using Linked Data, graph databases, and semantic algorithms
(Section 2): we describe a use case for discovering similar genes (Section 3) and derive
lessons learned (Section 4). We mention related work (Section 5) and conclude
(Section 6).
2</p>
    </sec>
    <sec id="sec-2">
      <title>Semantic Clinical Data Warehouse</title>
      <p>See Figure 1 for the architecture. Information in the semantic clinical data
warehouse is presented to the user by a visualisation and analysis tool. To store, query,
and visualise arbitrary information we use a graph database and the following
intuitive data model (property graph): Relevant objects such as patients,
interactions, and observations are represented as vertices in the graph. Such objects
have properties with values of primitive datatypes such as String and Integer,
e.g., the surname of a patient. Objects are related to each other via edges in the
graph, e.g., a patient is diagnosed with a disease. Such relationships also can have
properties, e.g., provenance information about the algorithm or human expert
that has generated the relationship. The integrator and reasoner component 1)
translates an RDF graph to a property graph, 2) derives implicit information
useful for data integration and decision support of users, and 3) imports the
graph to the graph database. The RDF graph is crawled based on the Linked
Data principles.</p>
      <p>HTTP Get
User HTTP Result</p>
      <p>Visualisation &amp;
Analysis Tool</p>
      <p>Graph Query</p>
      <p>Results</p>
      <p>Graph
Database Results</p>
      <p>Semantic Clinical Data Warehouse
Graph Query
Bulk Load IRneteags&amp;roanteorr RDF</p>
      <p>Crawler</p>
      <p>Data Sources
RDF Linked (Open) Data
RDF LD-Wrapper
RDF LD-Wrapper</p>
      <p>CSV
XML</p>
      <p>
        This architecture has the following advantages: Already, there are large
amounts of life science data { directly or using LD wrappers { published using
such widely-adopted access mechanisms and standard vocabularies [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. A graph
database is more schema- exible than a relational database, i.e., if new data
sources introduce new vertices, edges, and properties in the graph, no database
administrator has to modify the schema. Linked Data allows to easily add new
data sources to the data warehouse by following new links to further objects on
the Web. Implicit information can be derived in two ways: 1) by evaluating OWL
axioms represented in RDF; for instance, semantics from the OWL 2 RL pro le
such as equality can be evaluated using rule engines, and 2) by ML algorithms
that make use of ontological information, e.g., to discover similar genes. Also,
graph databases usually are designed to e ciently process analytical operations
over large graphs, i.e., can be used to e ciently compute and write-back results
from ML algorithms.
      </p>
    </sec>
    <sec id="sec-3">
      <title>Case Study of Discovering Similar Genes</title>
      <p>
        In this section, we apply our approach to a use case for discovering similar genes
from a plant [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. Similarity is an important basis for other relationships. For
instance, the e ect of a drug depends on the genes it targets. If drugs target
similar genes, they likely have similar e ects.
      </p>
      <p>Relevant data sources for our prototype { HANA Linked Data AnnSim
(HLA) { are descriptions of genes4, gene annotations from experts5, and the
Gene Ontology (GO) with a concept hierarchy6.</p>
      <p>Using OpenRe ne with RDF extension, we translate the former two sources
to RDF and reuse links from the GO RDF representation. Crawling such data
results in one RDF graph with genes, concepts, and annotations between genes
and concepts.</p>
      <p>
        HLA uses as a graph database HANA Graph, an extension to the HANA
in-memory database for storing and querying of property graphs [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. Graphs are
logically stored in HANA using two (virtual) tables: one table for vertices and one
table for edges each with columns for an id and every possible property. Graph
queries over HANA Graph are issued using the so-called GEM language and are
translated to SQL queries over the two tables. Based on a column-oriented and
in-memory database, HANA Graph allows fast query processing.
      </p>
      <p>An importer program then maps the crawled RDF graph to a property graph
and bulk loads the property graph to HANA Graph. Intuitively, the importer
generates for every triple two vertices for the subject and object (if not
existing), and an edge for the predicate. HANA Graph then contains genes (e.g.,
AT5G23810) and concepts (e.g., Amino Acid Transport) as vertices, and
relationships between genes and concepts as edges. For instance, there are annotation
relationships between genes and concepts as well as is-a relationships between
concepts. Vertices and edges can have properties, e.g., a concept has a textual
description. The graph is then extended with edges between genes describing
their similarity, and edges between concepts describing their distance in the is-a
concept hierarchy.</p>
      <p>
        Such information we compute based on an existing algorithm, AnnSim [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
AnnSim makes use of the distances between the concepts of two genes.
Intuitively, the shorter the average path between any two concepts of two genes the
more similar the two genes. Both computed distances and similarities are written
back to HANA Graph as edges between concepts and genes, respectively. Every
edge has as a property a numeric value between 0 and 1 for the similarity and
distance, and { in case several di erent algorithms are used { the name of the
algorithm.
4 ftp://ftp.arabidopsis.org/home/tair/Genes/TAIR10_genome_release/TAIR10_
functional_descriptions
5 ftp://ftp.arabidopsis.org/home/tair/Ontologies/Gene_Ontology/ATH_GO_
      </p>
      <p>GOSLIM.txt
6 http://purl.obolibrary.org/obo/go.owl</p>
      <p>As also visible in a screencast on our paper website7, the user of HLA gets
an overview of genes and their similarities to other genes; can zoom into
single genes to see textual descriptions of concepts; can visit concepts along the
concept hierarchy (Figure 2). Also, the user can ask for a graph view showing
similarities between genes based on distances between concepts (Figure 3). For
the visualisation, we used a visualisation engine called Symbiosis that can be
con gured with a JSON-based template language to visualise a graph.
Symbiosis uses HANA Graph and GEM for querying the data.</p>
      <p>Correct Computation of Similarities. There is a mean squared error
between the results of HLA and AnnSim of 0.09. This di erence we expect is
due to newer, possibly more elaborate versions of annotations and GO (version
1.2) used by HLA. We compared the results of both approaches with the gold
standard, a similarity metric based on the DNA sequence of genes (SeqSim).
7 http://people.aifb.kit.edu/bka/hla/
The mean squared error between HLA and SeqSim (0.19) is lower than between
AnnSim 1.0 and SeqSim (0.36), indicating that AnnSim similarities improve with
newer data sources; yet, further experiments are needed to con rm this claim.</p>
      <p>Scalable Computation of Similarities. Table 2 gives an overview of the
time for the di erent steps in the execution. Loading of data is estimated with
a connection of 6.7 Mbps download speed. Although HLA takes considerably
longer than AnnSim 1.0, we argue that HLA's bottlenecks can be resolved and
that HLA is more promising for larger datasets.</p>
      <p>AnnSim 1.0 uses a proprietary graph data format with reduced information
that is probably fast to generate (Prepare), download, and load. HLA uses a
more verbose but also more expressive graph model (RDF) and has to
generate (Prepare), download, transform to property graph (Map) and load 15 times
more vertices, 22 times more edges and comprehensive properties such as
textual descriptions. Loading graph data to HANA Graph showed fast and the
preprocessing steps we believe can be optimised by parallelisation.</p>
      <p>AnnSim 1.0 uses program logic in C/C++ over arrays to compute the 400
similarities and displays the results to users. HLA loads the relevant data to
a graph database and uses program logic in Java to issue database queries to
e ciently compute the similarities and to write back the results to the data
warehouse. The query language GEM was useful and intuitive for graph-traversal
queries. For instance, the following GEM read query is issued using a special-type
function WIPE() to the SQL interface of HANA, recursively visiting one or more
edges of type rdfs:subClassOf, and returns a vertex table with all ancestors
of a GO concept: RESULT uri:myResult FROM f GO:0005634 g-[@core:type
= 'rdfs:subClassOf']-&gt;(1,*);.</p>
      <p>The program logic in HLA spent more than 90% of the time to compute a
speci c part of AnnSim, a min-weight perfect matching (Blossom IV). We believe
we can optimise the Blossom IV execution, e.g., by running part of it directly in
HANA Graph via built-in and user-de ned functions. Writing back of the results
to the data warehouse took a lot of time since done using single write queries
instead of a bulk load. In this case, since read queries to HANA Graph showed
fast, HLA should also scale with larger datasets, in contrast to AnnSim 1.0 that
does not outsource bulk loading, reading, and writing to an external database.
Computing additional information can be done o ine. Interactive visualisation
over HANA Graph were possible using the Symbiosis engine.</p>
      <p>Flexible Computation and Visualisation of Similarities. Whereas
AnnSim 1.0 was implemented speci cally for the problem of e ciently computing
similarities of objects described in a proprietary format, HLA uses Linked Data
as a uni ed data model and standard access mechanism.</p>
      <p>
        New data sources can be added to HLA by providing more links to crawleable
Linked Data. We believe that e orts such as by Bio2RDF [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] to release life
science Linked Data will allow to semi-automatically resolve semantic con icts
using OWL semantics and rules.
      </p>
      <p>Other objects such as patients can be compared in HLA; AnnSim only
requires objects to be annotated with concepts and concepts to be described in an
is-a hierarchy. Algorithms that use other relationships and derive other
information can be added to HLA. The Symbiosis engine showed that { given su cient
understanding of the domain experts' problem { it is easily possible (5{10h of
manual work) to provide exible visualisations over a graph-based data model.</p>
    </sec>
    <sec id="sec-4">
      <title>Related Work</title>
      <p>
        According to Haussler et al. [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] a Million Genome Warehouse has to pro-actively
process relevant data in data analysis pipelines to draw valid and useful medical
inferences. HLA accesses the Gene Ontology and computes AnnSim [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ], yet,
can be extended with other biomedical ontologies and other semantic similarity
measures [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. HLA uses the HANA Graph in-memory database [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ] but may
also use other graph databases such as Graphium. Callahan and Dumontier [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]
present an approach to represent and evaluate scienti c hypotheses based on
RDF and SPARQL.
6
      </p>
    </sec>
    <sec id="sec-5">
      <title>Conclusions</title>
      <p>In this work-in-progress paper, in a small case study of discovering similar genes
we illustrated the potential of modular access mechanisms with Linked Data,
queries over a schema- exible graph database, and semantic algorithms to derive
new information. Continuously adding new data sources and data items, new
algorithms, and new visualisations leave exciting future work.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Callahan</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Cruz-Toledo</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ansell</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dumontier</surname>
            ,
            <given-names>M.:</given-names>
          </string-name>
          <article-title>Bio2RDF Release 2: Improved Coverage, Interoperability and Provenance of Life Science Linked Data</article-title>
          .
          <source>In: The Semantic Web: Semantics and Big Data</source>
          (
          <year>2013</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Callahan</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dumontier</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          :
          <article-title>Evaluating Scienti c Hypotheses Using the SPARQL Inferencing Notation</article-title>
          .
          <source>In: The Semantic Web: Research and Applications</source>
          (
          <year>2012</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Haussler</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Patterson</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Diekhans</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fox</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jordan</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Joseph</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ma</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Paten</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shenker</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sittler</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Stoika</surname>
            ,
            <given-names>I.: A Million</given-names>
          </string-name>
          <string-name>
            <surname>Cancer Genome Warehouse</surname>
          </string-name>
          .
          <source>Tech. rep.</source>
          , University of California at Berkeley (
          <year>2012</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Palma</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Vidal</surname>
            ,
            <given-names>M.E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Haag</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Raschid</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thor</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Measuring Relatedness Between Scienti c Entities in Annotation Datasets</article-title>
          . In: International Conference on Bioinformatics, Computational Biology and Biomedical
          <string-name>
            <surname>Informatics</surname>
          </string-name>
          (
          <year>2013</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Pesquita</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Faria</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Falco</surname>
            ,
            <given-names>A.O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lord</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Couto</surname>
            ,
            <given-names>F.M.</given-names>
          </string-name>
          :
          <article-title>Semantic Similarity in Biomedical Ontologies</article-title>
          .
          <source>PLOS Computational Biology</source>
          <volume>5</volume>
          (
          <issue>7</issue>
          ) (
          <year>2009</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Vasilyeva</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thiele</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          , Bornhovd,
          <string-name>
            <given-names>C.</given-names>
            ,
            <surname>Lehner</surname>
          </string-name>
          ,
          <string-name>
            <surname>W.</surname>
          </string-name>
          :
          <article-title>Leveraging Flexible Data Management with Graph Databases</article-title>
          .
          <source>First International Workshop on Graph Data Management Experiences and Systems</source>
          (
          <year>2013</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>