GitPlanet

Cheap and reliable Node.js hosting starts at $3/month, and $1/month static HTML hosting

Created with love in Canada, visit hostnodejs.com today

Feel like to post an Ad? Learn Details
All Projects → SUwonglab → PECA

SUwonglab / PECA

Licence: other
PECA is a software for inferring context specific gene regulatory network from paired gene expression and chromatin accessibility data

Programming Languages

matlab
3953 projects
shell
77523 projects

Labels

dnase-seqrna-seqgene-expressionatac-seqgene-regulatory-networkchromatin-accessibiity

Projects that are alternatives of or similar to PECA

GGR-cwl
CWL tools and workflows for GGR
Stars: ✭ 20 (-35.48%)
Mutual labels:  dnase-seq, rna-seq, atac-seq
Awesome Single Cell
Community-curated list of software packages and data resources for single-cell, including RNA-seq, ATAC-seq, etc.
Stars: ✭ 1,937 (+6148.39%)
Mutual labels:  gene-expression, atac-seq
GREIN
GREIN : GEO RNA-seq Experiments Interactive Navigator
Stars: ✭ 40 (+29.03%)
Mutual labels:  rna-seq, gene-expression
DGCA
Differential Gene Correlation Analysis
Stars: ✭ 32 (+3.23%)
Mutual labels:  rna-seq, gene-expression
scATAC-pro
A comprehensive tool for processing, analyzing and visulizing single cell chromatin accessibility sequencing data
Stars: ✭ 63 (+103.23%)
Mutual labels:  atac-seq, chromatin-accessibiity
dolphinnext
A graphical user interface for distributed data processing of high throughput genomics
Stars: ✭ 92 (+196.77%)
Mutual labels:  rna-seq, atac-seq
Protocols-4pub
Multi-omics analysis protocols by Lyu.
Stars: ✭ 37 (+19.35%)
Mutual labels:  rna-seq, atac-seq
reg-gen
Regulatory Genomics Toolbox: Python library and set of tools for the integrative analysis of high throughput regulatory genomics data.
Stars: ✭ 64 (+106.45%)
Mutual labels:  dnase-seq, atac-seq
chromVAR
chromatin Variability Across Regions (of the genome!)
Stars: ✭ 116 (+274.19%)
Mutual labels:  dnase-seq, atac-seq
SCope
Fast visualization tool for large-scale and high dimensional single-cell data
Stars: ✭ 62 (+100%)
Mutual labels:  gene-expression, gene-regulatory-network
haystack bio
Haystack: Epigenetic Variability and Transcription Factor Motifs Analysis Pipeline
Stars: ✭ 42 (+35.48%)
Mutual labels:  rna-seq, gene-expression
atacr
Analysing Capture Seq Count Data
Stars: ✭ 14 (-54.84%)
Mutual labels:  atac-seq
dorothea
R package to access DoRothEA's regulons
Stars: ✭ 98 (+216.13%)
Mutual labels:  gene-expression
kana
Single cell analysis in the browser
Stars: ✭ 81 (+161.29%)
Mutual labels:  rna-seq
grape-nf
An automated RNA-seq pipeline using Nextflow
Stars: ✭ 30 (-3.23%)
Mutual labels:  rna-seq
netZooM
netZooM is a network biology package implemented in MATLAB.
Stars: ✭ 25 (-19.35%)
Mutual labels:  gene-regulatory-network
ALPS
AnaLysis routines for ePigenomicS data - 🏫 Bioconductor project
Stars: ✭ 13 (-58.06%)
Mutual labels:  atac-seq
TCC-GUI
📊 Graphical User Interface for TCC package
Stars: ✭ 35 (+12.9%)
Mutual labels:  rna-seq
pipeline-pinfish-analysis
Pipeline for annotating genomes using long read transcriptomics data with pinfish
Stars: ✭ 27 (-12.9%)
Mutual labels:  rna-seq
GNE
This repository contains the tensorflow implementation of "GNE: A deep learning framework for gene network inference by aggregating biological information"
Stars: ✭ 27 (-12.9%)
Mutual labels:  gene-expression
View All Similar Projects ➔

PECA

Introduction:

PECA is a software for inferring context specific gene regulatory network from paired gene expression and chromatin accessibility data. please cite PECA and PECA2 papers:

Duren, Zhana, et al. "Modeling gene regulation from paired expression and chromatin accessibility data." Proceedings of the National Academy of Sciences 114.25 (2017): E4914-E4923.

Duren, Zhana, et al. "Time course regulatory analysis based on paired expression and chromatin accessibility data." Genome research 30.4 (2020): 622-634.

Quickly start:

wget https://github.com/SUwonglab/PECA/archive/master.zip

unzip master.zip

cd PECA-master/

bash install.sh

bash PECA.sh sampleName genome

Install:

bash install.sh

Run PECA:

Run PECA by following two steps:

Step 1: Input

