Quick start for perform cross-species analysis on gene regulatory networks

Jiang junyao

2024-05-13

introduction

This tutorial is made for users who want to quickly use CACIMAR to identify conserved gene regulatory networks from different modules. Input of this part is a table contains Source, SourceGroup, Target, TargetGroup.

library(CACIMAR)
load(system.file("extdata", "gene_network.rda", package = "CACIMAR"))
head(mm_gene_network)
##               Source SourceGroup             Target TargetGroup
## 1 ENSMUSG00000003949           1 ENSMUSG00000021250           1
## 2 ENSMUSG00000003949           1 ENSMUSG00000022707           1
## 3 ENSMUSG00000003949           1 ENSMUSG00000025892           1
## 4 ENSMUSG00000003949           1 ENSMUSG00000026989           1
## 5 ENSMUSG00000003949           1 ENSMUSG00000036523           1
## 6 ENSMUSG00000003949           1 ENSMUSG00000038248           1

Run the code

ConservedNetworks_ct <- Identify_ConservedNetworks(OrthG_Mm_Zf, mm_gene_network, zf_gene_network, 'mm', 'zf')
## Using NCS as value column: use value.var to override.
Heatmap_Cor(ConservedNetworks_ct[[2]], cluster_cols=F, cluster_rows=F, Color1 = c(rgb(102/255,46/255,115/255), rgb(31/255,153/255,139/255), rgb(251/255,232/255,48/255)))