• Explore Vox
  • Culture
  • Entertainment
  • Life
  • Music
  • News & Politics
  • Technology
  • Join Vox
  • Take a Tour
  • Already a Member? Sign in
Allyson

Systems Biology & Bioinformatics

Ontologies, Data Integration, Workshops, News, Papers, etc..

  • Allyson’s Blog
  • Profile
  • Neighbors
  • Photos
  • More 
    • Audio
    • Videos
    • Books
    • Links
    • Collections

ResearchBlogging.org now available

  • Jan 22, 2008
  • 6 comments

An article from the BPR3.org website announcing the ResearchBlogging.org website appeared today in my RSS reader. I'd had a look on the BPR3 website the other day and thought it a good idea to have some way of marking those posts of mine specifically dealing with peer-reviewed research. Well, now ResearchBlogging.org provides a citation-composer on their website to make adding the BPR3 logo to your blog even easier: just put in your DOI, and then the code to embed in your blog for both the appropriately-formatted citation and the logo itself magically appear. I'll post my first BPR3-identified post shortly, and you'll see that Vox doesn't render the embedded code exactly correctly, but it's a start. There's no RSS feed on ResearchBlogging.org yet, but it's a must, and they do say that one is in the works.

For those that haven't heard of this, here's a summary from their website:

"Do you like to read about new developments in science and other fields? Are you tired of "science by press release"? BPR3 is your place. BPR3 allows readers to easily find blog posts about serious peer-reviewed research, instead of just news reports and press releases. We provide bloggers with an icon they can use to show when they're talking about a peer-reviewed work that they've read and analyzed closely."

How do others feel about this idea? There are already a number of posts (here, here, here, and here, just to name a few). I like it :) The fact that everyone writing about a paper will be using a similar citation display might mean its easier to find other people that have opinions on a particular research topic. But perhaps you think it is too much classification, or that we shouldn't try to judge the "quality" of various posts. After all, though you can report people misusing the logo, just having it isn't a mark of the seriousness or quality of a paper, or of the seriousness or quality of the person writing about the paper ;)

6 comments Tags: standards, bioinformatics, systems biology, bpr3

Whole-Genome Reference Networks for the Community

  • Jan 22, 2008
  • Post a comment
Srinivasan et al use this paper as a call to the community to begin the development of whole-genome reference networks for key model organisms. This paper is a combination of a review (in that it summarizes methods of network generation and analysis) and a call to arms, stating that reference networks are needed. It begins by describing systems biology as "the science of quantitatively defining and analyzing" functional modules, or components of biological systems.

There are many different definitions of systems biology (see here, here, here, here and here, just to name a few), but generally it seems the twin pillars of data integration and study - at various levels of granularity - of biological systems are present in most of them. A focus on integration and top-down research rather than the more traditional reductionist point of view is also often mentioned.

The authors then divide systems biology into three broad categories: high-level networks of the interactome or metabolome, deterministic models of kinetics and diffusion, and finally stochastic models of variation in cell lines. This division would be slightly clearer if they specified continuous deterministic models and discrete stochastic models. I realize that these adjectives are generally assumed for these model types, but as it is their discrete- or continuous-ness that increases the complexity of the models, it is something that would be useful to include.

They collapse many different types of network data into a single global interaction network, stating that it would be prohibitively expensive to try to prise out all of the sub-graphs, as variables such as time or sub-cellular location are often not simple to pull out on their own. This "lowest common denominator" method of network generation is not ideal, but does provide more information than, they attest, a simple genome sequence. In their networks, nodes represent proteins and edge weights are probabilities of association between proteins.

Noise is a real problem in most of these high-throughput data sets, and such data sets are not all created equal: one group may make a very good gene expression data set, and another may not. How can variable quality of data be dealt with? Early efforts focused on integrating multiple networks and only taking those nodes and edges that were present in more than one network. After that, methods of network generation that used "gold standards" created better integrated networks.

Descriptions of network analyses (rather than network creation) focus on network alignment and experiment prioritization. The latter is a general term for pulling out elements of the network that haven't been experimentally verified, such as likely additions to known pathways or important disease genes. The former is an interesting extrapolation to networks of sequence alignments for genomes. In network alignments, conserved modules of nodes are identified if they have "both conserved primary sequences and conserved pair-wise interactions between species". They specifically mention Graemlin, which is a tool they have developed that can identify conserved functional modules across multiple networks.

Finally, they suggest that the reference networks should show only those reactions present in the "‘average cell’ of a given organism near the median of the norm of reaction".

