Cell-cycle markers
Usage
CellCycleMarkers(object, geneToSymbol)
importCellCycleMarkers(file, organism, release, ignoreVersion = TRUE)
Arguments
- object
Object.
- geneToSymbol
GeneToSymbol
. Gene-to-symbol mappings. Must containgeneId
andgeneName
columns. SeeGeneToSymbol
for 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 intxIdVersion
andgeneIdVersion
, 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/l1/8y8sjzmn15v49jgrqglghcfr0000gn/T/RtmpJ7yUtj/temp_libpath2a0db37ca07/AcidSingleCell/extdata/markers/cell-cycle/drosophila-melanogaster.csv using base::`read.table()`.
#> → Making <GRanges> from Ensembl.
#> → Getting <EnsDb> from AnnotationHub 3.10.0 (2023-10-23).
#> ℹ "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()`.
#> ! Dropping 1 identifier without defined gene symbol.
print(object)
#> CellCycleMarkers
#> Drosophila melanogaster (Ensembl 101)
#> g2m(37): Bub1 BubR1 ... twe vih
#> s(29): Blm CDC45L ... spn-A tos