/Open_Orthobench

Primary LanguagePython

Orthobench Revisited

This is the GitHub repository for a complete reanalysis of the Orthobench benchmarks for orthogroup inference accuracy. It contains 70 curated Bilaterian orthogroups based on the orthogroups from the original Orthobench study.

There have been a number of major corrections to the curated orthogroups, which together improve the accuracy of the benchmarks considerably. The updated orthogroups together with the complete set of data supporting the analysis are provided in this repository. The issues page is open, so anyone can identify any further corrections that are required and can submit the data supporting the correction.

Citation

If using this work please cite:

D M Emms, S Kelly, Benchmarking Orthogroup Inference Accuracy: Revisiting Orthobench, Genome Biology and Evolution, evaa211, https://doi.org/10.1093/gbe/evaa211

As well as the original study:

Trachana, K., Larsson, T.A., Powell, S., Chen, W.‐H., Doerks, T., Muller, J. and Bork, P. (2011), Orthology prediction methods: A quality assessment using curated protein families. Bioessays, 33: 769-780. doi:10.1002/bies.201100062

Benchmarking an Orthogroup Inference Method

  1. Download BENCHMARKS.tar.gz from https://github.com/davidemms/Open_Orthobench/releases
  2. Run the method on the proteomes in BENCHMARKS/Input/
  3. Write the predicted orthogroups to a file, one orthogroup per line
  4. Run the script on the results file: python benchmark.py orthogroups_results_file.txt

The predicted orthogroups file should have one orthogroup per line, the script uses regular expressions to extract the genes from each line for a wide variety of formats. Lines starting with '#' are ignored.

Description of Files

There are two main directories

  • BENCHMARKS: Contains all the files you need to benchmark an orthogroup inference method.
  • Supporting_Data: Contains all the data supporting the curation of the 70 Bilaterian orthogroups (RefOGs). If you want to critique the curation of the orthogroups that form the benchmark then these files explain the reasoning behind the exclusion/inclusion of genes in each RefOG.

The contents of this directory are as follows.

BENCHMARKS - data & script for benchmarking orthogroup inference methods

  • Input - the 12 input proteomes for the orthogroup inference method to be tested
  • RefOGs - the genes in each RefOG as determined by this study
  • benchmark.py - python script to calculate the benchmarks for a user-supplied set of inferred orthogroups.

Supporting_Data - Data supporting the inferred RefOGs

This contains two subdirectories

Data_for_RefOGs - The data associated with the inference of each RefOG

  • sequences - the sets of sequences used for tree inference
  • alignments - the MSA for each sequence file
  • trees - the trees inferred from the MSAs
  • trees_tight - a tree for each RefOG cropped tightly around the RefOG
  • evidence - the evidence considered for each RefOG in determining which genes belong to each RefOG
  • low_certainty_assignments - those genes that have been included/excluded from a RefOG but for which this could only be done with low certainty

Additional_Files

  • proteomes - the proteomes for the 12 study species plus additional in-group and out-group species used to help determine the gene membership for each RefOG
  • original_study_trees - gene trees from the original study
  • hmm_profiles - hmm profiles from the original study for the RefOGs
  • hmmer_results - HMMER results for searching the hmm profiles against the proteomes used for this study
  • additional_identified_orthogroups - Additionally identified orthogroups determined in the process of determining the membership of the target orthogroups