While they acknowledge that, like the reference human genome sequence, such a creation is a "useful fiction", it is my opinion that finding the average cell will be much more difficult, and perhaps less illuminating, than its equivalent in the sequencing world. Further, describing what is "normal" is something that is truly difficult, and will vary from species to species. The PATO / quality ontology people (http://obofoundry.org/cgi-bin/detail.cgi?quality) have known about the problems facing the "average" phenotype for a while now. I do, however, like their idea of storing the reference networks using RDF, as that seems a fitting format for networks. Overall, a laudable goal but one which will need some more thinking about. I've tried to run Graemlin using one of their example searches, and it didn't run (at least today), and the main author's website won't load for me to today, though one of the other author's pages did work.

All-in-all, a useful review of recent network methods in bioinformatics, and an interesting goal described. Low-noise reference networks for key model organisms, together with the annotation tracks that would describe deviations from the norm is a good idea.

Topics for discussion (aka leading questions): More fine-grained reference implementations are available, such as Reactome. Reactome provides a curated database of human biological pathways, with inferred orthologous events for 22 other organisms. Do we need reference networks when we're gradually growing our knowledge of reference pathways? Are reference networks of "normal" organisms states helpful? How do we define average? Would the median of the norm of a reaction be different under different environmental conditions? What if what one group considers an average cell differs from another group's average cell? Having reference networks would mean easier comparisons of different network analysis programs. Would this end up being a major use of the networks? Would such comparisons just lead to network analysis programs that fit the reference network, but not work in a generic manner? What do others think?

Srinivasan, B.S., Shah, N.H., Flannick, J.A., Abeliuk, E., Novak, A.F., Batzoglou, S. (2007). Current progress in network research: toward reference networks for key model organisms. Briefings in Bioinformatics, 8(5), 318-332. DOI: 10.1093/bib/bbm038

Post a comment Tags: bioinformatics, systems biology, network analysis

Google Notebook and the Google Cloud in Bioinformatics

  • Jan 10, 2008
  • 8 comments

Ever since my friend Frank posted an article about how he has one foot in the internet "cloud", I have been a little worried about how much of my cloud is composed of Google apps. I currently have a Google cloud composed of Google Mail, Calendar, Reader, Docs, Groups, and most recently Notebook. It's convenient, the apps are useful, and it means I can access much of the information relevant to my work in bioinformatics wherever I am. I'm by no means a Google-only person (I use Vox rather than Blogger, for instance), and my set of a dozen or so "home" pages that I like to load whenever I open Firefox doesn't even have Google in a majority of tabs. But I can't help feeling slightly guilty and a little be worried about how dependent I've become on Google products. Then on top of that, I can't help feeling silly to be feeling guilty. It's my choice, after all!

So, moving past all that: today I've added Google Notebook to my cloud. I wanted to have somewhere where I could stash all those interesting papers that are sent to me each week via my personal tame searches on PubMed (aka MyNCBI: perform a regular search, click on "Save Search" and then ask email updates to be sent to you). I tried Connotea but never went and actually looked at it. I tried reading them when the emails arrived, but the time I have to read papers is never regular. So I deliberately had a look at the Google App list, and found Notebook. You can make notes on anything, but there is also a Firefox plugin that allows you to just highlight any portion of a web page you want and then place it directly in its own section of one of your notebooks. Further, Google Notebook is private unless you make it public. Here I can put all of my thoughts on these papers, whenever I have a minute to do it. Each clipping can be labelled, and comments can be added. The plugin is nice, but I've already made it break on one of the PubMed result pages. I suppose I'm talented?

Perhaps I'll finally get both the time to read all those interesting-sounding papers as well as have a single location that I'll actually visit. Only time will tell.

8 comments Tags: tools, bioinformatics, google notebook

3 Bioinformatics Research Associate Positions: Newcastle University

  • Dec 18, 2007
  • Post a comment

There are three bioinformatics jobs (one in pure bioinformatics, one in network analysis, and another in modelling/mathematical biology) currently available within CISBAN, an interdisciplinary centre studying the systems biology of ageing and nutrition. The full particulars are posted both on Nature Jobs and on the Newcastle University Job Vacancies web pages.

Below are links to the various job advertisements, as well as summaries of the jobs themselves. This is a summary of the three Nature Jobs postings, put together on a single page for easy perusal. The closing date for all of these positions is 11 January 2008. This is a great opportunity, though I may be speaking from a biased perspective as I work at CISBAN and find it an interesting and challenging workplace.

  1. Centre for Integrated Systems Biology of Ageing and Nutrition, Institute for Ageing and Health

    Research Positions

    Level F £25,134 – £32,796 p.a.
    Level G: £33,779 – £40,335 p.a.

    We seek scientists to join CISBAN, an exciting new research centre established following a major award (£6.4m) from BBSRC and EPSRC, to participate in studies of the mechanisms responsible for ageing and how they are affected by nutrition. Ageing is recognised internationally as a ‘grand challenge’ and is a field prioritised for growth. This post offer opportunities to work in an intensely multidisciplinary, world-class centre and contribute to the development and application of systems science.

    Research Associate (Bioinformation/Computing Scientist – Applications)

    To develop and maintain the computing software and hardware infrastructure for systems biology, including a central web portal integrating applications for data capture, storage and visualisation and high performance computing systems and databases, including a large Linux cluster.

    Job reference: A1091R

    Posts are tenable until 30 September 2010.

    Enquiries for the post may be directed to Dr Anil Wipat, School of Computing Science (email: anil.wipat@ncl.ac.uk) Further particulars for this post can be found on the University’s web page at http://www.ncl.ac.uk/vacancies/list.phtml?category=Research.

    Applications should be submitted by 11 January 2008 to Professor Tom Kirkwood, CISBAN Director, Institute for Ageing and Health, Henry Wellcome Laboratory for Biogerontology Research, Newcastle University, Newcastle upon Tyne NE4 6BE (email: tom.kirkwood@ncl.ac.uk).
    Committed to Equal Opportunities

  2. Centre for Integrated Systems Biology of Ageing and Nutrition, Institute for Ageing and Health

    Research Positions

    Level F £25,134 – £32,796 p.a.
    Level G: £33,779 – £40,335 p.a.

    We seek scientists to join CISBAN, an exciting new research centre established following a major award (£6.4m) from BBSRC and EPSRC, to participate in studies of the mechanisms responsible for ageing and how they are affected by nutrition. Ageing is recognised internationally as a ‘grand challenge’ and is a field prioritised for growth. This post offer opportunities to work in an intensely multidisciplinary, world-class centre and contribute to the development and application of systems science.

    Research Associate (Bioinformatician – Network Analysis)

    To research and develop novel methods of representing and integrating molecular and cellular data as networks and apply this methodology to identify novel proteins and elucidate novel pathways involved in the process of cellular ageing and senescence.

    Job reference: A1090R

    Posts are tenable until 30 September 2010.

    Enquiries for the post may be directed to Dr Anil Wipat, School of Computing Science (email: anil.wipat@ncl.ac.uk) Further particulars for this post can be found on the University’s web page at http://www.ncl.ac.uk/vacancies/list.phtml?category=Research.

    Applications should be submitted by 11 January 2008 to Professor Tom Kirkwood, CISBAN Director, Institute for Ageing and Health, Henry Wellcome Laboratory for Biogerontology Research, Newcastle University, Newcastle upon Tyne NE4 6BE (email: tom.kirkwood@ncl.ac.uk).

    Committed to Equal Opportunities

  3. Centre for Integrated Systems Biology of Ageing and Nutrition, Institute for Ageing and Health

    Research Positions

    Level F £25,134 – £32,796 p.a.
    Level G: £33,779 – £40,335 p.a.

    We seek scientists to join CISBAN, an exciting new research centre established following a major award (£6.4m) from BBSRC and EPSRC, to participate in studies of the mechanisms responsible for ageing and how they are affected by nutrition. Ageing is recognised internationally as a ‘grand challenge’ and is a field prioritised for growth. This post offer opportunities to work in an intensely multidisciplinary, world-class centre and contribute to the development and application of systems science.

    Research Associate (Modeller/Mathematical Biologist)

    To develop models of molecular and cellular mechanisms of ageing and to explore links between ageing, development and evolution from a life-course perspective. This post will also involve collaboration within the EU Network of Excellence LifeSpan, linking development and ageing.

    Job Ref: A1092R

    Posts are tenable until 30 September 2010.

    Enquiries for the post may be directed to to Professor Tom Kirkwood, Institute for Ageing and Health (email: tom.kirkwood@ncl.ac.uk) Further particulars for this post can be found on the University’s web page.

    Applications should be submitted by 11 January 2008 to Professor Tom Kirkwood, CISBAN Director, Institute for Ageing and Health, Henry Wellcome Laboratory for Biogerontology Research, Newcastle University, Newcastle upon Tyne NE4 6BE (email:* tom.kirkwood@ncl.ac.uk).

    Committed to Equal Opportunities

Post a comment Tags: jobs, modelling, bioinformatics, systems biology, cisban, newcastle university, research associate, network analysis …

FuGE Users' Workshop: 13-14 December, 2007

  • Dec 14, 2007
  • 15 comments

The two-day FuGE Users' workshop was organized by Norman Paton and held at the University of Manchester. It was great fun, and if you just want the short summary of my time there, then just know that there was loads of enthusiasm for FuGE as well as interesting talks, both by communities who were already extending FuGE, and by  developers who were already building tools and databases based on it. There were only a dozen or so people, which kept the discussions lively but neither too long nor too divergent. The workshop dinner was great, though the trip to the restaurant was correctly described by one of the attendees as an Odyssey. For more information on the social aspect of the FuGE workshop, please have a look at Phil Lord's humorous posting on the matter. For another post on the workshop, see the peanutbutter Bioinformatics blog by Frank Gibson.

If you wish to read the longer notes rather than the short summary, then please read on!

Please note that these are my own notes, and are in no way considered to be an "official" FuGE report on the workshop. As such, any errors or inconsistencies are entirely my own. However, if you see a problem with this post, then please let me know, and I'll fix it!

The objectives of the workshop were to share and document experiences in the use of FuGE, to identify good-practices, to document guidelines, and to make known these experiences and guidelines. Hopefully, the result will be a paper that documents the current users' experiences and increases communities' understanding of FuGE. It will hopefully help people who who have read the Nature Biotechnology paper and want to use FuGE, but aren't completely sure what to do next.

Attendees were:

Peter Wilkinson, from Montreal, who was interested in FuGE for flow cytometry.

Khalid Belhajjame: works with Norman Paton in Manchester, and who may soon be a full-time developer of FuGE

Javier Santoyo: University of Edinburgh, part of a consortium trying to develop standards for RNAi work

Andy Jones: one of the original developers of FuGE, from Liverpool, developed GelML with Frank Gibson.

Heiko Rosenfelder: German Cancer Centre at Heidelberg, here as part of MIACA, and wants to use FuGE for the cellular assay format.

Martin Eisenacher: Proteome Centre (mzML and analysisXML) and wants to use FuGE

Phil Lord, Frank Gibson: via CARMEN, wants to use FuGE. Frank also developed GelML with Andy Jones.

Neil Wipat, Matt Pocock, Allyson Lister: We use FuGE in our internal application for storing HT data. Matt and Allyson also involved in OBI.

Leandro Hermida: SIB, they're part of a group that is making SystemsX. Want to use FuGE to store and manage the data. Also want to make an extension of FuGE for deep sequencing technologies.

Norman Paton: originally from proteomics field, but developer of FuGE and organizer of the workshop.

Session 1: Experiences Using and Extending FuGE

GelML – Frank Gibson and Andy Jones

GelML is a FuGE extension that has passed the PSI approval process. PSI defines community standards for data representation in proteomics. There are a variety of working groups, including gel electrophoresis, mass spectrometry, protein modifications, etc. Within the Gel WG there are three specifications: MIAPE-GE (minimum checklist for reporting gel elecrotphoresis experiments), sepCV (controlled vocabulary), and GelML (data transfer format, based on FuGE).

GelML covers the model of a gel, 1-D and 2-D GE, other (N-dimensional) GE's, sample loading, electrophoresis, detection, image acquisition, the excision of locations on gels, and SubstanceMixtureProtocol and SubstanceAction.

The first extended FuGE class described was the Material abstract class. The first of such classes is the Gel class. A Gel has Dimensions, MeasuredMaterial, and others. You use the “Measurement” package to describe the characteristics of the Gel. Measurements include PercentageAcrylamide, while information about the gel (i.e. if purchased, from where), information on the Dimensional Units and CrossLinkerRatio are all FuGE OntologyTerms). MeasuredMaterial was not originally in FuGE because it was planned that such substances could be captured by ontology terms. Rather than using named associations to GenericParameter, they tended to use either GenericParameter (with a CV term) or extended the Parameter class. This was just a design decision, and he would like to see how others do it.

