| Name | SyntenyQC JSON |
| Version |
1.1
JSON |
| download |
| home_page | None |
| Summary | A command line app for the pre-processing and quality control of synteny plot neighbourhoods |
| upload_time | 2024-08-10 09:54:03 |
| maintainer | None |
| docs_url | None |
| author | None |
| requires_python | >=3.10 |
| license | MIT License Copyright (c) 2024 Timothy David John Kirkwood Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. |
| keywords |
bgc
biosynthetic gene cluster
genome neighbourhood
neighbourhood
synteny
synteny plot
|
| VCS |
 |
| bugtrack_url |
|
| requirements |
No requirements were recorded.
|
| Travis-CI |
No Travis.
|
| coveralls test coverage |
No coveralls.
|
# NOTE - this app has only been tested on the Windows OS
# SyntenyQC
## Motivation:
Synteny plots are widely used for the comparison of genomic neighbourhoods. Whilst synteny plots are often included as part of larger software suites (e.g. the `antiSMASH` ClusterBlast module), various low-code, stand-alone tools are now available that allow users to source candidate neighbourhoods and build their own synteny plots. However, a gap remains between:
**(i) tools that source these candidate neighbourhoods** (e.g. `cblaster`, which can find hundreds of candidates),
**(ii) tools that build the synteny plots** (e.g. `clinker`, which struggles as the number of neighbourhoods exceeds 30-50) and
**(iii) the synteny plots themselves**, which become much harder to analyse/present as the number of neighbourhoods they include increases.
## Description:
`SyntenyQC` is a python app for the curation of neighbourhoods immediately prior to synteny plot creation. `SyntenyQC collect` supports the systematic definition and annotation of candidate neighbourhoods based on a direct integration to `cblaster`. `SytenyQC sieve` offers a flexible method for objectively removing redundant neighbourhoods (sourced using `cblaster` or any other tool) prior to synteny plot creation. This is in some cases an absolute requirement (e.g. `cblaster` called via the `CAGECAT` webserver places a limit of 50 neighbourhoods).
## Installation
```
pip install SyntenyQC
```
Note - `SyntenyQC` depends on [BLAST+](https://blast.ncbi.nlm.nih.gov/doc/blast-help/downloadblastdata.html), which must be [installed](https://www.ncbi.nlm.nih.gov/books/NBK569861/) by the user (tested with v2.12.0 - but should work with other versions unless there are parameter changes). If this is installed correctly, you should be able to see help messages after typing `blastp -h` and `makeblastdb -h` in the command line.
## Tests
Tests are performed using [pytest](https://pypi.org/project/pytest/), but are not distributed with `SyntenyQC`. To run tests:
1) Install pytest
2) Clone the `SyntenyQC` github repository.
3) Update `path/to/cloned/repository/tests/email.txt` with your email (if left blank, two webscraping tests will fail).
4) Navigate to the cloned repository via command line (on Windows, type `cd path/to/cloned/repository`).
5) Type `pytest` in the command line.
## Usage
### General help:
```
>SyntenyQC -h
usage: SyntenyQC [-h] {collect,sieve} ...
options:
-h, --help show this help message and exit
subcommands:
{collect,sieve} Synteny quality control options
```
### Collect subcommand:
```
>SyntenyQC collect -h
usage: SyntenyQC collect [-h] -bp -ns -em [-fn] [-sp] [-wg]
Write genbank files corresponding to cblaster neighbourhoods from a specified CSV-format binary file located at
BINARY_PATH. For each cblaster hit accession in the binary file:
1) A record is downloaded from NCBI using the accession. NCBI requires a user EMAIL to search for this record
programatically. If WRITE_GENOMES is specified, this record is written to a local file according to FILENAMES
(see final bulletpoint).
2) A neighbourhood of size NEIGHBOURHOOD_SIZE bp is defined, centered on the cblaster hits defined in the binary
file for the target accession.
3) (If STRICT_SPAN is specified:) If the accession's record is too small to contain a neighbourhood of the
desired size, it is discarded. For example, if an accession record is a 25kb contig and NEIGHBOURHOOD_SIZE
is 50000, the record is discarded.
4) If FILENAMES is "organism", the nighbourhood is written to file called *organism*.gbk. If FILENAMES is
"accession", the neighbourhood is written to *accession*.gbk. Synteny softwares such as clinker can use these
filenames to label synetny plot samples.
Once COLLECT has been run, a new folder with the same name as the binary file will be seen in the directory
that holds the binary file (i.e. the file "path/to/binary/file.txt" will generate the folder "path/to/binary/file").
This folder will have a subdirectory called "neighbourhood", containing all of the neighbourhood genbank files
(i.e. "path/to/binary/file/neighbourhood"). If WRITE_GENOMES is specified, a second direcory ("genome") will also
be present, containing the entire record associated with each cblaster accession (i.e. "path/to/binary/file/genome").
Finally, a log file will be present in the folder "path/to/binary/file", containing all run details.
options:
-h, --help show this help message and exit
-bp, --binary_path
Full filepath to the CSV-format cblaster binary file containing neighbourhoods that should
be extracted
-n, --neighbourhood_size
Size (basepairs) of neighbourhood to be extracted (centered on middle of CBLASTER-defined
neighbourhood)
-em, --email Email - required for NCBI entrez querying
-fn, --filenames
If "organism", all collected files will be named according to organism. If "accession", all
files will be named by NCBI accession. (default:
organism)
-sp, --strict_span
If set, will discard all neighbourhoods that are smaller than neighbourhood_size bp. For
example, if you set a neighbourhood_size of 50000, a 50kb neighbourhood will be extracted
from the NCBI record associateed with each cblaster hit. If the record is too small for this
to be done (i.e. the record is smaller then 50kb) it is discarded.
-wg, --write_genomes
If set, will write entire NCBI record containing a cblaster hit to file (as well as just the
neighbourhood)
```
### Sieve subcommand:
```
>SyntenyQC sieve -h
usage: SyntenyQC sieve [-h] -gf [-ev] [-mi] [-mts] [-mev] -sf
Filter redundant genomic neighbourhoods based on neighbourhood similarity:
- First, an all-vs-all BLASTP is performed with user-specified BLASTP settings and the neighbourhoods in GENBANK_FOLDER.
- Secondly, these are parsed to define reciprocal best hits between every pair of neighbourhoods.
- Thirdly, these reciprocal best hits are used to derive a neighbourhood similarity network. Nodes are neighbourhood
filenames and edges indicate two neighbourhood nodes that have a similarity > SIMILARITY_FILTER.
Similarity = Number of RBHs / Number of proteins in smallest neighbourhood in pair.
- Finally, this network is pruned to remove neighbourhoods that exceed the user's SIMILARITY_FILTER threshold.
Nodes that remain are copied to the newly created folder 'genbank_folder/sieve_results/genbank'.
options:
-h, --help show this help message and exit
-g, --genbank_folder
Full path to folder containing neighbourhood genbank files requiring de-duplication
-ev, --e_value BLASTP evalue threshold. (default: 1e-05)
-mi, --min_percent_identity
BLASTP percent identity threshold. (default: 50)
-mts, --max_target_seqs
BLASTP -max_target_seqs. Maximum number of aligned sequences to keep. (default: 200)
-mev, --min_edge_view
Minimum similarity between two neighbourhoods for an edge to be drawn betweeen them in the RBH
graph. Purely for visualisation of the graph HTML file - has no impact on the graph pruning
results. (default: None)
-sf, --similarity_filter
Similarity threshold above which two neighbourhoods are considered redundant
```
### Sieve pruning algorithm:
```
Data: RBH graph
Result: Nodes from pruned RBH graph
Procedure:
while max(node degrees in RBH graph) > 0:
delete nodes = []
for node in RBH graph:
if node degree == max(node degrees in RBH graph):
delete nodes + node
delete node = random node from delete nodes
RBH graph = RBH graph - delete node
return nodes in RBH graph
```
# Example use:
## (1) Preliminaries
Take the [actinorhodin BGC from MIBIG](https://mibig.secondarymetabolites.org/repository/BGC0000194/index.html#r1c1), and analyse with `cblaster` via the command line or `CAGECAT` web server. This will generate a set of files, one of which will be a 'binary file', stored at `folder/with/binary.csv`. To show that even stringent searches can generate many hits, all core biosynthetic genes (`-r CAC44200.1;CAC44201.1`) and >= 11 hits (`-mh 11`) were required for each hit neighbourhood.
**⚠️WARNING: To use a binary file in `SyntenyQC collect`, `cblaster` must be run with a ',' binary delimeter and no intermediate genes (`-bde ','` and no `-ig` flag)⚠️**
## (2) Collect neighbourhoods
### Starting directory structure:
```
folder/with/binary.csv
```
### Command (neighbourhood size 42566kb = 2 x BGC length):
```
#Note - you should add your email
SyntenyQC collect -bp path/to/BGC0000194_binary.txt -ns 42566 -em my_email@domain.com -fn organism -sp -wg
```
### Finishing directory structure:
```
folder/with/binary/neighbourhood/organism1.gbk, organism2.gbk...organism157.gbk
/genomes /organism1.gbk, organism2.gbk...organism157.gbk ###only if -wg!
/log.txt
```
#### 🔴 157 neighbourhoods is a lot for a synteny plot 🔴
## (3) Sieve neighbourhoods
### Starting directory structure
From `SyntenyQC collect` in this example, but any folder with at least one genbank file can be used:
```
folder/with/binary/neighbourhood/organism1.gbk, ...
```
### Command:
```
SyntenyQC sieve -gf folder/with/binary/neighbourhood -sf 0.7
```
### Finishing directory structure:
```
folder/with/binary/neighbourhood/organism1.gbk, ...
/sieve_results/blastp /results.xml, db.txt, db.pin, ... #call all be deleted
/genbank /organism1.gbk, ...organism38.gbk #use as e.g. clinker input
/visualisations/RBH_graph.html, RBH_histogram.html #see what is being pruned
/log.txt
```
#### :green_heart: 38 neighbourhoods is OK for a synteny plot :green_heart:
## Notes
- Given filenames are purely to show number of files - `neighbourhood/sieve_results/genbank/organism1.gbk` is one of the genbanks in the `neighbourhood` folder, but may be different to `neighbourhood/organism1.gbk`.
- `RBH_graph` is an interactive html picture of the similarity graph created by `SyntenyQC sieve` (before pruning), only showing edges with a similarity > `min_edge_view`.
Edges are black (< `similairty_filter`) or red (>= `similarity_filter`).
- `RBH_histogram` shows the distribution of edge weights.
- Most neighbourhoods that meet the user-defined similarity threshold `-sf` will be removed in a single sieve run. However, the `sieve -mts` setting can impact final results. If a protein has homologs in 251 neighboughoods and `-mts` is 250, then one of the homologs will be missed by BLASTP, and the host neighbourhood may appear less similar to the neighbourhood with the query. Whilst a high `-mts` setting could be used for `sieve` runs pruning many neighbourhoods, this will generate large blast files that may take up a lot of space and slow down the run. Thus, if users wish to remove the (typically 2-3) redundant neighbourhoods remaining after a single `sieve` call, they can run sieve again on the pruned results (`syntenyqc sieve -g path/to/sieve_results/genbank -sf 0.7`).
## References
### `cblaster`
**Paper:** Cameron L M Gilchrist, Thomas J Booth, Bram van Wersch, Liana van Grieken, Marnix H Medema, Yit-Heng Chooi, cblaster: a remote search tool for rapid identification and visualization of homologous gene clusters, Bioinformatics Advances, Volume 1, Issue 1, 2021, vbab016, https://doi.org/10.1093/bioadv/vbab016
**Docs:** https://cblaster.readthedocs.io/en/latest/
### `clinker`
**Paper:** Cameron L M Gilchrist, Yit-Heng Chooi, clinker & clustermap.js: automatic generation of gene cluster comparison figures, Bioinformatics, Volume 37, Issue 16, August 2021, Pages 2473–2475, https://doi.org/10.1093/bioinformatics/btab007
**Docs:** https://github.com/gamcil/clinker
### `CAGECAT`
**Paper:** van den Belt, M., Gilchrist, C., Booth, T.J. et al. CAGECAT: The CompArative GEne Cluster Analysis Toolbox for rapid search and visualisation of homologous gene clusters. BMC Bioinformatics 24, 181 (2023). https://doi.org/10.1186/s12859-023-05311-2
**Website:** https://cagecat.bioinformatics.nl/
### `antiSMASH`
**Paper (ClusterBlast was introduced in version 1):** Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R. antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011 Jul;39(Web Server issue):W339-46. doi: 10.1093/nar/gkr466. Epub 2011 Jun 14. PMID: 21672958; PMCID: PMC3125804
**Website (latest version):** https://antismash.secondarymetabolites.org/#!/start
Raw data
{
"_id": null,
"home_page": null,
"name": "SyntenyQC",
"maintainer": null,
"docs_url": null,
"requires_python": ">=3.10",
"maintainer_email": "Timothy Kirkwood <tdjkirkwood@hotmail.com>",
"keywords": "bgc, biosynthetic gene cluster, genome neighbourhood, neighbourhood, synteny, synteny plot",
"author": null,
"author_email": "Timothy Kirkwood <tdjkirkwood@hotmail.com>",
"download_url": "https://files.pythonhosted.org/packages/df/57/c0bcf625a35514c5f596fdb39568f8323bcdd348d5f9bd82a18f6d7affab/syntenyqc-1.1.tar.gz",
"platform": null,
"description": "# NOTE - this app has only been tested on the Windows OS\n\n# SyntenyQC\n## Motivation: \nSynteny plots are widely used for the comparison of genomic neighbourhoods. Whilst synteny plots are often included as part of larger software suites (e.g. the `antiSMASH` ClusterBlast module), various low-code, stand-alone tools are now available that allow users to source candidate neighbourhoods and build their own synteny plots. However, a gap remains between: \n\n**(i) tools that source these candidate neighbourhoods** (e.g. `cblaster`, which can find hundreds of candidates), \n\n**(ii) tools that build the synteny plots** (e.g. `clinker`, which struggles as the number of neighbourhoods exceeds 30-50) and \n\n**(iii) the synteny plots themselves**, which become much harder to analyse/present as the number of neighbourhoods they include increases.\n\n## Description: \n`SyntenyQC` is a python app for the curation of neighbourhoods immediately prior to synteny plot creation. `SyntenyQC collect` supports the systematic definition and annotation of candidate neighbourhoods based on a direct integration to `cblaster`. `SytenyQC sieve` offers a flexible method for objectively removing redundant neighbourhoods (sourced using `cblaster` or any other tool) prior to synteny plot creation. This is in some cases an absolute requirement (e.g. `cblaster` called via the `CAGECAT` webserver places a limit of 50 neighbourhoods). \n\n## Installation \n```\npip install SyntenyQC\n```\nNote - `SyntenyQC` depends on [BLAST+](https://blast.ncbi.nlm.nih.gov/doc/blast-help/downloadblastdata.html), which must be [installed](https://www.ncbi.nlm.nih.gov/books/NBK569861/) by the user (tested with v2.12.0 - but should work with other versions unless there are parameter changes). If this is installed correctly, you should be able to see help messages after typing `blastp -h` and `makeblastdb -h` in the command line. \n\n## Tests\nTests are performed using [pytest](https://pypi.org/project/pytest/), but are not distributed with `SyntenyQC`. To run tests:\n\n1) Install pytest\n2) Clone the `SyntenyQC` github repository.\n3) Update `path/to/cloned/repository/tests/email.txt` with your email (if left blank, two webscraping tests will fail).\n4) Navigate to the cloned repository via command line (on Windows, type `cd path/to/cloned/repository`).\n5) Type `pytest` in the command line.\n \n## Usage\n### General help:\n```\n>SyntenyQC -h\nusage: SyntenyQC [-h] {collect,sieve} ...\n\noptions:\n -h, --help show this help message and exit\n\nsubcommands:\n {collect,sieve} Synteny quality control options\n```\n### Collect subcommand:\n```\n>SyntenyQC collect -h\nusage: SyntenyQC collect [-h] -bp -ns -em [-fn] [-sp] [-wg]\n\nWrite genbank files corresponding to cblaster neighbourhoods from a specified CSV-format binary file located at\nBINARY_PATH. For each cblaster hit accession in the binary file:\n\n1) A record is downloaded from NCBI using the accession. NCBI requires a user EMAIL to search for this record\n programatically. If WRITE_GENOMES is specified, this record is written to a local file according to FILENAMES\n (see final bulletpoint).\n2) A neighbourhood of size NEIGHBOURHOOD_SIZE bp is defined, centered on the cblaster hits defined in the binary\n file for the target accession.\n3) (If STRICT_SPAN is specified:) If the accession's record is too small to contain a neighbourhood of the\n desired size, it is discarded. For example, if an accession record is a 25kb contig and NEIGHBOURHOOD_SIZE\n is 50000, the record is discarded.\n4) If FILENAMES is \"organism\", the nighbourhood is written to file called *organism*.gbk. If FILENAMES is\n \"accession\", the neighbourhood is written to *accession*.gbk. Synteny softwares such as clinker can use these\n filenames to label synetny plot samples.\n\nOnce COLLECT has been run, a new folder with the same name as the binary file will be seen in the directory\nthat holds the binary file (i.e. the file \"path/to/binary/file.txt\" will generate the folder \"path/to/binary/file\").\nThis folder will have a subdirectory called \"neighbourhood\", containing all of the neighbourhood genbank files\n(i.e. \"path/to/binary/file/neighbourhood\"). If WRITE_GENOMES is specified, a second direcory (\"genome\") will also\nbe present, containing the entire record associated with each cblaster accession (i.e. \"path/to/binary/file/genome\").\nFinally, a log file will be present in the folder \"path/to/binary/file\", containing all run details.\n\noptions:\n -h, --help show this help message and exit\n -bp, --binary_path\n Full filepath to the CSV-format cblaster binary file containing neighbourhoods that should\n be extracted\n -n, --neighbourhood_size\n Size (basepairs) of neighbourhood to be extracted (centered on middle of CBLASTER-defined\n neighbourhood)\n -em, --email Email - required for NCBI entrez querying\n -fn, --filenames\n If \"organism\", all collected files will be named according to organism. If \"accession\", all\n files will be named by NCBI accession. (default:\n organism)\n -sp, --strict_span\n If set, will discard all neighbourhoods that are smaller than neighbourhood_size bp. For\n example, if you set a neighbourhood_size of 50000, a 50kb neighbourhood will be extracted\n from the NCBI record associateed with each cblaster hit. If the record is too small for this\n to be done (i.e. the record is smaller then 50kb) it is discarded.\n -wg, --write_genomes\n If set, will write entire NCBI record containing a cblaster hit to file (as well as just the\n neighbourhood)\n```\n### Sieve subcommand:\n```\n>SyntenyQC sieve -h\nusage: SyntenyQC sieve [-h] -gf [-ev] [-mi] [-mts] [-mev] -sf\n\nFilter redundant genomic neighbourhoods based on neighbourhood similarity:\n- First, an all-vs-all BLASTP is performed with user-specified BLASTP settings and the neighbourhoods in GENBANK_FOLDER.\n- Secondly, these are parsed to define reciprocal best hits between every pair of neighbourhoods.\n- Thirdly, these reciprocal best hits are used to derive a neighbourhood similarity network. Nodes are neighbourhood\n filenames and edges indicate two neighbourhood nodes that have a similarity > SIMILARITY_FILTER.\n Similarity = Number of RBHs / Number of proteins in smallest neighbourhood in pair.\n- Finally, this network is pruned to remove neighbourhoods that exceed the user's SIMILARITY_FILTER threshold.\n Nodes that remain are copied to the newly created folder 'genbank_folder/sieve_results/genbank'.\n\noptions:\n -h, --help show this help message and exit\n -g, --genbank_folder\n Full path to folder containing neighbourhood genbank files requiring de-duplication\n -ev, --e_value BLASTP evalue threshold. (default: 1e-05)\n -mi, --min_percent_identity\n BLASTP percent identity threshold. (default: 50)\n -mts, --max_target_seqs\n BLASTP -max_target_seqs. Maximum number of aligned sequences to keep. (default: 200)\n -mev, --min_edge_view\n Minimum similarity between two neighbourhoods for an edge to be drawn betweeen them in the RBH\n graph. Purely for visualisation of the graph HTML file - has no impact on the graph pruning\n results. (default: None)\n -sf, --similarity_filter\n Similarity threshold above which two neighbourhoods are considered redundant\n```\n### Sieve pruning algorithm:\n```\nData: RBH graph\nResult: Nodes from pruned RBH graph\nProcedure:\n while max(node degrees in RBH graph) > 0:\n delete nodes = []\n for node in RBH graph:\n if node degree == max(node degrees in RBH graph):\n delete nodes + node\n delete node = random node from delete nodes\n RBH graph = RBH graph - delete node \n return nodes in RBH graph \n```\n\n# Example use:\n## (1) Preliminaries\nTake the [actinorhodin BGC from MIBIG](https://mibig.secondarymetabolites.org/repository/BGC0000194/index.html#r1c1), and analyse with `cblaster` via the command line or `CAGECAT` web server. This will generate a set of files, one of which will be a 'binary file', stored at `folder/with/binary.csv`. To show that even stringent searches can generate many hits, all core biosynthetic genes (`-r CAC44200.1;CAC44201.1`) and >= 11 hits (`-mh 11`) were required for each hit neighbourhood.\n\n**\u26a0\ufe0fWARNING: To use a binary file in `SyntenyQC collect`, `cblaster` must be run with a ',' binary delimeter and no intermediate genes (`-bde ','` and no `-ig` flag)\u26a0\ufe0f** \n## (2) Collect neighbourhoods \n### Starting directory structure:\n```\nfolder/with/binary.csv\n```\n \n### Command (neighbourhood size 42566kb = 2 x BGC length):\n \n```\n#Note - you should add your email\nSyntenyQC collect -bp path/to/BGC0000194_binary.txt -ns 42566 -em my_email@domain.com -fn organism -sp -wg\n```\n \n### Finishing directory structure: \n```\nfolder/with/binary/neighbourhood/organism1.gbk, organism2.gbk...organism157.gbk\n /genomes /organism1.gbk, organism2.gbk...organism157.gbk ###only if -wg!\n /log.txt\n```\n#### \ud83d\udd34 157 neighbourhoods is a lot for a synteny plot \ud83d\udd34\n\n\n## (3) Sieve neighbourhoods\n### Starting directory structure \nFrom `SyntenyQC collect` in this example, but any folder with at least one genbank file can be used:\n```\nfolder/with/binary/neighbourhood/organism1.gbk, ...\n```\n### Command:\n```\nSyntenyQC sieve -gf folder/with/binary/neighbourhood -sf 0.7\n```\n### Finishing directory structure: \n```\nfolder/with/binary/neighbourhood/organism1.gbk, ...\n /sieve_results/blastp /results.xml, db.txt, db.pin, ... #call all be deleted\n /genbank /organism1.gbk, ...organism38.gbk #use as e.g. clinker input\n /visualisations/RBH_graph.html, RBH_histogram.html #see what is being pruned\n /log.txt\n```\n#### :green_heart: 38 neighbourhoods is OK for a synteny plot :green_heart:\n## Notes \n- Given filenames are purely to show number of files - `neighbourhood/sieve_results/genbank/organism1.gbk` is one of the genbanks in the `neighbourhood` folder, but may be different to `neighbourhood/organism1.gbk`.\n- `RBH_graph` is an interactive html picture of the similarity graph created by `SyntenyQC sieve` (before pruning), only showing edges with a similarity > `min_edge_view`. \n Edges are black (< `similairty_filter`) or red (>= `similarity_filter`).\n- `RBH_histogram` shows the distribution of edge weights.\n- Most neighbourhoods that meet the user-defined similarity threshold `-sf` will be removed in a single sieve run. However, the `sieve -mts` setting can impact final results. If a protein has homologs in 251 neighboughoods and `-mts` is 250, then one of the homologs will be missed by BLASTP, and the host neighbourhood may appear less similar to the neighbourhood with the query. Whilst a high `-mts` setting could be used for `sieve` runs pruning many neighbourhoods, this will generate large blast files that may take up a lot of space and slow down the run. Thus, if users wish to remove the (typically 2-3) redundant neighbourhoods remaining after a single `sieve` call, they can run sieve again on the pruned results (`syntenyqc sieve -g path/to/sieve_results/genbank -sf 0.7`). \n\n## References\n### `cblaster`\n\n**Paper:** Cameron L M Gilchrist, Thomas J Booth, Bram van Wersch, Liana van Grieken, Marnix H Medema, Yit-Heng Chooi, cblaster: a remote search tool for rapid identification and visualization of homologous gene clusters, Bioinformatics Advances, Volume 1, Issue 1, 2021, vbab016, https://doi.org/10.1093/bioadv/vbab016 \n\n**Docs:** https://cblaster.readthedocs.io/en/latest/ \n\n### `clinker` \n\n**Paper:** Cameron L M Gilchrist, Yit-Heng Chooi, clinker & clustermap.js: automatic generation of gene cluster comparison figures, Bioinformatics, Volume 37, Issue 16, August 2021, Pages 2473\u20132475, https://doi.org/10.1093/bioinformatics/btab007 \n\n**Docs:** https://github.com/gamcil/clinker \n\n### `CAGECAT`\n\n**Paper:** van den Belt, M., Gilchrist, C., Booth, T.J. et al. CAGECAT: The CompArative GEne Cluster Analysis Toolbox for rapid search and visualisation of homologous gene clusters. BMC Bioinformatics 24, 181 (2023). https://doi.org/10.1186/s12859-023-05311-2\n\n**Website:** https://cagecat.bioinformatics.nl/ \n\n### `antiSMASH` \n\n**Paper (ClusterBlast was introduced in version 1):** Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R. antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011 Jul;39(Web Server issue):W339-46. doi: 10.1093/nar/gkr466. Epub 2011 Jun 14. PMID: 21672958; PMCID: PMC3125804\n\n**Website (latest version):** https://antismash.secondarymetabolites.org/#!/start \n",
"bugtrack_url": null,
"license": "MIT License Copyright (c) 2024 Timothy David John Kirkwood Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the \"Software\"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. THE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.",
"summary": "A command line app for the pre-processing and quality control of synteny plot neighbourhoods",
"version": "1.1",
"project_urls": {
"Homepage": "https://github.com/Tim-Kirkwood/SyntenyQC",
"Issues": "https://github.com/Tim-Kirkwood/SyntenyQC/issues"
},
"split_keywords": [
"bgc",
" biosynthetic gene cluster",
" genome neighbourhood",
" neighbourhood",
" synteny",
" synteny plot"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "c6d11a9e5d0b0a5863e4bb5d5092e49efbf50bf54dc6bc6313703ccb0e8b23c8",
"md5": "83453c6b6b80946bd4d1d27e2f3de924",
"sha256": "87d43e2383e815177c3e690ca24c84dbdfc06c08135612a5a8f4385dbd292467"
},
"downloads": -1,
"filename": "syntenyqc-1.1-py3-none-any.whl",
"has_sig": false,
"md5_digest": "83453c6b6b80946bd4d1d27e2f3de924",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.10",
"size": 34929,
"upload_time": "2024-08-10T09:53:55",
"upload_time_iso_8601": "2024-08-10T09:53:55.624002Z",
"url": "https://files.pythonhosted.org/packages/c6/d1/1a9e5d0b0a5863e4bb5d5092e49efbf50bf54dc6bc6313703ccb0e8b23c8/syntenyqc-1.1-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "df57c0bcf625a35514c5f596fdb39568f8323bcdd348d5f9bd82a18f6d7affab",
"md5": "35c0cabd5f3b3511fcc9f0f78a70114a",
"sha256": "1dd53b9422250f7ee3bdbd3d51d7a3d1bd1739f44822f60b41c58d59957b6882"
},
"downloads": -1,
"filename": "syntenyqc-1.1.tar.gz",
"has_sig": false,
"md5_digest": "35c0cabd5f3b3511fcc9f0f78a70114a",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.10",
"size": 5640140,
"upload_time": "2024-08-10T09:54:03",
"upload_time_iso_8601": "2024-08-10T09:54:03.978435Z",
"url": "https://files.pythonhosted.org/packages/df/57/c0bcf625a35514c5f596fdb39568f8323bcdd348d5f9bd82a18f6d7affab/syntenyqc-1.1.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-08-10 09:54:03",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "Tim-Kirkwood",
"github_project": "SyntenyQC",
"github_not_found": true,
"lcname": "syntenyqc"
}
