view hyphy_absrel.xml @ 17:0b58ef9e84b2 draft

"planemo upload for repository commit f6cabbd9eb2eea9f5d44fe9836cc5f6d6087a839"
author iuc
date Fri, 30 Oct 2020 16:26:35 +0000
parents 4cc4dbca2f02
children 9275e81b5590
line wrap: on
line source

<?xml version="1.0"?>
<tool id="hyphy_absrel" name="HyPhy-aBSREL" version="@VERSION@+galaxy0" profile="19.09">
    <description>adaptive Branch Site Random Effects Likelihood</description>
    <expand macro="requirements"/>
    <command detect_errors="exit_code"><![CDATA[
        ln -s '$input_file' absrel_input.fa &&
        ln -s '$input_nhx' absrel_input.nhx &&
        @HYPHYMPI@ absrel
            --alignment ./absrel_input.fa
            --tree ./absrel_input.nhx
            --code '$gencodeid'
            --branches '$branches'
            --output '$absrel_output'
            --multiple-hits $multiple_hits
            #if $log:
                > '$absrel_log'
            #end if
        <expand macro="inputs"/>
        <param name="log" type="boolean" truevalue="" falsevalue="" label="Output aBSREL log to history"/>
        <expand macro="gencode"/>
        <param name="multiple_hits" type="select" argument="--multiple-hits" label="Include support for multiple nucleotide substitutions">
            <option value="None">Use standard models which permit only single nucleotide changes to occur instantly</option>
            <option value="Double">Include branch-specific rates for double nucleotide substitutions</option>
            <option value="Double+Triple">Include branch-specific rates for double and triple nucleotide substitutions</option>
        <expand macro="branches"/>
        <data name="absrel_log" format="txt">
        <data name="absrel_output" format="hyphy_results.json" />
            <param name="input_file" ftype="fasta" value="absrel-in1.fa"/>
            <param name="input_nhx" ftype="nhx" value="absrel-in1.nhx"/>
            <param name="multiple_hits" value="None" />
            <param name="log" value="false"/>
            <output name="absrel_output" file="absrel-out1.json" compare="sim_size"/>
            <param name="input_file" ftype="fasta" value="absrel-in1.fa"/>
            <param name="input_nhx" ftype="nhx" value="absrel-in1.nhx"/>
            <param name="multiple_hits" value="Double" />
            <param name="log" value="false"/>
            <output name="absrel_output" file="absrel-out1.json" compare="sim_size"/>

aBSREL : Adaptive Branch-Site Random Effects Likelihood

What question does this method answer?

Is some proportion of sites is subject to positive selection along specific branches or lineages of a phylogeny?

Recommended Applications

1. Exploratory testing for evidence of lineage-specific positive diversifying selection in small to medium sized alignments (up to 100 sequences).
2. Targeted testing of branches selected a priori for positive diversifying selection, including alignments with prohibitive runtimes under older branch-site models (up to 1,000 sequences)

Brief description

aBSREL (Adaptive branch-site random effects likelihood) uses an adaptive
random effects branch-site model framework to test whether each branch
has evolved under positive selection, using a procedure which infers an
optimal number of rate categories per branch.


1. A *FASTA* sequence alignment.
2. A phylogenetic tree in the *Newick* format

Note: the names of sequences in the alignment must match the names of the sequences in the tree.


A JSON file with analysis results (

For each tested branch the analysis will infer the appropriate number of selective regimes,
and whether or not there is statistical evidence of positive selection on that branch.

A custom visualization module for viewing these results is available (see for an example)

Further reading

Tool options

    --code              Which genetic code to use

    --branches          Which branches should be tested for selection?
                            All [default] : test all branches

                            Internal : test only internal branches (suitable for
                            intra-host pathogen evolution for example, where terminal branches
                            may contain polymorphism data)

                            Leaves: test only terminal (leaf) branches

                            Unlabeled: if the Newick string is labeled using the {} notation,
                            test only branches without explicit labels

    <expand macro="citations">
        <citation type="doi">10.1093/molbev/msv022</citation>