A R/MATLAB package to perform virtual knockout experiments on single-cell gene regulatory networks. **scTenifoldKnk** is a machine learning workflow that performs virtual knockout experiments using single-cell RNA sequencing (scRNAseq) data from wild-type (WT) control samples as input. Constructs a single-cell gene regulatory network (scGRN) and knocks out a target gene from the adjacency matrix of the WT scGRN by setting the gene’s outdegree edges to zero. **scTenifoldKnk** then compares the knocked out scGRN with the WT scGRN to identify differentially regulated genes, called virtual-knockout perturbed genes, which are used to assess the impact of the gene knockout and reveal the gene’s function in the analyzed cells.

This package is under active development, you can install **scTenifoldKnk**, using the following command:

```
library(remotes)
install_github('cailab-tamu/scTenifoldKnk')
library(scTenifoldKnk)
```

Code | Function |
---|---|

scTenifoldKnk | Perform virtual knockout experiments on single-cell gene regulatory networks |

The required input for **scTenifoldKnk** is an expression matrix with genes in the rows and cells (barcodes) in the columns. Data is expected to be *not normalized*.

The running time of scTenifoldKnk is largely dependent on how long it takes to construct scGRNs from subsampled expression matrices. Time increases proportional to the number of cells and genes in the dataset used as input. Below is a table of running times under different scenarios:

Number of Cells | Number of Genes | Running Time |
---|---|---|

300 | 1000 | 3.45 min |

1000 | 1000 | 4.25 min |

1000 | 5000 | 171.88 min (2 h 51.6 min) |

2500 | 5000 | 175.29 min (2 h 55.3 min) |

5000 | 5000 | 188.88 min (3 h 8.9 min) |

5000 | 7500 | 189.51 min (3 h 9.5 min) |

7500 | 5000 | 615.45 min (10 h 15.5 min) |

7500 | 7500 | 616.12 min (10 h 16.1 min) |

The output of **scTenifoldKnk** is a list with 3 slots as follows: * **tensorNetworks**: The computed weight-averaged denoised gene regulatory networks after CANDECOMP/PARAFAC (CP) tensor decomposition. It includes two slots with: * **X**: The constructed network for the *X* sample. * **Y**: The constructed network for the *Y* sample. * **manifoldAlignment**: The generated low-dimensional features result of the non-linear manifold alignment. It is a data frame with *2 times the number of genes* in the rows and *d* (default= 2) dimensions in the columns * **diffRegulation**: The results of the differential regulation analysis. It is a data frame with 6 columns as follows: * **gene**: A character vector with the gene id identified from the manifoldAlignment output. * **distance**: A numeric vector of the Euclidean distance computed between the coordinates of the same gene in both conditions. * **Z**: A numeric vector of the Z-scores computed after Box-Cox power transformation. * **FC**: A numeric vector of the FC computed with respect to the expectation. * **p.value**: A numeric vector of the p-values associated to the fold-changes, probabilities are asigned as P[X > x] using the Chi-square distribution with one degree of freedom. * **p.adj**: A numeric vector of adjusted p-values using Benjamini & Hochberg (1995) FDR correction.