Add profile for ARM compatibility

docs/usage.md

@@ -116,6 +116,10 @@ If you would like to reduce the number of reads used in the analysis, for exampl
 
 ## Alignment options
 
+:::note
+Not currently available using ARM architecture ('-profile arm')
+:::
+
 By default, the pipeline uses [STAR](https://github.com/alexdobin/STAR) (i.e. `--aligner star_salmon`) to map the raw FastQ reads to the reference genome, project the alignments onto the transcriptome and to perform the downstream BAM-level quantification with [Salmon](https://salmon.readthedocs.io/en/latest/salmon.html). STAR is fast but requires a lot of memory to run, typically around 38GB for the Human GRCh37 reference genome. Since the [RSEM](https://github.com/deweylab/RSEM) (i.e. `--aligner star_rsem`) workflow in the pipeline also uses STAR you should use the [HISAT2](https://ccb.jhu.edu/software/hisat2/index.shtml) aligner (i.e. `--aligner hisat2`) if you have memory limitations.
 
 You also have the option to pseudoalign and quantify your data directly with [Salmon](https://salmon.readthedocs.io/en/latest/salmon.html) or [Kallisto](https://pachterlab.github.io/kallisto/) by specifying `salmon` or `kallisto` to the `--pseudo_aligner` parameter. The selected pseudoaligner will then be run in addition to the standard alignment workflow defined by `--aligner`, mainly because it allows you to obtain QC metrics with respect to the genomic alignments. However, you can provide the `--skip_alignment` parameter if you would like to run Salmon or Kallisto in isolation. By default, the pipeline will use the genome fasta and gtf file to generate the transcripts fasta file, and then to build the Salmon index. You can override these parameters using the `--transcript_fasta` and `--salmon_index` parameters, respectively.
@@ -298,6 +302,10 @@ By default, the input GTF file will be filtered to ensure that sequence names co
 
 ## Contamination screening options
 
+:::note
+Not currently available using ARM architecture ('-profile arm')
+:::
+
 The pipeline provides the option to scan unaligned reads for contamination from other species using [Kraken2](https://ccb.jhu.edu/software/kraken2/), with the possibility of applying corrections from [Bracken](https://ccb.jhu.edu/software/bracken/). Since running Bracken is not computationally expensive, we recommend always using it to refine the abundance estimates generated by Kraken2.
 
 It is important to note that the accuracy of Kraken2 is [highly dependent on the database](https://doi.org/10.1099/mgen.0.000949) used. Specifically, it is [crucial](https://doi.org/10.1128/mbio.01607-23) to ensure that the host genome is included in the database. If you are particularly concerned about certain contaminants, it may be beneficial to use a smaller, more focused database containing primarily those contaminants instead of the full standard database. Various pre-built databases [are available for download](https://benlangmead.github.io/aws-indexes/k2), and instructions for building a custom database can be found in the [Kraken2 documentation](https://github.com/DerrickWood/kraken2/blob/master/docs/MANUAL.markdown). Additionally, genomes of contaminants detected in previous sequencing experiments are available on the [OpenContami website](https://openlooper.hgc.jp/opencontami/help/help_oct.php).
@@ -420,6 +428,7 @@ If `-profile` is not specified, the pipeline will run locally and expect all sof
   - A generic configuration profile to enable [Wave](https://seqera.io/wave/) containers. Use together with one of the above (requires Nextflow ` 24.03.0-edge` or later).
 - `conda`
   - A generic configuration profile to be used with [Conda](https://conda.io/docs/). Please only use Conda as a last resort i.e. when it's not possible to run the pipeline with Docker, Singularity, Podman, Shifter, Charliecloud, or Apptainer.
+  - A configuration profile that will set `docker.runOptions` appropriately for ARM architectures, and apply overrides supplying ARM-compatible containers and Conda environments.
 
 ### `-resume`
 

modules.json

@@ -245,7 +245,7 @@
                     },
                     "trimgalore": {
                         "branch": "master",
-                        "git_sha": "49f4e50534fe4b64101e62ea41d5dc43b1324358",
+                        "git_sha": "8c5eeedd45e295fc9a4f164631da6a8b37e6b9c6",
                         "installed_by": ["fastq_fastqc_umitools_trimgalore"]
                     },
                     "tximeta/tximport": {
@@ -333,7 +333,7 @@
                     },
                     "fastq_fastqc_umitools_trimgalore": {
                         "branch": "master",
-                        "git_sha": "49f4e50534fe4b64101e62ea41d5dc43b1324358",
+                        "git_sha": "8c5eeedd45e295fc9a4f164631da6a8b37e6b9c6",
                         "installed_by": ["fastq_qc_trim_filter_setstrandedness", "subworkflows"]
                     },
                     "fastq_qc_trim_filter_setstrandedness": {

nextflow.config

@@ -186,6 +186,7 @@ profiles {
     }
     arm {
         docker.runOptions       = '-u $(id -u):$(id -g) --platform=linux/amd64'
+        includeConfig 'conf/arm.config'
     }
     singularity {
         singularity.enabled     = true

subworkflows/local/align_star/nextflow.config

@@ -1,36 +1,55 @@
+def generateStarAlignArgs(save_unaligned, contaminant_screening, extra_star_align_args) {
+    def argsToMap = { String args ->
+        args.split(/\s(?=--)/).collectEntries {
+            def parts = it.trim().split(/\s+/, 2)
+            [(parts[0]): parts.size() > 1 ? parts[1] : '']
+        }
+    }
+
+    def base_args = """
+        --quantMode TranscriptomeSAM
+        --twopassMode Basic
+        --outSAMtype BAM Unsorted
+        --readFilesCommand zcat
+        --runRNGseed 0
+        --outFilterMultimapNmax 20
+        --alignSJDBoverhangMin 1
+        --outSAMattributes NH HI AS NM MD
+        --outSAMstrandField intronMotif
+    """.trim()
+
+    if (save_unaligned || contaminant_screening) {
+        base_args += "\n--outReadsUnmapped Fastx"
+    }
+
+    def final_args_map = argsToMap(base_args) + (extra_star_align_args ? argsToMap(extra_star_align_args) : [:])
+    final_args_map.collect { key, value -> "${key} ${value}".trim() }.join(' ')
+}
+
 if (!params.skip_alignment && params.aligner == 'star_salmon') {
     process {
-        withName: '.*:ALIGN_STAR:STAR_ALIGN|.*:ALIGN_STAR:STAR_ALIGN_IGENOMES' {
-            ext.args   = {
-                    // Function to convert argument strings into a map
-                    def argsToMap = { String args ->
-                        args.split("\\s(?=--)").collectEntries {
-                            def parts = it.trim().split(/\s+/, 2)
-                            [(parts.first()): parts.last()]
-                        }
-                    }
-
-                    // Initialize the map with preconfigured values
-                    def preset_args_map = argsToMap("""
-                        --quantMode TranscriptomeSAM
-                        --twopassMode Basic
-                        --outSAMtype BAM Unsorted
-                        --readFilesCommand zcat
-                        --runRNGseed 0
-                        --outFilterMultimapNmax 20
-                        --alignSJDBoverhangMin 1
-                        --outSAMattributes NH HI AS NM MD
-                        --quantTranscriptomeSAMoutput BanSingleEnd
-                        --outSAMstrandField intronMotif
-                        ${params.save_unaligned || params.contaminant_screening ? '--outReadsUnmapped Fastx' : ''}
-                    """.trim())
 
-                    // Consolidate the extra arguments
-                    def final_args_map = preset_args_map + (params.extra_star_align_args ? argsToMap(params.extra_star_align_args) : [:])
+        // We have to condition this, because the args are slightly different between the latest STAR and the one compatible with iGenomes
 
-                    // Convert the map back to a list and then to a single string
-                    final_args_map.collect { key, value -> "${key} ${value}" }.join(' ').trim()
+        withName: '.*:ALIGN_STAR:STAR_ALIGN' {
+            ext.args   = {
+                generateStarAlignArgs(
+                    params.save_unaligned,
+                    params.contaminant_screening,
+                    (params.extra_star_align_args ?: '') + ' --quantTranscriptomeSAMoutput BanSingleEnd'
+                )
+            }
+        }
+        withName: '.*:ALIGN_STAR:STAR_ALIGN_IGENOMES' {
+            ext.args   = {
+                generateStarAlignArgs(
+                    params.save_unaligned,
+                    params.contaminant_screening,
+                    (params.extra_star_align_args ?: '') + ' --quantTranscriptomeBan Singleend'
+                )
             }
+        }
+        withName: '.*:ALIGN_STAR:STAR_ALIGN|.*:ALIGN_STAR:STAR_ALIGN_IGENOMES' {
 
             publishDir = [
                 [