diff --git a/Degradome-seq/Degradome-seq.snakefile b/Degradome-seq/Degradome-seq.snakefile
index ccd8c3859f4c6a8edf021dafce767d05516bdd73..3d1b6e54982c637d8caf258891e88289a03e34bc 100644
--- a/Degradome-seq/Degradome-seq.snakefile
+++ b/Degradome-seq/Degradome-seq.snakefile
@@ -130,7 +130,7 @@ rule all:
         #     lib=LIBS, biotype=["alltypes"], orientation=ORIENTATIONS),
         expand(
             OPJ(counts_dir, "diff_%s" % eff_name, "{contrast}", "{orientation}_{biotype}", "{contrast}_diff_%s.txt" % eff_name),
-            contrast=DD_CONTRASTS, orientation=["fwd"], biotype=["alltypes"]),
+            contrast=DD_CONTRASTS, orientation=["fwd"], biotype=["alltypes", "protein_coding"]),
 
 
 include: ensure_relative(irules["link_raw_data"], workflow.basedir)
diff --git a/PRO-seq/PRO-seq.snakefile b/PRO-seq/PRO-seq.snakefile
index f8717af88fc945416e4a494bcf025fda8556a7ab..9c290c9accfa796df3a58c768892394e6ad0fd53 100644
--- a/PRO-seq/PRO-seq.snakefile
+++ b/PRO-seq/PRO-seq.snakefile
@@ -357,7 +357,7 @@ rule all:
         expand(OPJ(
             "{trimmer}", aligner, "mapped_C_elegans", "{counter}",
             "all_on_C_elegans", "{biotype}_{orientation}_TPM.txt"),
-            trimmer=TRIMMERS, counter=COUNTERS, biotype=["alltypes"], orientation=ORIENTATIONS),
+            trimmer=TRIMMERS, counter=COUNTERS, biotype=["alltypes", "protein_coding"], orientation=ORIENTATIONS),
         #expand(OPJ("{trimmer}", "figures", aligner, "{lib}_mean", "{orientation}_on_merged_isolated_%d_{biotype}_min_%d_meta_profile.pdf" % (MIN_DIST, META_MIN_LEN)), trimmer=TRIMMERS, lib=LIBS, orientation=["all"], biotype=["protein_coding"]),
         #expand(OPJ("{trimmer}", "figures", aligner, "{lib}_mean", "{orientation}_on_merged_isolated_%d_{biotype}_min_%d_meta_profile.pdf" % (MIN_DIST, META_MIN_LEN)), trimmer=TRIMMERS, lib=LIBS, orientation=["all"], biotype=METAGENE_BIOTYPES),
         # TODO: Add metagene profiles similar to small RNA-seq
diff --git a/RNA_Seq_Cecere/RNA-seq.snakefile b/RNA_Seq_Cecere/RNA-seq.snakefile
index a420f593126c87d8e5527e0b4b06cee8d7c11e2b..0853fde408276cf43ba60d89bc855d2d60bad830 100644
--- a/RNA_Seq_Cecere/RNA-seq.snakefile
+++ b/RNA_Seq_Cecere/RNA-seq.snakefile
@@ -318,7 +318,7 @@ counts_files = [
         OPJ(aligner, f"mapped_{genome}", "{counter}",
             "all_{mapped_type}", "{biotype}_{orientation}_TPM.txt"),
         counter=COUNTERS, mapped_type=[f"on_{genome}", f"unique_on_{genome}"],
-        biotype=["alltypes"], orientation=ORIENTATIONS),
+        biotype=["alltypes", "protein_coding"], orientation=ORIENTATIONS),
     expand(
         OPJ(aligner, f"mapped_{genome}", "{counter}",
             "all_{mapped_type}", "{biotype}_{orientation}_RPK.txt"),
diff --git a/Ribo-seq/Ribo-seq.snakefile b/Ribo-seq/Ribo-seq.snakefile
index f48f972c35cbf3ba301dd4c7c657b8a0b63fd2fb..7e60128343f0bd70388fcbaf291c7932d848911b 100644
--- a/Ribo-seq/Ribo-seq.snakefile
+++ b/Ribo-seq/Ribo-seq.snakefile
@@ -579,7 +579,7 @@ rule all:
         # translation_efficiency
         expand(
             OPJ(counts_dir, "{lib}_mean_{read_type}_on_%s" % genome, "{lib}_{biotype}_{orientation}_%s.txt" % eff_name),
-            lib=EFF_LIBS, read_type=["RPF"], biotype=["alltypes"], orientation=["fwd"]),
+            lib=EFF_LIBS, read_type=["RPF"], biotype=["alltypes", "protein_coding"], orientation=["fwd"]),
         # DESeq2 results
         expand(
             OPJ(
@@ -599,7 +599,7 @@ rule all:
                 counts_dir,
                 "diff_%s_{read_type}" % eff_name,
                 "{contrast}", "{orientation}_{biotype}", "{contrast}_diff_%s.txt" % eff_name),
-            read_type=["RPF"], contrast=DT_CONTRASTS, orientation=["fwd"], biotype=["alltypes"]),
+            read_type=["RPF"], contrast=DT_CONTRASTS, orientation=["fwd"], biotype=["alltypes", "protein_coding"]),
 
 # rule future_all:
 #         meta_profiles,