Cell-type markers
Usage
CellTypeMarkers(object, geneToSymbol)
importCellTypeMarkers(file, organism, release, ignoreVersion = TRUE)Arguments
- object
Object.
- geneToSymbol
GeneToSymbol. Gene-to-symbol mappings. Must containgeneIdandgeneNamecolumns. SeeGeneToSymbolfor more information.- file
character(1). File path.- organism
character(1). Full Latin organism name (e.g."Homo sapiens").- release
integer(1). Ensembl release version (e.g.100). We recommend setting this value if possible, for improved reproducibility. When left unset, the latest release available via AnnotationHub/ensembldb is used. Note that the latest version available can vary, depending on the versions of AnnotationHub and ensembldb in use.- ignoreVersion
logical(1). Ignore identifier (e.g. transcript, gene) versions. When applicable, the identifier containing version numbers will be stored intxIdVersionandgeneIdVersion, and the variants without versions will be stored intxId,txIdNoVersion,geneId, andgeneIdNoVersion.
Examples
markersDir <- system.file(
"extdata", "markers",
package = "AcidSingleCell"
)
cellTypeDir <- file.path(markersDir, "cell-type")
files <- sort(list.files(
path = cellTypeDir,
pattern = "*.csv",
full.names = TRUE
))
file <- files[[1L]]
organism <- syntactic::sentenceCase(
gsub(
pattern = "-",
replacement = " ",
x = AcidBase::basenameSansExt(file)
)
)
releaseFile <- file.path(markersDir, "ensembl-release.txt")
release <- as.integer(readLines(releaseFile))
object <- importCellTypeMarkers(
file = file,
organism = organism,
release = release
)
#> → Importing /private/var/folders/9b/4gh0pghx1b71jjd0wjh5mj880000gn/T/RtmpKjDQxC/temp_libpathb04a57984e76/AcidSingleCell/extdata/markers/cell-type/drosophila-melanogaster.csv using base::`read.table()`.
#> → Making <GRanges> from Ensembl.
#> → Getting <EnsDb> from AnnotationHub 3.15.0 (2025-04-08).
#> ℹ "AH83185": Ensembl 101 EnsDb for Drosophila melanogaster.
#> → Making <GRanges> from <EnsDb>.
#> Organism: Drosophila melanogaster
#> Genome build: BDGP6.28
#> Release: 101
#> Level: genes
#> → Defining names by `geneId` column in `mcols()`.
print(object)
#> CellTypeMarkers
#> Drosophila melanogaster (Ensembl 101)
#> crystal_cell(5): PPO1 Sox102F ... lz mam
#> lamellocyte(16): 14-3-3zeta Eip75B ... pnt tx
#> plasmatocyte(18): Col4a1 He ... pnt ttk
#> prohemocyte(8): Tep4 dome ... srp wg