{
  "@context": "https://w3id.org/ro/crate/1.1/context",
  "@graph": [
    {
      "@id": "ro-crate-metadata.json",
      "@type": "CreativeWork",
      "about": {
        "@id": "./"
      },
      "conformsTo": [
        {
          "@id": "https://w3id.org/ro/crate/1.1"
        },
        {
          "@id": "https://w3id.org/workflowhub/workflow-ro-crate/1.0"
        }
      ]
    },
    {
      "@id": "ro-crate-preview.html",
      "@type": "CreativeWork",
      "about": {
        "@id": "./"
      }
    },
    {
      "@id": "./",
      "@type": "Dataset",
      "mainEntity": {
        "@id": "Worlkflow/annotation_sw.cwl"
      },
      "datePublished": "2025-09-19 07:17:33 +0100",
      "author": [
        {
          "@id": "https://orcid.org/0009-0007-5830-6482"
        },
        {
          "@id": "https://orcid.org/0009-0004-1874-3117"
        }
      ],
      "license": "MIT",
      "identifier": "https://doi.org/10.48546/workflowhub.workflow.1955.1",
      "url": "https://workflowhub.eu/workflows/1955/ro_crate?version=1",
      "name": "Research Object Crate for Gene Expression Analysis Workflow in Complex Microbiomes",
      "description": "\n![GitHub last commit (branch)](https://img.shields.io/github/last-commit/RyoMameda/workflow_cwl/main)\n![Status](https://img.shields.io/badge/status-development-yellow)\n[![cwltool](https://img.shields.io/badge/cwltool-3.1.20250110105449-success)](https://github.com/common-workflow-language/cwltool/releases/tag/3.1.20250110105449)\n[![License](https://img.shields.io/badge/License-MIT-blue.svg)](./LICENSE)\n[![Version](https://img.shields.io/badge/version-1.0-brightgreen)](https://github.com/RyoMameda/ComplexMicrobiome_GeneExpression_CWL/releases/tag/v1.0)\n[![Open in Dev Containers](https://img.shields.io/static/v1?label=Dev%20Containers&amp;message=python3.11&amp;color=blue&amp;logo=docker)](https://github.com/yonesora56/plant2human/tree/main/.devcontainer)\n\n# Gene Expression Analysis Workflow in Complex Microbiomes\n\n&nbsp;\n\n## Workflow Schema \n- more details: [Optimization of Mapping Tools and Investigation of Ribosomal RNA Influence for Data-Driven Gene Expression Analysis in Complex Microbiomes](https://doi.org/10.3390/microorganisms13050995)\n\n![image](./image/microorganisms-13-00995-g001.png)\n\n&nbsp;\n\n### 1. Overview of the Workflow\n\nThis analysis focuses on transcriptional profiling of complex microbiomes. It requires both metagenomic and metatranscriptomic NGS short-read data, along with annotation reference information (e.g., ribosomal RNA sequences and referenced protein databases, listed below). The metagenomic and metatranscriptomic reads should be derived from the same microbiome samples. Assembled metagenomic contigs are then used as reference sequences to map both types of reads, enabling gene-level quantification.\n\n&nbsp;\n\n### 2. Minimum Requirements\n\n- `Docker`\n- `cwltool`\n\n&nbsp;\n\n### 3. Workflow Component\n\nThis analysis workflow is composed of three sub-workflows; metagenomic contig assembling, reads mapping and annotation. \n\n&nbsp;\n\n#### Metagenomic contig assembling\n\nIn this process, the following steps are performed:\n\n\n1. Assembly process using `Megahit`. \n2. Prediction Protein sequences using `Prodigal`.\n3. Statical analysis of contigs useing `SeqKit`.\n\n&nbsp;\n\n#### Reads mapping\n\nIn this process, the following steps are performed:\n\n1. Mapping process using `BWA MEM`.\n2. Statical analysis of mapping results using `SAMtools`\n\n&nbsp;\n\n#### Annotation\n\nIn this process, the following steps are performed:\n\n1. Searching contaminated ribosomal RNA sequences using `BLAST`.\n2. Searching referenced proteins using `DIAMOND`.\n3. Creation GTF formated file contained annotation informations.\n\n&nbsp;\n\n### 4. Test Dataset and Your Own Dataset\n\n- If you are testing with the following files, please place them in the `Data` directory!\n- You can also obtain metagenomic and metatranscriptomic FASTQ files either by downloading them from public databases or by using your own samples, and then place them in your `Data` directory.\n\n#### Metagenome data\n\n- [SRR27548858](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548858)\n\n#### Metatranscriptome data\n\n- [SRR27548863](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548863)\n- [SRR27548864](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548864)\n- [SRR27548865](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548865)\n\n&nbsp;\n\n### 5. Annotation References\n\nThese reference files are used in the BLAST and DIAMOND processes. The downloaded files are available in the `Data` directory (accessed on September 17, 2025). If you wish to use the latest versions of the references, please download them using the following scripts.\n\n```bash\n# rRNA data from SILVA website (release138.1; accessed on 17,September,2025)\ncurl -O https://ftp.arb-silva.de/release_138.1/Exports/SILVA_138.1_LSUParc_tax_silva.fasta.gz\ncurl -O https://ftp.arb-silva.de/release_138.1/Exports/SILVA_138.1_SSUParc_tax_silva.fasta.gz\n\n# Swiss-Prot data from UniProt for diamond makedb process (accessed on 17,September,2025)\ncurl -O https://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.fasta.gz\n\n# Pfam data from InterPro (accessed on 17,September,2025)) for hmmscan proess. Appling HMMER process in this workflow is on going, however this process takes time. This step will be optional.\n# curl -O https://ftp.ebi.ac.uk/pub/databases/Pfam/current_release/Pfam-A.hmm.gz\n```\n\n&nbsp;\n\n### 6. Command Execution\n\nWe recommend creating a `cache` directory to store cache and intermediate files. Since metagenomic and metatranscriptomic reads are mapped to contigs, the assembled results can be reused to reduce analytical costs. The `cwltool` properly recognizes caches when the `--cachedir` option is specified.\n\n```bash\n# main workflow\n\ncwltool --debug --cachedir  --outdir  ./Worlkflow/main_w.cwl ./config/main_w.yml\n\n```\n\n&nbsp;\n\n### 7. based shell script &amp; python script\n\nGitHub: https://github.com/RyoMameda/workflow\n\nThis workflow is developed at [DBCLS BioHackathon 2025](https://github.com/dbcls/bh25/), and the preprint of developing project is https://doi.org/10.37044/osf.io/qd5sz_v1.\n",
      "isBasedOn": "https://github.com/RyoMameda/ComplexMicrobiome_GeneExpression_CWL.git",
      "sdDatePublished": "2026-04-18 05:31:33 +0100",
      "creativeWorkStatus": "Work-in-progress",
      "hasPart": [
        {
          "@id": "Worlkflow/annotation_sw.cwl"
        },
        {
          "@id": "workflow-1955-95a536b953f51f802b0c80729fefc1cb765d10f1-diagram.svg"
        }
      ]
    },
    {
      "@id": "Worlkflow/annotation_sw.cwl",
      "@type": [
        "File",
        "SoftwareSourceCode",
        "ComputationalWorkflow"
      ],
      "contentSize": 4342,
      "programmingLanguage": {
        "@id": "#cwl"
      },
      "image": {
        "@id": "workflow-1955-95a536b953f51f802b0c80729fefc1cb765d10f1-diagram.svg"
      },
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/ComputationalWorkflow/1.0-RELEASE/"
      },
      "description": "\n![GitHub last commit (branch)](https://img.shields.io/github/last-commit/RyoMameda/workflow_cwl/main)\n![Status](https://img.shields.io/badge/status-development-yellow)\n[![cwltool](https://img.shields.io/badge/cwltool-3.1.20250110105449-success)](https://github.com/common-workflow-language/cwltool/releases/tag/3.1.20250110105449)\n[![License](https://img.shields.io/badge/License-MIT-blue.svg)](./LICENSE)\n[![Version](https://img.shields.io/badge/version-1.0-brightgreen)](https://github.com/RyoMameda/ComplexMicrobiome_GeneExpression_CWL/releases/tag/v1.0)\n[![Open in Dev Containers](https://img.shields.io/static/v1?label=Dev%20Containers&amp;message=python3.11&amp;color=blue&amp;logo=docker)](https://github.com/yonesora56/plant2human/tree/main/.devcontainer)\n\n# Gene Expression Analysis Workflow in Complex Microbiomes\n\n&nbsp;\n\n## Workflow Schema \n- more details: [Optimization of Mapping Tools and Investigation of Ribosomal RNA Influence for Data-Driven Gene Expression Analysis in Complex Microbiomes](https://doi.org/10.3390/microorganisms13050995)\n\n![image](./image/microorganisms-13-00995-g001.png)\n\n&nbsp;\n\n### 1. Overview of the Workflow\n\nThis analysis focuses on transcriptional profiling of complex microbiomes. It requires both metagenomic and metatranscriptomic NGS short-read data, along with annotation reference information (e.g., ribosomal RNA sequences and referenced protein databases, listed below). The metagenomic and metatranscriptomic reads should be derived from the same microbiome samples. Assembled metagenomic contigs are then used as reference sequences to map both types of reads, enabling gene-level quantification.\n\n&nbsp;\n\n### 2. Minimum Requirements\n\n- `Docker`\n- `cwltool`\n\n&nbsp;\n\n### 3. Workflow Component\n\nThis analysis workflow is composed of three sub-workflows; metagenomic contig assembling, reads mapping and annotation. \n\n&nbsp;\n\n#### Metagenomic contig assembling\n\nIn this process, the following steps are performed:\n\n\n1. Assembly process using `Megahit`. \n2. Prediction Protein sequences using `Prodigal`.\n3. Statical analysis of contigs useing `SeqKit`.\n\n&nbsp;\n\n#### Reads mapping\n\nIn this process, the following steps are performed:\n\n1. Mapping process using `BWA MEM`.\n2. Statical analysis of mapping results using `SAMtools`\n\n&nbsp;\n\n#### Annotation\n\nIn this process, the following steps are performed:\n\n1. Searching contaminated ribosomal RNA sequences using `BLAST`.\n2. Searching referenced proteins using `DIAMOND`.\n3. Creation GTF formated file contained annotation informations.\n\n&nbsp;\n\n### 4. Test Dataset and Your Own Dataset\n\n- If you are testing with the following files, please place them in the `Data` directory!\n- You can also obtain metagenomic and metatranscriptomic FASTQ files either by downloading them from public databases or by using your own samples, and then place them in your `Data` directory.\n\n#### Metagenome data\n\n- [SRR27548858](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548858)\n\n#### Metatranscriptome data\n\n- [SRR27548863](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548863)\n- [SRR27548864](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548864)\n- [SRR27548865](https://www.ncbi.nlm.nih.gov/sra/?term=SRR27548865)\n\n&nbsp;\n\n### 5. Annotation References\n\nThese reference files are used in the BLAST and DIAMOND processes. The downloaded files are available in the `Data` directory (accessed on September 17, 2025). If you wish to use the latest versions of the references, please download them using the following scripts.\n\n```bash\n# rRNA data from SILVA website (release138.1; accessed on 17,September,2025)\ncurl -O https://ftp.arb-silva.de/release_138.1/Exports/SILVA_138.1_LSUParc_tax_silva.fasta.gz\ncurl -O https://ftp.arb-silva.de/release_138.1/Exports/SILVA_138.1_SSUParc_tax_silva.fasta.gz\n\n# Swiss-Prot data from UniProt for diamond makedb process (accessed on 17,September,2025)\ncurl -O https://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.fasta.gz\n\n# Pfam data from InterPro (accessed on 17,September,2025)) for hmmscan proess. Appling HMMER process in this workflow is on going, however this process takes time. This step will be optional.\n# curl -O https://ftp.ebi.ac.uk/pub/databases/Pfam/current_release/Pfam-A.hmm.gz\n```\n\n&nbsp;\n\n### 6. Command Execution\n\nWe recommend creating a `cache` directory to store cache and intermediate files. Since metagenomic and metatranscriptomic reads are mapped to contigs, the assembled results can be reused to reduce analytical costs. The `cwltool` properly recognizes caches when the `--cachedir` option is specified.\n\n```bash\n# main workflow\n\ncwltool --debug --cachedir  --outdir  ./Worlkflow/main_w.cwl ./config/main_w.yml\n\n```\n\n&nbsp;\n\n### 7. based shell script &amp; python script\n\nGitHub: https://github.com/RyoMameda/workflow\n\nThis workflow is developed at [DBCLS BioHackathon 2025](https://github.com/dbcls/bh25/), and the preprint of developing project is https://doi.org/10.37044/osf.io/qd5sz_v1.\n",
      "name": "Gene Expression Analysis Workflow in Complex Microbiomes",
      "url": "https://workflowhub.eu/workflows/1955?version=1",
      "keywords": "Annotation, Metagenomics, metatranscriptomics",
      "version": 1,
      "license": "https://spdx.org/licenses/MIT",
      "creator": [
        {
          "@id": "https://workflowhub.eu/people/1195"
        },
        {
          "@id": "https://workflowhub.eu/people/805"
        }
      ],
      "producer": {
        "@id": "https://workflowhub.eu/projects/283"
      },
      "dateCreated": "2025-09-19T06:27:20Z",
      "dateModified": "2025-10-02T13:21:26Z",
      "isPartOf": [

      ],
      "identifier": "https://doi.org/10.48546/workflowhub.workflow.1955.1",
      "isBasedOn": "https://workflowhub.eu/workflows/1955?version=1",
      "input": [
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_THREADS"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_EVALUE"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_BLASTN_rRNA_FASTA_FILE1"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_BLASTN_rRNA_FASTA_FILE2"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_PRODIGAL_RESULT_DNA_FASTA_FILE"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_PRODIGAL_RESULT_PROTEIN_FASTA_FILE"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_DIAMOND_INDEX_FILE"
        },
        {
          "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_OUTPUT_GTF_FILE_NAME"
        }
      ],
      "output": {
        "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-outputs-OUTPUT_GTF_FILE"
      },
      "sdPublisher": {
        "@id": "https://about.workflowhub.eu/"
      }
    },
    {
      "@id": "workflow-1955-95a536b953f51f802b0c80729fefc1cb765d10f1-diagram.svg",
      "@type": [
        "File",
        "ImageObject"
      ],
      "contentSize": 13477
    },
    {
      "@id": "https://orcid.org/0009-0007-5830-6482",
      "@type": "Person",
      "name": "Ryo Mameda",
      "identifier": "https://orcid.org/0009-0007-5830-6482"
    },
    {
      "@id": "https://orcid.org/0009-0004-1874-3117",
      "@type": "Person",
      "name": "Sora Yonezawa",
      "identifier": "https://orcid.org/0009-0004-1874-3117"
    },
    {
      "@id": "https://workflowhub.eu/people/1195",
      "@type": "Person",
      "name": "Ryo Mameda"
    },
    {
      "@id": "https://workflowhub.eu/people/805",
      "@type": "Person",
      "name": "Sora Yonezawa"
    },
    {
      "@id": "https://workflowhub.eu/projects/283",
      "@type": [
        "Project",
        "Organization"
      ],
      "name": "bonohulab"
    },
    {
      "@id": "#cwl",
      "@type": "ComputerLanguage",
      "name": "Common Workflow Language",
      "alternateName": "CWL",
      "identifier": {
        "@id": "https://w3id.org/cwl/v1.0/"
      },
      "url": {
        "@id": "https://www.commonwl.org/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_THREADS",
      "@type": "FormalParameter",
      "name": "threads",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_EVALUE",
      "@type": "FormalParameter",
      "name": "evalue",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_BLASTN_rRNA_FASTA_FILE1",
      "@type": "FormalParameter",
      "name": "SILVA_138.1_LSUParc_tax_silva",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_BLASTN_rRNA_FASTA_FILE2",
      "@type": "FormalParameter",
      "name": "SILVA_138.1_SSUParc_tax_silva",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_PRODIGAL_RESULT_DNA_FASTA_FILE",
      "@type": "FormalParameter",
      "name": "input fasta file (nucleotide sequence generated by prodigal process)",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_PRODIGAL_RESULT_PROTEIN_FASTA_FILE",
      "@type": "FormalParameter",
      "name": "predicted protein coding sequences produced by Prodigal process",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_DIAMOND_INDEX_FILE",
      "@type": "FormalParameter",
      "name": "Protein fasta file for diamond index",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-inputs-SW_OUTPUT_GTF_FILE_NAME",
      "@type": "FormalParameter",
      "name": "Output GTF file name",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "#gene_expression_analysis_workflow_in_complex_microbiomes-outputs-OUTPUT_GTF_FILE",
      "@type": "FormalParameter",
      "name": "Output GTF file",
      "dct:conformsTo": {
        "@id": "https://bioschemas.org/profiles/FormalParameter/1.0-RELEASE/"
      }
    },
    {
      "@id": "https://about.workflowhub.eu/",
      "@type": "Organization",
      "name": "WorkflowHub",
      "url": "https://about.workflowhub.eu/"
    }
  ]
}