Package: edgehog
edgeHOGedgeHOGis a tool to infer the gene order of each ancestor in a species phylogeny. As such,edgeHOGenables both to explore ancestral microsyntenies (local scale) and to reconstruct ancestral chromosomes (global scale).edgeHOGrelies on objects called HOGs (Hierarchical Groups of Orthologs) to model gene lineages and ancestral gene content. Basically, genes that belong to the same HOG across extant genomes are inferred to have descended from the same common ancestral gene in the common ancestor of these genomes. Accordingly, adjacencies between extant genes can be converted to edges between HOGs, which enables parsimonious ancestral gene order inferences.InstallationFrom PyPi using pipedgeHOGcan be installed directly from pypi using pip. The command is the following:pipinstalledgehogFrom sourcesedgeHOGwas built and tested with python 3.6. To set upedgeHOGon your local machine, please follow the instructions belowgitclonehttps://github.com/dessimozlab/edgehog.gitcdedgehog
pythonsetup.pyinstallUsageusage: edgehog [-h] [--version] [--output_directory OUTPUT_DIRECTORY]
               --species_tree SPECIES_TREE --hogs HOGS
               [--gff_directory GFF_DIRECTORY] [--hdf5 HDF5] [--date_edges]
               [--max_gaps MAX_GAPS]

edgehog is a software tool that infers an ancestral synteny graph at each
internal node of an input species phylogenetic tree

optional arguments:
  -h, --help            show this help message and exit
  --version             print version number and exit
  --output_directory OUTPUT_DIRECTORY
                        path to output directory (default is current
                        directory)
  --species_tree SPECIES_TREE
                        path to species/genomes phylogenetic tree (newick
                        format)
  --hogs HOGS           path to the HierarchicalGroups.orthoxml file in which
                        HOGs are stored
  --gff_directory GFF_DIRECTORY
                        path to directory with the gffs of extant genomes
                        (each gff file must be named according to the name of
                        an extant genome / leaf on the species tree)
  --hdf5 HDF5           path to the hdf5 file (alternative to gff_directory to
                        run edgeHOG on the OMA database)
  --date_edges          whether the age of edges in extant species should be
                        predictedInput dataThree types of input data are needed foredgeHOGto run:a phylogenetic tree of species/genomes of interest (in newick format)the annotation of each of these genomes (e.g. in the form of a directory of gff files)an HierarchicalGroups.orthoxml file corresponding to the extant genomesSince these input data are intersected, they must comply with the following requirements:the prefix of a gff filename must correspond to a species/genome identifier in the phylogenetic treeall genome identifiers in the phylogenetic tree must correspond to a genome entry in the HierarchicalGroups.orthoxml filetheprotIdor the prefix of theprotIdfollowed by the' 'character of each entry of a given input genome in the HierarchicalGroups.orthoxml file must match theprotein_idof a CDS in the gff file of this genomeSpecies treeA phylogenetic tree of the input genomes/species must be provided in the newick format. If internal nodes are not named, they will be named based on the concatenation of the names of their descendant leavesThe species phylogenetic tree used in the OMA database can be downloadedhereThe high-quality GTDB archaeal and bacterial species trees, along with metadata can be foundhereTo use the tree or a subtree of the NCBI taxonomy database, the ete3 python package has someuseful build-in functionsTo prune a tree in order to obtain only the phylogeny of your genomes of interest, please refer to the correspondingete3 tutorialIf you don't have a species tree available for your genomes, you can followthis tutorialon how to use OMA Browser and OMA standalone for species tree inference.HierarchicalGroups.orthoxmlAn HierarchicalGroups.orthoxml file of HOGs defined based on the proteomes and the species tree of input genomes is required for genomic comparisons. HOGs xml files can be retrieved from several leading orthology databases such asOrthoDB,EggNOG,HieranoiDBorOMA.If you don't have a HierarchicalGroups.orthoxml for your genomes, HOGs can be inferred from your input dataset usingOMA_standalone.