edroaldo / cellrouter

Reconstruction of complex single-cell trajectories using CellRouter
45 stars 21 forks source link

regulate module #12

Open JianmeiZhong opened 6 years ago

JianmeiZhong commented 6 years ago

Hi, sorry that i still have several questions to bother you.=.= 1.in the part of regulators network,labels of some genes are printed as picture below shows while other equal-sized ,equal-colored don't.what's meaning when a gene is labeled? regulator_network

2.i'm a little confused about the definition of derivtive shown in heatmap of negtive and positive control.what's the difference and relationship with expression?

  1. your article mentioned that In the neutrophil trajectory, Mxd1 is upregulated in later stages of differentiation compared to Cebpe (Fig. 2c,d). According to GRN analyses, Mxd1 is predicted to regulate Mmp8, Mmp9, Retnlg, and Cd52 expression, all of which are known neutrophil markers (Fig. 2e, f).How to infer the predictive relationship of Mxd1 from the GRN results?

Thanks for your patience and your great work sincerely.

edroaldo commented 6 years ago

Hi, No problem.

  1. Only top "transcriptional regulators" are shown. The ones with similar size/color are then predicted targets of the "trascriptional regulators".

  2. The derivative captures the changes that happen in the expression kinects along the trajectory. Is is the same concept of the derivative that people usually have in math/calculus classes. So, the derivative should be the highest or lowest where the changes are happening in the kinetic profile of a given gene. Let me know if that clarifies.

  3. Ater you create the gene regulatory network with functions available in CellRouter, and after you calculate the GRN scores, using, for example: transitions <- c('6.1','6.2', '6.7', '6.8') #transitions to be analyzed grn.scores <- grnscores(cellrouter, ggrn = grn.data$GRN, tfs = grn.data$tfs, transitions = transitions, direction='both', flip=TRUE, q.up=0.75, q.down=0.25, dir.targets='up', columns=2, width=5, height=7, filename='results/GRN_scores')

You can check the targets of the top regulators by doing (using transition from cluster 6 to 1, as an example): grn.scores$6.1$targets$Mxd1. Let me know if that clarifies. I am working on some visualization functions for this, which will be available in a next release of CellRouter.

No problem! Thanks for your using CellRouter =D!

2018-05-17 5:33 GMT-04:00 JianmeiZhong notifications@github.com:

Hi, sorry that i still have several questions to bother you.=.= 1.in the part of regulators network,labels of some genes are printed as picture below shows while other equal-sized ,equal-colored don't.what's meaning when a gene is labeled? [image: regulator_network] https://user-images.githubusercontent.com/37535747/40165700-3eca236e-59ef-11e8-9bdd-d5cc20634239.png

2.i'm a little confused about the definition of derivtive shown in heatmap of negtive and positive control.what's the difference and relationship with expression?

  1. your article mentioned that In the neutrophil trajectory, Mxd1 is upregulated in later stages of differentiation compared to Cebpe (Fig. 2c,d). According to GRN analyses, Mxd1 is predicted to regulate Mmp8, Mmp9, Retnlg, and Cd52 expression, all of which are known neutrophil markers (Fig. 2e, f).How to infer the predictive relationship of Mxd1 from the GRN results?

Thanks for your patience and your great work sincerely.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/edroaldo/cellrouter/issues/12, or mute the thread https://github.com/notifications/unsubscribe-auth/AJqUR3n0tev1c-lVrNCBx6OiImJzdVy1ks5tzUP5gaJpZM4UCuW0 .

-- Edroaldo

JianmeiZhong commented 6 years ago

OK,I see them all.Thanks for your detailed explanation~

edroaldo commented 6 years ago

Great! Thanks! Please, let me know if you have further questions or suggestions to improve.

Thanks!

2018-05-26 7:27 GMT-04:00 JianmeiZhong notifications@github.com:

OK,I see them all.Thanks for your detailed explanation~

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/edroaldo/cellrouter/issues/12#issuecomment-392255186, or mute the thread https://github.com/notifications/unsubscribe-auth/AJqURx6AmPv5EKf1XGYtjIcKq6afPqABks5t2TwIgaJpZM4UCuW0 .

-- Edroaldo