Cancer-Testis (CT) genes, also called Cancer-Germline (CG), are a group of genes whose expression is normally restricted to the germline but that are found aberrantly activated in many types of cancers. These genes produce cancer specific antigens represent ideal targets for anti-cancer vaccines. Besides their interest in immunotherapy, they can also be used as cancer biomarkers and as target of anti-tumor therapies with limited side effect.
Many CT genes use DNA methylation as a primary mechanism of transcriptionnal regulation. This is another interesting point about CT genes as they represent suitable models to study DNA demethylation in cancer.
Currently the reference database gathering CT genes is the CTdatabase that was published in 2009, based on a literature screening(Almeida et al., 2009). This database is however not up to date. Recently identified CT genes are not referenced (in particular CT genes identified by omics methods that didn’t exist at the time) while some genes referred as CT genes appeared to be in reality expressed in many somatic tissues. Furthermore, the database is not in an easily importable format, some genes are not encoded properly (by synonyms names rather than by their official HGNC symbol names, or by a concatenation of both) resulting in poor interoperability for downstream analyses. More recent studies proposed other lists of CT genes like Wang’s CTatlas (Wang et al., 2016, Jamin et al., 2021, Carter et al., 2023). These lists were established using different criteria to define CT genes and hence differ substantially from each other. Moreover, these lists are usually provided as supplemental data files and are not strictly speaking databases. Finally, none of these studies describe the involvement of DNA methylation in the regulation of individual CT genes.
We therefore created CTexploreR
, a Bioconductor R
package, aimed to redefine rigorously the list of CT genes based on
publicly available RNAseq databases and to summarize their main
characteristics. We included methylation analyses to classify these
genes according to whether they are regulated or not by DNA methylation.
The package also offers tools to visualize CT genes expression and
promoter DNA methylation in normal and tumoral tissues. CTexploreR hence
represents an up-to-date reference database for CT genes and can be used
as a starting point for further investigations related to these
genes.
To install the package:
To install the package from GitHub:
The central element of CTexploreR
is the list of 280 CT
and CTP genes (see table below) selected based on their expression in
normal and tumoral tissues (selection details in the next section). The
table also summarises their main characteristics.
## Loading required package: CTdata
## # A tibble: 10 × 47
## ensembl_gene_id external_gene_name CT_gene_type testis_specificity
## <chr> <chr> <chr> <chr>
## 1 ENSG00000234593 KAZN-AS1 CT_gene testis_specific
## 2 ENSG00000131914 LIN28A CT_gene testis_specific
## 3 ENSG00000142698 C1orf94 CT_gene testis_specific
## 4 ENSG00000237853 NFIA-AS1 CT_gene testis_specific
## 5 ENSG00000137948 BRDT CT_gene testis_specific
## 6 ENSG00000198765 SYCP1 CT_gene testis_specific
## 7 ENSG00000237463 LRRC52-AS1 CT_gene testis_specific
## 8 ENSG00000215817 ZC3H11B CT_gene testis_specific
## 9 ENSG00000162843 WDR64 CT_gene testis_specific
## 10 ENSG00000231532 LINC01249 CT_gene testis_specific
## # ℹ 43 more variables: regulated_by_methylation <lgl>, X_linked <lgl>,
## # chr <chr>, strand <int>, transcript_start <int>, transcript_end <int>,
## # transcription_start_site <int>, GTEX_category <chr>, q75_TPM_somatic <dbl>,
## # max_TPM_somatic <dbl>, ratio_testis_somatic <dbl>, TPM_testis <dbl>,
## # lowly_expressed_in_GTEX <lgl>, multimapping_analysis <chr>,
## # HPA_scRNAseq_celltype_specific_nTPM <chr>, max_HPA_somatic <dbl>,
## # max_HPA_germcell <dbl>, not_detected_in_somatic_HPA <lgl>, …
CTdata
is the companion Package for
CTexploreR
and provides the omics data that was necessary
to select and characterize cancer testis genes as well as exploring
them. The data are served through the ExperimentHub
infrastructure. Currently available data are summarised in the table
below and details can be found in CTdata
vignette or
manuals.
## Title
## 1 GTEX_data
## 2 CCLE_data
## 3 normal_tissues_multimapping_data
## 4 DAC_treated_cells
## 5 DAC_treated_cells_multimapping
## 6 TCGA_TPM
## 7 CT_methylation_in_tissues
## 8 CT_mean_methylation_in_tissues
## 9 TCGA_CT_methylation
## 10 CT_genes
## 11 CCLE_correlation_matrix
## 12 testis_sce
## 13 scRNAseq_HPA
## 14 all_genes
## 15 FGC_sce
## 16 HPA_cell_type_specificities
## 17 mean_methylation_in_tissues
## 18 methylation_in_tissues
## 19 oocytes_sce
## 20 TCGA_methylation
## 21 embryo_sce_Zhu
## 22 embryo_sce_Petropoulos
## 23 hESC_data
## 24 methylation_in_hESC
## 25 mean_methylation_in_hESC
## 26 methylation_in_FGC
## 27 mean_methylation_in_FGC
## 28 methylation_in_embryo
## 29 mean_methylation_in_embryo
## Description
## 1 Gene expression data in normal tissues from GTEx database
## 2 Gene expression data in cancer cell lines from CCLE
## 3 Gene expression values in normal tissues with or without allowing multimapping
## 4 Gene expression values in a set of cell lines treated or not with 5-Aza-2'-Deoxycytidine (DAC), a demethylating agent
## 5 Gene expression values (multimapping allowed) in a set of cell lines treated or not with 5-Aza-2'-Deoxycytidine (DAC)
## 6 Gene expression data in TCGA samples
## 7 Methylation values of CpGs located within Cancer-Testis promoters in a set of normal tissues
## 8 Mean methylation values of all CpGs located within Cancer-Testis (CT) promoters in a set of normal tissues
## 9 Methylation of CT promoters in TCGA samples
## 10 Cancer-Testis (CT) genes description
## 11 Gene correlations in CCLE cancer cell lines
## 12 Testis scRNAseq data
## 13 Gene expression profiles in different human cell types based on scRNAseq data from The Human Protein Atlas
## 14 All genes description, according to the analysis done for CT genes
## 15 Fetal gonad scRNAseq
## 16 Cell type specificities (from HPA)
## 17 Mean methylation values of all CpGs located within all (CT) promoters in a set of normal tissues
## 18 Methylation values of CpGs located within all promoters in a set of normal tissues
## 19 Oocytes scRNAseq
## 20 Methylation of all genes promoters in TCGA samples
## 21 Gene expression values in early embryo (blastocyst)
## 22 Gene expression values in early embryo (blastocyst and Morula)
## 23 Gene expression data in human embryonic stem cells
## 24 Methylation values of CpGs located within all promoters in human embryonic stem cells
## 25 Mean methylation values of all CpGs located within all promoters in human embryonic stem cells
## 26 Methylation values of CpGs located within all promoters in fetal germ cells
## 27 Mean methylation values of all CpGs located within all promoters in fetal germ cells
## 28 Methylation values of CpGs located within all promoters in early embryo
## 29 Mean methylation values of all CpGs located within all promoters in early embryo
## BiocVersion Genome SourceType
## 1 3.19 hg38 tab
## 2 3.19 hg38 CSV
## 3 3.19 hg38 FASTQ
## 4 3.19 hg38 FASTQ
## 5 3.19 hg38 FASTQ
## 6 3.19 hg38 TSV
## 7 3.17 hg38 BED
## 8 3.17 hg38 BED
## 9 3.17 hg38 TXT
## 10 3.19 hg38 TSV
## 11 3.19 hg38 CSV
## 12 3.19 hg38 FASTQ
## 13 3.19 hg38 TSV
## 14 3.19 hg38 TSV
## 15 3.19 hg38 RDS
## 16 3.19 hg38 TSV
## 17 3.19 hg38 BED
## 18 3.19 hg38 BED
## 19 3.19 hg38 TXT
## 20 3.19 hg38 TSV
## 21 3.19 hg38 tab
## 22 3.19 hg38 TXT
## 23 3.19 hg38 FASTQ
## 24 3.19 hg38 BED
## 25 3.19 hg38 BED
## 26 3.19 hg19 TXT
## 27 3.19 hg19 TXT
## 28 3.19 hg19 BED
## 29 3.19 hg19 BED
## SourceUrl
## 1 https://storage.googleapis.com/gtex_analysis_v8/rna_seq_data/GTEx_Analysis_2017-06-05_v8_RNASeQCv1.1.9_gene_median_tpm.gct.gz
## 2 https://ndownloader.figshare.com/files/34989922, https://ndownloader.figshare.com/files/35020903
## 3 https://www.encodeproject.org/
## 4 https://www.encodeproject.org/
## 5 https://www.encodeproject.org/
## 6 https://portal.gdc.cancer.gov/
## 7 https://www.encodeproject.org/
## 8 https://www.encodeproject.org/
## 9 https://portal.gdc.cancer.gov/
## 10 https://zenodo.org/record/7537824/files/cancermine_collated.tsv?download=1
## 11 https://ndownloader.figshare.com/files/34989922, https://ndownloader.figshare.com/files/35020903
## 12 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE120508
## 13 https://www.proteinatlas.org/download/rna_single_cell_type.tsv.zip
## 14 https://zenodo.org/record/7537824/files/cancermine_collated.tsv?download=1
## 15 https://cellxgene.cziscience.com/collections/661a402a-2a5a-4c71-9b05-b346c57bc451Data
## 16 https://www.proteinatlas.org/
## 17 https://www.encodeproject.org/
## 18 https://www.encodeproject.org/
## 19 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE154762
## 20 https://portal.gdc.cancer.gov/
## 21 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81233
## 22 https://www.ebi.ac.uk/biostudies/files/E-MTAB-3929/
## 23 https://www.encodeproject.org/
## 24 https://www.encodeproject.org/
## 25 https://www.encodeproject.org/
## 26 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE107714
## 27 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE107714
## 28 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81233
## 29 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81233
## SourceVersion Species TaxonomyId Coordinate_1_based
## 1 8 Homo sapiens 9606 TRUE
## 2 22Q2 Homo sapiens 9606 TRUE
## 3 5 Homo sapiens 9606 TRUE
## 4 5 Homo sapiens 9606 TRUE
## 5 5 Homo sapiens 9606 TRUE
## 6 36.0 Homo sapiens 9606 TRUE
## 7 5 Homo sapiens 9606 TRUE
## 8 5 Homo sapiens 9606 TRUE
## 9 36.0 Homo sapiens 9606 TRUE
## 10 49 Homo sapiens 9606 TRUE
## 11 22Q2 Homo sapiens 9606 TRUE
## 12 1 Homo sapiens 9606 TRUE
## 13 1 Homo sapiens 9606 TRUE
## 14 49 Homo sapiens 9606 TRUE
## 15 1 Homo sapiens 9606 TRUE
## 16 23 Homo sapiens 9606 TRUE
## 17 5 Homo sapiens 9606 TRUE
## 18 5 Homo sapiens 9606 TRUE
## 19 1 Homo sapiens 9606 TRUE
## 20 36.0 Homo sapiens 9606 TRUE
## 21 1 Homo sapiens 9606 TRUE
## 22 1 Homo sapiens 9606 TRUE
## 23 5 Homo sapiens 9606 TRUE
## 24 5 Homo sapiens 9606 TRUE
## 25 5 Homo sapiens 9606 TRUE
## 26 1 Homo sapiens 9606 TRUE
## 27 1 Homo sapiens 9606 TRUE
## 28 1 Homo sapiens 9606 TRUE
## 29 1 Homo sapiens 9606 TRUE
## DataProvider
## 1 The Genotype-Tissue Expression
## 2 Cancer Cell Line Encyclopedia
## 3 ENCODE Project
## 4 ENCODE Project
## 5 ENCODE Project
## 6 The Cancer Genome Atlas
## 7 ENCODE Project
## 8 ENCODE Project
## 9 The Cancer Genome Atlas
## 10 CancerMine
## 11 Cancer Cell Line Encyclopedia
## 12 The adult human testis transcriptional cell atlas, Guo et al. 2018
## 13 The Human Protein Atlas
## 14 CancerMine
## 15 Single-cell roadmap of human gonadal development (Garcia-Alonso, Nature 2022)
## 16 The Human Protein Atlas
## 17 ENCODE Project
## 18 ENCODE Project
## 19 NCBI GEO : Decoding dynamic epigenetic landscapes in human oocytes\n using single-cell multi-omics sequencing (Yan et al. Cell Stem Cell 2021)
## 20 The Cancer Genome Atlas
## 21 Single Cell DNA Methylome Sequencing of Human Preimplantation Embryos (Zhu et al. 2018)
## 22 Single-Cell RNA-Seq Reveals Lineage and X Chromosome Dynamics in\n Human Preimplantation Embryos (Petropulous et al, 2014)
## 23 ENCODE Project
## 24 ENCODE Project
## 25 ENCODE Project
## 26 Data from Dissecting the epigenomic dynamics of human fetal germ cell\n development at single-cell resolution (Li et al. 2021)
## 27 Data from Dissecting the epigenomic dynamics of human fetal germ cell\n development at single-cell resolution (Li et al. 2021)
## 28 Single Cell DNA Methylome Sequencing of Human Preimplantation Embryos (Zhu et al. 2018)
## 29 Single Cell DNA Methylome Sequencing of Human Preimplantation Embryos (Zhu et al. 2018)
## Maintainer RDataClass
## 1 Axelle Loriot <[email protected]> SummarizedExperiment
## 2 Axelle Loriot <[email protected]> SummarizedExperiment
## 3 Axelle Loriot <[email protected]> SummarizedExperiment
## 4 Axelle Loriot <[email protected]> SummarizedExperiment
## 5 Axelle Loriot <[email protected]> SummarizedExperiment
## 6 Axelle Loriot <[email protected]> SummarizedExperiment
## 7 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 8 Axelle Loriot <[email protected]> SummarizedExperiment
## 9 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 10 Axelle Loriot <[email protected]> tibble
## 11 Axelle Loriot <[email protected]> matrix
## 12 Axelle Loriot <[email protected]> SingleCellExperiment
## 13 Axelle Loriot <[email protected]> SingleCellExperiment
## 14 Axelle Loriot <[email protected]> tibble
## 15 Axelle Loriot <[email protected]> SingleCellExperiment
## 16 Axelle Loriot <[email protected]> SingleCellExperiment
## 17 Axelle Loriot <[email protected]> SummarizedExperiment
## 18 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 19 Axelle Loriot <[email protected]> SingleCellExperiment
## 20 Axelle Loriot <[email protected]> SummarizedExperiment
## 21 Axelle Loriot <[email protected]> SingleCellExperiment
## 22 Axelle Loriot <[email protected]> SingleCellExperiment
## 23 Axelle Loriot <[email protected]> SummarizedExperiment
## 24 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 25 Axelle Loriot <[email protected]> SummarizedExperiment
## 26 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 27 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 28 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## 29 Axelle Loriot <[email protected]> RangedSummarizedExperiment
## DispatchClass RDataPath
## 1 Rda CTdata/eh_data/v2/GTEX_data.rda
## 2 Rda CTdata/eh_data/v2/CCLE_data.rda
## 3 Rda CTdata/eh_data/v2/normal_tissues_multimapping_data.rda
## 4 Rda CTdata/eh_data/v2/DAC_treated_cells.rda
## 5 Rda CTdata/eh_data/v2/DAC_treated_cells_multimapping.rda
## 6 Rda CTdata/eh_data/v2/TCGA_TPM.rda
## 7 Rda CTdata/eh_data/v2/CT_methylation_in_tissues.rda
## 8 Rda CTdata/eh_data/v2/CT_mean_methylation_in_tissues.rda
## 9 Rda CTdata/eh_data/v2/TCGA_CT_methylation.rda
## 10 Rda CTdata/eh_data/v3/CT_genes.rda
## 11 Rda CTdata/eh_data/v3/CCLE_correlation_matrix.rda
## 12 Rda CTdata/eh_data/v2/testis_sce.rda
## 13 Rda CTdata/eh_data/v2/scRNAseq_HPA.rda
## 14 Rda CTdata/eh_data/all_genes.rda
## 15 Rda CTdata/eh_data/FGC_sce.rda
## 16 Rda CTdata/eh_data/HPA_cell_type_specificities.rda
## 17 Rda CTdata/eh_data/mean_methylation_in_tissues.rda
## 18 Rda CTdata/eh_data/methylation_in_tissues.rda
## 19 Rda CTdata/eh_data/oocytes_sce.rda
## 20 Rda CTdata/eh_data/TCGA_methylation.rda
## 21 Rda CTdata/eh_data/embryo_sce_Zhu.rda
## 22 Rda CTdata/eh_data/embryo_sce_Petropoulos.rda
## 23 Rda CTdata/eh_data/hESC_data.rda
## 24 Rda CTdata/eh_data/methylation_in_hESC.rda
## 25 Rda CTdata/eh_data/mean_methylation_in_hESC.rda
## 26 Rda CTdata/eh_data/methylation_in_FGC.rda
## 27 Rda CTdata/eh_data/mean_methylation_in_FGC.rda
## 28 Rda CTdata/eh_data/methylation_in_embryo.rda
## 29 Rda CTdata/eh_data/mean_methylation_in_embryo.rda
## Tags
## 1 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:NormalTissues
## 2 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:CancerData:CellLines
## 3 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:NormalTissues:Multimapping
## 4 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:CellLines:DemethylatingAgent
## 5 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:CellLines:DemethylatingAgent:Multimapping
## 6 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:CancerData:TumorSamples
## 7 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:CancerTestis:NormalTissues
## 8 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:CancerTestis:NormalTissues
## 9 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:CancerTestis:CancerData:TumorSamples
## 10 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:CancerTestis:MethylationData:GeneExpression:CancerData:NormalTissues:TumorSamples
## 11 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:CancerTestis:GeneExpression:CancerData:CellLines
## 12 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:testis
## 13 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:normalCellTypes
## 14 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:GeneExpression:CancerData:NormalTissues:TumorSamples
## 15 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:FetalCells:Gonad
## 16 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneActivation:scRNAseq:testis:normalCellTypes
## 17 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:NormalTissues
## 18 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:NormalTissues
## 19 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:Oocyts:Gonad
## 20 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:CancerData:TumorSamples
## 21 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:embryo
## 22 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:scRNAseq:embryo
## 23 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:GeneExpression:hESC
## 24 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:hESC
## 25 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:hESC
## 26 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:FetalGermCells
## 27 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:FetalGermCells
## 28 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:embryo
## 29 ExperimentHub:ExperimentData:ReproducibleResearch:RepositoryData:Homo_sapiens_Data:MethylationData:embryo
In order to generate the list of CT genes, we followed a specific selection procedure (see figure below).
Testis-specific genes (expressed exclusively in testis) and testis-preferential genes (expressed in a few somatic tissues at a level lower than 10x testis expression) were first selected using the GTEx database (Aguet et al., 2020).
Note that some genes were undetectable in the GTEx database due to
multimapping issues (these were flagged as “lowly expressed” in
GTEX_category
column). A careful inspection of these genes
showed that many of them are well-known Cancer-Testis genes, belonging
to gene families (MAGEA, SSX, CT45A, GAGE, …) from which members have
identical or nearly identical sequences. This is likely the reason why
these genes are not detected in the GTEx database, as GTEx processing
pipeline specifies that overlapping intervals between genes are excluded
from all genes for counting.For these genes, as testis-specificity could
not be assessed using GTEx database, RNAseq data from a set of normal
tissues were reprocessed in order to allow multimapping. Expression of
several genes became detectable in some tissues. Genes showing a
testis-specific expression (expression at least 10x higher in testis
than in any somatic tissues when multimapping was allowed) were
selected, and flagged as testis-specific in
multimapping_analysis
column.
Additionally, as our selection procedure is based on bulk RNAseq data, we wanted to ensure that the selected genes are not expressed in rare populations of somatic cells. We used the Single Cell Type Atlas classification from the Human Protein Atlas (Uhlén et al., 2015) to exclude the ones that were flagged as specific of any somatic cell type.
To assess activation in cancers, RNAseq data from cancer cell lines from CCLE (Barretina et al., 2012) and from TCGA cancer samples (Cancer Genome Atlas Research Network et al., 2013) were used. This allowed to select among testis-specific and testis-preferential genes those that are activated in cancers.
In the CCLE_category
and TCGA_category
columns, genes are tagged as “activated” when they are highly expressed
in at one percent of cancer cell line/sample (TPM >= 1). However
genes that were found to be expressed in all -or almost all-cancer cell
lines/samples were removed, as this probably reflects a constitutive
expression rather than a true activation. We filtered out genes that
were not completely repressed in at least 20 % of cancer cell
lines/samples (TPM <= 0.5). We also made use of the normal
peritumoral samples available in TCGA data to remove from our selection
genes that were already detected in a significant fraction of these
cells.
All selected CT genes were visualised on IGV (Thorvaldsdóttir et al.,
2013) using a RNA-seq alignment from testis, to ensure that expression
in testis really corresponded to the canonical transcript. The aim was
initially to identify precisely the transcription start site of each
gene, but unexpectedly we observed that for some genes, the reads were
not properly aligned on exons, but were instead spread across a wide
genomic region spanning the genes. These genes, flagged as “unclear” in
IGV_backbone
column, were removed from the CT_gene
category, as their expression values in GTEX, TCGA and CCLE might
reflect a poorly defined transcription in these regions and are hence
likely unreliable.
Genes flagged as TRUE
in
regulated_by_methylation
column correspond to
Genes that are significantly induced by a demethylating agent (RNAseq analysis of cell lines reated with DAC (5-Aza-2′-Deoxycytidine)).
Genes that have a highly methylated promoter in normal somatic tissues (WGBS analysis of a set of normal tissues).
For some genes showing a strong activation in cells treated with 5-Aza-2′-Deoxycytidine, methylation analysis was not possible due to multimapping issues. In this case, genes were still considered as regulated by methylation unless their promoter appeared unmethylated in somatic tissues or methylated in germ cells.
For details about functions, see their respective manual pages. For
all functions, an option values_only
can be set to
TRUE
in order to get the values instead of the
visualisation.
All visualisation functions can be used on all GTEx genes, not only
on Cancer-Testis genes, as the data they refer to contains all genes. By
default, if no genes are specified in a function, only strict CT genes
will be used. If one wants all 280 CT genes to be used, an option
include_CTP
can be set to TRUE
to also use CT
preferential genes.
GTEX_expression()
Allows to visualise gene expression in GTEx tissues. We can for example see the difference of expression between testis-specific and testis-preferential genes. Testis-specific genes have been determined with a stricter specificity to the testis : they are lowly expressed in all somatic tissues and at least 10 times more in the testis. Whereas testis-preferential accepts a little expression outside the testis : they are lowly expressed in at least 75% of somatic tissues, but still 10 times more in the testis.
testis_specific <- dplyr::filter(
CT_genes,
CT_gene_type == "CT_gene")
GTEX_expression(testis_specific$external_gene_name, units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
testis_preferential <- dplyr::filter(
CT_genes, CT_gene_type == "CTP_gene")
GTEX_expression(testis_preferential$external_gene_name, units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
normal_tissue_expression_multimapping()
Allows to visualise the expression values obtained by counting or not multi-mapped reads in normal tissues. We can apply that to genes that are lowly expressed in GTEx (as seen above) in order to visualise their expression. First heatmap shows expression without multimapping while multimapping was allowed in the second one, where expression can be observed. This allowed their determination as testis-specific genes.
testis_specific_in_multimapping_analysis <-
dplyr::filter(CT_genes, lowly_expressed_in_GTEX)
normal_tissue_expression_multimapping(
testis_specific_in_multimapping_analysis$external_gene_name,
multimapping = FALSE, units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
normal_tissue_expression_multimapping(
testis_specific_in_multimapping_analysis$external_gene_name,
multimapping = TRUE, units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
testis_expression()
Allows to visualise gene expression in all different testis cell type, somatic or germ, using data from the adult human testis transcriptional cell atlas.
Using CT genes localisation on the X chromosome, we can see that they tend to be expressed in th early stages of spermatogenesis when located on the X. We also visualise clearly that genes mainly expressed in an early stage of spermatogenesis aren’t expressed later and vice-versa.
X_CT <-
dplyr::filter(CT_genes, X_linked)
testis_expression(X_CT$external_gene_name,
cells = "germ_cells")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 12 out of 82 names invalid: PPP4R3C, CENPVL2, CENPVL1, CT47A1, CT45A6,
## CT45A7, CT45A8, CT45A9, MAGEA4-AS1, MAGEA2, GAGE2E, GAGE12C.
## See the manual page for valid types.
## `use_raster` is automatically set to TRUE for a matrix with more than
## 2000 columns You can control `use_raster` argument by explicitly
## setting TRUE/FALSE to it.
##
## Set `ht_opt$message = FALSE` to turn off this message.
## 'magick' package is suggested to install to give better rasterization.
##
## Set `ht_opt$message = FALSE` to turn off this message.
notX_CT <-
dplyr::filter(CT_genes, !X_linked)
testis_expression(notX_CT$external_gene_name,
cells = "germ_cells")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 38 out of 198 names invalid: KAZN-AS1, NFIA-AS1, LRRC52-AS1, ZC3H11B,
## THORLNC, LINC02377, SLC7A11-AS1, LINC02241, LINC02228, LINC03004,
## SEC61G-DT, ZNF722, GTF2I-AS1, SPATA31C2, TEX48, LINC02663, LINC02750,
## FAM230C, GOLGA6L7, LINC02864, ZNF723, VCY, VCY1B, RBMY1B, LINC01783,
## DHCR24-DT, LINC03102, C2orf74-AS1, ZNF775-AS1, MIR3150BHG, CHKA-DT,
## UVRAG-DT, CCDC196, POLG-DT, APPBP2-DT, LINC01970, TUBB8B, FAM230I.
## See the manual page for valid types.
## `use_raster` is automatically set to TRUE for a matrix with more than
## 2000 columns You can control `use_raster` argument by explicitly
## setting TRUE/FALSE to it.
##
## Set `ht_opt$message = FALSE` to turn off this message.
## 'magick' package is suggested to install to give better rasterization.
##
## Set `ht_opt$message = FALSE` to turn off this message.
oocytes_expression()
Allows to visualise gene expression oocytes at different stage, using scRNA-seq data form “Decoding dynamic epigenetic landscapes in human oocytes using single-cell multi-omics sequencing” (Yan et al. Cell Stem Cell 2021)
We can here again compare the expression of CT genes that are located or not on the X chromosome. Showing a lesser expression of X-linked genes in oocytes.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 4 out of 82 names invalid: PPP4R3C, XAGE1A, TEX13C, MAGEA4-AS1.
## See the manual page for valid types.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 52 out of 198 names invalid: KAZN-AS1, LRRC52-AS1, ZC3H11B, THORLNC,
## LINC02377, LINC02241, LINC02228, LINC03004, SEC61G-DT, ZNF722,
## GTF2I-AS1, LINC02663, LINC02750, FAM230C, CBY2, GARIN2, GOLGA6L7,
## CDRT15, CHCT1, LINC02864, TSPY9, DHCR24-DT, LINC03102, SPMIP3,
## C2orf74-AS1, NECTIN3-AS1, CFAP92, LINC02901, TMEM270, ERVW-1, GARIN1B,
## LLCFC1, ZNF775-AS1, PPP1R3B-DT, DNAJC5B, MIR3150BHG, C10orf55, CHKA-DT,
## UVRAG-DT, RB1-DT, CCDC196, INSYN1-AS1, POLG-DT, TEKT5, KCTD19,
## APPBP2-DT, MARCHF10, TUBB8B, SPMAP2, SAXO5, EFCAB8, FAM230I.
## See the manual page for valid types.
HPA_cell_type_expression()
Allows to visualise gene expression in all different healthy cell type, comparing different somatic cells to germ cells, using data from the human proteome atlas.
Visualising all CT_genes, the specificity to germ cells only is quite clear.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 29 out of 146 names invalid: KAZN-AS1, NFIA-AS1, LRRC52-AS1, LINC01249,
## THORLNC, LINC01206, LINC02475, LINC02377, SLC7A11-AS1, LINC01098,
## LINC01511, LINC02241, LINC02228, LINC03004, SEC61G-DT, GTF2I-AS1,
## LINC00200, LINC02663, LINC01518, LINC02750, TMEM132D-AS1, FAM230C,
## LINC01193, LINC01413, FLJ36000, LINC00470, LINC02864, DSCR8,
## MAGEA4-AS1.
## See the manual page for valid types.
No analysis was added to the all_genes
and
CT_genes
tables using these datasets.
embryo_expression()
Allows to visualise gene expression in human early embryos using two different scRNA-seq datasets.
This dataset contains different stages of blastocysts and morula.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 62 out of 146 names invalid: KAZN-AS1, LRRC52-AS1, ZC3H11B, THORLNC,
## LINC02475, LINC02377, LINC01511, LINC02241, LINC02228, LINC03004,
## SEC61G-DT, ZNF722, GTF2I-AS1, TEX48, LINC02663, LINC01518, LINC02750,
## MAJIN, TMEM132D-AS1, CBY2, GARIN2, GOLGA6L7, CHCT1, LINC02864, ZNF723,
## PPP4R3C, GAGE12E, GAGE2A, CENPVL2, CENPVL1, XAGE1A, XAGE1B, SSX2,
## SSX2B, CXorf49, CXorf49B, NXF2, CT47A12, CT47A11, CT47A1, TEX13C,
## CT45A6, CT45A7, CT45A8, CT45A9, CT45A10, SPANXB1, CXorf51B, MAGEA9B,
## MAGEA9, MAGEA4-AS1, MAGEA2, CTAG1A, CTAG1B, TSPY4, TSPY3, TSPY9, VCY,
## VCY1B, RBMY1B, RBMY1A1, RBMY1E.
## See the manual page for valid types.
This dataset only contains blastocyst stage.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 36 out of 146 names invalid: KAZN-AS1, LRRC52-AS1, ZC3H11B, THORLNC,
## LINC02475, LINC02377, LINC02241, LINC02228, LINC03004, SEC61G-DT,
## ZNF722, GTF2I-AS1, TEX48, LINC02663, LINC02750, MAJIN, TMEM132D-AS1,
## CBY2, GARIN2, GOLGA6L7, CHCT1, LINC02864, ZNF723, PPP4R3C, CENPVL2,
## CENPVL1, XAGE1A, CXorf49B, CT47A12, CT47A11, TEX13C, CT45A3, CT45A9,
## CXorf51B, MAGEA4-AS1, TSPY9.
## See the manual page for valid types.
fetal_germcells_expression()
Allows to visualise gene expression in fetal germ cells the scRNA-seq dataset from “Single-cell roadmap of human gonadal development” (Garcia-Alonso, Nature 2022).
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## Warning: 16 out of 146 names invalid: NFIA-AS1, SLC7A11-AS1, GTF2I-AS1,
## SPATA31C2, TEX48, LINC02864, CXorf49B, CT47A12, CT47A11, CT47A1,
## CT45A6, CT45A7, CT45A8, CT45A9, TSPY9, RBMY1B.
## See the manual page for valid types.
CCLE_expression()
Allows to visualise gene expression in different histological types of CCLE cancer cell lines. We can thus compare genes that are or not activated in tumoral cell lines.
frequently_activated <- dplyr::filter(
CT_genes,
percent_of_positive_CCLE_cell_lines >= 5)
CCLE_expression(
genes = frequently_activated$external_gene_name,
type = c(
"lung", "skin", "bile_duct", "bladder", "colorectal",
"lymphoma", "uterine", "myeloma", "kidney",
"pancreatic", "brain", "gastric", "breast", "bone",
"head_and_neck", "ovarian", "sarcoma", "leukemia",
"esophageal", "neuroblastoma"),
units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
not_frequently_activated <- dplyr::filter(
CT_genes,
percent_of_negative_CCLE_cell_lines >= 95)
CCLE_expression(
genes = not_frequently_activated$external_gene_name,
type = c(
"lung", "skin", "bile_duct", "bladder", "colorectal",
"lymphoma", "uterine", "myeloma", "kidney",
"pancreatic", "brain", "gastric", "breast", "bone",
"head_and_neck", "ovarian", "sarcoma", "leukemia",
"esophageal", "neuroblastoma"),
units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
TCGA_expression()
Allows to visualise gene expression in cancer samples from TCGA. We can thus for example compare with activation in CCLE cell lines.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## `use_raster` is automatically set to TRUE for a matrix with more than
## 2000 columns You can control `use_raster` argument by explicitly
## setting TRUE/FALSE to it.
##
## Set `ht_opt$message = FALSE` to turn off this message.
## 'magick' package is suggested to install to give better rasterization.
##
## Set `ht_opt$message = FALSE` to turn off this message.
TCGA_expression(
genes = not_frequently_activated$external_gene_name,
tumor = "all",
units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## `use_raster` is automatically set to TRUE for a matrix with more than
## 2000 columns You can control `use_raster` argument by explicitly
## setting TRUE/FALSE to it.
##
## Set `ht_opt$message = FALSE` to turn off this message.
## 'magick' package is suggested to install to give better rasterization.
##
## Set `ht_opt$message = FALSE` to turn off this message.
When visualising only one tumor, there is a separation between tumoral and peritumoral tissue.
TCGA_expression(
genes = frequently_activated$external_gene_name,
tumor = "LUAD",
units = "log_TPM")
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
DAC_induction()
Allows to visualise gene induction upon DAC treatment in a series of cell lines. We can see the difference between CT genes regulated by methylation or not
controlled_by_methylation <- dplyr::filter(CT_genes, regulated_by_methylation)
DAC_induction(genes = controlled_by_methylation$external_gene_name)
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
not_controlled_by_methylation <- dplyr::filter(
CT_genes,
!regulated_by_methylation)
DAC_induction(genes = not_controlled_by_methylation$external_gene_name)
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
normal_tissues_methylation()
Gives the methylation level of all CpGs located in each promoter region (by default defined as 1000 nt upstream TSS and 200 nt downstream TSS) in a set of normal tissues. Can be used to analyse more precisely the methylation of a particular promoter in normal tissues as methylation values are given CpG by CpG.
normal_tissues_mean_methylation()
Gives the mean methylation level of CpGs located in each promoter region (by default defined as 1000 nt upstream TSS and 200 nt downstream TSS) in a set of normal tissues. When comparing genes controlled or not by methylation, we can see a specific methylation pattern in the first group.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
embryos_mean_methylation()
Gives the mean methylation level of CpGs located in each promoter region of any genes in early embryos, using WGSB data from (“Single-cell DNA methylome sequencing of human preimplantation embryos”. Zhu et al. Nat genetics 2018). Methylation levels in tissues correspond to the mean methylation of CpGs located in range of 1000 pb upstream and 500 pb downstream from gene TSS.
Can be used to explore well known CT genes methylation in embryos as below.
embryos_mean_methylation(c("MAGEA1", "MAGEA3", "MAGEA4", "MAGEC2", "MAGEB16"),
stage = c( "MII Oocyte", "Sperm", "Zygote", "2-cell",
"4-cell", "8-cell", "Morula"))
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
fetal_germcells_mean_methylation()
Allows to visualise mean promoter methylation levels ofany genes in fetal germ cells, using WGSB data from “Dissecting the epigenomic dynamics of human fetal germ cell development at single-cell resolution” (Li et al. 2021). Methylation levels in tissues correspond to the mean methylation of CpGs located in range of 1000 pb upstream and 500 pb downstream from gene TSS.
Can be used to explore well known CT genes methylation in fetal germ cells, as below.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
hESC_mean_methylation()
Allows to visualise mean promoter methylation levels of any genes in human embryonic cell lines. WGBS methylation data was downloaded from Encode. Methylation levels in tissues correspond to the mean methylation of CpGs located in range of 1000 pb upstream and 200 pb downstream from gene TSS.
Here used on strict CT genes.
## see ?CTdata and browseVignettes('CTdata') for documentation
## loading from cache
TCGA_methylation_expression_correlation()
Shows the correlation between gene expression and promoter methylation in TCGA samples.
CTexploreR’s heatmaps are compatible with the
InteractiveComplexHeatmap
package. This allows the navigate
in the last heatmap generated, search a specific gene as well as
creating subsets to facilitate data exploration. More info about
InteractiveComplexHeatmap
available here.
Aguet, F., Anand, S., Ardlie, K. G., Gabriel, S., Getz, G. A., Graubert, A., Hadley, K., Handsaker, R. E., Huang, K. H., Kashin, S., Li, X., MacArthur, D. G., Meier, S. R., Nedzel, J. L., Nguyen, D. T., Segrè, A. V., Todres, E., Balliu, B., Barbeira, A. N., … Volpi, S. (2020). The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science. https://science.sciencemag.org/content/369/6509/1318.abstract
Almeida, L. G., Sakabe, N. J., deOliveira, A. R., Silva, M. C. C., Mundstein, A. S., Cohen, T., Chen, Y.-T., Chua, R., Gurung, S., Gnjatic, S., Jungbluth, A. A., Caballero, O. L., Bairoch, A., Kiesler, E., White, S. L., Simpson, A. J. G., Old, L. J., Camargo, A. A., & Vasconcelos, A. T. R. (2009). CTdatabase: a knowledge-base of high-throughput and curated data on cancer-testis antigens. Nucleic Acids Research, 37(Database issue), D816–D819.
Barretina, J., Caponigro, G., Stransky, N., Venkatesan, K., Margolin, A. A., Kim, S., Wilson, C. J., Lehár, J., Kryukov, G. V., Sonkin, D., Reddy, A., Liu, M., Murray, L., Berger, M. F., Monahan, J. E., Morais, P., Meltzer, J., Korejwa, A., Jané-Valbuena, J., … Garraway, L. A. (2012). The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature, 483(7391), 603–607.
Cancer Genome Atlas Research Network, Weinstein, J. N., Collisson, E. A., Mills, G. B., Shaw, K. R. M., Ozenberger, B. A., Ellrott, K., Shmulevich, I., Sander, C., & Stuart, J. M. (2013). The Cancer Genome Atlas Pan-Cancer analysis project. Nature Genetics, 45(10), 1113–1120.
Carter, J. A., Matta, B., Battaglia, J., Somerville, C., Harris, B. D., LaPan, M., Atwal, G. S., & Barnes, B. J. (2023). Identification of pan-cancer/testis genes and validation of therapeutic targeting in triple-negative breast cancer: Lin28a- and Siglece-based vaccination induces anti-tumor immunity and inhibits metastasis. In bioRxiv (p. 2023.05.09.539617).
Guo, J. et al. The adult human testis transcriptional cell atlas. Cell Res. 28, 1141–1157 (2018).
Uhlén, M. et al. Proteomics. Tissue-based map of the human proteome. Science 347, 1260419 (2015).
Thorvaldsdóttir, H., Robinson, J. T., & Mesirov, J. P. (2013). Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Briefings in Bioinformatics, 14(2), 178–192.
Wang, C., Gu, Y., Zhang, K., Xie, K., Zhu, M., Dai, N., Jiang, Y., Guo, X., Liu, M., Dai, J., Wu, L., Jin, G., Ma, H., Jiang, T., Yin, R., Xia, Y., Liu, L., Wang, S., Shen, B., … Hu, Z. (2016). Systematic identification of genes with a cancer-testis expression pattern in 19 cancer types. Nature Communications, 7, 10499.
## R version 4.4.2 (2024-10-31)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.1 LTS
##
## Matrix products: default
## BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
## LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: Etc/UTC
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] SingleCellExperiment_1.29.0 SummarizedExperiment_1.37.0
## [3] Biobase_2.67.0 GenomicRanges_1.59.0
## [5] GenomeInfoDb_1.43.0 IRanges_2.41.0
## [7] S4Vectors_0.45.0 BiocGenerics_0.53.1
## [9] generics_0.1.3 MatrixGenerics_1.19.0
## [11] matrixStats_1.4.1 CTexploreR_1.3.0
## [13] CTdata_1.7.0 BiocStyle_2.35.0
##
## loaded via a namespace (and not attached):
## [1] DBI_1.2.3 rlang_1.1.4 magrittr_2.0.3
## [4] clue_0.3-65 GetoptLong_1.0.5 compiler_4.4.2
## [7] RSQLite_2.3.7 png_0.1-8 vctrs_0.6.5
## [10] pkgconfig_2.0.3 shape_1.4.6.1 crayon_1.5.3
## [13] fastmap_1.2.0 dbplyr_2.5.0 XVector_0.47.0
## [16] labeling_0.4.3 utf8_1.2.4 rmarkdown_2.29
## [19] UCSC.utils_1.3.0 purrr_1.0.2 bit_4.5.0
## [22] xfun_0.49 zlibbioc_1.52.0 cachem_1.1.0
## [25] jsonlite_1.8.9 blob_1.2.4 highr_0.11
## [28] DelayedArray_0.33.1 parallel_4.4.2 cluster_2.1.6
## [31] R6_2.5.1 bslib_0.8.0 RColorBrewer_1.1-3
## [34] jquerylib_0.1.4 Rcpp_1.0.13-1 iterators_1.0.14
## [37] knitr_1.48 Matrix_1.7-1 tidyselect_1.2.1
## [40] abind_1.4-8 yaml_2.3.10 doParallel_1.0.17
## [43] codetools_0.2-20 curl_6.0.0 lattice_0.22-6
## [46] tibble_3.2.1 withr_3.0.2 KEGGREST_1.47.0
## [49] evaluate_1.0.1 BiocFileCache_2.15.0 circlize_0.4.16
## [52] ExperimentHub_2.15.0 Biostrings_2.75.0 pillar_1.9.0
## [55] BiocManager_1.30.25 filelock_1.0.3 foreach_1.5.2
## [58] BiocVersion_3.21.1 ggplot2_3.5.1 munsell_0.5.1
## [61] scales_1.3.0 glue_1.8.0 maketools_1.3.1
## [64] tools_4.4.2 AnnotationHub_3.15.0 sys_3.4.3
## [67] buildtools_1.0.0 grid_4.4.2 tidyr_1.3.1
## [70] AnnotationDbi_1.69.0 colorspace_2.1-1 GenomeInfoDbData_1.2.13
## [73] cli_3.6.3 rappdirs_0.3.3 fansi_1.0.6
## [76] S4Arrays_1.7.1 ComplexHeatmap_2.23.0 dplyr_1.1.4
## [79] gtable_0.3.6 sass_0.4.9 digest_0.6.37
## [82] SparseArray_1.7.1 ggrepel_0.9.6 farver_2.1.2
## [85] rjson_0.2.23 memoise_2.0.1 htmltools_0.5.8.1
## [88] lifecycle_1.0.4 httr_1.4.7 GlobalOptions_0.1.2
## [91] mime_0.12 bit64_4.5.2