Package: xebec
xebecSnakemake pipeline for microbiome diversity effect size benchmarkingNOTE: Please note that xebec is still under active development.InstallationTo use xebec, you will need several dependencies.
We recommend usingmambato install these packages when possible.mamba install -c conda-forge -c bioconda biom-format h5py==3.1.0 snakemake pandas unifrac scikit-bio bokeh unifrac-binaries jinja2

pip install evident>=0.4.0 gemelli>=0.0.8To install xebec, run the following command from the command line:pip install xebecUsageIf you runxebec --help, you should see the following:$ xebec --help
Usage: xebec [OPTIONS]

Options:
  --version                       Show the version and exit.
  -ft, --feature-table PATH       Feature table in BIOM format.  [required]
  -m, --metadata PATH             Sample metadata in TSV format.  [required]
  -t, --tree PATH                 Phylogenetic tree in Newick format.
                                  [required]
  -o, --output PATH               Output workflow directory.  [required]
  --max-category-levels INTEGER   Max number of levels in a category.
                                  [default: 5]
  --min-level-count INTEGER       Min number of samples per level per
                                  category.  [default: 3]
  --rarefy-percentile FLOAT       Percentile of sample depths at which to
                                  rarefy.  [default: 10]
  --n-pcoa-components INTEGER     Number of PCoA components to compuate.
                                  [default: 3]
  --validate-input / --no-validate-input
                                  Whether to validate input before creating
                                  workflow.  [default: validate-input]
  --help                          Show this message and exit.To create the workflow structure, pass in the filepaths for the feature table, sample metadata, and phylogenetic tree.
You must also pass in a path to a directory in which to create the workflow.
Additionally, you can provide parameters for determining how to process your sample metadata.After running this command, navigate inside the output directory you created.
There should be two subdirectories:workflow/andconfig/.To start the pipeline , run the following command:snakemake --cores 1You should see the Snakemake pipeline start running the jobs.
If this pipeline runs sucessfully, the processed results will be located atresults/.
Included in the results are the concatenated effect size values as well as interactive plots summarizing the effect sizes for each metadata column for each diversity metric.
These plots are generated usingBokehand can be visualized in any modern web browser.Workflow Overviewxebec performs four main steps, some of which have substeps.Process data (filter metadata, rarefaction)Run diversity analysesCalculate effect sizes (concatenate together)Generate visualizationsAn overview of the DAG is shown below:ConfigurationDiversity Metricsxebec allows configuration of what alpha and beta diversity metrics are included in the workflow.
To add or remove metrics, modify theconfig/alpha_div_metrics.ymlandconfig/beta_div_metrics.ymlfiles.
For alpha diversity, any metric that can be passed intoskbio.alpha_diversityshould work.
For beta diversity, any non-phylogenetic metric that can be passed intoskbio.beta_diversityshould work.
Valid phylogenetic beta diversity are those that can be passed intoStriped UniFrac.
Make sure that any additional diversity metrics are annotated withphyloornon_phyloso xebec knows how to process them.Snakemake OptionsThe xebec workflow can be decorated with many configuration options available in Snakemake, including resource usage and HPC scheduling.
We recommend reading through theSnakemake documentationfor details on these options.
Note that some of these options may require creating new configuration files.