changeset 1:28d4e824d3d3 draft

Deleted selected files
author sigven
date Tue, 27 Sep 2022 15:31:26 +0000
parents fb035154d720
children 3c61ef74a176
files oncoenrichr_wrapper.xml
diffstat 1 files changed, 0 insertions(+), 391 deletions(-) [+]
line wrap: on
line diff
--- a/oncoenrichr_wrapper.xml	Thu Sep 22 11:35:13 2022 +0000
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,391 +0,0 @@
-<tool id="oncoenrichr_wrapper" name="oncoEnrichR" version="1.3.1">
-  <description>Cancer-dedicated gene set interpretation</description>
-    <requirements>
-        <container type="docker">sigven/oncoenrichr:1.3.1</container>
-  </requirements>
-  <command detect_errors="aggressive"><![CDATA[
-      #if $query_set.query_choice.query_input == "text"
-        echo $query_set.query_choice.query_text | sed 's/__cn__/\n/g' > query_text.csv &&
-        #set input_file = './query_text.csv'
-      #else if $query_set.query_choice.query_input == "file"
-        ln -s $query_set.query_choice.query_file "$query_set.query_choice.query_file.element_identifier" &&
-        #set input_file = './' + str($query_set.query_choice.query_file.element_identifier)
-      #end if
-
-      #set background_file = ''
-      #if $fun_enrich.custom_bgset.def_background
-	      #if $fun_enrich.custom_bgset.bg_choice.bg_source == "text"
-	        echo $fun_enrich.custom_bgset.bg_choice.bg_enrich_text | sed 's/__cn__/\n/g' > custom_bgset.csv &&
-	        #set background_file = './custom_bgset.csv'
-	      #else if $fun_enrich.custom_bgset.bg_choice.bg_source == "file" and $fun_enrich.custom_bgset.bg_choice.bg_enrich_file
-	        ln -s $fun_enrich.custom_bgset.bg_choice.bg_enrich_file background_text.csv &&
-	        #set background_file = './custom_bgset.csv'
-	      #else
-	        #set background_file = ''
-	      #end if
-	 #end if
-
-      R -e 'suppressPackageStartupMessages(library(oncoEnrichR));
-      suppressWarnings(load(system.file("internal_db", "oedb.rda", package = "oncoEnrichR")));
-      gene_data <- read.csv("$input_file", stringsAsFactors = F, header = F);
-      oe_report <- oncoEnrichR::onco_enrich(
-        query = gene_data[[1]],
-        oeDB = oedb,
-      #if $query_set.query_id_type
-        query_id_type = "$query_set.query_id_type",
-      #end if
-      ignore_id_err = $query_set.ignore_id_err,
-
-      #if $report_metadata.project_title
-          project_title = "$report_metadata.project_title",
-      #end if
-      #if $report_metadata.project_owner
-          project_owner = "$report_metadata.project_owner",
-      #end if
-      #if $report_metadata.project_description
-          project_description = "$report_metadata.project_description",
-      #end if
-
-      show_enrichment = $modules.show_enrichment,
-      show_ppi = $modules.show_ppi,
-      show_disease = $modules.show_disease,
-      show_cancer_hallmarks = $modules.show_cancer_hallmarks,
-      show_drug = $modules.show_drug,
-      show_aberration = $modules.show_aberration,
-      show_coexpression = $modules.show_coexpression,
-      show_subcell_comp = $modules.show_subcell_comp,
-      show_complex = $modules.show_complex,
-	 show_domain = $modules.show_domain,
-      show_fitness = $modules.show_fitness,
-	 show_cell_tissue = $modules.show_cell_tissue,
-      show_ligand_receptor = $modules.show_ligand_receptor,
-      show_regulatory = $modules.show_regulatory,
-	 show_prognostic = $modules.show_prognostic,
-	 show_unknown_function = $modules.show_unknown_function,
-      show_synleth = $modules.show_synleth,
-
-      #if $background_file
-          bgset = read.csv("$background_file", stringsAsFactors = F, header = F)[[1]],
-		#if $fun_enrich.custom_bgset.bg_enrich_id_type
-            bgset_id_type = "$fun_enrich.custom_bgset.bg_enrich_id_type",
-          #end if
-          #if $fun_enrich.custom_bgset.bg_enrich_description
-            bgset_description = "$fun_enrich.custom_bgset.bg_enrich_description",
-          #end if
-      #else
-          bgset = NULL,
-      #end if
-
-      #if $fun_enrich.p_value_cutoff_enrichment
-          p_value_cutoff_enrichment = $fun_enrich.p_value_cutoff_enrichment,
-      #end if
-      #if $fun_enrich.p_value_adjustment_method
-          p_value_adjustment_method = "$fun_enrich.p_value_adjustment_method",
-      #end if
-      #if $fun_enrich.q_value_cutoff_enrichment
-          q_value_cutoff_enrichment = $fun_enrich.q_value_cutoff_enrichment,
-      #end if
-      #if $fun_enrich.min_geneset_size
-          min_geneset_size = $fun_enrich.min_geneset_size,
-      #end if
-      #if $fun_enrich.max_geneset_size
-          max_geneset_size = $fun_enrich.max_geneset_size,
-      #end if
-
-      #if $protein_interactions.ppi_add_nodes
-          ppi_add_nodes = $protein_interactions.ppi_add_nodes,
-      #end if
-      #if $protein_interactions.ppi_score_threshold
-          ppi_score_threshold = $protein_interactions.ppi_score_threshold,
-      #end if
-      show_drugs_in_ppi = $protein_interactions.show_drugs_in_ppi,
-	 ppi_node_shadow = $protein_interactions.ppi_node_shadow,
-
-	 #if $subcellular_compartments.min_subcellcomp_confidence
-          min_subcellcomp_confidence = $subcellular_compartments.min_subcellcomp_confidence,
-      #end if
-      #if $fitness.max_fitness_score
-          max_fitness_score = $fitness.max_fitness_score,
-      #end if
-      subcellcomp_show_cytosol = $subcellular_compartments.show_cytosol,
-      #if $disease.show_top_diseases_only
-          show_top_diseases_only = $disease.show_top_diseases_only,
-      #end if
-
-      min_confidence_reg_interaction = "$regulatory.min_confidence_reg_interaction",
-      num_terms_enrichment_plot = $fun_enrich.num_terms_enrichment_plot,
-      simplify_go = $fun_enrich.simplify_go,
-      html_floating_toc = $report_metadata.html_floating_toc,
-      html_report_theme = "$report_metadata.html_report_theme",
-      galaxy = TRUE
-      );
-
-      oncoEnrichR::write(report = oe_report, oeDB = oedb, file = "$report1", format = "html", selfcontained_html = F, extra_files_path = "$report1.extra_files_path", overwrite = T, ignore_file_extension = T);
-	 oncoEnrichR::write(report = oe_report, oeDB = oedb, file = "$report2", format = "excel", overwrite = T, ignore_file_extension = T)' 2>&1
-
-  ]]></command>
-  <inputs>
-      <section title="" name=""/>
-      <section name="query_set" title="Query gene set" expanded="true">
-          <conditional name="query_choice">
-                <param name="query_input" type="select" multiple="false" display="radio"
-                       label="Query gene set: do you want to upload a file OR paste into a text box?">
-                    <option value="text">Text field</option>
-				<option value="file">From file</option>
-                </param>
-                <when value="text">
-                    <param type="text" name="query_text" label="Query gene set identifiers (one per line)" area="true"/>
-                </when>
-			 <when value="file">
-                   <param name="query_file" type="data" format="txt" label="Query gene set identifiers" multiple="false"/>
-                </when>
-          </conditional>
-          <param name="query_id_type" type="select" label="Query identifier type" display="radio" multiple="false">
-              <option value="symbol">Primary gene symbol (HGNC) - e.g. KRAS</option>
-              <option value="uniprot_acc">UniProt accession - e.g. P01116</option>
-              <option value="entrezgene">NCBI Entrez gene identifier - e.g. 3845</option>
-              <option value="ensembl_gene">Ensembl gene identifier - e.g. ENSG00000133703</option>
-              <option value="ensembl_mrna">Ensembl transcript identifier - e.g. ENST00000311936</option>
-              <option value="ensembl_protein">Ensembl protein identifier - e.g. ENSP00000308495</option>
-              <option value="refseq_mrna">RefSeq mRNA identifier - e.g. NM_004985</option>
-              <option value="refseq_protein">RefSeq protein identifier - e.g. NP_004976</option>
-          </param>
-		<param name="ignore_id_err" type="boolean" label="Ignore erroneous idenfiers" truevalue="T" falsevalue="F" checked="true"/>
-      </section>
-
-      <section title="" name=""/>
-      <section name="report_metadata" title="Project metadata and output settings" expanded="true">
-          <param type="text" name="report_name" label="Output filename (prefix)" value="Report"/>
-          <param type="text" name="project_title" label="Project title" />
-          <param type="text" name="project_owner" label="Project owner" />
-          <param type="text" name="project_description" label="Project description" area="true"/>
-          <param name="html_floating_toc" type="boolean" label="HTML report - float the table of contents to the left of the main document content (always visible during scrolling)" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="html_report_theme" type="select" label="HTML report - bootswatch theme" expanded="true">
-              <option value="default">default</option>
-              <option value="cerulean">cerulean</option>
-              <option value="cosmo">cosmo</option>
-              <option value="journal">journal</option>
-              <option value="lumen">lumen</option>
-              <option value="paper">paper</option>
-              <option value="sandstone">sandstone</option>
-              <option value="simplex">simplex</option>
-              <option value="spacelab">spacelab</option>
-              <option value="united">united</option>
-              <option value="yeti">yeti</option>
-          </param>
-      </section>
-
-      <section title="" name=""/>
-      <section name="modules" title="Analysis modules included in the report" expanded="true">
-		<param name="show_disease" type="boolean" label="Gene-cancer associations" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_enrichment" type="boolean" label="Gene functional enrichment" truevalue="T" falsevalue="F" checked="true"/>
-		<param name="show_cell_tissue" type="boolean" label="Tissue/cell-type enrichment" truevalue="T" falsevalue="F" checked="false"/>
-          <param name="show_ppi" type="boolean" label="Protein-protein interaction network" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_regulatory" type="boolean" label="Regulatory (TF-target) interactions" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_ligand_receptor" type="boolean" label="Ligand-receptor interactions" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_cancer_hallmarks" type="boolean" label="Cancer hallmark associations" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_drug" type="boolean" label="Drug-target associations" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_aberration" type="boolean" label="Tumor aberration frequencies" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_coexpression" type="boolean" label="Tumor co-expression patterns" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_subcell_comp" type="boolean" label="Subcellular localizations" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_complex" type="boolean" label="Protein complex memberships" truevalue="T" falsevalue="F" checked="true"/>
-		<param name="show_domain" type="boolean" label="Protein domain frequencies" truevalue="T" falsevalue="F" checked="false"/>
-          <param name="show_fitness" type="boolean" label="Gene fitness effects" truevalue="T" falsevalue="F" checked="true"/>
-          <param name="show_synleth" type="boolean" label="Predicted synthetic lethality interactions" truevalue="T" falsevalue="F" checked="true"/>
-		<param name="show_unknown_function" type="boolean" label="Genes of poorly defined function" truevalue="T" falsevalue="F" checked="true"/>
-		<param name="show_prognostic" type="boolean" label="Prognostic cancer associations" truevalue="T" falsevalue="F" checked="true"/>
-      </section>
-
-      <section title="" name=""/>
-      <section name="fun_enrich" title="Options - gene functional enrichment" expanded="true">
-		 <conditional name="custom_bgset">
-			 <param name="def_background" type="boolean" label="Define custom background set (all annotated protein-coding genes by default)" truevalue="T" falsevalue="F" checked="false"/>
-			 <when value="T">
-		            <conditional name="bg_choice">
-		                <param name="bg_source" type="select" display="radio"
-		                       label="Custom background gene set: do you want to upload a file OR paste into a text box?">
-						<option value="text">Text field</option>
-		                    <option value="file">From file</option>
-
-		                </param>
-		                <when value="file">
-		                   <param type="data" format="txt" name="bg_enrich_file" label="Custom background gene set" optional="true" multiple="false"/>
-		                </when>
-		                <when value="text">
-		                    <param type="text" name="bg_enrich_text" label="Custom background gene set identifiers (one per line):" area="true"/>
-		                </when>
-		          </conditional>
-
-		          <param type="select" name="bg_enrich_id_type" label="Custom background identifier type" display="radio" multiple="false">
-		              <option value="symbol">Primary gene symbol (HGNC) - e.g. KRAS</option>
-                      <option value="uniprot_acc">UniProt accession - e.g. P01116</option>
-                      <option value="entrezgene">NCBI Entrez gene identifier - e.g. 3845</option>
-                      <option value="ensembl_gene">Ensembl gene identifier - e.g. ENSG00000133703</option>
-                      <option value="ensembl_mrna">Ensembl transcript identifier - e.g. ENST00000311936</option>
-                      <option value="ensembl_protein">Ensembl protein identifier - e.g. ENSP00000308495</option>
-                      <option value="refseq_mrna">RefSeq mRNA identifier - e.g. NM_004985</option>
-                      <option value="refseq_protein">RefSeq protein identifier - e.g. NP_004976</option>
-		          </param>
-		          <param type="text" name="bg_enrich_description" label="Custom background gene set description" value="Custom background description"/>
-			</when>
-		</conditional>
-
-          <param type="float" name="p_value_cutoff_enrichment" label="P-value cutoff for enrichment tests (clusterProfiler)" value="0.05"/>
-          <param type="select" name="p_value_adjustment_method" label="P-value adjustment method (clusterProfiler)">
-              <option value="holm">holm</option>
-              <option value="hochberg">hochberg</option>
-              <option value="hommel">hommel</option>
-              <option value="bonferroni">bonferroni</option>
-              <option value="BH">BH</option>
-              <option value="BY">BY</option>
-              <option value="fdr">fdr</option>
-              <option value="none">none</option>
-          </param>
-          <param type="float" name="q_value_cutoff_enrichment" label="Q-value cutoff for enrichment tests to report as significant (clusterProfiler)" value="0.2"/>
-          <param type="integer" name="min_geneset_size" label="Minimum number of genes annotated by ontology term for testing (clusterProfiler)" value="10"/>
-          <param type="integer" name="max_geneset_size" label="Maximum number of genes annotated by ontology term for testing (clusterProfiler)" value="500"/>
-          <param name="simplify_go" type="boolean" label="Simplify GO enrichment results by removal of redundant terms (recommended)" truevalue="T" falsevalue="F" checked="true"/>
-          <param type="integer" name="num_terms_enrichment_plot" label="Number of top enriched Gene Ontology terms (max) to show in enrichment barplot" min="10" max="30" value="20"/>
-      </section>
-
-      <section title="" name=""/>
-      <section name="fitness" title="Options - gene fitness scores" expanded="true">
-          <param type="float" name="max_fitness_score" label="Maximum loss-of-fitness score (Bayes Factor from BAGEL) for genes retrieved from Project Score" value="-2" min="-5" max="0"/>
-      </section>
-      <section title="" name=""/>
-      <section name="protein_interactions" title="Options - protein-protein interaction network" expanded="true">
-          <param type="integer" name="ppi_add_nodes" label="Addition of interacting non-queryset proteins to the protein-protein interaction network (maximum number)" value="50" min="0" max="50"/>
-          <param type="integer" name="ppi_score_threshold" label="Minimum confidence score for interactions to be included in the network (STRING confidence: 0-1000)" value="900" min="400" max="1000"/>
-          <param name="show_drugs_in_ppi" type="boolean" label="Show anti-cancer drugs in protein-protein interaction network" truevalue="T" falsevalue="F" checked="true"/>
-		<param name="ppi_node_shadow" type="boolean" label="Add shadow to nodes in protein-protein interaction network" truevalue="T" falsevalue="F" checked="true"/>
-      </section>
-      <section title="" name=""/>
-      <section name="regulatory" title="Options - regulatory interactions" expanded="true">
-         <param type="select" name="min_confidence_reg_interaction" label = "Minimum confidence level of regulatory interactions included (DoRothEA - A:highest, D:lowest)">
-            <option value="D">D</option>
-            <option value="C">C</option>
-            <option value="B">B</option>
-            <option value="A">A</option>
-        </param>
-      </section>
-     <section title="" name=""/>
-
-	 <section name="subcellular_compartments" title="Options - Subcellular localizations" expanded="true">
-           <param type="integer" name="min_subcellcomp_confidence" label="Minimum confidence level for subcellular localization annotations" value="1" min="1" max="6"/>
-          <param name="show_cytosol" type="boolean" label="Show cytosol annotations (very common localization) in subcellular heatmap " truevalue="T" falsevalue="F" checked="false"/>
-      </section>
-      <section title="" name=""/>
-
-      <section name="disease" title="Options - Disease associations" expanded="true">
-          <param type="boolean" name="show_top_diseases_only" label="Show top disease assocations only" truevalue="T" falsevalue="F" checked="true"/>
-      </section>
-
-  </inputs>
-    <outputs>
-        <data format="xlsx" name="report2" label="$report_metadata.report_name - xlsx"/>
-        <data format="html" name="report1" label="$report_metadata.report_name - html"/>
-    </outputs>
-
-
-  <help><![CDATA[
-.. class:: infomark
-
-The query gene set is limited to n = 500 identifiers. A limited query gene set (e.g. n < 5) will in general reduce the relevance and significance of many oncoEnrichR report modules.
-
------
-
-**Dataset formats**
-
-The input dataset is in tabular_ format. The two output datasets are html_ and xlsx.
-
-.. _tabular: ${static_path}/formatHelp.html#tab
-.. _html: ${static_path}/formatHelp.html#html
-
------
-
-**What it does**
-
-*OncoEnrichR* is intended for exploratory analysis and prioritization of a candidate hits (referred to as *query set* below) from high-throughput cancer biology experiments. The tool queries a number of high-quality data resources in order to interpret the query gene set along various dimensions, examples being cancer aberration frequencies, protein-protein interactions, pathway enrichment, subcellular compartment localization, target druggability, gene fitness scores, and tissue/cell-type specificity.
-
-The results from the various analysis modules are provided in an interactive HTML report where the user can interrogate the results further. A multisheet Excel workbook is also provided for convience. The following resources are currently utilized for annotation and analysis:
-
--  `Open Targets Platform <https://targetvalidation.org/>`_ - disease associations, drug-target associations, cancer hallmarks, and druggability/tractability rankings
-
--  `The Cancer Genome Atlas <https://portal.gdc.cancer.gov/>`_ - gene aberration frequencies and co-expression patterns in approximately 10,000 primary tumor samples
-
--  `The Human Protein Atlas <https://www.proteinatlas.org/>`_ - expression data for healthy human tissues (`GTex <https://gtexportal.org/home/>`_)/cell types, and prognostic gene expression associations in cancer (`The Pathology Atlas <https://www.proteinatlas.org/humanproteome/pathology/>`_)
-
--  `Molecular Signatures Database (MSigDB) <http://software.broadinstitute.org/gsea/msigdb/index.jsp/>`_ - collection of annotated (e.g. towards pathways) gene sets for enrichment/overrepresentation analysis. This includes gene sets from `Gene Ontology <http://geneontology.org/>`_, `Reactome <https://reactome.org/>`_, `KEGG <https://www.genome.jp/kegg/pathway.html/>`_, `WikiPathways <https://www.wikipathways.org/index.php/WikiPathways/>`_, `BIOCARTA <https://maayanlab.cloud/Harmonizome/dataset/Biocarta+Pathways/>`_, as well as curated `immunologic <https://www.gsea-msigdb.org/gsea/msigdb/collections.jsp#C7/>`_ and `cancer-specific <https://www.gsea-msigdb.org/gsea/msigdb/collections.jsp#C6/>`_ signatures.
-
--  `NetPath <http://www.netpath.org/>`_ - manually curated resource of signal transduction pathways in humans
-
--  `STRING <https://string-db.org/>`_ - protein-protein interaction database
-
--  `CellChatDB <http://www.cellchat.org/>`_ - database on ligand-receptor interactions
-
--  `DoRothEA <https://saezlab.github.io/dorothea/>`_ - gene set resource containing signed transcription factor (TF) - target interactions
-
--  `CORUM <https://mips.helmholtz-muenchen.de/corum/>`_ - protein complex database
-
--  `Compleat <https://fgr.hms.harvard.edu/compleat>`_ - protein complex resource
-
--  `ComplexPortal <https://www.ebi.ac.uk/complexportal/home/>`_ - manually curated, encyclopaedic resource of macromolecular complexes
-
--  `hu.MAP2 <http://humap2.proteincomplexes.org/>`_ - human protein complex map
-
--  `ComPPI <http://comppi.linkgroup.hu/>`_ - subcellular compartment database
-
--  `CancerMine <http://bionlp.bcgsc.ca/cancermine/>`_ - literature-mined resource on cancer drivers, oncogenes and tumor suppressor genes
-
--  `Network of Cancer Genes <http://ncg.kcl.ac.uk/>`_ - manually curated collection of cancer genes, healthy drivers and their properties
-
--  `Project Score <https://score.depmap.sanger.ac.uk/>`_ - database on the effects on cancer cell line viability elicited by CRISPR-Cas9 mediated gene activation
-
--  `Genetic determinants of survival in cancer <http://survival.cshl.edu/>`_ - resource on the prognostic impact of genetic aberrations (methylation, CNA, mutation, expression) in human cancers (TCGA)
-
--  `Predicted synthetic lethality interactions <https://pubmed.ncbi.nlm.nih.gov/34529928/>`_ - comprehensive prediction of synthetic lethality interactions in human cancer cell lines
-
-The contents of the gene set analysis report attempt to answer the following questions related to the query set:
-
--  Which diseases/tumor types are known to be associated with genes in the query set, and to what extent? Which genes are a classified as proto-oncogenes, tumor suppressors or cancer driver genes?
-
--  Which query genes have been linked (through literature) to the various hallmarks of cancer?
-
--  Which genes in the query set are poorly characterized or have an unknown function?
-
--  Which proteins in the query set can be targeted by inhibitors for diffferent cancer conditions (early and late clinical development phases)? What is the tractability/druggability status for other targets in the query set?
-
--  Which cancer-relevant protein complexes are involved for proteins in the query set?
-
--  Are there known cancer-relevant regulatory interactions (transcription factor (TF) - target) found in the query set?
-
--  Are there known ligand-receptor interactions in the query set?
-
--  Which subcellular compartments (nucleus, cytosol, plasma membrane etc.) are dominant localizations for members of the query set?
-
--  Are specific tissues or cell types enriched in the query set, considering healthy tissue/cell-type specific expression patterns (GTex/Human Protein Atlas) of query genes?
-
--  Which protein-protein interactions are known within the query set? Are there interactions between members of the query set and other cancer-relevant proteins (e.g. proto-oncogenes, tumor-suppressors or predicted cancer drivers)? Which proteins constitute hubs in the protein-protein interaction network?
-
--  Are there specific pathways, biological processes or molecular functions that are enriched within the query set, as compared to a reference/background set?
-
--  Which members of the query set are frequently mutated in tumor sample cohorts (TCGA - SNVs/InDels / homozygous deletions / copy number amplifications)? What are the most frequent recurrent somatic variants (SNVs/InDels) in the query set genes?
-
--  Which members of the query set are co-expressed (strong negative or positive correlations) with cancer-relevant genes (i.e. proto-oncogenes or tumor suppressors) in tumor sample cohorts (TCGA)?
-
--  Which members of the query set are associated with better/worse survival in different cancers, considering mutation, expression, methylation or copy number levels in tumors?
-
--  Which members of the query set are predicted as partners of synthetic lethality interactions?
-
--  Which members of the query set are associated with cellular loss-of-fitness in CRISPR/Cas9 whole-genome drop out screens of cancer cell lines (i.e. reduction of cell viability elicited by a gene inactivation)? Which genes should be prioritized considering genomic biomarkers and fitness scores in combination?
-
-
-]]>
-  </help>
-
- <citations>
-     <!-- Example of annotating a citation using a DOI. -->
-     <citation type="doi">10.48550/arXiv.2107.13247</citation>
-     <!-- Example of annotating a citation using a BibTex entry. -->     
-  </citations>
-</tool>