|
a |
|
b/settings.py |
|
|
1 |
import io |
|
|
2 |
import numpy as np |
|
|
3 |
import anndata as anndata |
|
|
4 |
import scipy as s |
|
|
5 |
import scanpy as sc |
|
|
6 |
import matplotlib.pyplot as plt |
|
|
7 |
import seaborn |
|
|
8 |
import pandas as pd |
|
|
9 |
import csv |
|
|
10 |
import os |
|
|
11 |
from dfply import * |
|
|
12 |
# from IPython.display import display |
|
|
13 |
from re import search |
|
|
14 |
from pathlib import Path |
|
|
15 |
import argparse |
|
|
16 |
|
|
|
17 |
######### |
|
|
18 |
## I/O ## |
|
|
19 |
######### |
|
|
20 |
|
|
|
21 |
io = {} |
|
|
22 |
|
|
|
23 |
host = os.uname()[1] |
|
|
24 |
if search("ricard", host): |
|
|
25 |
io["basedir"] = "/Users/ricard/data/gastrulation_multiome_10x/test" |
|
|
26 |
io["cellranger_output"] = "/Users/ricard/data/gastrulation_multiome_10x/original" |
|
|
27 |
io["gene_metadata"] = "/Users/ricard/data/ensembl/mouse/v87/BioMart/all_genes/Mmusculus_genes_BioMart.87.txt" |
|
|
28 |
io["gene_metadata_gtf"] = "/Users/ricard/data/ensembl/mouse/v98/GTF/Mus_musculus.GRCm38.98.gtf.gz" |
|
|
29 |
elif search("Workstation", host): |
|
|
30 |
io["basedir"]="/home/lijingyu/gastrulation/data/gastrulation_multiome_10x" |
|
|
31 |
io["gene_metadata"]='/home/lijingyu/gastrulation/data/ensembl/mouse/v87/BioMart/mRNA/Mmusculus_genes_BioMart.87.txt' |
|
|
32 |
elif search("BI2404M", host): |
|
|
33 |
io["basedir"] = "/Users/argelagr/data/gastrulation_multiome_10x/test" |
|
|
34 |
io["gene_metadata"] = "/Users/argelagr/data/ensembl/mouse/v87/BioMart/all_genes/Mmusculus_genes_BioMart.87.txt" |
|
|
35 |
elif search("pebble|headstone", host): |
|
|
36 |
io["basedir"] = "/bi/group/reik/ricard/data/gastrulation_multiome_10x/test" |
|
|
37 |
io["cellranger_output"] = "/bi/group/reik/ricard/data/gastrulation_multiome_10x/original" |
|
|
38 |
io["gene_metadata"] = "/bi/group/reik/ricard/data/ensembl/mouse/v87/BioMart/all_genes/Mmusculus_genes_BioMart.87.txt" |
|
|
39 |
io["gene_metadata_gtf"] = "/bi/group/reik/ricard/data/ensembl/mouse/v98/GTF/Mus_musculus.GRCm38.98.gtf.gz" |
|
|
40 |
|
|
|
41 |
else: |
|
|
42 |
print("Computer not recognised"); exit() |
|
|
43 |
|
|
|
44 |
io["metadata"] = io["basedir"] + "/sample_metadata.txt.gz" |
|
|
45 |
io["anndata"] = io["basedir"] + "/processed/rna/anndata.h5ad" |
|
|
46 |
io["anndata_scvelo"] = io["basedir"] + "/processed/rna/velocyto/anndata_scvelo.h5ad" |
|
|
47 |
|
|
|
48 |
############# |
|
|
49 |
## Options ## |
|
|
50 |
############# |
|
|
51 |
|
|
|
52 |
opts = {} |
|
|
53 |
|
|
|
54 |
opts["samples"] = [ |
|
|
55 |
"E7.5_rep1", |
|
|
56 |
"E7.5_rep2", |
|
|
57 |
"E7.75_rep1", |
|
|
58 |
"E8.0_rep1", |
|
|
59 |
"E8.0_rep2", |
|
|
60 |
"E8.5_rep1", |
|
|
61 |
"E8.5_rep2", |
|
|
62 |
"E8.75_rep1", |
|
|
63 |
"E8.75_rep2", |
|
|
64 |
"E8.5_CRISPR_T_KO", |
|
|
65 |
"E8.5_CRISPR_T_WT" |
|
|
66 |
] |
|
|
67 |
|
|
|
68 |
|
|
|
69 |
opts["stages"] = [ |
|
|
70 |
"E7.5", |
|
|
71 |
"E7.75", |
|
|
72 |
"E8.0", |
|
|
73 |
"E8.5", |
|
|
74 |
"E8.75" |
|
|
75 |
] |
|
|
76 |
|
|
|
77 |
opts["celltypes"] = [ |
|
|
78 |
"Epiblast", |
|
|
79 |
"Primitive_Streak", |
|
|
80 |
"Caudal_epiblast", |
|
|
81 |
"PGC", |
|
|
82 |
"Anterior_Primitive_Streak", |
|
|
83 |
"Notochord", |
|
|
84 |
"Def._endoderm", |
|
|
85 |
"Gut", |
|
|
86 |
"Nascent_mesoderm", |
|
|
87 |
"Mixed_mesoderm", |
|
|
88 |
"Intermediate_mesoderm", |
|
|
89 |
"Caudal_Mesoderm", |
|
|
90 |
"Paraxial_mesoderm", |
|
|
91 |
"Somitic_mesoderm", |
|
|
92 |
"Pharyngeal_mesoderm", |
|
|
93 |
"Cardiomyocytes", |
|
|
94 |
"Allantois", |
|
|
95 |
"ExE_mesoderm", |
|
|
96 |
"Mesenchyme", |
|
|
97 |
"Haematoendothelial_progenitors", |
|
|
98 |
"Endothelium", |
|
|
99 |
"Blood_progenitors_1", |
|
|
100 |
"Blood_progenitors_2", |
|
|
101 |
"Erythroid1", |
|
|
102 |
"Erythroid2", |
|
|
103 |
"Erythroid3", |
|
|
104 |
"NMP", |
|
|
105 |
"Rostral_neurectoderm", |
|
|
106 |
"Caudal_neurectoderm", |
|
|
107 |
"Neural_crest", |
|
|
108 |
"Forebrain_Midbrain_Hindbrain", |
|
|
109 |
"Spinal_cord", |
|
|
110 |
"Surface_ectoderm", |
|
|
111 |
"Visceral_endoderm", |
|
|
112 |
"ExE_endoderm", |
|
|
113 |
"ExE_ectoderm", |
|
|
114 |
"Parietal_endoderm" |
|
|
115 |
] |
|
|
116 |
|
|
|
117 |
opts["celltype_colors"] = { |
|
|
118 |
"Epiblast" : "#635547", |
|
|
119 |
"Primitive_Streak" : "#DABE99", |
|
|
120 |
"Caudal_epiblast" : "#9e6762", |
|
|
121 |
"PGC" : "#FACB12", |
|
|
122 |
"Anterior_Primitive_Streak" : "#c19f70", |
|
|
123 |
"Notochord" : "#0F4A9C", |
|
|
124 |
"Def._endoderm" : "#F397C0", |
|
|
125 |
"Gut" : "#EF5A9D", |
|
|
126 |
"Nascent_mesoderm" : "#C594BF", |
|
|
127 |
"Mixed_mesoderm" : "#DFCDE4", |
|
|
128 |
"Intermediate_mesoderm" : "#139992", |
|
|
129 |
"Caudal_Mesoderm" : "#3F84AA", |
|
|
130 |
"Paraxial_mesoderm" : "#8DB5CE", |
|
|
131 |
"Somitic_mesoderm" : "#005579", |
|
|
132 |
"Pharyngeal_mesoderm" : "#C9EBFB", |
|
|
133 |
"Cardiomyocytes" : "#B51D8D", |
|
|
134 |
"Allantois" : "#532C8A", |
|
|
135 |
"ExE_mesoderm" : "#8870ad", |
|
|
136 |
"Mesenchyme" : "#cc7818", |
|
|
137 |
"Haematoendothelial_progenitors" : "#FBBE92", |
|
|
138 |
"Endothelium" : "#ff891c", |
|
|
139 |
"Blood_progenitors" : "#c9a997", |
|
|
140 |
"Blood_progenitors_1" : "#f9decf", |
|
|
141 |
"Blood_progenitors_2" : "#c9a997", |
|
|
142 |
"Erythroid" : "#EF4E22", |
|
|
143 |
"Erythroid1" : "#C72228", |
|
|
144 |
"Erythroid2" : "#f79083", |
|
|
145 |
"Erythroid3" : "#EF4E22", |
|
|
146 |
"NMP" : "#8EC792", |
|
|
147 |
"Neurectoderm" : "#65A83E", |
|
|
148 |
"Rostral_neurectoderm" : "#65A83E", |
|
|
149 |
"Caudal_neurectoderm" : "#354E23", |
|
|
150 |
"Neural_crest" : "#C3C388", |
|
|
151 |
"Forebrain_Midbrain_Hindbrain" : "#647a4f", |
|
|
152 |
"Spinal_cord" : "#CDE088", |
|
|
153 |
"Surface_ectoderm" : "#f7f79e", |
|
|
154 |
"Visceral_endoderm" : "#F6BFCB", |
|
|
155 |
"ExE_endoderm" : "#7F6874", |
|
|
156 |
"ExE_ectoderm" : "#989898", |
|
|
157 |
"Parietal_endoderm" : "#1A1A1A" |
|
|
158 |
} |
|
|
159 |
|
|
|
160 |
# opts["stages_colors"] = { |
|
|
161 |
# 'E6.5':"#D53E4F", |
|
|
162 |
# 'E6.75':"#F46D43", |
|
|
163 |
# 'E7.0':"#FDAE61", |
|
|
164 |
# 'E7.5':"#FFFFBF", |
|
|
165 |
# 'E7.25':"#FEE08B", |
|
|
166 |
# 'E7.75':"#E6F598", |
|
|
167 |
# 'E8.0':"#ABDDA4", |
|
|
168 |
# 'E8.5':"#3288BD", |
|
|
169 |
# 'E8.25':"#66C2A5", |
|
|
170 |
# 'mixed_gastrulation': "#A9A9A9" |
|
|
171 |
# } |
|
|
172 |
|
|
|
173 |
opts["stage_colors"] = { |
|
|
174 |
'E7.5':"#440154FF", |
|
|
175 |
'E7.75':"#E6F598", |
|
|
176 |
'E8.0':"#31688EFF", |
|
|
177 |
'E8.5':"#35B779FF", |
|
|
178 |
'E8.75':"#FDE725FF" |
|
|
179 |
} |