Skip to contents

Cell-cycle markers

Usage

CellCycleMarkers(object, geneToSymbol)

importCellCycleMarkers(file, organism, release, ignoreVersion = TRUE)

Arguments

object

Object.

geneToSymbol

GeneToSymbol. Gene-to-symbol mappings. Must contain geneId and geneName columns. See GeneToSymbol 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 in txIdVersion and geneIdVersion, and the variants without versions will be stored in txId, txIdNoVersion, geneId, and geneIdNoVersion.

Value

CellCycleMarkers.

Note

Updated 2024-03-27.

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