Nextflow流程开发求助：基因组区域批量处理与按染色体合并

现有22个命名为chr1.txt至chr22.txt的文件，每个文件包含基因组区域数据，需通过Nextflow实现：

• 对每个区域独立运行process1
• 可单独跟踪每个区域的处理进度
• 待同染色体所有区域处理完成后，将结果合并为merged_chr1.txt、merged_chr2.txt等对应文件

原脚本因过早使用groupTuple将同染色体区域打包，导致process1仅按染色体批量处理，无法跟踪单个区域进度，需修正。

示例区域文件内容（chr10.txt） #

0       chr10   chr10:1-16637292
1       chr10   chr10:15435012-33506537
2       chr10   chr10:30788681-57663231
3       chr10   chr10:54643594-75924707
4       chr10   chr10:71749908-93488832
5       chr10   chr10:91114596-112850006

params.chunks_dir = "/workspace/gitpod/testing/chr*.txt"
params.results    = "/workspace/gitpod/testing/results"

workflow {
    // 读取区域文件，拆分出每个单独的区域记录
    ch_regions = Channel.fromPath(params.chunks_dir)
                    .splitCsv(header: false, sep: "\t")
                    .map{ chunk_no, chr, region -> tuple(chr, chunk_no, region) }

    // 对每个区域单独运行process1
    ch_processed = run_process1(ch_regions)

    // 按染色体分组，聚合同染色体的所有处理结果
    ch_processed
        .groupTuple()
        .set{ ch_chromosome_results }

    // 运行合并流程
    merge_results(ch_chromosome_results)
}

process run_process1 {
    tag "${chr}:${region}"  // 标记任务，方便跟踪单个区域的处理进度
    publishDir "${params.results}/processed", mode: 'copy'

    input:
    tuple val(chr), val(chunk_no), val(region)

    output:
    tuple val(chr), path("ABCD_${chr}.chunk_${chunk_no}.test.txt")

    script:
    """
    echo "${region}" > ABCD_${chr}.chunk_${chunk_no}.test.txt
    """
}

process merge_results {
    tag "${chr}"
    publishDir "${params.results}/merged", mode: 'copy'

    input:
    tuple val(chr), path(chunk_files)

    output:
    path "merged_${chr}.txt"

    script:
    """
    cat ${chunk_files.join(' ')} > merged_${chr}.txt
    """
}

• 调整通道逻辑：原脚本过早用groupTuple打包同染色体区域，修改后先将每个区域拆分为独立的(chr, chunk_no, region)元组，确保run_process1为每个区域生成独立任务，Nextflow控制台会显示每个区域的处理状态，实现进度跟踪。
• 添加任务标记：给run_process1设置tag "${chr}:${region}"，执行时能直观看到每个区域的任务启动、完成情况。
• 分阶段流程：拆分为三个清晰阶段：读取拆分区域 → 单独处理每个区域 → 按染色体合并结果，完全匹配需求。
• 规范输出目录：将中间处理文件和最终合并结果分别发布到processed和merged子目录，避免文件混乱。

内容的提问来源于stack exchange，提问作者kinolodeon