Put the input files in folder named ./Input. Three files: ${SampleName}.txt, ${SampleName}.bam, ${SampleName}.bam.bai.

${SampleName}.txt is gene expression file containing two columns (tab delimited), gene Symbol and FPKM (or TPM).

${SampleName}.bam is chromatin accessibility data, DNase-seq or ATAC-seq.

${SampleName}.bam.bai is the index file of bam file.

Note that all the three files should have same before-dot-file-name ${SampleName},only difference is after dot ".txt", ".bam" or ".bam.bai". Please see the example of RAd4 in the ./Input directory.

Step 2: Run

sh PECA.sh ${SampleName} ${genome}

Example: sh PECA.sh RAd4 mm9

To make sure the code run smoothly, please provide at least 64GB memory.

The results will be ./Results/${SampleName}/ . ${SampleName}_network.txt is the tissue specific network.

TFTG_score.txt is regulation strength for the all TF to TG. Each row represent one TF and each column represents one target gene. Higher value represents higher possibility of regulation.

CRB_pval.txt is the Chromatin regulators' (CR) binding site matrix, each column represent one CR, each row represent one region, the values are p-values.

Run PECA without ENCODE data information

PECA model uses prior information from ENCODE data. One can learn this prior information using their own data without using the ENCODE data if the number of paired samples are greater than 5.

sh PECA_withoutENCODE.sh FullPath_to_sampleNameFile ${genome}

Example: sh PECA_withoutENCODE.sh /home/user/sampleName.txt hg19 Here /home/user/sampleName.txt is a txt file that contain sample names (contain one sample name per line). For example

ES_day0

ES_day2

ES_day4

ES_day6

ES_day10

ES_day20

Under Input folder you should have ES_day0.txt, ES_day0.bam, and ES_day0.bam.bai, and the same for other samples. The reults of ES_day0 will be stored in ./Results__withoutENCODE/ES_day0/.

Run PECA_net_dif:

If you have two samples and want to compare the two samples at network level, please do it by following steps:

1, Prepare two networks: Run PECA on two samples one by one by "sh PECA.sh ${sampleName} ${genome}"

2, Run: sh PECA_compare_dif.sh ${Sample1} ${Sample2} ${Organism}

Example: sh PECA_compare_dif.sh K562 GM12878 human ; sh PECA_compare_dif.sh mESC RAd4 mouse

The results will be ./Results/Compare_${Sample1}_${Sample2}. Containing six files:

specific network of two samples: ${Sample1}_specific_network.txt and ${Sample2}_specific_network.txt

common network of two samples: ${Sample1}_${Sample2}_common_network.txt

specific module of two networks: ${Sample1}_specific_module.txt and ${Sample2}_specific_module.txt

common module of two samples: ${Sample1}_${Sample2}_common_module.txt

Files PooledNetwork.txt or PooledModule.txt can be used to visualize the network by cytoscype, and the node lable is given in file Node_lable.txt. "1" and "-1" in PooledNetwork.txt or PooledModuole.txt represent "Activation" and "Repression" respectively. "1" and "2" in Node_lable.txt represent the gene is Sample1 specific or Sample2 specific.

Run PECA_net_dif_multiple:

If you have two conditions (multiple samples in each conditions) and want to compare the two conditions at network level, please do it by following steps:

1, Prepare networks: Run PECA on all the samples from two conditions one by one by "sh PECA.sh ${sampleName} ${genome}"

2, Construct lables: Write the sample names of Group1 and Group2 into text files named $Group1 and $Group2, respectively. (eg. create one text file named "Control" and put the sample names of one condition to this file, create other text file named "Case" and put the names of the other condition to this file. Note that the sample name files contain one sample name per line )

3, Run: sh PECA_compare_dif_multiple.sh $Group1 $Group2 ${Organism} Example： sh PECA_compare_dif_multiple.sh Control Case human

The results will be ./Results/CompareGroup_${Group1}_${Group2}. Containing six files:

specific network of two conditions: ${Group1}_specific_network.txt and ${Group2}_specific_network.txt

common network of two conditions: ${Group1}_${Group2}_common_network.txt

specific module of two conditions: ${Group1}_specific_module.txt and ${Group2}_specific_module.txt

common module of two conditions: ${Group1}_${Group2}_common_module.txt

Files PooledNetwork.txt or PooledModuole.txt can be used to visualize the network by cytoscype, and the node lable is given in file Node_lable.txt. "1" and "-1" in PooledNetwork.txt or PooledModuole.txt represent "Activation" and "Repression" respectively. "1" and "2" in Node_lable.txt represent the gene is Group1 specific or Group2 specific.

Requirements:

  • Matlab (Optimization Toolbox)

  • macs2

  • homer

  • samtools

  • bedtools

Contact:

If you have any issues, please contact Zhana Duren by [email protected]

Note that the project description data, including the texts, logos, images, and/or trademarks, for each open source project belongs to its rightful owner. If you wish to add or remove any projects, please contact us at [email protected].