Skip to content

Commit 06a46f1

Browse files
authored
Merge pull request #204 from bnovak32/DEV_NF_MAAgilent_1ch
- Refactor custom annotation for Agilent catalog arrays - pulls ProbeID -> Ensembl transcript ID mapping from Agilent AllAnnotations file - uses BioMart to map from transcript to Ensembl probe ID - uses standard GeneLab annotations for all other annotation assignment - Updated outdated links (nextflow, and Genome Reference) - Updated syntax for notes/tips - Removed software in pipeline doc that is only used in the workflow (it remains in the licenses section) - automated TOC generation - fixed typos in documentation
2 parents d6d864a + 5fec3dd commit 06a46f1

6 files changed

Lines changed: 338 additions & 271 deletions

File tree

Microarray/Agilent_1-channel/Pipeline_GL-DPPD-7112_Versions/GL-DPPD-7112-A.md

Lines changed: 208 additions & 176 deletions
Large diffs are not rendered by default.

Microarray/Agilent_1-channel/README.md

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -26,5 +26,6 @@
2626
Jonathan Oribello
2727

2828
**Maintained by:**
29-
Crystal Han (crystal.han@nasa.gov)
29+
Jihan Yehia (jihan.yehia@nasa.gov)
3030
Alexis Torres (alexis.torres@nasa.gov)
31+
Crystal Han (crystal.han@nasa.gov)

Microarray/Agilent_1-channel/Workflow_Documentation/NF_MAAgilent1ch/README.md

Lines changed: 20 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@
44

55
### Implementation Tools <!-- omit in toc -->
66

