Rsubread处理RNA Bulk Seq配对数据：featureCounts无结果问题排查

核心问题分析 #

你遇到的featureCounts计数为0的情况，大概率是比对路径错误或参考序列与注释ID不匹配导致的，以下是具体排查和修复步骤：

1. 紧急修复比对代码的路径错误 #

你的比对代码中拼接read路径时遗漏了FASTQ_Raw目录：
原错误代码片段：

read1_file <- paste0("D:/", prefix, "_1.fq")
read2_file <- paste0("D:/", prefix, "_2.fq")

这会导致程序去D:/根目录查找fq文件，但你的原始文件实际在D:/FASTQ_Raw/下，修正后：

read1_file <- paste0("D:/FASTQ_Raw/", prefix, "_1.fq")
read2_file <- paste0("D:/FASTQ_Raw/", prefix, "_2.fq")

重新运行比对生成正确BAM文件后，再尝试featureCounts计数。

2. 验证BAM文件的有效性 #

先确认BAM中是否存在成功比对的reads：

• 在R中用Rsamtools检查：

library(Rsamtools)
# 替换为你的目标BAM文件路径
bam_stats <- scanBamFlagstat("aligned_reads/sample_aligned.bam")
print(bam_stats)

查看mapped对应的数值：如果为0，说明比对本身未成功，需优先解决比对问题；如果有映射reads，再排查注释匹配问题。

3. 解决注释与参考序列的ID匹配问题 #

你用NCBI转录组序列构建的索引，序列ID是NCBI转录本编号（如NM_001301717），但featureCounts内置的hg38注释采用Ensembl格式ID（如ENST00000373020），两者无法匹配导致计数为0。解决方法二选一：

• 方案一：使用NCBI官方的hg38注释文件（GTF格式），确保注释中的转录本ID与索引的转录组序列ID完全一致，在featureCounts中指定该文件：

fc <- featureCounts(files = "aligned_reads/sample_aligned.bam", 
                    annot.ext = "path/to/ncbi_hg38.gtf",
                    isGTFAnnotationFile = TRUE,
                    nthreads = 4)

• 方案二：重新用hg38基因组序列构建索引（而非转录组），重新比对后再使用内置hg38注释计数，此时参考序列与注释的染色体ID、特征ID可完全匹配。

4. 适配转录组比对的featureCounts参数 #

如果是基于转录组索引的比对，需调整参数适配转录本水平计数：

fc <- featureCounts(files = "aligned_reads/sample_aligned.bam",
                    annot.ext = "path/to/ncbi_hg38.gtf",
                    useMetaFeatures = FALSE, # 直接对转录本计数而非基因
                    nthreads = 4)

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