.. image:: https://badge.fury.io/py/sequana-rnaseq.svg
:target: https://pypi.python.org/pypi/sequana_rnaseq
.. image:: http://joss.theoj.org/papers/10.21105/joss.00352/status.svg
:target: http://joss.theoj.org/papers/10.21105/joss.00352
:alt: JOSS (journal of open source software) DOI
.. image:: https://github.com/sequana/rnaseq/actions/workflows/main.yml/badge.svg
:target: https://github.com/sequana/rnaseq/actions/workflows/main.yaml
This is is the **RNA-seq** pipeline from the `Sequana <https://sequana.readthedocs.org>`_ project
:Overview: RNASeq analysis from raw data to feature counts
:Input: A set of Fastq Files and genome reference and annotation.
:Output: MultiQC and HTML reports, BAM and bigwig files, feature Counts, script to launch differential analysis
:Status: Production.
:Citation(sequana): Cokelaer et al, (2017), ‘Sequana’: a Set of Snakemake NGS pipelines, Journal of Open Source Software, 2(16), 352, JOSS DOI doi:10.21105/joss.00352
:Citation(pipeline):
.. image:: https://zenodo.org/badge/DOI/10.5281/zenodo.4047837.svg
:target: https://doi.org/10.5281/zenodo.4047837
Installation
~~~~~~~~~~~~
**sequana_rnaseq** is based on Python3, just install the package as follows:
pip install sequana_rnaseq --upgrade
You will need third-party software such as bowtie2/star. However, if you choose to use aptainer/singularity,
then nothing to install except singularity itself ! See below for details.
Usage
~~~~~
::
sequana_rnaseq --help
sequana_rnaseq --input-directory DATAPATH --genome-directory genome --aligner-choice star
This creates a directory with the pipeline and configuration file. You will then need
to execute the pipeline::
cd rnaseq
sh rnaseq.sh # for a local run
This launch a snakemake pipeline. If you are familiar with snakemake, you can
retrieve the pipeline itself and its configuration files and then execute the pipeline yourself with specific parameters::
snakemake -s rnaseq.rules -c config.yaml --cores 4 --stats stats.txt
Or use `sequanix <https://sequana.readthedocs.io/en/main/sequanix.html>`_ interface.
Usage with apptainer:
~~~~~~~~~~~~~~~~~~~~~~~~~
With apptainer, initiate the working directory as follows::
sequana_rnaseq --use-apptainer
Images are downloaded in the working directory but you can store then in a directory globally (e.g.)::
sequana_rnaseq --use-apptainer --apptainer-prefix ~/.sequana/apptainers
and then::
cd rnaseq
sh rnaseq.sh
if you decide to use snakemake manually, do not forget to add apptainer options::
snakemake -s rnaseq.rules -c config.yaml --cores 4 --use-apptainer --apptainer-prefix ~/.sequana/apptainers --apptainer-args "-B /home:/home"
Requirements
~~~~~~~~~~~~
This pipelines requires lots of third-party executable(s). Here is a list that
may change. A Message will inform you would you be missing an executable:
- bowtie
- bowtie2>=2.4.2
- STAR
- featureCounts (subread package)
- picard
- multiqc
- samtools
Note that bowtie>=2.4.2 is set to ensure the pipeline can be used with python 3.7-3.8-3.9 and the sequana-wrappers that supports bowtie2 with option --threads only (not previous versions). See environment.yaml or conda.yaml for latest list of required third-party tools.
You can install most of the tools using `damona <https://damona.readthedocs.io>`_::
damona create --name sequana_tools
damona activate sequana_tools
damona install sequana_tools
Or use the conda.yaml file available in this repository. If you start a new
environment from scratch, those commands will create the environment and install
all dependencies for you::
conda create --name sequana_env python 3.7.3
conda activate sequana_env
conda install -c anaconda qt pyqt>5
pip install sequana
pip install sequana_rnaseq
conda install --file https://raw.githubusercontent.com/sequana/rnaseq/main/conda.yaml
For Linux users, we provide singularity images available through within the damona project (https://damona.readthedocs.io).
.. image:: https://raw.githubusercontent.com/sequana/sequana_rnaseq/main/sequana_pipelines/rnaseq/dag.png
Details
~~~~~~~~~
This pipeline runs a **RNA-seq** analysis of sequencing data. It runs in
parallel on a set of input FastQ files (paired or not).
A brief HTML report is produced together with a MultiQC report.
This pipeline is complex and requires some expertise for the interpretation.
Many online-resources are available and should help you deciphering the output.
Yet, it should be quite straigtforward to execute it as shown above. The
pipeline uses bowtie1 to look for ribosomal contamination (rRNA). Then,
it cleans the data with cutapdat if you say so (your data may already be
pre-processed). If no adapters are provided (default), reads are
trimmed for low quality bases only. Then, mapping is performed with standard mappers such as
star or bowtie2 (--aligner option). Finally,
feature counts are extracted from the previously generated BAM files. We guess
the strand and save the feature counts into the directoy
./rnadiff/feature_counts.
The pipelines stops there. However, RNA-seq analysis are followed by a different
analysis (DGE hereafter). Although the DGE is not part of the pipeline, you can
performed it with standard tools using the data in ./rnadiff directory. One such
tool is provided within our framework (based on the well known DEseq2 software).
Using our framework::
cd rnadiff
sequana rnadiff --design design.csv --features all_features.out --annotation ANNOT \
--feature-name FEAT --attribute-name ATTR
where ANNOT is the annotation file of your analysis, FEAT and ATTR the attribute
and feature used in your analysis (coming from the annotation file).
This produces a HTML repot summarizing you differential analysis.
Rules and configuration details
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Here is the `latest documented configuration file <https://raw.githubusercontent.com/sequana/sequana_rnaseq/main/sequana_pipelines/rnaseq/config.yaml>`_
to be used with the pipeline. Each rule used in the pipeline may have a section in the configuration file.
.. warning:: the RNAseQC rule is switch off and is not currently functional in
version 0.9.X
Issues
~~~~~~
In the context of eukaryotes, you will need 32G of memory most probably. If this is too much,
you can try to restrict the memory. Check out the config.yaml file in the star section.
Changelog
~~~~~~~~~
========= ====================================================================
Version Description
========= ====================================================================
0.20.0 * Fix regression due to new sequana version
* Update summary html to use new sequana plots
0.19.3 * fix regression with click to set the default rRNA to 'rRNA' again.
0.19.2 * fix bowtie1 regression in the log file, paired end case in
multiqc and rnadiff script (regression)
* set genome directory default to None to enforce its usage
0.19.1 * add rnaseqc container.
* Update rseqc rules (redirection)
* cleanup onsuccess rule
0.19.0 * Refactorisation to use click
0.18.1 * fastp multiqc regression. Fixed missing sample names by updating
multiqc_config and adding sample names in the output filename
0.18.0 * New plots in the HTML reports. Includes version of executables.
0.17.2 * CHANGES: in star section, added --limitBAMsortRAM and set to 30G
* BUG: Fix missing params (options) in star_mapping rule not taken
into account
0.17.1 * use new rulegraph / graphviz apptainer
0.17.0 * fastp step changed to use sequana-wrappers. Slight change in
config file. The reverse and forward adapter options called
rev and fwd have been dropped in favor of a single adapters option.
v0.17.0 config and schema are not compatible with previous
versions.
* Update singularity containers and add new one for fastp
0.16.1 * fix bug in feature counts automatic strand balance detection. Was
always using the stranded case (2).
* add singularity workflow for testing
* fix documentation in config.yaml
0.16.0 * star, salmon, bam_coverage are now in sequana wrappers, updated
the pipeline accordingly
* updated config file and schema to include resources inside the
config file (so as to use new --profile option)
* set singularity images in all rules
* star wrappers has changed significantly to use star
recommandation. To keep using previous way, a legacy option
is available and set to True in this version.
* bamCoverage renamed in bam_coverage in the config file
* multiqc_config removed redundant information and ordered
the output in a coherent way (QC and then analysis)
0.15.2 * Fix bowtie2 rule to use new wrappers. Use wrappers in
add_read_group and mark_duplicates
0.15.1 * Adapt to new bowtie2 align wrapper
0.15.0 * fix typo reported in https://github.com/sequana/rnaseq/issues/12
* fix feature counts plots not showing anymore
* Script for differential analysis is now in the main pipeline
0.14.2 * fix feature counts plots missing in multiqc results
0.14.1 * fix regression bug introduced in snakemake 6.9.0
0.14.0 * Allow the aligners to have dedicated index for each version in the
same genome directory.
* Ribosomal is now estimated on the first 100,000 reads to speed up
analysis
* --indexing and --force-indexing options not required anymore.
Indexing will be done automatically and not redone if present.
* Use of the new sequana-wrappers repository
0.13.0 * Major update to use the new sequana version and the RNADiff tools.
* remove fastq_screen. One can use sequana_multitax for taxonomic
content and contamination.
* cutadapt is now replaced by fastp, although it can still be used.
* full integration of salmon for prokaryotes and eukaryotes
* user interface has now a --skip-gff-check option. Better handling of
input gff with more meaningful messages
* integration of rseqc tool
0.12.1 * indexing was always set to True in the config after 0.9.16 update.
0.12.0 * BUG fix: Switch mark_duplicates correctly beore feature counts
0.11.0 * rnadiff one factor is simplified
* When initiating the pipeline, provide information about the GFF
* mark duplicates off by default
* feature_counts has more options in the help. split options into
feature/attribute/extra_attributes.
* HTML reports better strand picture and information about rRNA
* refactorising the main standalone and config file to split feature
counts optiions into feature and attribute. Sanoty checks are ow
provided (--feature-counts-attribute, --feature-counts-feature-type)
* can provide a custom GFF not in the genome directory
* can provide several feature from the GFF. Then, a custom GFF is
created and used
* fix the --do-igvtools and --do-bam-coverage with better doc
0.10.0 * 9/12/2020
* Fixed bug in sequana/star_indexing for small genomes (v0.9.7).
Changed the rnaseq requirements to benefit from this bug-fix that
could lead to seg fault with star aligner for small genomes.
* Report improved with strand guess and plot
0.9.20 * 7/12/2020
* BUG in sequana/star rules v0.9.6. Fixed in this release.
* In config file, bowtie section 'do' option is removed. This is now
set automatically if rRNA_feature or rRNA_file is provided. This
allows us to skip the rRNA mapping entirely if needed.
* fastq_screen should be functional. Default behaviour is off. If
set only phiX174 will be search for. Users should build their own
configuration file.
* star/bowtie1/bowtie2 have now their own sub-directories in the
genome directory.
* added --run option to start pipeline automatically (if you know
what you are doing)
* rnadiff option has now a default value (one_factor)
* add strandness plot in the HTML summary page
0.9.19 * Remove the try/except around tolerance (guess of strandness) to
make sure this is provided by the user. Final onsuccess benefits
from faster GFF function (sequana 0.9.4)
0.9.18 * Fix typo (regression bug) + add tolerance in schema + generic
title in multiqc_config. (oct 2020)
0.9.17 * add the *tolerance* parameter in the feature_counts rule as a user
parameter (config and pipeline).
0.9.16 * Best feature_counts is now saved into rnadiff/feature_counts
directory and rnadiff scripts have been updated accordingly
* the most probable feature count option is now computed more
effectivily and incorporated inside the Snakemake pipeline (not in
the onsuccess) so that multiqc picks the best one (not the 3
results)
* the target.txt file can be generated inside the pipeline if user
fill the rnadiff/conditions section in the config file
* indexing options are filled automatically when calling
sequana_rnaseq based on the presence/absence of the index
of the aligner being used.
* salmon now integrated and feature counts created (still WIP in
sequana)
0.9.15 * FastQC on raw data skipped by default (FastQC
for processed data is still available)
* Added paired options (-p) for featureCounts
* Switch back markduplicates to False for now.
0.9.14 * Use only R1 with bowtie1
* set the memory requirements for mark_duplicates in cluster_config
file
* Set temporary directory for mark_duplicates to be local ./tmp
0.9.13 * set mark_duplicate to true by default
* use new sequana pipeline manager
* export all features counts in a single file
* custom HTML report
* faster --help calls
* --from-project option added
0.9.12 * include salmon tool as an alternative to star/bowtie2
* include rnadiff directory with required input for Differential
analysis
0.9.11 * Automatic guessing of the strandness of the experiment
0.9.10 * Fix multiqc for RNAseQC rule
0.9.9 * Fix RNAseQC rule, which is now available.
* Fix ability to use existing rRNA file as input
0.9.8 * Fix indexing for bowtie1 to not be done if aligner is different
* add new options: --feature-counts-options and --do-rnaseq-qc,
--rRNA-feature
* Based on the input GFF, we now check the validity of the rRNA
feature and feature counts options to check whether the feature
exists in the GFF
* schema is now used to check the config file values
* add a data test for testing and documentation
0.9.7 * fix typo found in version 0.9.6
0.9.6 * Fixed empty read tag in the configuration file
* Possiblity to switch off cutadapt section
* Fixing bowtie2 rule in sequana and update the pipeline accordingly
* Include a schema file
* output-directory parameter renamed into output_directory (multiqc
section)
* handle stdout correctly in fastqc, bowtie1, bowtie2 rules
0.9.5 * Fixed https://github.com/sequana/sequana/issues/571
* More cutadapt commands and sanity checks
* Fixed bowtie2 options import in rnaseq.rules
0.9.4
0.9.3 if a fastq_screen.conf is provided, we switch the fastqc_screen
section ON automatically
0.9.0 **Major refactorisation.**
* remove sartools, kraken rules.
* Indexing is now optional and can be set in the configuration.
* Configuration file is simplified with a general section to enter
the genome location and aligner.
* Fixed rules in sequana (0.8.0) that were not up-to-date with
several executables used in the pipeline including picard,
fastq_screen, etc. See Sequana Changelog for details with respect
to rules changes.
* Copying the feature counts in main directory ready to use for
a differential analysis.
========= ====================================================================
Raw data
{
"_id": null,
"home_page": "https://github.com/sequana/rnaseq",
"name": "sequana-rnaseq",
"maintainer": null,
"docs_url": null,
"requires_python": "<4.0,>=3.8",
"maintainer_email": null,
"keywords": "snakemake, sequana, RNAseq, RNADiff, differential analysis",
"author": "Sequana Team",
"author_email": null,
"download_url": "https://files.pythonhosted.org/packages/1a/58/59ffe3f081bf37b01993949b347b5c45b2a58804a22d3099fe8b89c0f959/sequana_rnaseq-0.20.0.tar.gz",
"platform": null,
"description": "\n\n.. image:: https://badge.fury.io/py/sequana-rnaseq.svg\n :target: https://pypi.python.org/pypi/sequana_rnaseq\n\n.. image:: http://joss.theoj.org/papers/10.21105/joss.00352/status.svg\n :target: http://joss.theoj.org/papers/10.21105/joss.00352\n :alt: JOSS (journal of open source software) DOI\n\n.. image:: https://github.com/sequana/rnaseq/actions/workflows/main.yml/badge.svg\n :target: https://github.com/sequana/rnaseq/actions/workflows/main.yaml\n\n\n\nThis is is the **RNA-seq** pipeline from the `Sequana <https://sequana.readthedocs.org>`_ project\n\n:Overview: RNASeq analysis from raw data to feature counts\n:Input: A set of Fastq Files and genome reference and annotation.\n:Output: MultiQC and HTML reports, BAM and bigwig files, feature Counts, script to launch differential analysis\n:Status: Production.\n:Citation(sequana): Cokelaer et al, (2017), \u2018Sequana\u2019: a Set of Snakemake NGS pipelines, Journal of Open Source Software, 2(16), 352, JOSS DOI doi:10.21105/joss.00352\n:Citation(pipeline):\n .. image:: https://zenodo.org/badge/DOI/10.5281/zenodo.4047837.svg\n :target: https://doi.org/10.5281/zenodo.4047837\n\nInstallation\n~~~~~~~~~~~~\n\n**sequana_rnaseq** is based on Python3, just install the package as follows:\n\n pip install sequana_rnaseq --upgrade\n\nYou will need third-party software such as bowtie2/star. However, if you choose to use aptainer/singularity,\nthen nothing to install except singularity itself ! See below for details.\n\n\nUsage\n~~~~~\n\n::\n\n sequana_rnaseq --help\n sequana_rnaseq --input-directory DATAPATH --genome-directory genome --aligner-choice star\n\nThis creates a directory with the pipeline and configuration file. You will then need\nto execute the pipeline::\n\n cd rnaseq\n sh rnaseq.sh # for a local run\n\nThis launch a snakemake pipeline. If you are familiar with snakemake, you can\nretrieve the pipeline itself and its configuration files and then execute the pipeline yourself with specific parameters::\n\n snakemake -s rnaseq.rules -c config.yaml --cores 4 --stats stats.txt\n\nOr use `sequanix <https://sequana.readthedocs.io/en/main/sequanix.html>`_ interface.\n\n\nUsage with apptainer:\n~~~~~~~~~~~~~~~~~~~~~~~~~\n\nWith apptainer, initiate the working directory as follows::\n\n sequana_rnaseq --use-apptainer\n\nImages are downloaded in the working directory but you can store then in a directory globally (e.g.)::\n\n sequana_rnaseq --use-apptainer --apptainer-prefix ~/.sequana/apptainers\n\nand then::\n\n cd rnaseq\n sh rnaseq.sh\n\nif you decide to use snakemake manually, do not forget to add apptainer options::\n\n snakemake -s rnaseq.rules -c config.yaml --cores 4 --use-apptainer --apptainer-prefix ~/.sequana/apptainers --apptainer-args \"-B /home:/home\"\n\n\nRequirements\n~~~~~~~~~~~~\n\nThis pipelines requires lots of third-party executable(s). Here is a list that\nmay change. A Message will inform you would you be missing an executable:\n\n- bowtie\n- bowtie2>=2.4.2\n- STAR\n- featureCounts (subread package)\n- picard\n- multiqc\n- samtools\n\nNote that bowtie>=2.4.2 is set to ensure the pipeline can be used with python 3.7-3.8-3.9 and the sequana-wrappers that supports bowtie2 with option --threads only (not previous versions). See environment.yaml or conda.yaml for latest list of required third-party tools.\n\nYou can install most of the tools using `damona <https://damona.readthedocs.io>`_::\n\n damona create --name sequana_tools\n damona activate sequana_tools\n damona install sequana_tools\n\nOr use the conda.yaml file available in this repository. If you start a new\nenvironment from scratch, those commands will create the environment and install\nall dependencies for you::\n\n conda create --name sequana_env python 3.7.3\n conda activate sequana_env\n conda install -c anaconda qt pyqt>5\n pip install sequana\n pip install sequana_rnaseq\n conda install --file https://raw.githubusercontent.com/sequana/rnaseq/main/conda.yaml\n\nFor Linux users, we provide singularity images available through within the damona project (https://damona.readthedocs.io).\n\n\n.. image:: https://raw.githubusercontent.com/sequana/sequana_rnaseq/main/sequana_pipelines/rnaseq/dag.png\n\n\nDetails\n~~~~~~~~~\n\nThis pipeline runs a **RNA-seq** analysis of sequencing data. It runs in\nparallel on a set of input FastQ files (paired or not).\nA brief HTML report is produced together with a MultiQC report.\n\nThis pipeline is complex and requires some expertise for the interpretation.\nMany online-resources are available and should help you deciphering the output.\n\nYet, it should be quite straigtforward to execute it as shown above. The\npipeline uses bowtie1 to look for ribosomal contamination (rRNA). Then,\nit cleans the data with cutapdat if you say so (your data may already be\npre-processed). If no adapters are provided (default), reads are\ntrimmed for low quality bases only. Then, mapping is performed with standard mappers such as\nstar or bowtie2 (--aligner option). Finally,\nfeature counts are extracted from the previously generated BAM files. We guess\nthe strand and save the feature counts into the directoy\n./rnadiff/feature_counts.\n\nThe pipelines stops there. However, RNA-seq analysis are followed by a different\nanalysis (DGE hereafter). Although the DGE is not part of the pipeline, you can\nperformed it with standard tools using the data in ./rnadiff directory. One such\ntool is provided within our framework (based on the well known DEseq2 software).\n\nUsing our framework::\n\n cd rnadiff\n sequana rnadiff --design design.csv --features all_features.out --annotation ANNOT \\\n --feature-name FEAT --attribute-name ATTR\n\nwhere ANNOT is the annotation file of your analysis, FEAT and ATTR the attribute\nand feature used in your analysis (coming from the annotation file).\n\nThis produces a HTML repot summarizing you differential analysis.\n\n\nRules and configuration details\n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\n\nHere is the `latest documented configuration file <https://raw.githubusercontent.com/sequana/sequana_rnaseq/main/sequana_pipelines/rnaseq/config.yaml>`_\nto be used with the pipeline. Each rule used in the pipeline may have a section in the configuration file.\n\n\n.. warning:: the RNAseQC rule is switch off and is not currently functional in\n version 0.9.X\n\nIssues\n~~~~~~\n\nIn the context of eukaryotes, you will need 32G of memory most probably. If this is too much,\nyou can try to restrict the memory. Check out the config.yaml file in the star section.\n\n\n\nChangelog\n~~~~~~~~~\n\n========= ====================================================================\nVersion Description\n========= ====================================================================\n0.20.0 * Fix regression due to new sequana version\n * Update summary html to use new sequana plots\n0.19.3 * fix regression with click to set the default rRNA to 'rRNA' again.\n0.19.2 * fix bowtie1 regression in the log file, paired end case in\n multiqc and rnadiff script (regression)\n * set genome directory default to None to enforce its usage\n0.19.1 * add rnaseqc container.\n * Update rseqc rules (redirection)\n * cleanup onsuccess rule\n0.19.0 * Refactorisation to use click\n0.18.1 * fastp multiqc regression. Fixed missing sample names by updating\n multiqc_config and adding sample names in the output filename\n0.18.0 * New plots in the HTML reports. Includes version of executables.\n0.17.2 * CHANGES: in star section, added --limitBAMsortRAM and set to 30G\n * BUG: Fix missing params (options) in star_mapping rule not taken\n into account\n0.17.1 * use new rulegraph / graphviz apptainer\n0.17.0 * fastp step changed to use sequana-wrappers. Slight change in\n config file. The reverse and forward adapter options called\n rev and fwd have been dropped in favor of a single adapters option.\n v0.17.0 config and schema are not compatible with previous\n versions.\n * Update singularity containers and add new one for fastp\n0.16.1 * fix bug in feature counts automatic strand balance detection. Was\n always using the stranded case (2).\n * add singularity workflow for testing\n * fix documentation in config.yaml\n0.16.0 * star, salmon, bam_coverage are now in sequana wrappers, updated\n the pipeline accordingly\n * updated config file and schema to include resources inside the\n config file (so as to use new --profile option)\n * set singularity images in all rules\n * star wrappers has changed significantly to use star\n recommandation. To keep using previous way, a legacy option\n is available and set to True in this version.\n * bamCoverage renamed in bam_coverage in the config file\n * multiqc_config removed redundant information and ordered\n the output in a coherent way (QC and then analysis)\n0.15.2 * Fix bowtie2 rule to use new wrappers. Use wrappers in\n add_read_group and mark_duplicates\n0.15.1 * Adapt to new bowtie2 align wrapper\n0.15.0 * fix typo reported in https://github.com/sequana/rnaseq/issues/12\n * fix feature counts plots not showing anymore\n * Script for differential analysis is now in the main pipeline\n0.14.2 * fix feature counts plots missing in multiqc results\n0.14.1 * fix regression bug introduced in snakemake 6.9.0\n0.14.0 * Allow the aligners to have dedicated index for each version in the\n same genome directory.\n * Ribosomal is now estimated on the first 100,000 reads to speed up\n analysis\n * --indexing and --force-indexing options not required anymore.\n Indexing will be done automatically and not redone if present.\n * Use of the new sequana-wrappers repository\n0.13.0 * Major update to use the new sequana version and the RNADiff tools.\n * remove fastq_screen. One can use sequana_multitax for taxonomic\n content and contamination.\n * cutadapt is now replaced by fastp, although it can still be used.\n * full integration of salmon for prokaryotes and eukaryotes\n * user interface has now a --skip-gff-check option. Better handling of\n input gff with more meaningful messages\n * integration of rseqc tool\n0.12.1 * indexing was always set to True in the config after 0.9.16 update.\n0.12.0 * BUG fix: Switch mark_duplicates correctly beore feature counts\n0.11.0 * rnadiff one factor is simplified\n * When initiating the pipeline, provide information about the GFF\n * mark duplicates off by default\n * feature_counts has more options in the help. split options into\n feature/attribute/extra_attributes.\n * HTML reports better strand picture and information about rRNA\n * refactorising the main standalone and config file to split feature\n counts optiions into feature and attribute. Sanoty checks are ow\n provided (--feature-counts-attribute, --feature-counts-feature-type)\n * can provide a custom GFF not in the genome directory\n * can provide several feature from the GFF. Then, a custom GFF is\n created and used\n * fix the --do-igvtools and --do-bam-coverage with better doc\n0.10.0 * 9/12/2020\n * Fixed bug in sequana/star_indexing for small genomes (v0.9.7).\n Changed the rnaseq requirements to benefit from this bug-fix that\n could lead to seg fault with star aligner for small genomes.\n * Report improved with strand guess and plot\n0.9.20 * 7/12/2020\n * BUG in sequana/star rules v0.9.6. Fixed in this release.\n * In config file, bowtie section 'do' option is removed. This is now\n set automatically if rRNA_feature or rRNA_file is provided. This\n allows us to skip the rRNA mapping entirely if needed.\n * fastq_screen should be functional. Default behaviour is off. If\n set only phiX174 will be search for. Users should build their own\n configuration file.\n * star/bowtie1/bowtie2 have now their own sub-directories in the\n genome directory.\n * added --run option to start pipeline automatically (if you know\n what you are doing)\n * rnadiff option has now a default value (one_factor)\n * add strandness plot in the HTML summary page\n0.9.19 * Remove the try/except around tolerance (guess of strandness) to\n make sure this is provided by the user. Final onsuccess benefits\n from faster GFF function (sequana 0.9.4)\n0.9.18 * Fix typo (regression bug) + add tolerance in schema + generic\n title in multiqc_config. (oct 2020)\n0.9.17 * add the *tolerance* parameter in the feature_counts rule as a user\n parameter (config and pipeline).\n0.9.16 * Best feature_counts is now saved into rnadiff/feature_counts\n directory and rnadiff scripts have been updated accordingly\n * the most probable feature count option is now computed more\n effectivily and incorporated inside the Snakemake pipeline (not in\n the onsuccess) so that multiqc picks the best one (not the 3\n results)\n * the target.txt file can be generated inside the pipeline if user\n fill the rnadiff/conditions section in the config file\n * indexing options are filled automatically when calling\n sequana_rnaseq based on the presence/absence of the index\n of the aligner being used.\n * salmon now integrated and feature counts created (still WIP in\n sequana)\n0.9.15 * FastQC on raw data skipped by default (FastQC\n for processed data is still available)\n * Added paired options (-p) for featureCounts\n * Switch back markduplicates to False for now.\n0.9.14 * Use only R1 with bowtie1\n * set the memory requirements for mark_duplicates in cluster_config\n file\n * Set temporary directory for mark_duplicates to be local ./tmp\n0.9.13 * set mark_duplicate to true by default\n * use new sequana pipeline manager\n * export all features counts in a single file\n * custom HTML report\n * faster --help calls\n * --from-project option added\n0.9.12 * include salmon tool as an alternative to star/bowtie2\n * include rnadiff directory with required input for Differential\n analysis\n0.9.11 * Automatic guessing of the strandness of the experiment\n0.9.10 * Fix multiqc for RNAseQC rule\n0.9.9 * Fix RNAseQC rule, which is now available.\n * Fix ability to use existing rRNA file as input\n0.9.8 * Fix indexing for bowtie1 to not be done if aligner is different\n * add new options: --feature-counts-options and --do-rnaseq-qc,\n --rRNA-feature\n * Based on the input GFF, we now check the validity of the rRNA\n feature and feature counts options to check whether the feature\n exists in the GFF\n * schema is now used to check the config file values\n * add a data test for testing and documentation\n0.9.7 * fix typo found in version 0.9.6\n0.9.6 * Fixed empty read tag in the configuration file\n * Possiblity to switch off cutadapt section\n * Fixing bowtie2 rule in sequana and update the pipeline accordingly\n * Include a schema file\n * output-directory parameter renamed into output_directory (multiqc\n section)\n * handle stdout correctly in fastqc, bowtie1, bowtie2 rules\n0.9.5 * Fixed https://github.com/sequana/sequana/issues/571\n * More cutadapt commands and sanity checks\n * Fixed bowtie2 options import in rnaseq.rules\n0.9.4\n0.9.3 if a fastq_screen.conf is provided, we switch the fastqc_screen\n section ON automatically\n0.9.0 **Major refactorisation.**\n\n * remove sartools, kraken rules.\n * Indexing is now optional and can be set in the configuration.\n * Configuration file is simplified with a general section to enter\n the genome location and aligner.\n * Fixed rules in sequana (0.8.0) that were not up-to-date with\n several executables used in the pipeline including picard,\n fastq_screen, etc. See Sequana Changelog for details with respect\n to rules changes.\n * Copying the feature counts in main directory ready to use for\n a differential analysis.\n========= ====================================================================\n",
"bugtrack_url": null,
"license": "BSD-3",
"summary": "A RNAseq pipeline from raw reads to feature counts",
"version": "0.20.0",
"project_urls": {
"Homepage": "https://github.com/sequana/rnaseq",
"Repository": "https://github.com/sequana/rnaseq"
},
"split_keywords": [
"snakemake",
" sequana",
" rnaseq",
" rnadiff",
" differential analysis"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "4b098ab7e654376a648294f69a31e8ac1a6982140659b6a2485a5869b364b1fa",
"md5": "e417fb388a32306c838611fbd9fb5f3e",
"sha256": "433bf548073b069881ffc2db23d85e0c0a2619adaa635c01447c26460e9560d6"
},
"downloads": -1,
"filename": "sequana_rnaseq-0.20.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "e417fb388a32306c838611fbd9fb5f3e",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": "<4.0,>=3.8",
"size": 80712,
"upload_time": "2024-08-19T18:26:08",
"upload_time_iso_8601": "2024-08-19T18:26:08.195706Z",
"url": "https://files.pythonhosted.org/packages/4b/09/8ab7e654376a648294f69a31e8ac1a6982140659b6a2485a5869b364b1fa/sequana_rnaseq-0.20.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "1a5859ffe3f081bf37b01993949b347b5c45b2a58804a22d3099fe8b89c0f959",
"md5": "25f98ceb0431a58373b2c3616d3ef89a",
"sha256": "a689f868705a0c61e34c0a99f432a9dacf060166e8bc37d938d80a9e7cf92126"
},
"downloads": -1,
"filename": "sequana_rnaseq-0.20.0.tar.gz",
"has_sig": false,
"md5_digest": "25f98ceb0431a58373b2c3616d3ef89a",
"packagetype": "sdist",
"python_version": "source",
"requires_python": "<4.0,>=3.8",
"size": 81432,
"upload_time": "2024-08-19T18:26:09",
"upload_time_iso_8601": "2024-08-19T18:26:09.896504Z",
"url": "https://files.pythonhosted.org/packages/1a/58/59ffe3f081bf37b01993949b347b5c45b2a58804a22d3099fe8b89c0f959/sequana_rnaseq-0.20.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-08-19 18:26:09",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "sequana",
"github_project": "rnaseq",
"travis_ci": false,
"coveralls": false,
"github_actions": true,
"lcname": "sequana-rnaseq"
}
JSON
