if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
# orthogene is only available on Bioconductor>=3.14
if(BiocManager::version()<"3.14")
BiocManager::install(update = TRUE, ask = FALSE)
BiocManager::install("orthogene")
library(orthogene)
data("exp_mouse")
# Setting to "homologene" for the purposes of quick demonstration.
# We generally recommend using method="gprofiler" (default).
method <- "homologene"
It’s not always clear whether a dataset is using the original species gene names, human gene names, or some other species’ gene names.
infer_species takes a list/matrix/data.frame with genes and
infers the species that they best match to!
For the sake of speed, the genes extracted from gene_df
are tested against genomes from only the following 6 test_species by default:
- human
- monkey
- rat
- mouse
- zebrafish
- fly
However, you can supply your own list of test_species, which will
be automatically be mapped and standardised using map_species.
matches <- orthogene::infer_species(gene_df = exp_mouse,
method = method)
## Preparing gene_df.
## sparseMatrix format detected.
## Extracting genes from rownames.
## 15,259 genes extracted.
## Testing for gene overlap with: human
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: human
## Common name mapping found for human
## 1 organism identified from search: 9606
## Gene table with 19,129 rows retrieved.
## Returning all 19,129 genes from human.
## Testing for gene overlap with: monkey
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: monkey
## Common name mapping found for monkey
## 1 organism identified from search: 9544
## Gene table with 16,843 rows retrieved.
## Returning all 16,843 genes from monkey.
## Testing for gene overlap with: rat
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: rat
## Common name mapping found for rat
## 1 organism identified from search: 10116
## Gene table with 20,616 rows retrieved.
## Returning all 20,616 genes from rat.
## Testing for gene overlap with: mouse
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: mouse
## Common name mapping found for mouse
## 1 organism identified from search: 10090
## Gene table with 21,207 rows retrieved.
## Returning all 21,207 genes from mouse.
## Testing for gene overlap with: zebrafish
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: zebrafish
## Common name mapping found for zebrafish
## 1 organism identified from search: 7955
## Gene table with 20,897 rows retrieved.
## Returning all 20,897 genes from zebrafish.
## Testing for gene overlap with: fly
## Retrieving all genes using: homologene.
## Retrieving all organisms available in homologene.
## Mapping species name: fly
## Common name mapping found for fly
## 1 organism identified from search: 7227
## Gene table with 8,438 rows retrieved.
## Returning all 8,438 genes from fly.
## Top match:
## - species: mouse
## - percent_match: 92%
To create an example dataset, turn the gene names into rat genes.
exp_rat <- orthogene::convert_orthologs(gene_df = exp_mouse,
input_species = "mouse",
output_species = "rat",
method = method)
matches <- orthogene::infer_species(gene_df = exp_rat,
method = method)
To create an example dataset, turn the gene names into human genes.
exp_human <- orthogene::convert_orthologs(gene_df = exp_mouse,
input_species = "mouse",
output_species = "human",
method = method)
matches <- orthogene::infer_species(gene_df = exp_human,
method = method)
test_speciesYou can even supply test_species with the name of one of the R packages that
orthogene gets orthologs from. This will test against all species available
in that particular R package.
For example, by setting test_species="homologene" we automatically test for
% gene matches in each of the 20+ species available in homologene.
matches <- orthogene::infer_species(gene_df = exp_human,
test_species = method,
method = method)
utils::sessionInfo()
R version 4.5.1 Patched (2025-09-10 r88807)
Platform: x86_64-apple-darwin20
Running under: macOS Monterey 12.7.6
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/4.5-x86_64/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.5-x86_64/Resources/lib/libRlapack.dylib; LAPACK version 3.12.1
locale:
[1] C/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
time zone: America/New_York
tzcode source: internal
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] orthogene_1.15.02 BiocStyle_2.37.1
loaded via a namespace (and not attached):
[1] gtable_0.3.6 babelgene_22.9
[3] xfun_0.53 bslib_0.9.0
[5] ggplot2_4.0.0 htmlwidgets_1.6.4
[7] rstatix_0.7.2 lattice_0.22-7
[9] vctrs_0.6.5 tools_4.5.1
[11] generics_0.1.4 yulab.utils_0.2.1
[13] parallel_4.5.1 tibble_3.3.0
[15] pkgconfig_2.0.3 Matrix_1.7-4
[17] data.table_1.17.8 homologene_1.4.68.19.3.27
[19] ggplotify_0.1.3 RColorBrewer_1.1-3
[21] S7_0.2.0 lifecycle_1.0.4
[23] compiler_4.5.1 farver_2.1.2
[25] treeio_1.33.0 tinytex_0.57
[27] carData_3.0-5 ggtree_3.99.0
[29] gprofiler2_0.2.3 fontLiberation_0.1.0
[31] fontquiver_0.2.1 ggfun_0.2.0
[33] htmltools_0.5.8.1 sass_0.4.10
[35] yaml_2.3.10 lazyeval_0.2.2
[37] plotly_4.11.0 Formula_1.2-5
[39] pillar_1.11.1 car_3.1-3
[41] ggpubr_0.6.1 jquerylib_0.1.4
[43] tidyr_1.3.1 cachem_1.1.0
[45] grr_0.9.5 magick_2.9.0
[47] abind_1.4-8 fontBitstreamVera_0.1.1
[49] nlme_3.1-168 tidyselect_1.2.1
[51] aplot_0.2.9 digest_0.6.37
[53] dplyr_1.1.4 purrr_1.1.0
[55] bookdown_0.45 labeling_0.4.3
[57] fastmap_1.2.0 grid_4.5.1
[59] cli_3.6.5 magrittr_2.0.4
[61] patchwork_1.3.2 dichromat_2.0-0.1
[63] broom_1.0.10 ape_5.8-1
[65] withr_3.0.2 gdtools_0.4.4
[67] scales_1.4.0 backports_1.5.0
[69] rappdirs_0.3.3 httr_1.4.7
[71] rmarkdown_2.30 ggsignif_0.6.4
[73] evaluate_1.0.5 knitr_1.50
[75] viridisLite_0.4.2 gridGraphics_0.5-1
[77] rlang_1.1.6 ggiraph_0.9.2
[79] Rcpp_1.1.0 glue_1.8.0
[81] tidytree_0.4.6 BiocManager_1.30.26
[83] jsonlite_2.0.0 R6_2.6.1
[85] systemfonts_1.3.1 fs_1.6.6