@@ -29,6 +29,9 @@ if os.environ.get("CI", "false").lower() == "true":
29
29
cwd = os .getcwd ()
30
30
old_path = os .environ ["PATH" ]
31
31
os .environ ["PATH" ] = f"{ cwd } /test-mode-bin:{ old_path } "
32
+ RUNS_IN_CI = True
33
+ else :
34
+ RUNS_IN_CI = False
32
35
33
36
34
37
# ===============================================================================================
@@ -71,7 +74,6 @@ rule all:
71
74
f"work/genes/entrez/{ DV .today } /gene_info.jsonl" ,
72
75
f"work/genes/gnomad/{ DV .gnomad_constraints } /gnomad_constraints.tsv" ,
73
76
f"work/genes/hgnc/{ DV .today } /hgnc_info.jsonl" ,
74
- f"work/genes/hgnc/{ DV .today } /hgnc_xlink.tsv" ,
75
77
f"work/genes/mim2gene/{ DV .today } /mim2gene.tsv" ,
76
78
# reference-specific annotations
77
79
# -- background/population sequence variants and annotations thereof
@@ -96,38 +98,16 @@ rule all:
96
98
f"work/annos/grch38/seqvars/gnomad_mtdna/{ DV .gnomad_mtdna } /gnomad_mtdna.vcf.gz" ,
97
99
f"work/download/annos/grch38/seqvars/gnomad_exomes/{ DV .gnomad_v2 } /.done" ,
98
100
f"work/download/annos/grch38/seqvars/gnomad_genomes/{ DV .gnomad_v3 } /.done" ,
99
- # -- background/population structural variants and annoations thereof
100
- # ---- GRCh37
101
- f"work/annos/grch37/strucvars/dbvar/{ DV .dbvar } /dbvar.bed.gz" ,
102
- f"work/annos/grch37/strucvars/dgv/{ DV .dgv } /dgv.bed.gz" ,
103
- f"work/annos/grch37/strucvars/dgv_gs/{ DV .dgv_gs } /dgv_gs.bed.gz" ,
104
- f"work/annos/grch37/strucvars/exac/{ DV .exac_cnv } /exac.bed.gz" ,
105
- f"work/annos/grch37/strucvars/g1k/{ DV .g1k_svs } /g1k.bed.gz" ,
106
- f"work/annos/grch37/strucvars/gnomad/{ DV .gnomad_sv } /gnomad_sv.bed.gz" ,
107
- # ---- GRCh38
108
- f"work/annos/grch38/strucvars/dbvar/{ DV .dbvar } /dbvar.bed.gz" ,
109
- f"work/annos/grch38/strucvars/dgv/{ DV .dgv } /dgv.bed.gz" ,
110
- f"work/annos/grch38/strucvars/dgv_gs/{ DV .dgv_gs } /dgv_gs.bed.gz" ,
111
101
# NB: gnomAD-SV GRCh38 was announced end of 2020 but not released yet
112
102
# -- genome browser "features" (position-specific)
113
103
# ---- GRCh37
114
104
f"work/annos/grch37/features/cons/{ DV .ucsc_cons_37 } /ucsc_conservation.tsv" ,
115
105
f"work/annos/grch37/features/ensembl/{ DV .ensembl_37 } /ensembl_genes.bed.gz" ,
116
106
f"work/annos/grch37/features/refseq/{ DV .refseq_37 } /refseq_genes.bed.gz" ,
117
- "work/annos/grch37/features/tads/dixon2015/hesc.bed" ,
118
- f"work/annos/grch37/features/ucsc/{ DV .ucsc_genomic_super_dups_37 } /genomicSuperDups.bed.gz" ,
119
- f"work/annos/grch37/features/ucsc/{ DV .ucsc_rmsk_37 } /rmsk.bed.gz" ,
120
- f"work/annos/grch37/features/ucsc/{ DV .ucsc_alt_seq_liftover_37 } /altSeqLiftOverPsl.bed.gz" ,
121
- f"work/annos/grch37/features/ucsc/{ DV .ucsc_fix_seq_liftover_37 } /fixSeqLiftOverPsl.bed.gz" ,
122
107
# ---- GRCh38
123
108
f"work/annos/grch38/features/cons/{ DV .ucsc_cons_38 } /ucsc_conservation.tsv" ,
124
109
f"work/annos/grch38/features/ensembl/{ DV .ensembl_38 } /ensembl_genes.bed.gz" ,
125
110
f"work/annos/grch38/features/refseq/{ DV .refseq_38 } /refseq_genes.bed.gz" ,
126
- "work/annos/grch38/features/tads/dixon2015/hesc.bed" ,
127
- f"work/annos/grch38/features/ucsc/{ DV .ucsc_genomic_super_dups_38 } /genomicSuperDups.bed.gz" ,
128
- f"work/annos/grch38/features/ucsc/{ DV .ucsc_rmsk_38 } /rmsk.bed.gz" ,
129
- f"work/annos/grch38/features/ucsc/{ DV .ucsc_alt_seq_liftover_38 } /altSeqLiftOverPsl.bed.gz" ,
130
- f"work/annos/grch38/features/ucsc/{ DV .ucsc_fix_seq_liftover_38 } /fixSeqLiftOverPsl.bed.gz" ,
131
111
#
132
112
# == output directory ===================================================================
133
113
#
@@ -165,6 +145,41 @@ rule all:
165
145
# ---- UCSC conservation
166
146
f"output/worker/annos/seqvars/cons-grch37-{ DV .ucsc_cons_37 } +{ PV .annonars } /rocksdb/IDENTITY" ,
167
147
f"output/worker/annos/seqvars/cons-grch38-{ DV .ucsc_cons_38 } +{ PV .annonars } /rocksdb/IDENTITY" ,
148
+ # ----- Genes
149
+ f"output/worker/genes-{ DV .acmg_sf } +{ DV .gnomad_constraints } +{ DV .dbnsfp } +{ DV .today } +{ PV .worker } /rocksdb/IDENTITY" ,
150
+ f"output/worker/genes-xlink-{ DV .today } /genes-xlink.tsv" ,
151
+ f"output/worker/genes-txs-grch37-{ DV .mehari_tx } /mehari-data-txs-grch37-{ DV .mehari_tx } .bin.zst" ,
152
+ f"output/worker/genes-txs-grch38-{ DV .mehari_tx } /mehari-data-txs-grch38-{ DV .mehari_tx } .bin.zst" ,
153
+ # ----- HPO
154
+ f"output/viguno/hpo-{ DV .hpo } +{ PV .viguno } /hp.obo" ,
155
+ f"output/viguno/hpo-{ DV .hpo } +{ PV .viguno } /phenotype.hpoa" ,
156
+ f"output/viguno/hpo-{ DV .hpo } +{ PV .viguno } /phenotype_to_genes.txt" ,
157
+ f"output/viguno/hpo-{ DV .hpo } +{ PV .viguno } /hpo.bin" ,
158
+ # ----- background/population structural variants and annotations thereof
159
+ f"output/worker/annos/strucvars/dbvar-grch37-{ DV .dbvar } /dbvar.bed.gz" ,
160
+ f"output/worker/annos/strucvars/dbvar-grch38-{ DV .dbvar } /dbvar.bed.gz" ,
161
+ f"output/worker/annos/strucvars/dgv-grch37-{ DV .dgv } /dgv.bed.gz" ,
162
+ f"output/worker/annos/strucvars/dgv-grch38-{ DV .dgv } /dgv.bed.gz" ,
163
+ f"output/worker/annos/strucvars/dgv-gs-grch37-{ DV .dgv_gs } /dgv-gs.bed.gz" ,
164
+ f"output/worker/annos/strucvars/dgv-gs-grch38-{ DV .dgv_gs } /dgv-gs.bed.gz" ,
165
+ f"output/worker/annos/strucvars/exac-grch37-{ DV .exac_cnv } /exac.bed.gz" ,
166
+ f"output/worker/annos/strucvars/g1k-grch37-{ DV .g1k_svs } /g1k.bed.gz" ,
167
+ f"output/worker/annos/strucvars/gnomad-grch37-{ DV .gnomad_sv } /gnomad.bed.gz" ,
168
+ # ----- known pathogenic MMS
169
+ f"output/worker/annos/strucvars/patho-mms-grch37-{ DV .patho_mms } /patho-mms.bed" ,
170
+ f"output/worker/annos/strucvars/patho-mms-grch38-{ DV .patho_mms } /patho-mms.bed" ,
171
+ # ----- problematic regions (rmsk, genomicSuperDups, altSeqLiftOverPsl, fixSeqLiftOverPsl)
172
+ f"output/worker/annos/features/ucsc-genomicsuperdups-grch37-{ DV .ucsc_genomic_super_dups_37 } /genomicSuperDups.bed.gz" ,
173
+ f"output/worker/annos/features/ucsc-genomicsuperdups-grch38-{ DV .ucsc_genomic_super_dups_38 } /genomicSuperDups.bed.gz" ,
174
+ f"output/worker/annos/features/ucsc-rmsk-grch37-{ DV .ucsc_rmsk_37 } /rmsk.bed.gz" ,
175
+ f"output/worker/annos/features/ucsc-rmsk-grch38-{ DV .ucsc_rmsk_38 } /rmsk.bed.gz" ,
176
+ f"output/worker/annos/features/ucsc-altseqliftoverpsl-grch37-{ DV .ucsc_alt_seq_liftover_37 } /altSeqLiftOverPsl.bed.gz" ,
177
+ f"output/worker/annos/features/ucsc-altseqliftoverpsl-grch38-{ DV .ucsc_alt_seq_liftover_38 } /altSeqLiftOverPsl.bed.gz" ,
178
+ f"output/worker/annos/features/ucsc-fixseqliftoverpsl-grch37-{ DV .ucsc_fix_seq_liftover_37 } /fixSeqLiftOverPsl.bed.gz" ,
179
+ f"output/worker/annos/features/ucsc-fixseqliftoverpsl-grch38-{ DV .ucsc_fix_seq_liftover_38 } /fixSeqLiftOverPsl.bed.gz" ,
180
+ # ----- tads
181
+ "output/worker/annos/strucvars/tads-grch37-dixon2015/hesc.bed" ,
182
+ "output/worker/annos/strucvars/tads-grch38-dixon2015/hesc.bed" ,
168
183
169
184
170
185
# ===============================================================================================
@@ -173,11 +188,14 @@ rule all:
173
188
174
189
175
190
# -- work directory -----------------------------------------------------------------------------
191
+ # Misc rules.
192
+ include : "rules/work/misc/hpo.smk"
176
193
# Gene-related rules.
177
194
include : "rules/work/genes/dbnsfp.smk"
178
195
include : "rules/work/genes/ensembl.smk"
179
196
include : "rules/work/genes/gnomad.smk"
180
197
include : "rules/work/genes/hgnc.smk"
198
+ include : "rules/work/genes/mehari_data_tx.smk"
181
199
include : "rules/work/genes/ncbi.smk"
182
200
# Reference sequence--related rules.
183
201
include : "rules/work/reference/human.smk"
@@ -202,13 +220,21 @@ include: "rules/work/annos/strucvars/exac.smk"
202
220
include : "rules/work/annos/strucvars/g1k.smk"
203
221
include : "rules/work/annos/strucvars/gnomad.smk"
204
222
# -- output directory ---------------------------------------------------------------------------
223
+ # ---- mehari
205
224
include : "rules/output/mehari/freqs.smk"
206
- include : "rules/output/worker/cadd.smk"
207
- include : "rules/output/worker/dbsnp.smk"
208
- include : "rules/output/worker/dbnsfp.smk"
209
- include : "rules/output/worker/dbscsnv.smk"
210
- include : "rules/output/worker/gnomad_mtdna.smk"
211
- include : "rules/output/worker/gnomad_exomes.smk"
212
- include : "rules/output/worker/gnomad_genomes.smk"
213
- include : "rules/output/worker/helix.smk"
214
- include : "rules/output/worker/cons.smk"
225
+ # ---- viguno
226
+ include : "rules/output/viguno/hpo.smk"
227
+ # ------ annonars
228
+ include : "rules/output/annonars/cadd.smk"
229
+ include : "rules/output/annonars/cons.smk"
230
+ include : "rules/output/annonars/dbnsfp.smk"
231
+ include : "rules/output/annonars/dbscsnv.smk"
232
+ include : "rules/output/annonars/dbsnp.smk"
233
+ include : "rules/output/annonars/gnomad_exomes.smk"
234
+ include : "rules/output/annonars/gnomad_genomes.smk"
235
+ include : "rules/output/annonars/gnomad_mtdna.smk"
236
+ include : "rules/output/annonars/helix.smk"
237
+ # ---- worker
238
+ # ------ global
239
+ include : "rules/output/worker/genes.smk"
240
+ include : "rules/output/worker/patho_mms.smk"
0 commit comments