Another extended FuGE class is the Protocol abstract class. The GelML SampleLoadingProtocol has an AddBufferAction which points to a SubstanceMixtureProtocol. 2DGelProtocol has a SampleLoadingAction, a FirstDimensionAction, a SecondDimensionAction (both Electrophoresis protocols), and an InterDimensionAction (for when something happens between the first and second dimension actions), and the DetectionAction.

Within the Electrophoresis protocols there is the ElectrophoresisStop (an Action) which contains a StopTime, which is a TimeParameter, with has Duration and TimePoint. They'd be really interested to see how others have/would like to model time.

It was also a design decision to guide people with the structure of the XML to help them know what to fill out, e.g. you must have a 2dGelProtocol.

For each case, should we extend the FuGE model or add experiment-specific semantics through the use of ontologies? I think this is a case of using both, depending on the circumstances.

They have used standard XML references within the documents. But, for instance, do we still need internal document identifiers when the ontologyURI is a globally-unique identifer anyway? Maybe required if the terms are created ad hoc within the group making the XML file. What is the best way to use ontologies?

AnalysisXML – Martin Eisenacher

http://www.fp6-prodac.eu

He is a member of the ProDaC Consortium. ProDaC is a funded consortium that is meant as a “coordination action” within the 6th EU Framework Programme. Its aims are the development of international standards, standardized data submission pipelines, systematic data collection in public standards-compliant repositories, and data access for community and publication. There was a kick-off meeting of ProDaC in Long Beach in October 2006, and there have been two workshops since.

