Querying phenotype-genotype relationships on patient datasets using semantic web technology: the example of cerebrotendinous xanthomatosis
- María Taboada†1Email author,
- Diego Martínez†2,
- Belén Pilo†3,
- Adriano Jiménez-Escrig†4,
- Peter N Robinson†5 and
- María J Sobrido†6
© Taboada et al.; licensee BioMed Central Ltd. 2012
Received: 7 December 2011
Accepted: 26 July 2012
Published: 31 July 2012
Semantic Web technology can considerably catalyze translational genetics and genomics research in medicine, where the interchange of information between basic research and clinical levels becomes crucial. This exchange involves mapping abstract phenotype descriptions from research resources, such as knowledge databases and catalogs, to unstructured datasets produced through experimental methods and clinical practice. This is especially true for the construction of mutation databases. This paper presents a way of harmonizing abstract phenotype descriptions with patient data from clinical practice, and querying this dataset about relationships between phenotypes and genetic variants, at different levels of abstraction.
Due to the current availability of ontological and terminological resources that have already reached some consensus in biomedicine, a reuse-based ontology engineering approach was followed. The proposed approach uses the Ontology Web Language (OWL) to represent the phenotype ontology and the patient model, the Semantic Web Rule Language (SWRL) to bridge the gap between phenotype descriptions and clinical data, and the Semantic Query Web Rule Language (SQWRL) to query relevant phenotype-genotype bidirectional relationships. The work tests the use of semantic web technology in the biomedical research domain named cerebrotendinous xanthomatosis (CTX), using a real dataset and ontologies.
A framework to query relevant phenotype-genotype bidirectional relationships is provided. Phenotype descriptions and patient data were harmonized by defining 28 Horn-like rules in terms of the OWL concepts. In total, 24 patterns of SWQRL queries were designed following the initial list of competency questions. As the approach is based on OWL, the semantic of the framework adapts the standard logical model of an open world assumption.
This work demonstrates how semantic web technologies can be used to support flexible representation and computational inference mechanisms required to query patient datasets at different levels of abstraction. The open world assumption is especially good for describing only partially known phenotype-genotype relationships, in a way that is easily extensible. In future, this type of approach could offer researchers a valuable resource to infer new data from patient data for statistical analysis in translational research. In conclusion, phenotype description formalization and mapping to clinical data are two key elements for interchanging knowledge between basic and clinical research.
Cerebrotendinous xanthomatosis (CTX) is a rare lipid-storage disease that leads to a complex combination of neurologic dysfunctions including cerebellar, pyramidal and extrapyramidal signs, neuropathy, dementia and psychiatric disturbances, as well as extra-neurological manifestations (chronic diarrhea, cataracts, tendon xanthomas, premature arteriosclerosis) . CTX is caused by mutations in the gene CYP27A1, which encodes the mitochondrial enzyme sterol 27-hydroxylase, deficiency of which causes an excess of intermediate metabolites such as cholestanol to accumulate in virtually every tissue. Like many neurodegenerative disorders, CTX is characterized by an insidious onset, progressive course and variable combination of clinical manifestations in each patient, which, together with the rarity of the disease, hampers correct and early diagnosis. Therapeutic delay is especially catastrophic in CTX, since there is a specific treatment (chenodeoxycolic acid), which is effective in reducing the plasma levels of cholestanol but has not been demonstrated to improve established neurological deficits. Mutation analysis of the CYP27A1 gene is a key step in the diagnosis of CTX and is routinely performed. The availability of comprehensive genotype to phenotype data sets will be crucial in order to promote early recognition and optimize the diagnostic process.
In any disease but most especially in rare diseases the possibility of accessing detailed patient datasets from research and clinical studies, including genetic variants and phenotypic manifestations, would significantly improve diagnosis and treatment. Electronic patient records are able to gather diverse types and growing amounts of phenotypic data, while the use of genome-scale DNA sequencing techniques allows the collection of an increasing number of genetic variants per individual. Thus, integrating complex phenotype descriptions with genetic testing records has become one of the main challenges of biomedicine . As the number of openly accessible datasets continues to rise, the integration of research repositories and patient clinical data will be more viable. However, bioinformatics tools are needed to help explore complex genotype-phenotype relationships. Geneticists would request software tools able to retrieve and analyze the data produced in diverse clinical settings and associated to a new given genetic variant; that is, answering questions like what are the phenotype traits that have been identified in patients with this genetic variation? Clinicians, on the other hand, would see their work greatly facilitated by being able to answer queries like what genes or genetic variants are associated with this particular combination of observable features?
The development of locus-specific mutation databases (LSDBs) and tools to build them such as the Leiden Open Variation Database (LOVD) , and the Universal Mutation Database (UMD)  started to pave the way to solve the problem of collecting genetic datasets produced by diverse experimental methods in different laboratories. However, the phenotype description in most LSDBs is very scarce. The Human Variome Project (HVP)  is an international initiative aiming ultimately at the worldwide collection and harmonization of all human genetic variations and associated phenotypic data. The GEN2PHEN project also represents an international attempt to undertake the logistical and technical challenges to join disparate genotype-phenotype resources in a shared mode . In order to achieve that goal, communication standards are needed to allow interoperability between clinical and genetic datasets. Standards to represent genetic findings are already available, such as those produced by the HUGO committee (http://www.genenames.org/aboutHGNC.html), gene relationships provided by Gene Ontology  or the nomenclature for description of sequence variants proposed by the Human Genome Variation Society (HGVS, http://www.hgvs.org/mutnomen). However, such a level of consensus on the best descriptors for phenotypic information is far more complex and has not been reached in clinical medicine.
Although the term phenotype covers an extensive range of information varying from molecular to organism level observable characteristics , in this work phenotype is meant only as any observable human trait, such as an anatomical abnormality (e.g., juvenile cataracts) or a clinical feature ( e.g., tendon xanthomas). Currently, the most useful catalog of human Mendelian disorders is OMIM, the Online Mendelian Inheritance in Man , a text-based knowledge source of human phenotypes and related genes. OMIM describes phenotypes using narrative sentences (e.g. normal to slightly elevated plasma cholesterol). Although these textual descriptions are highly expressive, capturing phenotype information using free-text fields in databases hampers computational processing and inference . The use of a standard terminology provides a more appropriate method of expressing unambiguous, computable, and interoperable phenotype descriptions. Standard terminologies organize the concepts of a particular domain into a taxonomy (e.g., epilepsy is a seizure disorder, which is an abnormality of the central nervous system), assigning them identifiers which do not change with new versions. They also address the issues of different synonyms for the same concept (e.g., convulsions vs. epileptic seizures).
Patient data from clinical and research settings are usually stored in different formats, from simple spreadsheets to relational databases, being extremely difficult to integrate genotype-phenotype data across multiple formats. The semantic web technology provides an adequate instrument of recording phenotypes in a standardized fashion and with a high degree of expressivity. Using this technology to represent data will ensure the compatibility of them with the future knowledge and data resources. Additionally, one of the main challenges of articulating queries on phenotype-genotype relationships is discrepancy in the level of abstraction between phenotype descriptions and patient clinical data. The semantic web technology provides a layer of abstraction that makes it simple to use. Moreover, this technology is based on open world assumption: everything we do not know is undefined. Hence, unknown relationships will be interpreted as not computable instead of false. This approach naturally deals with incomplete information, which is very usual in biomedicine, and it is able to refine knowledge when new information comes along.
One option to deal with the phenotype complexity can be to define a minimum set of phenotype template fields . In contrast, an ontology-based technology would provide a more open and flexible representation mechanism [10, 12–14], thus facilitating the continuing incorporation and interpretation of new phenotype characteristics. An ontology is a data model that represents a set of entities in some domain and the relationships among those entities. One of the benefits of using ontologies is the potential to apply reasoners (logical inference tools), which can infer new data to subsequently facilitate query answering and statistical analysis. In the present work, we used patient data from a specific rare genetic disease (CTX) to formally represent phenotype descriptions using the ontological paradigm . We then engineered the patient data in an ontology-based patient model and finally executed queries on genotype-phenotype relationships with a Semantic Web approach.
Overview of the query process
Additionally, queries such as ‘What genetic variants have been identified in patients without childhood-onset chronic diarrhea?’ will return those genetic variants for patients with asserted absence of childhood-onset chronic diarrhea, that is, the genetic variant p.R395C (patient 564385–1) in Figure 1. But, it will not include those patients for which there is unknown or missing information (open-world assumption); that is, it will say “Do not know” for the genetic variants c.844 + 1 G- > T y p.Q525X, unless a blank is interpreted as ‘not (childhood-onset chronic diarrhea)’.
Data and knowledge sources
Patient Datasets: The complete and detailed patient dataset published on a PhD dissertation written in Spanish  was used in our approach. Currently, these data are also available in English in a more summarized form . These patient data originated from a thorough collaborative clinical and genetic study on CTX carried out at the Hospital Ramón y Cajal (Madrid) and the Fundación Pública Galega de Medicina Xenómica (Santiago de Compostela) in Spain.
In brief, 25 patients from 19 families were thoroughly studied through personal examination of the patients and their medical records by the authors (BP, AJE, MJS). Detailed clinical history data, neurological signs, neurophysiologic, biochemical and neuroimaging data were collected in extensive table and text format.
Text-based resources: Two text-based resources from highly reliable web sites were searched in our approach, the OMIM clinical synopsis (http://www.ncbi.nlm.nih.gov/omim/213700?dopt=Synopsis; ) and the Summary, Diagnosis and Clinical Description sections of GeneReviews (http://www.ncbi.nlm.nih.gov/bookshelf/br.fcgi?book=gene&part=ctx), both of them containing known phenotype manifestations of CTX.
Terminology systems: The Unified Medical Language System (UMLS¸ http://umlsks.nlm.nih.gov ) Metathesaurus was chosen to codify the relevant CTX terminology extracted from the text-based resources.
Domain ontologies: Two web sites actively used in biomedical communities, the NCBO Bioportal (http://bioportal.bioontology.org; ) and the OBO Foundry (http://www.obofoundry.org; ), were accessed to search ontologies and terminologies in the specific domain of CTX.
CTX phenotype-genotype relationships are represented in the Semantic Web Layer via a patient data model expressed in terms of a CTX phenotype ontology. We developed this ontology following the reuse-based methodology NeOn , thanks to the current availability of ontological and terminological resources that have reached some consensus in biomedicine. Although NeOn is focused on collaboratively building ontology networks, it is currently the only one covering complex scenarios with reuse of ontological and non-ontological resources. Initially, a set of use cases describing phenotype-genotype bidirectional relationships and query support were specified in order to determine the content of our ontology. Thus, the ontology should represent CTX phenotypes, genetic variants, and bidirectional relationships between them. Next, we developed the CTX ontology carrying out several activities proposed by the NeOn methodology: extracting the CTX terminology, assessing the publicly available resources, reusing and reengineering the selected resources. The ontology was developed using the Web Ontology Language (OWL)  and Protégé-OWL 3.4 .
Extracting the CTX terminology
One of the critical steps when developing ontologies is to identify the particular knowledge to be represented in the ontology. Usually, this activity is carried out at the beginning of the ontology project through a set of interviews with domain experts. Since this activity is labor-intensive, we decided to semi-automatically extract the glossary of CTX terms from the different resources employed in this work. Then, we manually checked the extracted terminology, and we eliminated erroneously extracted concepts, very general concepts (e.g., Findings) and some terms that did not represent concepts relevant to CTX.
Assessing terminological and ontological resources
We searched the term Cerebrotendinous Xanthomatosis in the NCBO Bioportal and in the OBO Foundry resources. Next, we assessed the candidate ontologies from these resources by calculating their coverage with respect to the glossary of CTX terms.
Reusing terminological and ontological resources
We built the CTX ontology by reusing concepts and relationships from the selected resources in the last step.
Construction of patient data model in OWL
We developed an OWL model describing the patient dataset and the relationships between phenotypic and genotypic data, where the former are expressed in terms of the CTX phenotype ontology.
Querying phenotype-genotype bidirectional relationships
In order to have the ability to execute queries at different levels of abstraction, firstly a set of rules was defined to successfully interpret the meaning of phenotype descriptions (step 1 in Figure 1). This set of rules was implemented by means of the Semantic Web Rule Language (SWRL) , a language to express Horn-like rules in terms of OWL concepts. Next, in order to formulate queries traversing bidirectional phenotype-genotype relationships on patient data (step 2 in Figure 1), we used SQWRL , a language that expands SWRL in order to provide querying of OWL ontologies.
The CTX ontology
Extracting the CTX terminology
Percentage of match of CTX terminology extracted from different sources
Patient Data (ExactMatch)
Patient Data (ExactMatch + MetaMap)
Assessing terminological and ontological resources
Coverage of relevant CTX terminology in Snomed CT, HPO and PATO
Disorders & Abnormalities (Concept #)
Anatomical Structures (Concept #)
Diagnostic Studies (Concept #)
Qualifier Values (Concept #)
Reusing terminological and ontological resources
Although the coverage of relevant CTX terminology in Snomed CT was slightly higher than in HPO (Table 2), we decided to reuse HPO and extend it with the required Snomed CT concepts and relationships. The justification is the following: the organization of HPO follows the structure of OMIM, which is very close to experts in CTX, HPO is less complex than Snomed CT and it can be easily translated to OWL.
Descriptive statistics on the origin of CTX ontology concepts
HPO Concept #
SNOMED CT Concept #
New Concept #
In addition, it was necessary to distinguish among the presence and the absence of a phenotypic trait in a given patient. Therefore, we subdivided each phenotype into two OWL instances covering the two possibilities. For example, for the phenotype Epilepsy, there are two possible instances: presence of epilepsy and absence of epilepsy. Hence, there are three types of patients: those ones with asserted presence of epilepsy, those ones with asserted absence of epilepsy, and finally those ones for which there is no information about the presence or absence of epilepsy.
The patient data model
We created object properties to represent patient manifestations: gene mutations, extra-neurological symptoms, neurological symptoms, results from diagnostic studies and other manifestations. Then, the patient model was populated manually. In total, 1 class, 17 properties and 25 individuals represent the Patient Model in our ontology.
Querying phenotype-genotype bidirectional relationships
Before formulating queries traversing phenotype-genotype relationships, we defined a set of 28 rules to deduce abstract phenotype descriptions from the patient data. Supplementary data provides all rules defined in our system. The execution of the rules generates the corresponding links (e.g., has extra-neurological manifestations’) between the corresponding patient instances and the abstract phenotype instances. Next, we formulate some of the queries defined in the initial list of competency questions, showing how the approach can help us to analyze the CTX dataset. (Supplementary data provides the complete list of queries, showing the corresponding SQWRL implementation and the answers in every case).
Example 1: In the CTX study, the presence of childhood-onset chronic diarrhea was confirmed in 44 % of patients. The reply to the query ‘What genetic variants are associated with childhood-onset chronic diarrhea?’ substantiates that only 7 from 14 genetic variants (50 %) are associated with this manifestation. Moreover, the answer to ‘What are the genetic variants that have always been associated with childhood-onset chronic diarrhea?’ verifies that only one genetic variant (p.R395C) is associated with both the presence and absence of the manifestation, whereas the remaining 6 genetic variants are always associated with the presence of this manifestation. Our approach also provides queries like ‘What are the genetic variants that have never been associated with childhood-onset chronic diarrhea?. In the CTX study, only the generic variant p.Q525X was never associated.
Four examples of queries to the patient data
Pattern of CQ
Examples of CQ
Example of SQWRL
Example of answers
What are the genetic variants that have been associated with a combination of traits?
What are the genetic variants that have been associated with epilepsy and dementia?
Patient(?p) ^ hasNervousSystemDisorder(?p, ?x) ^ Epilepsy(?x) ^ hasPresence(?x, ?y) ^ Yes(?y) ^ hasNervousSystemDisorder(?p, ?d) ^ Dementia(?d) ^ hasPresence(?d, ?y) ^ hasGeneMutation(?p, ?g) ^ GeneticMutation(?g) ^mutation(?g, ?m) → sqwrl:columnNames("GeneMutation") ^ sqwrl:selectDistinct(?m)
c.844 + 1 G- > T
What are the abnormalities that have been associated with a specific genetic variant?
What are the Abnormalities of the Central Nervous System that have been associated with p.R395C?
Patient(?p1) ^ hasNervousSystemDisorder(?p1, ?z) ^ hasPresence(?z, ?y) ^ Yes(?y) ^ AbnormalityoftheCerebellum(?z) ^ hasGeneMutation(?p1, ?g) ˚ sqwrl:makeSet(?s1, ?z) ^ Patient(?p2) ^ hasOtherManifestations(?p2, ?x) ^ hasPresence(?x, ?y) ^ AbnormalityoftheCerebellum(?x) ^ hasGeneMutation(?p2, ?g) ^ sqwrl:makeSet(?s2, ?x) ^ GeneMutation(?g) ^ mutation(?g, ?m) ^ swrlb:equal(?m, "p.R395C") ˚ sqwrl:append(?s3, ?s1, ?s2) ^ sqwrl:element(?e, ?s3) → sqwrl:select(?e)
How often has a specific genetic variant been associated with a specific trait?
How often has p.R395C been associated with Ataxia?
Patient(?p1) ^ hasNervousSystemDisorder(?p1, ?x) ^ Ataxia(?x) ^ hasPresence(?x, ?y) ^ Yes(?y) ^ hasGeneMutation(?p1, ?g) ^ GeneMutation(?g) ^ mutation(?g, ?m) ^ swrlb:equal(?m, "p.R395C") ˚ sqwrl:makeSet(?s1, ?p1) ˚ sqwrl:size(?size1, ?s1) ^ Patient(?p2) ^ hasGeneMutation(?p2, ?g) ^ sqwrl:makeSet(?s2, ?p2) ^ sqwrl:size(?size2, ?s2) ^ swrlb:multiply(?mu, ?size1, 100.0) ^ swrlb:divide(?d, ?mu, ?size2) → sqwrl:select(?d)
What is the average number of years from the onset of a symptom to the onset of another symptom (or diagnosis/death) in patients with a given genetic variant?
What is the average number of years from the onset of diarrhea to the first neurological symptom in patients with the genetic variant p.R395C?
Patient(?p1) ^ hasGeneMutation(?p1, ?g) ^ GeneticMutation(?g) ^ mutation(?g, "p.R395C") ^ hasDiarrheaAge(?p1, ?d) ^ AgeatFirstSymptom(?d) ^ age(?d, ?da) ^ hasNeurologicalSymptomsOnsetAge(?p, ?a) ^ AgeatFirstSymptom(?a) ^ age(?a, ?ca) ^ swrlb:subtract(?di, ?ca, ?da) → sqwrl:columnNames("Average age from diarrhea to neurological symptoms onset") ^ sqwrl:avg(?di)
Example 3: The study compiled the set of main neurological manifestations and their frequency. The request ‘What are the central nervous system manifestations that have been identified in patients with p.R395C?’ supplies us with the specific manifestations (Ataxia and Chiari Type I) for a particular genetic variant (see Table 4).
Example 4: The CTX study linked the onset age of neurological symptoms with the presence of xanthomas. In patients with xanthoma s, the average onset age of neurological symptoms was 23 years in contrast with 16 years for patients without xanthomas. Inquiring about What is the average onset age of nervous system symptoms in patients with xanthomas and p.R405W? returns 31 years in contrast with 23,4 years (almost 8 years before) when the same query is made for the genetic variant p.R395C. Similarly, the query What is the average onset age of nervous system symptoms in patients without xanthomas and with p.R395C? confirms 12,7 years (more than 3 years before the general average of patients). The system provides a specific query to ask this type of questions directly: How long after the onset of xanthomas the first neurological symptom appears in patients with the genetic variant p.R395C compared to all patients?
Table 4 shows two other examples of queries to the patient data. Supplementary data provide the complete set of 24 patterns of SWQRL queries, which were designed following the competency questions. The results obtained for all queries in the ontology were tested with manual analysis during implementation. If a query did not return the expected result (from manual analysis), we assumed that the query was formulated wrongly and so, it was changed until to get the manual result.
Our approach shows how the ontological paradigm and the semantic web languages OWL, SWRL and SQWRL can be used in combination to develop tools to explore phenotype-genotype bidirectional associations in the particular clinical domain of CTX, a rare, autosomal-recessive neurometabolic disease. Two quite similar approaches were previously demonstrated in familial hypercholesterolemia  and autism cases . In the first one, an ontology was used to guide the expert with the choice of meaningful subsets of a large mass of genomic and post-genomic data. In the second one, a specific domain ontology was developed and complemented with an information model and a rule-based inference function, to infer high-level phenotypic abstractions. An interesting alternative to define some of these rules is the use of OWL 2, a new version of OWL, which considerably improves the OWL data types. Using OWL2, concepts such as Childhood-onset chronic diarrhea will be directly expressible with OWL2 using data type restrictions to define ranges for the concepts. Another option to define some SQWRL queries is to create defined OWL classes (through asserting necessary and sufficient conditions) and to use the reasoner to automatically compute the inferred types.
Although the present work is limited to a very specific domain and with data collected from only a limited number of patients through a Spanish research study, it shows the potential use in other rare diseases and larger datasets. One of the main problems for the correct diagnosis and handling of rare genetic disorders is the difficulty in recognizing their variable clinical expression, as well as the early symptoms of the disease. We propose here an example of the type of tools that can aid clinical practice by querying about potential genetic causes of a specific symptom combination. The power of this type of tools will greatly depend on both the amount and quality of data that nurture the system. Thus, two sets of advancements are extremely necessary in order to improve diagnosis and treatment of genetic diseases, especially for rare diseases: 1) Promoting international initiatives to gather large sets of accessible patient data and 2) Developing phenotype ontologies and bioinformatics tools, in order to query phenotype-genotype relationships. Achieving this necessary level of detailed genotype and phenotype description while handling patient datasets appropriately to protect individual confidentiality, is a question that is being addressed by the Human Variome Project consortium ethics committee .
The approach proposed here is a pilot model that must be taken with due caution, especially with small datasets. Variable penetrance and expressivity are common in neurogenetic disorders and phenotype trait description of a given patient will evolve with time. Thus, phenotype-genotype correlations are certainly not absolute. That is, we cannot assert that a patient with a specific mutation will necessarily have the same phenotype as another patient stored in the database with the same mutation. Still, if used with awareness of its limitations, this type of genotype-phenotype exploration will be of clinical utility. Furthermore, if the database is big enough it might be interesting to derive risk figures, although this is not at all easy to do directly in OWL/SQWRL.
Although the procedure used to populate the CTX ontology in this work was manual, this is not the most suitable method to incorporate new patients to the patient model. Managing this step is important for further use of the ontology or work inspired in this approach. Incorporating an ontology-based workflow to annotate patient data resources automatically, as proposed by , will provide the way of combining the new resources with the existing CTX ontology. Additionally, temporal dimensions to phenotypic data are important for this clinical domain and competency questions. Modeling of timelines and temporal occurrence of events is one of the most challenging problems in ontology development. The patient data used in this study only included temporal events referred to the patient age when a given manifestation had taken place (e.g., age onset of cataract) or temporal intervals referred to the duration of some manifestations. Even with this limited information, some relevant queries on temporal aspects of phenotype-genotype relationships could be designed. For example, ‘On average, how long was it from the onset of chronic diarrhea to the first neurological symptom in patients with CTX?’, ‘and in patients with CTX and the mutation p.395 C?’ Another aspect to be improved in the future is to explore the possibility of designing more complex SWRL rules dealing with probabilistic classification in order to deduce abstract phenotypes.
Our approach could also have limitations if the number of patient data increases extremely. Scalability issues or the slowness of OWL stores and reasoners are some of the disadvantages of using such a young technology . Even so, there are some ways of speeding up searches in OWL. One interesting way is the use of the OWL 2 EL (i.e. existential logic) subsets  to enable more efficient reasoning and complex queries to support scientific analyses. Although OWL EL is focused on large-scale ontologies, exploring the OWL EL expressivity power, which is more limited than OWL, for representing and processing our approach, could be an interesting area for future research. Another current limitation to the use of this technology is the scarcity of publicly available data sources. Although we cannot anticipate the format of data sources in the future, even in the case they were available in traditional databases and not in OWL, this technology could be used by implementing a bridge between the relational patient databases and the phenotype ontology.
Currently two communities are making a big effort to supply a full and comprehensive representation of the clinical domain: The International Health Terminology Standards Development Organization (IHTSDO), with the Snomed CT common vocabulary, and the Open Biomedical Ontologies (OBO) Foundry initiative , providing an ontology repository covering different domains of biomedicine. A recent study analyzing both communities concluded that it is premature to know whether one or the other will supply the solution to the breadth of coverage . Thus, instead of choosing between these two strategies to represent the particular CTX domain we selected the ontology HPO from the OBO Foundry initiative, and extended it with Snomed CT from IHTSDO and a limited number of newly created terms. It should be emphasized that the goal of this work was not to build an exhaustive and definitive ontology for CTX, but to develop semantic web tools to query phenotype-genotype relationships. An extra benefit of this approach is that the ontology can then be reused for other applications, such as differential diagnosis .
Our study shows that the Semantic Web paradigm provides the technology required to represent phenotype-genotype relationships in diseases with complex and variable manifestations such as CTX. We developed a tool to query this type of relationships on patient data in both directions and at different levels of abstraction. While ontologies have been generally used in medicine to describe unambiguous and standard terminologies agreed by consensus, our approach makes use of the ontological paradigm and semantic web technologies to provide a structured framework to query about individuals presenting a combination of phenotype traits or carrying specific genetic variants. The proposed approach implies querying the patient data by designing a patient data model in OWL and accessing them via the phenotype ontology. A semantic web rule language allowed us to infer phenotypic abstractions from patient data and provided the required bridge between phenotypic abstractions and clinical data. Through the use of this language, queries about phenotype-genotype relationships can now be formulated on the abstraction level that is common in genetic databases. We believe that our strategy is a promising approach for translational medical research, which will help improve diagnosis and thus early and effective treatment of genetic disorders. This is especially true for rare diseases, where the number of affected individuals is small and therefore easy data access and query is essential to the health care community. Coordinated international initiatives such as the HVP are crucial to promote the development of the necessary tools as well as to provide openly accessible patient data.
Availability and requirements
The CTX ontology in OWL and all information relevant to the paper is provided in following link: http://www.usc.es/keam/CTX/TheCTXOntology.html, which contains
The complete list of Competency Questions (QC).
Several diagrams showing the main hierarchies of the CTX phenotype management ontology.
The complete list of SWRL rules implemented in our approach.
The complete list of patterns of queries implemented in our approach. These groups of queries were designed following the list of CQ.
The patient model and an example of a fictitious patient.
The work presented in this paper has been developed in the funded national project Gestión de Terminologías Médicas para Arquetipos (TIN2009-14159-C05-05) by the Ministerio de Educación y Ciencia. This work was partly supported by the network REGENPSI (2009/019) from the Program of Consolidation and Structure of Competitive Units, Consellería de Educación e Ordenación Universitaria, Xunta de Galicia, and by FEDER funds for regional development. PNR was supported by a grant from the Deutsche Forschungsgemeinschaft (DFGRO2005/4-2).
- Federico A, Dotti MT: Cerebrotendinous Xanthomatosis: Clinical Manifestations, Diagnostic Criteria, Pathogenesis, and Therapy. J Child Neurol. 2003, 18: 633-638. 10.1177/08830738030180091001.View ArticlePubMedGoogle Scholar
- Lindblom A, Robinson PN: Bioinformatics for Human Genetics: Promises and Challenges. Hum Mutat. 2011, 32: 495-500. 10.1002/humu.21468.View ArticlePubMedGoogle Scholar
- Fokkema I, den Dunnen J, Taschner P: LOVD: easy creation of a locus specific sequence variation database using an LSDB-in-a-box approach. Hum Mutat. 2005, 26: 63-68. 10.1002/humu.20201.View ArticlePubMedGoogle Scholar
- Beroud C, Collod-Beroud G, Boileau C, Soussi T, Junien C: UMD (Universal Mutation Database): a generic software to build and analyze locus-specific databases. Hum Mutat. 2000, 15: 86-94. 10.1002/(SICI)1098-1004(200001)15:1<86::AID-HUMU16>3.0.CO;2-4.View ArticlePubMedGoogle Scholar
- Kaput J, Cotton R, Hardman L, Watson M, Al Aqeel AI, Al-Aama JY, Al-Mulla F, Alonso S, Aretz S, Auerbach AD, Bapat B, Bernstein IT, Bhak J, Bleoo SL, Blöcker H, Brenner SE, Burn J, Bustamante M, Calzone R, Cambon-Thomsen A, Cargill M, Carrera P, Cavedon L, Cho YS, Chung YJ, Claustres M, Cutting G, Dalgleish R, den Dunnen JT, Díaz C, Dobrowolski S, Dos Santos MR, Ekong R, Flanagan SB, Flicek P, Furukawa Y, Genuardi M, Ghang H, Golubenko MV, Greenblatt MS, Hamosh A, Hancock JM, Hardison R, Harrison TM, Hoffmann R, Horaitis R, Howard HJ, Barash CI, Izagirre N, Jung J, Kojima T, Laradi S, Lee YS, Lee JY, Gil-da-Silva-Lopes VL, Macrae FA, Maglott D, Marafie MJ, Marsh SG, Matsubara Y, Messiaen LM, Möslein G, Netea MG, Norton ML, Oefner PJ, Oetting WS, O'Leary JC, de Ramirez AM, Paalman MH, Parboosingh J, Patrinos GP, Perozzi G, Phillips IR, Povey S, Prasad S, Qi M, Quin DJ, Ramesar RS, Richards CS, Savige J, Scheible DG, Scott RJ, Seminara D, Shephard EA, Sijmons RH, Smith TD, Sobrido MJ, Tanaka T, Tavtigian SV, Taylor GR, Teague J, Töpel T, Ullman-Cullere M, Utsunomiya J, van Kranen HJ, Vihinen M, Webb E, Weber TK, Yeager M, Yeom YI, Yim SH, Yoo HS: Contributors to the Human Variome Project Planning Meeting: Planning the Human Variome Project: the Spain report. Hum Mutat. 2009, 30: 496-510. 10.1002/humu.20972.View ArticlePubMedGoogle Scholar
- Webb AJ, Thorisson GA, Brookes AJ: An informatics project and online “Knowledge Centre” supporting modern genotype-to-phenotype research. Hum Mutat. 2011, 32: 543-550. 10.1002/humu.21469.View ArticlePubMedGoogle Scholar
- Bada M, Stevens R, Goble C, Gil Y, Ashburner M, Blake JA, Cherry JM, Harris M, Lewis S: A short study on the success of the Gene Ontology. Web Semantics. 2004, 1: 235-240. 10.1016/j.websem.2003.12.003.View ArticleGoogle Scholar
- Schulze TG, McMahon FJ: Defining the Phenotype in Human Genetic Studies: Forward Genetics and Reverse Phenotyping. Hum Hered. 2004, 58: 131-138. 10.1159/000083539.View ArticlePubMedGoogle Scholar
- McKusick VA: Mendelian Inheritance in Man and Its Online Version, OMIM. Am J Hum Genet. 2007, 80: 588-604. 10.1086/514346.View ArticlePubMedPubMed CentralGoogle Scholar
- Mungall C, Gkoutos G, Smith C, Haendel MA, Lewis SE, Ashburner M: Integrating phenotype ontologies across multiple species. Genome Biol. 2010, 11: R9-10.1186/gb-2010-11-1-r9.View ArticleGoogle Scholar
- Cheung K, Frost HR, Marshall MS: Prud'hommeaux E, Samwald M, Zhao J. Paschke A: A journey to Semantic Web query federation in the life sciences. BMC Bioinformatics. 2009, 10: S10-PubMedGoogle Scholar
- Hoehndorf R, Loebe F, Kelso J, Herre H: Representing default knowledge in biomedical ontologies: Application to the integration of anatomy and phenotype ontologies. BMC Bioinformatics. 2007, 8: 377-10.1186/1471-2105-8-377.View ArticlePubMedPubMed CentralGoogle Scholar
- Ruttenberg A, Clark T, Bug W, Samwald M, Bodenreider O, Chen H, Doherty D, Forsberg K, Gao Y, Kashyap V, Kinoshita J, Luciano J, Marshall MS, Ogbuji C, Rees J, Stephens S, Wong GT, Wu E, Zaccagnini D, Hongsermeier T, Neumann E, Herman I, Cheung KH: Advancing translational research with the semantic web. BMC Bioinformatics. 2007, 8: S2-View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson P, Köhler S, Bauer S, Seelow D: HornD, Mundlos S: The human phenotype ontology: A tool for annotating and analyzing human hereditary disease. Am J Hum Genet. 2008, 83: 610-15. 10.1016/j.ajhg.2008.09.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Gómez-Pérez A, Fernández-López M, Corcho O: Ontological Engineering. 2004, Berlin, Springer VerlagGoogle Scholar
- Gómez-Pérez A, Suárez-Figueroa MC: Scenarios for building ontology networks within the NeOn methodology. In Proceedings of the K-CAP: 1–4: Redondo Beach. New York: ACM. September 2009, 2009: 183-184.Google Scholar
- Pilo B: Xantomatosis Cerebrotendinosa en España: mutaciones, aspectos clínicos y terapéuticos. 2009, PhD thesis, University of Alcalá de Henares, Faculty of MedicineGoogle Scholar
- Pilo B, Jimenez-Escrig A, Lorenzo JR, Pardo J, Arias M, Ares-Luque A, Duarte J, Muñiz-Pérez S, Sobrido MJ: Cerebrotendinous xanthomatosis in Spain: clinical, prognostic, and genetic survey. Eur J of Neurol. 2011, 18 (10): 1203-1211. 10.1111/j.1468-1331.2011.03439.x.View ArticleGoogle Scholar
- Cohen R, Gefen A, Elhadad M, Birk OS: CSI-OMIM–Clinical Synopsis Search in OMIM. BMC Bioinformatics. 2011, 12: 65-10.1186/1471-2105-12-65.View ArticlePubMedPubMed CentralGoogle Scholar
- Noy N, Shah N, Whetzel P, Dai B, Dorf M, Griffith N, Jonquet C, Rubin D, Storey M, Chute C, Musen M: BioPortal: ontologies and integrated data resources at the click of a mouse. Nucleic Acids Res. 2009, 37: W170-W173. 10.1093/nar/gkp440.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith B, Ashburner M, Rosse C, Bard J, Bug W, Ceusters W, Goldberg LJ, Eilbeck K, Ireland A: Mungall CJ; OBI Consortium, Leontis N, Rocca-Serra P, Ruttenberg A, Sansone SA, Scheuermann RH, Shah N, Whetzel PL, Lewis S: The OBO foundry: Coordinated evolution of ontologies to support biomedical data integration. Nat Biotechnol. 2007, 25: 1251-1255. 10.1038/nbt1346.View ArticlePubMedPubMed CentralGoogle Scholar
- McGuinness D, van Harmelen F: OWL web ontology language overview. 2012, [http://www.w3.org/TR/owl-features/]Google Scholar
- Knublauch H, Fergerson R, Noy N, Musen M: The Protege OWL plugin: an open development environment for Semantic Web applications. In Proceedings of the Third ISWC: 7–11: Hiroshima: Edited by McIlraith S, Plexousakis D, van Harmelen F. Hiroshima: Springer, LNCS. November 2004, 3298: 229-243.Google Scholar
- O'Connor M, Knublauch H, Tu SW, Musen MA: In Proceedings of the 8th International Protégé Conference. 2005, Madrid, Protégé with Rules WorkshopGoogle Scholar
- O'Connor M, Das A: SQWRL: a Query Language for OWL. 2009, In Proceedings of the Fifth International Workshop on OWL, Chantilly, 23-24.Google Scholar
- Aronson A, Lang F: An overview of MetaMap: historical perspective and recent advances. JAMIA. 2010, 17: 229-236.PubMedPubMed CentralGoogle Scholar
- Robinson PN, Mundlos S: The Human Phenotype Ontology. Clin Genet. 2010, 77: 525-534. 10.1111/j.1399-0004.2010.01436.x.View ArticlePubMedGoogle Scholar
- Coulet A, Smaïl-Tabbone M, Benlian P, Napoli A, Devignes MD: Ontology-guided data preparation for discovering genotype-phenotype relationships. BMC Bioinformatics. 2008, 9 (4): S3-10.1186/1471-2105-9-S4-S3.View ArticlePubMedPubMed CentralGoogle Scholar
- Tu SW, Tennakoon L, O'Connor M, Shankar R, Das A: Using an Integrated Ontology and Information Model for Querying and Reasoning about Phenotypes. 2008, In Proceedings of AMIA Annu Symp: 8–12 November 2008, Washington DC, 727-731.Google Scholar
- Povey S: Al Aqeel AI, Cambon-Thomsen A, Dalgleish R, den Dunnen JT, Firth HV, Greenblatt MS, Barash CI, Parker M, Patrinos GP, Savige J, Sobrido MJ, Winship I, Cotton RG; Ethics Committee of the Human Genome Organization (HUGO): Practical guidelines addressing ethical issues pertaining to the curation of human locus-specific variation databases (LSDBs). Hum Mutat. 2010, 31: 1179-1184. 10.1002/humu.21339.View ArticlePubMedPubMed CentralGoogle Scholar
- Jonquet C, LePendu P, Falconer S, Coulet A, Noy NF, Musen MA, Shah NH: NCBO Resource Index: Ontology-Based Search and Mining of Biomedical Resources. Web Semantics. 2011, 9: 316-324. 10.1016/j.websem.2011.06.005.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoehndorf R, Dumontier M, Oellrich A, Wimalaratne S, Rebholz-Schuhmann D, Schofield P, Gkoutos GV: A common layer of interoperability for biomedical ontologies based on OWL EL. Bioinformatics. 2011, 27: 1001-1008. 10.1093/bioinformatics/btr058.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith B, Brochhausen M: Putting biomedical ontologies to work. Methods Inf Med. 2010, 49: 135-40. 10.3414/ME9302.View ArticlePubMedPubMed CentralGoogle Scholar
- Köhler S, Schulz MH, Krawitz P, Bauer S, Dölken S, Ott CE, Mundlos C, Horn D, Mundlos S, Robinson PN: Clinical diagnostics in human genetics with semantic similarity searches in ontologies. Am J Hum Genet. 2009, 85: 457-64. 10.1016/j.ajhg.2009.09.003.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6947/12/78/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.