Convert Seurat objects to 10x Genomics Loupe files.
● How To Use ● Installation ● Tutorials ● Loupe Browser Compatibility ● Troubleshooting
loupeR
creates 10x Genomics Loupe files from Seurat objects and other 10x Genomics data in R. 10x Genomics Loupe Browser can visualize single-cell and spatial data from 10x Genomics. Only single-cell gene expression datasets are supported in LoupeR.
Converting a Seurat object to a Loupe file is as simple as the following:
# import the library
library("loupeR")
# convert the SeuratObject named `seurat_obj` to a Loupe file
create_loupe_from_seurat(seurat_obj)
Use the function create_loupe
if you need more control in the clusters and projections that included in the Loupe file.
# import the library
library("loupeR")
# Gene Expression RNA assay
assay <- seurat_obj[["RNA"]]
# get counts matrix from either the old or newer formats of assay
counts <- counts_matrix_from_assay(assay)
# convert the count matrix, clusters, and projections into a Loupe file
create_loupe(
counts,
clusters = select_clusters(seurat_obj),
projections = select_projections(seurat_obj)
)
Additionally, use the utility function read_feature_ids_from_tsv
to read the Ensemble ids from the 10x dataset. A Seurat object will only have imported the feature names or ids and attached these as rownames to the count matrix. In order for the Ensemble id links to work correctly within Loupe Browser, one must manually import them and include them.
# import the library
library("loupeR")
# Gene Expression RNA assay
assay <- seurat_obj[["RNA"]]
# A character vector of ensemble ids.
feature_ids <- read_feature_ids_from_tsv("PATH_TO_10X_DATA/features.tsv.gz")
create_loupe_from_seurat(seurat_obj, feature_ids = feature_ids)
Before using loupeR
, make sure that your system has installed HDF5. The HDF5 organization requires registration before being able to download the installer. Below are some other more convenient methods for installing HDF5 if you happen to have these package managers installed.
In RStudio, or your R shell, run the following to install this package.
# install dependencies
if (!require("hdf5r")) install.packages("hdf5r")
if (!require("Seurat")) install.packages("Seurat")
# install platform specific source package
os <- sub("Darwin", "macOS", Sys.info()["sysname"])
url <- paste0("https://github.com/10XGenomics/loupeR/releases/latest/download/loupeR_", os, ".tar.gz")
install.packages(url, repos = NULL, type = "source")
Another installation option is to use the remotes
package to directly install loupeR
and its dependencies. The installed package won't include the prebundled louper executable, so you must invoke the loupeR::setup()
function which will go and download it.
remotes::install_github("10XGenomics/loupeR")
loupeR::setup()
If you are interested in automating LoupeR installation and execution (and are blocked by interactive license acceptance), please write to support@10xgenomics.com for further assistance.
With new versions of the Loupe Browser, new version of LoupeR need to be released. The table below shows version requirements between the two.
LoupeR Version | Loupe Browser Version |
---|---|
v1.0.x | Loupe Browser >= 7.0 |
v1.1.1 | Loupe Browser >= 8.0 |
v1.1.2 | Loupe Browser >= 8.1 |
- Demo notebook with basic processing of an example 10x dataset
Barcodes must be from 10x Genomics exeriments to work with LoupeR. Valid 10x Genomics Single Cell Gene Expression barcodes have the characters ACGT repeated 16 times, followed by an optional GEM well suffix, for example:
AAACCCAAGAAATTGC
AAACCCAAGAAATTGC-1
Barcodes can also have an additional optional prefix or suffix. These optional prefixes and suffixes must be delineated either by a :
or a _
:
prefix_AAACCCAAGAAATTGC
AAACCCAAGAAATTGC_suffix
prefix_AAACCCAAGAAATTGC_suffix
prefix:AAACATACAAACAG
AAACATACAAACAG:suffix
prefix:AAACATACAAACAG:suffix
prefix_AAACCCAAGAAATTGC-1
AAACCCAAGAAATTGC-1_suffix
prefix_AAACCCAAGAAATTGC-1_suffix
prefix:AAACCCAAGAAATTGC-1
AAACCCAAGAAATTGC-1:suffix
prefix:AAACCCAAGAAATTGC-1:suffix
Note: Visium and Xenium barcodes are formatted differently. Visium and Xenium data are currently enabled for use with LoupeR, but not fully supported. Expression data for these assays can be processed by loupeR, but not image data.
See test-validate.R
for further examples of both valid and invalid barcode formatting, as well as validater.R
for the exact formatting requirements as code.
For more in depth documentation and support please head to our support page or send an email to support@10xgenomics.com
Additionally, we have provided utility functions to help gather useful information when contacting support or creating a Github issue.
# import the library
library("loupeR")
# print extra debug information
create_bugreport_from_seurat(seurat_obj)