Proteomics data includes spectra (peak lists), and peptide lists. He works specifically with the MS (for peak lists and instruments, mzML) and Proteomics Informatics (for “results”, analysisXML) PSI WG's. mzML is a merger of mzData and mzXML. Perhaps this merger is one of the reasons that it is not currently FuGE-based. AnalysisXML includes annotation of search databases, search, algorithms, search parameters, instrument characteristics, peptides (peptide-spectrum link, peptide scores), proteins (protein-peptide link, protein scores, significance values, false-discovery estimation) and quantisation. In September 2007 they added comments into the UML that are passed into the XML.

They use the MagicDraw Community Edition, which is available for free. The Analysis package is subdivided into process, quantisation, and search. Process contains things that aren't directly related to the search protocol applications, but other steps such as ProteinDetermination and PolypeptideProcess. Some of the classes they have made that inherit from the Data class inside the search package include AnalysisResultSet (a set of spectra), AnalysisResult (a spectrum), and AnalysisResultItem (all peptides found for that spectrum). These are all abstract classes, whose concrete subclasses include PolypeptideResultSet, PolypeptideSearchResult, and PolypeptideResultItem.

At the moment they are assembling their own CV (to include search parameters that are most commonly used in search engines like MASCOT), but they can also use Pride CV. They use the ontology classes directly from FuGE, without extending it. This means that it fits what they need without modifications.

