Cell-cycle markers
Usage
CellCycleMarkers(object, geneToSymbol)
importCellCycleMarkers(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"
)
cellCycleDir <- file.path(markersDir, "cell-cycle")
files <- sort(list.files(
path = cellCycleDir,
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 <- importCellCycleMarkers(
file = file,
organism = organism,
release = release
)
#> → Importing /private/var/folders/9b/4gh0pghx1b71jjd0wjh5mj880000gn/T/RtmpKjDQxC/temp_libpathb04a57984e76/AcidSingleCell/extdata/markers/cell-cycle/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)
#> CellCycleMarkers
#> Drosophila melanogaster (Ensembl 101)
#> g2m(37): Bub1 BubR1 ... twe vih
#> s(29): Blm CDC45L ... spn-A tos