Get Single Nucleux Gex

Retrieve Single Nucleus Gene Expression Data for a given Gene.

GTEx Portal API documentation

Usage

get_single_nucleus_gex(
  gencodeIds,
  datasetId = "gtex_snrnaseq_pilot",
  tissueSiteDetailIds = NULL,
  excludeDataArray = TRUE,
  page = 0,
  itemsPerPage = 250
)

Arguments

gencodeIds

A character vector of Versioned GENCODE IDs, e.g. c("ENSG00000132693.12", "ENSG00000203782.5").

datasetId

String. Unique identifier of a dataset. Usually includes a data source and data release. Options: "gtex_v8", "gtex_snrnaseq_pilot".

tissueSiteDetailIds

Character vector of IDs for tissues of interest. Can be GTEx specific IDs (e.g. "Whole_Blood"; use get_tissue_site_detail() to see valid values) or Ontology IDs.

excludeDataArray

String. Options are TRUE or FALSE

page

Integer (default = 0).

itemsPerPage

Integer (default = 250).

Value

A tibble

See also

Other Expression Data Endpoints: get_clustered_median_exon_expression(), get_clustered_median_gene_expression(), get_clustered_median_junction_expression(), get_clustered_median_transcript_expression(), get_expression_pca(), get_gene_expression(), get_median_exon_expression(), get_median_gene_expression(), get_median_junction_expression(), get_median_transcript_expression(), get_single_nucleus_gex_summary(), get_top_expressed_genes()

Examples

# \dontrun{
# Search for one or more genes - returns a tibble with one row per tissue.
# Column "cellTypes" now contains a tibble of expression summary data, with
# one row for each cell type
get_single_nucleus_gex(gencodeIds = c("ENSG00000203782.5",
                                      "ENSG00000132693.12"))
#> 
#> ── Paging info ─────────────────────────────────────────────────────────────────
#> • numberOfPages = 1
#> • page = 0
#> • maxItemsPerPage = 250
#> • totalNumberOfItems = 14
#> # A tibble: 14 × 7
#>    tissueSiteDetailId  ontologyId datasetId gencodeId geneSymbol cellTypes unit 
#>    <chr>               <chr>      <chr>     <chr>     <chr>      <list>    <chr>
#>  1 Muscle_Skeletal     UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  2 Breast_Mammary_Tis… UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  3 Esophagus_Mucosa    UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  4 Esophagus_Muscular… UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  5 Lung                UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  6 Prostate            UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  7 Skin_Sun_Exposed_L… UBERON:00… gtex_snr… ENSG0000… LOR        <tibble>  log(…
#>  8 Muscle_Skeletal     UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#>  9 Breast_Mammary_Tis… UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 10 Esophagus_Mucosa    UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 11 Esophagus_Muscular… UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 12 Heart_Left_Ventric… UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 13 Lung                UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 14 Prostate            UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…

# `excludeDataArray = FALSE` - expression values are stored under "celltypes"
# in an additional column called "data"
response <- get_single_nucleus_gex(gencodeIds = "ENSG00000132693.12",
                                   excludeDataArray = FALSE,
                                   itemsPerPage = 2)
#> Warning: ! Total number of items (7) exceeds maximum page size (2).
#> ℹ Try increasing `itemsPerPage`.
#> 
#> ── Paging info ─────────────────────────────────────────────────────────────────
#> • numberOfPages = 4
#> • page = 0
#> • maxItemsPerPage = 2
#> • totalNumberOfItems = 7

response
#> # A tibble: 2 × 7
#>   tissueSiteDetailId   ontologyId datasetId gencodeId geneSymbol cellTypes unit 
#>   <chr>                <chr>      <chr>     <chr>     <chr>      <list>    <chr>
#> 1 Muscle_Skeletal      UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…
#> 2 Breast_Mammary_Tiss… UBERON:00… gtex_snr… ENSG0000… CRP        <tibble>  log(…

# "cellTypes" contains a tibble of data with one row for each
# cell type e.g. for Breast_Mammary_Tissue
response$cellTypes[[2]]
#> # A tibble: 2 × 8
#>   cellType                  count meanWithZeros meanWithoutZeros medianWithZeros
#>   <chr>                     <int>         <dbl>            <dbl>           <dbl>
#> 1 Epithelial cell (luminal)     2      0.000902             2.19               0
#> 2 Endothelial cell (vascul…     1      0.00401              3.05               0
#> # ℹ 3 more variables: medianWithoutZeros <dbl>, numZeros <int>, data <list>

# when `excludeDataArray = FALSE`, expression values are stored in "data"
# e.g. for Breast_Mammary_Tissue, Epithelial cell (luminal):
response$cellTypes[[2]]$data[[1]]
#> [1] 1.593315 2.782725
# }