In analysisXML, peptides and sequences are listed only once. Different types of analyses in one file or in separate files with external cross-references. Also, the AnalysisProtocol could be used as parameter input for search engines. However, there are many cross-references and unique identifiers that are not validated by the FuGE Schema. Further, there are external cross-references to mzML, which can be difficult if you have only local files and not public URI's. Also, sequences (just the letters) are not polypeptides (“real” molecules with modifications). Therefore, the ConceptualMolecule FuGE class is not appropriate for polypeptides, though it is suitable for sequences (though they are still able to use that class). Additionally, the ResultSet-ResultItem hierarchy does not fit all analysis types. Finally, many FuGE elements seem to have very long names that aren't always useful (but you shouldn't be typing XML manually!).

All items of the collections have unique identifiers. References to them are attributes called “..._ref”. Schema validation does not consider whether _ref links to an allowed section (or that used CV's are allowed). In mzML, for example, “semantic validation” of CV's is possible (suggested/implemented by the EBI). Are identifier checks possible? ProDaC has an online validator for mzML, analysisXML, mzData and prideXML that performs semantic validation, though the extra ontology/CV checks are only supported for mzML.

Still to do is the finalization of analysisXML, which is a deliverable for last October! They also want to provide “Quality Determination” as a process. They also want to make some use-cases and instance documents. They will have some from Matrix Science, MPC. Also, they need to finalize the CV they are using.

SyMBA – Allyson Lister

I gave this talk, so I didn't write anything about it! Instead, have a look at the SourceForge website (http://symba.sf.net):

The Integrative Bioinformatics Group, headed by Neil Wipat and part of The Centre for Integrated Systems Biology of Ageing and Nutrition (CISBAN), has developed a data archive and integrator (SyMBA) based on Milestone 3 of the Functional Genomics Experiment (FuGE) Object Model (FuGE-OM), and which archives, stores, and retrieves raw high-throughput data. Until now, few published systems have successfully integrated multiple omics data types and information about experiments in a single database. SyMBA is the first published implementation of FuGE that includes a database back-end, expert and standard interfaces, and a Life Science Identifier (LSID) Resolution and Assigning service to identify objects and provide programmatic access to the database. Having a central data repository prevents deletion, loss, or accidental modification of primary data, while giving convenient access to the data for publication and analysis. It also provides a central location for storage of metadata for the high-throughput data sets, and will facilitate subsequent data integration strategies.

Developing Flow Cytometry FuGE extensions – Peter Wilkinson

Developing MIFlowCyt. Originally, they stored the metadata and data in a single file, but their latest format (ACS) will separate these two types. They are considering having some of their data formats be in RDF as well as XML, even for those formats that will be built on FuGE – is there a good XML to RDF converter? I suppose so, as I've been able to save OWL/RDF as OWL/XML in Protege 4.

One example of their extension is Cytometer, which is a subclass of equipment. How descriptive should they get with their samples? Should it be at the entity or attribute level? For instance, there is a conceptual difference between prepared samples and “generic” materials. But why not draw an association to material and call it “sample”? They can't do that because sample has a lot of associations itself that aren't present in Material. For things like buffers and solutions, spML doesn't seem to view them as things that exist – you just talk about them in the protocol. This way you don't have to list them 1000s of times. In FC, you have to know exactly which thing is used in the protocol (e.g. they must record batch numbers). However, you could have a single buffer instance, and then in the ProtocolApplication you have a specific parameter that is modified in that particular application of the Protocol, such as the batch number.

Open issues include: FuGE should reference a stable version of AndroMDA, there should be a best-practice for deciding when a Generic* class is replaced by a specific omics-type class, how is the OntologyTerm abstract class intended to be used for specific controlled lists, fitting the organism into FuGE::Bio somewhere, and versioning.

He's also trying to write a FuGE database by hand, rather than using what is generated by AndroMDA, as he needs to squeeze as much performance out of the system as he can. Much more difficult, but could conceivably be much more efficient.

Generic and Custom Extensions – Andy Jones

spML is for sample processing. SubstanceMixtureProtocol is for describing a mixture of substances, e.g. buffers and solutions and the method of their creation. Actions relate to constituents. Timings relate to constituents and volume, concentration, or mass. SetPropertyAction is a generic model to be used in conjunction with protocols where parameters may be set with associated ActionText. Their chromatography extension comprises extensions of Protocol, Equipment, and ProtocolApplication. The ChromatographyProtocol contains extensions of Parameter, has a child protocol for sample injection, and various uses of GenericActions. With ChromatographyEquipment, there is column-associated sub-components. All extensions of Chromatography equipment can have additional parameters, including specific named parameters where they are always required. Uses Equipment:make. The mobile phase of LCProtocol is described using the SubstanceMixtureProtocol. Inputs are defined with GenericMaterialMeasurement, and the outputs are either Chromatogram (ExternalData) and SeparationFraction. You can also have two-dimensional chromatography.

GenericSeparation is a protocol that uses generic models for defining substance used to create a separation gradient and the parameters applied. In this case, the equipment defines the type of separation and criteria using ontology terms – but how do you communicate how this should be used to all of the developers? In contrast, we don't want to have huge models. Inputs also defined using GenericMaterialMeasurement, and the outputs are either SeparationLog (ExternalData) and SeparationFraction.

TreatmentProtocol is a simple model for treatments, intended for labelling, mixing, splitting, and washing, for example. The treatment IO in TreatmentProtocolApplication is restricted to having material inputs and material outputs only.

There seems to be three sorts of models: column-oriented, category-oriented, and completely generic protocols. Much of what is in spML might be useful for the “library of models” we've been discussing.

The generic model is very flexible for different types of separation, and could be used for LC, GC, capillary electrophoresis, rotofors etc. It is also unlikely to break if new type of experiment is defined, and the Treatments model could potentially be useful in the context of any experiment type. Also, the generic model is much smaller, and can be used in various ways. However, this last one could be a “con” as well, because different users/implementers are likely to encode the same information in different ways. Further, a specific model can guide the user to provide specific details, e.g. for MIAPE compliance.

spML units are derived from the OBO Unit Ontology. Should FuGE extensions be allowed to have user-defined terms? It would be useful for the creation of in-house lists to populate drop-down menus.

Issues

spML units Below are a list of questions and suggestions that we came up with while the initial talks progressed in the first couple of sessions. Many were discussed, and some were answered, in breakout sessions later. Notes from the discussions I was a part of are included below. The unanswered points in the list may have been discussed at other breakout sessions, or may still be untouched.


Discussion on Semantic Validation and Identifiers: Khalid Belhajjame, Norman Paton, Allyson Lister, Martin Eisenacher

Identifiers and Auxillary/Semantic Validation: Types of Validation and How Simple Support can be done.


Unique in Document

Not Dangling

Globally Unique

Type Correct

Notes


Property

Checked by XML Tooling

Property

Checked by XML Tooling

Property

Checked by XML Tooling

Property

Checked by XML Tooling

Instances of Identifiable

yes

GP

yes

GP (yes?)

(+)

no

yes

no

See (1)

Ontology Terms

(#)

n/a

yes

no

yes

no

yes (not in UML)

no (^)


External Data

($)

n/a

yes

no

yes

no

yes

(*)

(*)want to know it's a file of a particular kind

GP: Can be checked with a generic program.

All things marked GP or X could be attacked by people wanting to write a semantic validation tool.

(+) Only for some types of Globally-unique identifiers would we be able to check that they were truly unique and well-formed.

(#): Should OntologyTerm elements be unique (irrespective of their identifiers, which must be unique)? If people compare OT identifiers they may think two terms are different when in fact they are the same, and someone was sloppy when making OT elements. However, if they have linked their OT to an OntologySource then it can be checked if it is both unique in document and globally unique (if it is a logical/physical uri)? In that case, why should OS be optional at all, if custom CV's can be included in the OS.

(^): This is where the ontology mapping files come in.

($): The same argument for uniqueness of ED applies as that in OT.

(1) Will we suggest a type of identifier to use with FuGE as a best-practice?

Do we still need internal document identifiers when the ontologyURI is a globally-unique identifier anyway?

Should identifiers be human readable?

Do community extensions automatically have their own namespace/prefix? That is, if “sample” is used in the FC community and also in another extension, will it be problematic if you try to create a multi-omics FuGE-ML file? This is all about linkage between different FuGE-based instances (unique identifiers, both within a single document and between documents.What is the identifier an identifier of? Is every Identifiable object a “first-class citizen”? We shouldn't force all (any!) identifiers to be URI's.

Should you use a logical or physical naming scheme?

  • Physical naming schemes:

    • Are fragile

    • May not work for all users (i.e. if the URI points to a laptop that isn't publicly accessible)

  • Logical naming schemes:

    • Are robust, but require a greater investment of time, as they need tools that provide resolving facilities.

  • If locally-unique identifiers are used:

    • it means that you may get into trouble in the long run

  • If globally-unique identifiers are used:

    • clashes between different FuGE-ML files will be avoided

  • People should look this over and discover which is the best setup for their situation.

If we use URI's, should URI's be resolvable? What is the scope?

  • Martin has a URI that points to a data file, and a (possibly locally unique) identifier that points to a spectrum within the data file. How to deal with this? Do we have a best-practice for it?

Schema Validation

  • Schema validation does not consider whether _ref links to an allowed section (or that used CV's are allowed). Native XML validation does not do this, but you could make a tool. In theory, the prefix before the _ref is always the name of the class. FuGE needs semantic validation.

  • How should user-defined ontology terms be validated in the XML?

Discussion on Versioning: Khalid Belhajjame, Norman Paton, Allyson Lister, Phil Lord, Matt Pocock, Leandro Hermida

  • Is there a best way to implement versioning?

Characteristics of (SyMBA) Versioning:

  1. Complete History of Atomic Changes

  2. Low Cost of Updates – No Cascades

  3. Higher Cost of Retrieval

This is actually a transaction-time database with tuple-level timestamps. In a transaction-time database, the time is in the world of the database and not the time in the real world (vs valid-time database, where the time you insert does not match the time that you actually wanted to input). If you don't put the timestamp in the tuple, you put it in the attribute. In this context, people have looked at the properties of update operations.

Can't just use LSID versioning because there is no specification of how the version should be updated.

SyMBA Versioning Requirements:

  1. Preserving the semantics of the LSID

  2. Getting exactly the version requested, and getting all versions

  3. Nothing should disappear


This isn't necessarily versioning, or what versioning in FuGE should be.

Leandro's Versioning Requirements:

  1. Getting exactly the version requested, and getting all versions
  2. Nothing should disappear

Should this be done in FuGE, or in the FuGE-OM specifically? Perhaps just in the Maven build? We could put hooks in FuGE that would allow fine-grained logging. The current Audit setup does not allow linking back to previous versions unless you put the delta in free text somehow. The Audit classes may be suitable for XML, and you could make a log of such changes and roll-back (in a non-RDBMS way) to whatever version you want.

While it is clear we could make an STK that could have versioning of some type, whether or not this should be a (optionally-used) change to the OM is a much bigger thing. It is certainly a worry that versioning has to be dealt with at the application level. However, versioning at the file or XML level means multiple files, otherwise you'd have to apply a diff to a very large file.

We haven't really had the time to scan the space of options here. We could circulate a general document, and then outline what's actually been done so far. A paper would, in any case, be centered around pros/cons, and a bit less on current implementations, but definitely not say which is the “right” way to do things, as there is no single right way.

There are different technical solutions, and not all of these solutions should necessarily be provided in the model.


Discussion on Tools – Leandro Hermida, Heiko Rosenfelder, Neil Wipat, Phil Lord, Allyson Lister

  • What about trying to get some automatic mapping between the XML classes and the Hibernate/Spring classes?

    • There is a disconnect between the XSD that is generated from the XML schema cartridge and the code generated from your persistence cartridge.

    • This means you have to write your mapping manually.

    • There is a possibility that we could get hyperjaxb3 to work for this (Allyson had tried with an earlier version but it didn't work properly). Hyperjaxb3 generates both Entity POJO's and the jaxb2 classes. So, in theory you could only use the Andromda XSD cartridge and hyperjaxb3 for the rest. However, then you loose all the information that is present in the UML but not in the XSD.

    • Hyperjaxb3 uses both hibernate and ejb3 natively (you can choose). Leandro wants to work on a merged persistence/hyperjaxb3 extended cartridge, or perhaps its own cartridge. So perhaps the generation of a hyperjaxb3 cartridge.

  • Is there an XSLT that could be made to have a “standard” way of viewing a fuge experiment?

    • Khalid mentioned that it is important to allow input from the programmer in such a tool, so they can see as little or as much of the FuGE structure as you wish to present to the user.

    • Leandro is working on an ejb3 cartridge from the androMDA plugins project (not part of the AndroMDA core yet), and have used FuGE as a test-case. What this cartridge does is generate a mapping file and load it into any application server running hibernate and it will generate your database. Whereas the Hibernate+Spring cartridge generates 1) Entity POJO's + mapping files 2) Spring DAO + DAOException + DAOImpl. With the ejb3 cartridge you get 1) ejb3-annotated Entity POJO's + DAO*.
      You can use Spring, if you wish, to build your web framework. Leandro decided to instead use Seam, which is the business layer of a web framework that builds on top of ejb3. Seam then uses the JSF (Facelets) and Jboss RichFaces for the actual web UI.
      To get the Seam classes, you model <> classes and then draw dependencies, which then auto-generate Seam-enabled ejb service beans. However, the Facelets and RichFaces have their xhtml files manually, though AndroMDA creates the entire web/ structure and base Seam classes for you.
      This doesn't answer our simple UI question.
      The ejb3 cartridge has a web service (jax-ws, via soap) to your DAO's and Entity POJO's.

    • With MAGE, someone wrote a regular Java Swing program where you download the jar which opens a little tabbed client that views MAGE. We could do something similar. (A J2SE app to write/read FuGE-ML of nice wizard interface)

    • The GSC has a lightweight XSD-to-web-form software app.

    • An XSLT, which is a style-sheet that is richer than CSS, but it is a tough language to use (convert XML to “HTML”). XSLT's don't have first-class functions, so you can't do anything generic.

    • Also would like to have simple jar that has input XML, output HTML.

    • This means three tool types here: 1) heavyweight (already existing in SyMBA and SystemsX) 2) midweight (J2SE app to read/write with a wizard-like interface) 3) lightweight (input XML, output HTML with some simple options).

  • Tool support for FuGE STK version 1, including a validator.

    • The MAGE STK includes a validator.

    • XML validation can be done with JAXB2 as is with the Milestone 3 STK, but longer-term we need the semantic validation tool.

    • Perhaps have some ontology lookup helper classes (OLS from the EBI?) to help users and developers add terms from (a certain set of?) ontologies. This may help people populate their databases, choose a term from a list on a front-end tool, etc.

  • Tool support for database schema / AndroMDA / Alternatives.

    • Dealt with in the other sections

Discussion on Challenging Constructs, including Investigation Package, Abstract Associations, and the Ontology Package - the entire group

  • What is the real meaning of the Investigation package? It's one of the few parts of FuGE that isn't meant to be extended.

  • How is the OntologyTerm abstract class intended to be used for specific controlled lists? One example is taxonomies as opposed to ontologies.

The intention is that this package would not be changed or extended by communities. Each technology would be reported in the InvestigationComponent. The Factor class actually is meant as a summary report of the factors used in the experiment. There is currently no direct link between the Factors and the protocol workflows – the detail can be recorded in other places in FuGE. It's a summary and duplication of the factor information.

So, if you want to say that your Investigation compared two different values for a single factor, the Investigation has the factor type, but not the data for the factor or the values themselves. However, you can connect to the data made from the various omics technologies via the DataPartition class. There could be a problem where it is a set of factors that only together make a particular set of data useful. Example: if your important aggregate of factors is time1.mouse1.foodstuff1. However, you would have to have each of these factors would be named separately, and you would get a different slice of data for the time1 than you would for the mouse1. How to you join them up? Perhaps allow multiple FactorValues (and OntologyTerms) for a single Factor. Not a very nice solution, though. Perhaps you don't need to change it at all, as you would only add Factors that are relevant to a particular InvestigationComponent. How do you describe which combinations of Factors are the combinations you're interested in? Norman did this by seeing an IC as a particular run of an experiment.

Dimensions are used in FuGE as a way of naming coordinates in a matrix. This does not mean that the data has to be stored here. You can store the data internally via the InternalData class, or you can reference it externally via the ExternalData class (or, of course, create your own subclass of Data).

There are 21 <>'s in FuGE, and all but 6 have identical concrete associations. Some auto-generated AndroMDA code mistakenly ignores the “abstract” parts and incorrectly generates the methods etc. In this case, you can just delete the abstract association in your copy of the UML and re-generate the code. It should be fixed within AndroMDA, though.

For multi-dimensional data, DataPartition are meant as a mechanism to relate back to the data from the Investigation, but many groups will choose not to use DataPartition. Very big, regularly-shaped data sets will be good things to use DataPartition with (e.g. Flow Cytometry). In the case of proteomics data, this may be more of a challenge. A best-practices documents should contain information on which data types are best-suited to this system, and which aren't. It should also include any alternatives to this system. One alternative solution to using DataPartition and its associated coordination system for dimensional data is to build an association from their data of interest back to FactorValue.

What is PartitionPair? In the case where you have two data files, and you wish to associate a particular row of one data file (for example) to a particular spectrum in another (to continue the example). So, it is a “shortcut” to linking particular data sets.

How should users of FuGE build CV lists using the Ontology Package? An OntologyTerm has an OntologyProperty, which contains both a DataProperty and ObjectProperty (these are the relationships within an ontology). Also inside OntologyTerm is OntologyIndividual. OI is the individual itself. Why not just provide the term – why try to recreate the structure of an ontology into UML? However, in OWL, every single class, relationship etc has a URI, so why not use those in UML? An example use: you have in an ontology the concept of age, which has an initial time point and a unit. How do you pull that concept into the UML? We're essentially creating a cut-down version of the ontology to allow extensibility in FuGE. But why would you want this? To create an individual of an ontology within the UML. It also allows restrictions of the name-values (left and right-hand side of a relationships) to those that are allowed within the ontology. One opinion is that there shouldn't be a purpose-built extensibility point in UML, as the entire purpose of UML is that it is extensible everywhere. It also means that users of your FuGE file don't need to parse both that file and the ontology file. However, the users of your file must understand your extensibility point that you've made, which isn't useful. The extra knowledge should be stored in that ontology, in the same way as analysisXML links to mzML. One solution is to have a Property class with term “height” and a Value class with term “meters”, and a PV class with associations to both Property and Value that provides the link. In the end, this is optional. In the guidelines, these concerns should be expressed.


Other questions not fully addressed:

  • How do we find out when Generic* classes should be replaced by a specific omics-type class? Rather than using named associations to GenericParameter, GelML tended to use either GenericParameter (with a CV term) or extended the Parameter class. Is there a best way to use Parameter/GenericParameter? If it is the same shape as the Generic class, and you are just renaming it, that is a good argument for using an ontology term. However, there is less of a learning curve for users if you subclass GenericParameter with your own name. Subclassing can lock you in, and may make life more difficult further down the line if your requirements change. Remember though, hardly anyone will write XML by hand, and we shouldn't worry too much about tool implementation. Still want to make it easier for tool developers, though!

  • How should we model time?

  • For experiment-specific semantics, when should we extend the FuGE model rather than add information through the use of ontologies?

  • How descriptive should extending communities get with their samples? Should it be at the entity or attribute level? Is there a best-practice that should be documented?

  • How do we find out if two classes from two different communities are actually the same? Recurring model requirements, e.g. a library of model fragments e.g time and sample.

  • Could organism be fitted into FuGE::Bio somewhere?

  • There is no date of the Action in the ActionApplication. You could have a time parameter that comes in when you add it to your own subclass of Action/ActionApplication, and then provide a different value for that parameter in ActionApplication.

  • Somewhere, the distinction between Action and Protocol should be defined.

  • In general, we should describe a modelling best-practice to tell what is considered “standard” procedure.

  • Data package: internal versus external data

  • There may be an issue with describing physical materials within Protocols versus ProtocolApplications (theoretical materials vs physical materials, SubstanceMixtureProtocol was designed to account for this problem)


15 comments Tags: workshops, bioinformatics, systems biology, cisban, newcastle university, fuge

1st RSBI Workshop, 6-8 December 2007

  • Dec 12, 2007
  • Post a comment

Last week I attended the first RSBI (Reporting Structure for Biological Investigations) Workshop, carrying with me a multitude of hats. RSBI is a working group committed to the progression of standardization in multi-omics investigations. The purpose of the workshop was to examine and offer suggestions on the initial draft of ISA-TAB (more on that in a moment).

My first hat was a FuGE-user's hat, as the triumvirate of standards upon which RSBI is built is the Functional Genomics Experiment Model (FuGE), the Minimum Information for Biological and Biomedical Investigations (MIBBI) Project, and the Ontology for Biomedical Investigations (OBI). I was asked to give a current status update on FuGE itself, and on any communities that have already built extensions to FuGE. Andy Jones from Liverpool provided me with all of the hot-off-the-press information (my FuGE slides) - thanks Andy!

My second hat was a SyMBA-developer's hat. SyMBA uses FuGE to build a database and web front-end for storing data and experimental metadata. We use it in-house to store all of our large, high-throughput 'omics data. The use of FuGE in the system made it relevant for the workshop (my SyMBA slides, more SyMBA slides).

My final hat was a CISBAN-employee's hat. I work in the Wipat group there, and CISBAN is one of the "leading groups" involved in RSBI. As such, I was CISBAN's representative to the workshop.

The reason for the workshop, as stated earlier, was the evaluation of ISA-TAB, a proposed tabular format whose purpose is to provide a standard format for data and metadata submission into the formative BioMAP database at the EBI. ISA-TAB would have two uses:

  1. Humans: As a tabular format, it is quite easy for people to view and manipulate such templates within spreadsheet software such as Excel.
  2. Computers: As an interim solution only, ISA-TAB would be used as a computational exchange format until such time as each of the FuGE-based community extensions are complete for Metabolomics, Proteomics, and Transcriptomics. At this time, ISA-TAB would remain available for human use, but there would be a conversion step into "FuGE-ML".
The scope for ISA-TAB is large, and this was reflected in the attendees of the meeting. Representatives from ArrayExpress, Pride, and BioMAP were of course present, but also attending were people from the Metabolomics community, the MIACA project, toxico- and environmental genomics, and the FDA's NCTR.

A full write-up of the results of the workshop will soon be available online at the project's RSBI Google Group, so I'll leave it there. It was an exciting meeting, with fantastic food and even better discussions on getting public databases organized quickly for simple, straightforward multi-omics investigation data and metadata submission.

You can contact the RSBI via rsbi@googlegroups.com.

Post a comment Tags: workshops, bioinformatics, obi, systems biology, fuge, mibbi, multi-omics, rsbi …

CISB07 - Day 2

  • Nov 23, 2007
  • Post a comment

Tuesday, 20 November was the final day of the joint-CISB meeting. Again, due to the possible sensitivity of the talks, unfortunately I cannot post my notes as I normally do. For the first day, I was able to present my notes on one of the keynote speakers, but I cannot do so today.

However, while there were no more gos on the Crazy Motion Fun Ride (theme: Haunted House) as there were the night before, I do think everyone enjoyed the event. The most important thing I drew away from it was a true idea of the size of people doing systems biology research in the UK. There may be many more who do not profess a specialty in SB (but who are, in any case, working in that field), and indeed the 6 CISBs and the new SABR grantees are hardly the total sum of SB researchers in the country, but the fact remains that the close to 150 people who attended a 2-day internal meeting gives you an idea of the interest such a field has to the biological community.

Not only that, but it seems that everyone had a clear understanding of why data, ontology, and minimal information reporting standards (put them together and what do you have? RSBI!) are so relevant to their work. Systems biology produces data in so many different forms that it just makes sense to ensure that others can understand all of it.

Post a comment Tags: conferences, bioinformatics, systems biology, cisban

CISB07 - Day 1

  • Nov 19, 2007
  • Post a comment

Monday (19 November, 2007) saw the start of the two-day CISB07 conference, this year's internal conference for the 6 Centres for Integrated Systems Biology (for more information, see http://www.cisbs.org). There were also a number of invited keynote speakers. As I am unsure of the private/public status of many of the talks, I will just present some notes on those talks that I know are in the public domain.

CISB07 is being held at the Centre for Life in Newcastle. The Newcastle University CISB (CISBAN) is the host of the meeting, and I was one of about 8 people from CISBAN that helped organize it. It seems to have come together really well - it looked like everyone had a good time today. The Centre for Life is a science museum mainly for kids, but it has conferencing facilities too. I have to say that the people running the CfL today were very helpful and friendly. I'd definitely recommend the venue. After the afternoon session we had drinks on the mezzanine level of the museum, where many of the interactive (ok, and meant for children) games were. We had free reign on these games (reaction time games, soccer goal-scoring games, etc) and it was really fantastic to see all these scientists playing them with what can truly only be described as glee. Further, the CfL had their Motion Ride on for us - it's like a small cinema that plays a movie that's about 5 minutes long, and all of the chairs are on their own hydraulic pumps, and the chairs move in time with the movie. The theme tonight was "Dracula's Haunted House", and I have to admit I went on it - twice. Then it was in for dinner. However, it is the games that were really fun!

I gave a talk on SyMBA today as well, which I won't go into any more other than to say I'll put up the slides on the SourceForge site later this week, when I get a minute. As I can't be sure of the private/public status of any of the talks other than mine and the Data Integration Keynote (the keynote for the session I was in), I'll just pop up my notes from that keynote, from Dr. Chris Taylor of the EBI. It should be an interesting day tomorrow! Other than Chris' talk (notes below), I found that Michael Wilson's (from CPIB) and Mark Sansom's talks were also some of the more interesting talks of the day.

Standards development: A Two-Way Street (or, I Saw Six CISBs Come Sailing By)
(Chris Taylor, Keynote, Session 2)

There are three big omics standards bodies in biological science: MGED, PSI, MSI. HUPO's parent organization is PSI, which deals with molecular interactions, post-translational modifications, mass spec, separations, gels, etc. They also have a number of formal groups e.g. Steering Committee, working groups, and a document release process.

Many community formats are already built on FuGE-OM (GelML, mzML, analysisXML, etc). These standards provide increased efficiency (methods remain properly associated with the results generated, no need to repeatedly construct sets of contextualizing information, for industry specifically (in the light of 21 CFR Part 11)), enhanced confidence in data (enables fully-informed assessment of results, supports assessment of results that may have been generated months ago, facilitates better-informed comparisons of data sets, supports the discovery of sources of systematic or random error by correlating errors with metadata features such as the date of the operator concerned, follow-through with experiments performed), added value & tool development (re-using existing data sets for a purpose significantly different to that for which the data were generated, building aggregate data sets containing similar data from different sources, integrating data from different domains, design requirements become both explicit and stable).

Generic features of experiments; technologically and biologically-delineated sections as well. But y