forked from sankaranlab/SCAVENGE
-
Notifications
You must be signed in to change notification settings - Fork 0
/
README.Rmd
75 lines (50 loc) · 4.57 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
---
title: "SCAVENGE"
author: "`r rworkflows::use_badges()`"
date: "<h4>README updated: <i>`r format( Sys.Date(), '%b-%d-%Y')`</i></h4>"
output:
github_document
---
<!-- To modify Package/Title/Description/Authors fields, edit the DESCRIPTION file -->
```{r, echo=FALSE, include=FALSE}
pkg <- read.dcf("DESCRIPTION", fields = "Package")[1]
title <- read.dcf("DESCRIPTION", fields = "Title")[1]
description <- read.dcf("DESCRIPTION", fields = "Description")[1]
URL <- read.dcf('DESCRIPTION', fields = 'URL')[1]
owner <- tolower(strsplit(URL,"/")[[1]][4])
```
## ``r pkg``: `r title`
### Overview:
Co-localization approaches using genetic variants and single-cell epigenomic data are unfortunately uninformative for many cells given the extensive sparsity across single-cell profiles. Therefore, only a few cells from the truly relevant population demonstrate reliable phenotypic relevance. The global high-dimensional features of individual single cells are sufficient to represent the underlying cell identities or states, which enables the relationships among such cells to be readily inferred. By taking advantage of these attributes, SCAVENGE identifies the most phenotypically-enriched cells by co-localization and explores the transitive associations across the cell-to-cell network to assign each cell a probability representing the cell’s relevance to those phenotype-enriched cells via network propagation.
We developed a novel enrichment method (**SCAVENGE**) (Single Cell Analysis of Variant Enrichment through Network propagation of GEnomic data) that can discriminate between closely related cell types/states and score single cells for GWAS enrichment.
<div align=center> <img src="image/schematic-view_1.png" width="680" height="278"> </div>
<p align="center">Schematic view of SCAVENGE</p>
We've implemented **SCAVENGE** as an `R` package for computing single-cell based GWAS enrichments from fine-mapped posterior probabilities and quantitative epigenomic data (i.e. scATAC-seq and potentially other single-cell epigenome profiling methods).
As single-cell genomic datasets grow in volume, we expect SCAVENGE will have great promise for efficiently uncovering relevant cell populations for more phenotypes or functions in different scenarios, which may expand beyond the complex trait genetic variants we have examined here. We welcome you to use SCAVENGE to discover more phenotype relevant cells!
### Installation:
Once all of the dependencies for `SCAVENGE` are installed, the package can be installed
directly from GitHub by typing the following into an `R` console:
```
if(!require("remotes")) install.packages("remotes")
remotes::install_github("`r URL`")
library(`r pkg`)
```
## Documentation
### Tutorial:
This web resource and vignette compiliation shows how to reproduce results of SCAVENGE analysis with monocyte count on a 10X PBMC dataset [[**Vignette-pdf**]](doc/SCAVENGE-vignette.pdf), [[**Vignette-R markdown code**]](doc/SCAVENGE-vignette.Rmd).
See the [[**Wiki page**]](https://github.com/sankaranlab/SCAVENGE/wiki) for some information on preparing your GWAS data for SCAVENGE, including finemapping.
### [Website](https://`r owner`.github.io/`r pkg`)
### [Getting started](https://`r owner`.github.io/`r pkg`/articles/`r pkg`)
### FAQs:
- *What input data are accepted for SCAVENGE analysis?*
A: The count matrix of scATAC-seq data and fine-mapped variants from GWAS summary statistics (we provided a tutorial for fine-mapping analysis from GWAS [[**Wiki page**]](https://github.com/sankaranlab/SCAVENGE/wiki)). Theoretically, GWAS summary statistics can be used as input but we do not recommend it because LD can obscure causal cell type identification.
- *Can I use scRNA-seq instead of scATAC-seq?*
A: It is not feasible for SCAVENGE analysis from scRNA-seq currently. We are actively developing this tool to be scalable to scRNA-seq, please stay tuned.
- *How can I request new feature?*
A: We open [[**Discussions**]](https://github.com/sankaranlab/SCAVENGE/discussions) page, please feel free to discuss and post your ideas.
### Citation:
If you used or adapted SCAVENGE in your study, please cite our paper [[*Nat Biotechnol*]](https://www.nature.com/articles/s41587-022-01341-y) || [[*PubMed*]](https://pubmed.ncbi.nlm.nih.gov/35668323/).
*Variant to function mapping at single-cell resolution through network propagation.*
### Contact:
If you run into issues and would like to report them, you can use the "Issues" tab on the left hand side.
Alternatively, you can contact authors: fyu{at}broadinstitute.org, lcato{at}broadinstitute.org and/or sankaran{at}broadinstitute.org.