Inference_of_Gene_regulatory_networks

The aim of this project was to infer a gene regulatory network from single cell data using Copulas to find differentially coexpressed genes.
The data used a example data consists of 3k PBMCs from a Healthy Donor and is freely available from 10x Genomics. It can be downloaded from this web site or using

wget http://cf.10xgenomics.com/samples/cell-exp/1.1.0/pbmc3k/pbmc3k_filtered_gene_bc_matrices.tar.gz -O data/pbmc3k_filtered_gene_bc_matrices.tar.gz

on linux.

Method

After doing the preprocessing on should end up with 2 clusters of expression data. In order to infer a gene regulatory network from the single cell data a copula is computed for a pair of transcription factor and gene. The distributions of each genes expression data are used as the marginal distributions. This process is done twice, once in each cluster. Afterwards the KS-Distance between the 2 Copulas is computed.

Documentation

Script name	functionality
Preprocessing_singlecell.ipynb	Does the preprocessing of the single cell data as described in the scanpy tutorial
Extract_clusters.ipynb	Extracts clusters from the preprocessed data
Copula_Networkinference.py	Computes KS-Distances between a Tf and a gene Can generate random KS-Distances
Merge_bootstrapps.py	Merges multiple Bootstrapps into one file and computes the mean squared error
Filtering.py	Filters the result based on error and assigns p-values to the KS-Distances (if a file with random KS-Distances is provided)
shufflegeneexpressiondata.py	shuffles the entire data shuffles the entrie data except for the transcription factors
PlottingFunctions.py	plot all KS-Distances for edges starting from a Tf plot the gene expression data of a gene plot the distribution of errors plot the distribution of p-values plot the KS-Distances for edges starting from a Tf against the same edges from a differnt file plot the copula for a Tf/gene pair in each cluster
PlottingVenn.py	Compare the edges most likely to be in the network against a gold standard to validate the result

benjamin-kroeger / inference_of_gene_regulatory_networks Goto Github PK