7-
The current GeneLab Agilent 1 Channel Microarray consensus processing pipeline (NF_MAAgilent1ch), [GL-DPPD-7112-A](../../Pipeline_GL-DPPD-7112_Versions/GL-DPPD-7112-A.md), is implemented as a [Nextflow](https://nextflow.io/) DSL2 workflow and utilizes [Singularity](https://docs.sylabs.io/guides/3.10/user-guide/introduction.html) to run all tools in containers. This workflow (NF_MAAgilent1ch) is run using the command line interface (CLI) of any unix-based system. While knowledge of creating workflows in Nextflow is not required to run the workflow as is, [the Nextflow documentation](https://nextflow.io/docs/latest/index.html) is a useful resource for users who want to modify and/or extend this workflow.
7+
The current GeneLab Agilent 1 Channel Microarray consensus processing pipeline (NF_MAAgilent1ch), [GL-DPPD-7112-A](../../Pipeline_GL-DPPD-7112_Versions/GL-DPPD-7112-A.md), is implemented as a [Nextflow](https://nextflow.io/) DSL2 workflow and utilizes [Singularity](https://docs.sylabs.io/guides/3.10/user-guide/introduction.html) or [Docker](https://docs.docker.com/get-started/) containers to run all tools. This workflow (NF_MAAgilent1ch) is run using the command line interface (CLI) of any unix-based system. While knowledge of creating workflows in Nextflow is not required to run the workflow as is, [the Nextflow documentation](https://docs.seqera.io/nextflow/) is a useful resource for users who want to modify and/or extend this workflow.
88

99
### Workflow & Subworkflows <!-- omit in toc -->
1010

@@ -33,13 +33,13 @@ Below is a description of each subworkflow and the additional output files gener
3333
- This subworkflow performs validation and verification (V&V) on the raw and processed data files. It performs a series of checks on the output files generated and flags the results, using the flag codes indicated in the table below, which are outputted into a log file.
3434
**V&V Flags**:
3535

36-
|Flag Codes|Flag Name|Interpretation|
37-
|:---------|:--------|:-------------|
38-
| 2 | MANUAL | Special flag that indicates a manual check that is advised. Often used to advise what assess in QA plots. |
39-
| 20 | GREEN | Indicates the check passed all validation conditions |
40-
| 30 | YELLOW | Indicates the check was flagged for minor issues (e.g. slight outliers) |
41-
| 50 | RED | Indicates the check was flagged for moderate issues (e.g. major outliers) |
42-
| 80 | HALT | Indicates the check was flagged for severe issues that trigger a processing halt (e.g. missing data) |
36+
| Flag Codes | Flag Name | Interpretation |
37+
| :--------- | :-------- | :-------------------------------------------------------------------------------------------------------- |
38+
| 2 | MANUAL | Special flag that indicates a manual check that is advised. Often used to advise what assess in QA plots. |
39+
| 20 | GREEN | Indicates the check passed all validation conditions |
40+
| 30 | YELLOW | Indicates the check was flagged for minor issues (e.g. slight outliers) |
41+
| 50 | RED | Indicates the check was flagged for moderate issues (e.g. major outliers) |
42+
| 80 | HALT | Indicates the check was flagged for severe issues that trigger a processing halt (e.g. missing data) |
4343

4444
<br>
4545

@@ -64,9 +64,11 @@ Below is a description of each subworkflow and the additional output files gener
6464

6565
#### 1a. Install Nextflow
6666

67-
Nextflow can be installed either through [Anaconda](https://anaconda.org/bioconda/nextflow) or as documented on the [Nextflow documentation page](https://www.nextflow.io/docs/latest/getstarted.html).
67+
Nextflow can be installed either through the [Anaconda bioconda channel](https://anaconda.org/bioconda/nextflow) or as documented in the [Nextflow installation documentation](https://docs.seqera.io/nextflow/install).
6868

69-
> Note: If you want to install Anaconda, we recommend installing a Miniconda, Python3 version appropriate for your system, as instructed by [Happy Belly Bioinformatics](https://astrobiomike.github.io/unix/conda-intro#getting-and-installing-conda).
69+
70+
> [!TIP]
71+
> If you wish to install Conda, we recommend installing a Miniforge version appropriate for your system, as documented on the [conda-forge website](https://conda-forge.org/download/), where you can find basic binaries for most systems. More detailed miniforge documentation is available in the [miniforge github repository](https://github.com/conda-forge/miniforge).
7072
>
7173
> Once conda is installed on your system, you can install the latest version of Nextflow by running the following commands:
7274
>
@@ -83,7 +85,9 @@ Singularity is a container platform that allows usage of containerized software.
8385
8486
We recommend installing Singularity on a system wide level as per the associated [documentation](https://docs.sylabs.io/guides/3.10/admin-guide/admin_quickstart.html).
8587
86-
> Note: Singularity is also available through [Anaconda](https://anaconda.org/conda-forge/singularity).
88+
> [!TIP]
89+
> - Singularity is also available through the [Anaconda conda-forge channel](https://anaconda.org/conda-forge/singularity).
90+
> - Alternatively, Docker can be used in place of Singularity. To get started with Docker, see the [Docker CE installation documentation](https://docs.docker.com/engine/install/).
8791
8892
<br>
8993
@@ -106,7 +110,8 @@ unzip NF_MAAgilent1ch_1.0.5.zip
106110
### 3. Run the Workflow
107111
108112
While in the location containing the `NF_MAAgilent1ch_1.0.5` directory that was downloaded in [step 2](#2-download-the-workflow-files), you are now able to run the workflow. Below are three examples of how to run the NF_MAAgilent1ch workflow:
109-
> Note: Nextflow commands use both single hyphen arguments (e.g. -help) that denote general nextflow arguments and double hyphen arguments (e.g. --ensemblVersion) that denote workflow specific parameters. Take care to use the proper number of hyphens for each argument.
113+
> [!NOTE]
114+
> Nextflow commands use both single hyphen arguments (e.g. -help) that denote general nextflow arguments and double hyphen arguments (e.g. --ensemblVersion) that denote workflow specific parameters. Take care to use the proper number of hyphens for each argument.
110115
111116
<br>
112117
@@ -196,7 +201,7 @@ All parameters listed above and additional optional arguments for the NF_MAAgile
196201
nextflow run NF_MAAgilent1ch_1.0.5/main.nf --help
197202
```
198203
199-
See `nextflow run -h` and [Nextflow's CLI run command documentation](https://nextflow.io/docs/latest/cli.html#run) for more options and details common to all nextflow workflows.
204+
See `nextflow run -h` and [Nextflow's CLI run command documentation](https://docs.seqera.io/nextflow/cli) for more options and details common to all nextflow workflows.
200205
201206
<br>
202207
@@ -211,7 +216,7 @@ All R code steps and output are rendered within a Quarto document yielding the f
211216
212217
213218
The outputs from the Analysis Staging and V&V Pipeline Subworkflows are described below:
214-
> Note: The outputs from the Agilent 1 Channel Microarray Processing Subworkflow are documented in the [GL-DPPD-7112-A.md](../../../Pipeline_GL-DPPD-7112_Versions/GL-DPPD-7112-A.md) processing protocol.
219+
> Note: The outputs from the Agilent 1 Channel Microarray Processing Subworkflow are documented in the [GL-DPPD-7112-A.md](../../Pipeline_GL-DPPD-7112_Versions/GL-DPPD-7112-A.md) processing protocol.
215220
216221
**Analysis Staging Subworkflow**
217222
> Note: only applicable for [Approach 1](#3a-approach-1-run-the-workflow-on-a-genelab-agilent-1-channel-microarray-dataset) and [Approach 3](#3c-approach-3-run-the-workflow-using-an-isa-archive)
@@ -229,7 +234,7 @@ The outputs from the Analysis Staging and V&V Pipeline Subworkflows are describe
229234
<br>
230235
231236
Standard Nextflow resource usage logs are also produced as follows:
232-
> Further details about these logs can also found within [this Nextflow documentation page](https://www.nextflow.io/docs/latest/tracing.html#execution-report).
237+
> Further details about these logs can also found in the [Nextflow Report Documentation](https://docs.seqera.io/nextflow/reports).
233238
234239
**Nextflow Resource Usage Logs**
235240
- Output:

Microarray/Agilent_1-channel/Workflow_Documentation/NF_MAAgilent1ch/examples/annotations/README.md

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -16,8 +16,8 @@
1616
| Column Name | Type | Description | Example |
1717
|:------------|:-----|:------------|:--------|
1818
| array_design | string | A bioMart attribute identifier denoting the microarray probe/probeset attribute used for annotation mapping. | AGILENT SurePrint G3 GE 8x60k v3 |
19-
| annot_type | string | Used to determine how the custom annotations are parsed before merging to the data. Currently, only the below are supported: <ul><li>`agilent`: Annotations file is expected to be in the AA (All Annotations) format by [Agilent](https://earray.chem.agilent.com/earray/)</li><li>`custom`: Annotations file is merged as is, expected to have the following columns: `ProbesetID`, `ENTREZID`, `SYMBOL`, `GENENAME`, `ENSEMBL`, `REFSEQ`, `GOSLIM_IDS`, `STRING_id`, `count_gene_mappings`, `gene_mapping_source`</li></ul> | agilent |
20-
| annot_filename | string | Name of the custom annotations file. This is the AllAnnotations file downloaded from Agilent's eArray web portal. | 072363_D_AA_20240521.txt |
19+
| annot_type | string | Used to determine how the custom annotations are parsed before merging to the data. Currently, only the below are supported: <ul><li>`agilent`: Annotations file is expected to be in the AA (All Annotations) format created by [Agilent](https://earray.chem.agilent.com/earray/) for catalog arrays. The Ensembl Transcript Stable IDs in the file are used to map to Ensembl Gene Stable IDs after which the standard GeneLab annotation process is used.</li><li>`custom`: Annotations file is merged as-is and expected to have the following columns: `ProbeName`, `ENTREZID`, `SYMBOL`, `GENENAME`, `ENSEMBL`, `REFSEQ`, `GOSLIM_IDS`, `STRING_id`, `count_gene_mappings`, `gene_mapping_source`, where `ENSEMBL` is the Ensembl Gene Stable ID, not the Transcript ID</li></ul> | agilent |
20+
| annot_filename | string | Name of the custom annotations file. | 072363_D_AA_20240521.txt |
2121

2222
## Optional columns
2323
If the file was downloaded from a website, provide the download link used and date

0 commit comments

Comments
 (0)