. . . . "WorkflowHub" . "https://about.workflowhub.eu/" . . "Workflow RO-Crate Profile" . "0.2.0" . . "Małgorzata Wolniewicz" . . "https://doi.org/10.1093/bioinformatics/bts480" . "Snakemake" . "https://snakemake.readthedocs.io/" . . "https://orcid.org/0000-0003-3876-6581" . "Agata Kilar" . . "CWL" . "https://w3id.org/cwl/v1.0/" . "Common Workflow Language" . "https://www.commonwl.org/" . . "analysis" . . "3.0936454849498327" . "3.7" . "NCBI database" . . "6.079027355623101" . "4.0" . "directory from the file" . . "2.8875379939209727" . "1.9" . "name" . . "5.1003344481605355" . "6.1" . "GERONIMO is a bioinformatics pipeline designed to conduct high-throughput homology searches of structural genes using covariance models." . . "27.215189873417724" . "4.3" . "pipeline" . . "8.646003262642742" . "5.3" . "genome" . . "10.114192495921698" . "6.2" . "computer science" . . "58.58585858585859" . "29.0" . "documentation file" . . "2.8875379939209727" . "1.9" . "covariance" . . "8.646003262642742" . "5.3" . "covariance model" . . "2.7355623100303954" . "1.8" . "database" . . "7.34094616639478" . "4.5" . "installation of Geronimo dependency" . . "2.8875379939209727" . "1.9" . "Geronimo" . . "4.893964110929853" . "3.0" . "atmospheric sciences" . . "100.0" . "0.9349703192710876" . "The" . . "6.035889070146819" . "3.7" . "IT-computer sciences" . . "Science and technology/Technology and engineering/IT-computer sciences" . "gene family" . . "2.759197324414716" . "3.3" . "software" . . "4.1806020066889635" . "5.0" . "mathematical and computer sciences" . . "100.0" . "0.41168108582496643" . "RNA" . . "8.156606851549755" . "5.0" . "GERONIMO installation" . . "3.0395136778115504" . "2.0" . "Genetics" . . "Science and technology/Natural science/Biology/Genetics" . "GERONIMO pipeline" . . "6.231003039513677" . "4.1" . "Oil and gas - upstream activities" . . "Economy, business and finance/Economic sector/Energy and resource/Oil and gas - upstream activities" . "GERONIMO repository" . . "3.0395136778115504" . "2.0" . "summary file" . . "3.4954407294832825" . "2.3" . "storage space" . . "3.5953177257525084" . "4.3" . "license" . . "4.431438127090301" . "5.3" . "database" . . "2.9264214046822743" . "3.5" . "question" . . "4.933110367892977" . "5.9" . "directory structure" . . "15.653495440729484" . "10.3" . "You can safely remove the following:\r\n- `GERONIMO/results`\r\n- `GERONIMO/database`\r\n- `GERONIMO/taxonomy`\r\n- `GERONIMO/temp`\r\n- `.create_genome_list.touch`\r\n- `list_of_genomes.txt`" . . "30.37974683544304" . "4.8" . "By default, the GERONIMO pipeline conducts high-throughput searches of homology sequences in downloaded genomes utilizing covariance models." . . "21.518987341772153" . "3.4" . "service-account-enrichment" . . . . "1715468"^^ . "https://api.rohub.org/api/ros/892d05bf-da84-4be9-9ab7-60eba5f73c70/crate/download/" . . "2023-09-08 12:04:04.643799+00:00" . "2024-03-05 12:23:14.455859+00:00" . "2023-09-08 12:04:04.643799+00:00" . "# GERONIMO\r\n\r\n## Introduction\r\nGERONIMO is a bioinformatics pipeline designed to conduct high-throughput homology searches of structural genes using covariance models. These models are based on the alignment of sequences and the consensus of secondary structures. The pipeline is built using Snakemake, a workflow management tool that allows for the reproducible execution of analyses on various computational platforms. \r\n\r\nThe idea for developing GERONIMO emerged from a comprehensive search for [telomerase RNA in lower plants] and was subsequently refined through an [expanded search of telomerase RNA across Insecta]. GERONIMO can test hundreds of genomes and ensures the stability and reproducibility of the analyses performed.\r\n\r\n\r\n[telomerase RNA in lower plants]: https://doi.org/10.1093/nar/gkab545\r\n[expanded search of telomerase RNA across Insecta]: https://doi.org/10.1093/nar/gkac1202\r\n\r\n## Scope\r\nThe GERONIMO tool utilises covariance models (CMs) to conduct homology searches of RNA sequences across a wide range of gene families in a broad evolutionary context. Specifically, it can be utilised to:\r\n\r\n* Detect RNA sequences that share a common evolutionary ancestor\r\n* Identify and align orthologous RNA sequences among closely related species, as well as paralogous sequences within a single species\r\n* Identify conserved non-coding RNAs in a genome, and extract upstream genomic regions to characterise potential promoter regions. \r\nIt is important to note that GERONIMO is a computational tool, and as such, it is intended to be run on a computer with a small amount of data. Appropriate computational infrastructure is necessary for analysing hundreds of genomes.\r\n\r\nAlthough GERONIMO was primarily designed for Telomerase RNA identification, its functionality extends to include the detection and alignment of other RNA gene families, including **rRNA**, **tRNA**, **snRNA**, **miRNA**, and **lncRNA**. This can aid in identifying paralogs and orthologs across different species that may carry specific functions, making it useful for phylogenetic analyses. \r\n\r\nIt is crucial to remember that some gene families may exhibit similar characteristics but different functions. Therefore, analysing the data and functional annotation after conducting the search is essential to characterise the sequences properly.\r\n\r\n## Pipeline overview\r\n\r\n\r\nBy default, the GERONIMO pipeline conducts high-throughput searches of homology sequences in downloaded genomes utilizing covariance models. If a significant similarity is detected between the model and genome sequence, the pipeline extracts the upstream region, making it convenient to identify the promoter of the discovered gene. In brief, the pipeline:\r\n- Compiles a list of genomes using the NCBI's [Entrez] database based on a specified query, *e.g. \"Rhodophyta\"[Organism]*\r\n- Downloads and decompresses the requested genomes using *rsync* and *gunzip*, respectively\r\n- *Optionally*, generates a covariance model based on a provided alignment using [Infernal]\r\n- Conducts searches among the genomes using the covariance model [Infernal]\r\n- Supplements genome information with taxonomy data using [rentrez]\r\n- Expands the significant hits sequence by extracting upstream genomic regions using [*blastcmd*]\r\n- Compiles the results, organizes them into a tabular format, and generates a visual summary of the performed analysis.\r\n\r\n[Entrez]: https://www.ncbi.nlm.nih.gov/books/NBK179288/\r\n[Infernal]: http://eddylab.org/infernal/\r\n[rentrez]: https://github.com/ropensci/rentrez\r\n[*blastcmd*]: https://www.ncbi.nlm.nih.gov/books/NBK569853/\r\n\r\n## Quick start\r\nThe GERONIMO is available as a `snakemake pipeline` running on Linux and Windows operating systems.\r\n\r\n### Windows 10\r\nInstal Linux on Windows 10 (WSL) according to [instructions], which bottling down to opening PowerShell or Windows Command Prompt in *administrator mode* and pasting the following:\r\n```shell\r\nwsl --install\r\nwsl.exe --install UBUNTU\r\n```\r\nThen restart the machine and follow the instructions for setting up the Linux environment.\r\n\r\n[instructions]: https://learn.microsoft.com/en-us/windows/wsl/install\r\n\r\n### Linux:\r\n#### Check whether the conda is installed:\r\n```shell\r\nconda -V\r\n```\r\n> GERONIMO was tested on conda 23.3.1\r\n#### 1) If you do not have installed `conda`, please install `miniconda`\r\nPlease follow the instructions for installing [miniconda]\r\n\r\n[miniconda]: https://conda.io/projects/conda/en/stable/user-guide/install/linux.html\r\n\r\n#### 2) Continue with installing `mamba` (recommended but optional)\r\n```shell\r\nconda install -n base -c conda-forge mamba\r\n```\r\n#### 3) Install `snakemake`\r\n```shell\r\nconda activate base\r\nmamba create -p env_snakemake -c conda-forge -c bioconda snakemake\r\nmamba activate env_snakemake\r\nsnakemake --help\r\n```\r\nIn case of complications, please check the section `Questions & Answers` below or follow the [official documentation] for troubleshooting.\r\n\r\n[official documentation]: https://snakemake.readthedocs.io/en/stable/getting_started/installation.html\r\n\r\n### Clone the GERONIMO repository\r\nGo to the path in which you want to run the analysis and clone the repository:\r\n```shell\r\ncd \r\ngit clone https://github.com/amkilar/GERONIMO.git\r\n```\r\n\r\n### Run sample analysis to ensure GERONIMO installation was successful\r\nAll files are prepared for the sample analysis as a default. Please execute the line below:\r\n```shell\r\nsnakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n```\r\n\r\nThis will prompt GERONIMO to quickly scan all modules, verifying the correct setup of the pipeline without executing any analysis.\r\nYou should see the message `Building DAG of jobs...`, followed by `Nothing to be done (all requested files are present and up to date).`, when successfully completed.\r\n\r\nIf you want to run the sample analysis fully, please remove the folder `results` from the GERONIMO directory and execute GERONIMO again with:\r\n\r\n`snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx`\r\n\r\n> You might consider allowing more cores to speed up the analysis, which might take up to several hours.\r\n\r\n#### You might want to clean `GERONIMO/` directory from the files produced by the example analysis. You can safely remove the following:\r\n- `GERONIMO/results`\r\n- `GERONIMO/database`\r\n- `GERONIMO/taxonomy`\r\n- `GERONIMO/temp`\r\n- `.create_genome_list.touch`\r\n- `list_of_genomes.txt`\r\n\r\n## Setup the inputs\r\n\r\n### 1) Prepare the `covariance models`:\r\n\r\n#### Browse the collection of available `covariance models` at [Rfam] (*You can find the covariance model in the tab `Curation`.*) \r\nPaste the covariance model to the folder `GERONIMO/models` and ensure its name follows the convention: `cov_model_`\r\n\r\n[Rfam]: https://rfam.org/\r\n\r\n#### **OR**\r\n\r\n#### Prepare your own `covariance model` using [LocARNA]\r\n1. Paste or upload your sequences to the web server and download the `.stk` file with the alignment result. \r\n \r\n > *Please note that the `.stk` file format is crucial for the analysis, containing sequence alignment and secondary structure consensus.*\r\n \r\n > The LocARNA web service allows you to align 30 sequences at once - if you need to align more sequences, please use the standalone version available [here] \r\n > After installation run: \r\n ```shell\r\n mlocarna my_fasta_sequences.fasta\r\n ```\r\n \r\n2. Paste the `.stk` alignment file to the folder `GERONIMO/model_to_build` and ensure its name follows the convention: `.stk`\r\n\r\n > Please check the example `heterotrichea.stk` format in `GERONIMO/models_to_built` for reference\r\n \r\n\r\n[LocARNA]: http://rna.informatik.uni-freiburg.de/LocARNA/Input.jsp\r\n[here]: http://www.bioinf.uni-freiburg.de/Software/LocARNA/\r\n\r\n\r\n### 2) Adjust the `config.yaml` file\r\nPlease adjust the analysis specifications, as in the following example:\r\n\r\n> - database: ' [Organism]' (in case of difficulties with defining the database query, please follow the instructions below)\r\n> - extract_genomic_region-length: (here you can determine how long the upstream genomic region should be extracted; tested for 200)\r\n> - models: [\"\", \"\"] (here specify the names of models that should be used to perform analysis)\r\n> \r\n> *Here you can also insert the name of the covariance model you want to build with GERONIMO - just be sure you placed `.stk` file in `GERONIMO/models_to_build` before starting analysis*\r\n> - CPU_for_model_building: (specify the number of available CPUs devoted to the process of building model (cannot exceed the CPU number allowed to snakemake with `--cores`)\r\n>\r\n> *You might ignore this parameter when you do not need to create a new covariance model*\r\n\r\n\r\nKeep in mind that the covariance models and alignments must be present in the respective GERONIMO folders.\r\n \r\n### 3) Remove folder `results`, which contains example analysis output\r\n### 4) **Please ensure you have enough storage capacity to download all the requested genomes (in the `GERONIMO/` directory)**\r\n\r\n## Run GERONIMO\r\n```shell\r\nmamba activate env_snakemake\r\ncd ~/GERONIMO\r\nsnakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n```\r\n \r\n## Example results\r\n\r\n### Outputs characterisation\r\n\r\n#### A) Summary table\r\nThe Excel table contains the results arranged by taxonomy information and hit significance. The specific columns include:\r\n* family, organism_name, class, order, phylum (taxonomy context)\r\n* GCA_id - corresponds to the genome assembly in the *NCBI database*\r\n* model - describes which covariance model identified the result\r\n* label - follows the *Infernal* convention of categorizing hits\r\n* number - the counter of the result\r\n* e_value - indicates the significance level of the hit\r\n* HIT_sequence - the exact HIT sequence found by *Infernal*, which corresponds to the covariance model\r\n* HIT_ID - describes in which part of the genome assembly the hit was found, which may help publish novel sequences\r\n* extended_genomic_region - upstream sequence, which may contain a possible promoter sequence\r\n* secondary_structure - the secondary structure consensus of the covariance model\r\n\r\n\r\n#### B) Significant Hits Distribution Across Taxonomy Families\r\nThe plot provides an overview of the number of genomes in which at least one significant hit was identified, grouped by family. The bold black line corresponds to the number of genomes present in each family, helping to minimize bias regarding unequal data representation across the taxonomy.\r\n\r\n\r\n#### C) Hits Distribution in Genomes Across Families\r\nThe heatmap provides information about the most significant hits from the genome, identified by a specific covariance model. Genomes are grouped by families (on the right). Hits are classified into three categories based on their e-values. Generally, these categories correspond to hit classifications (\"HIT,\" \"MAYBE,\" \"NO HIT\"). The \"HIT\" category is further divided to distinguish between highly significant hits and moderately significant ones.\r\n\r\n\r\n\r\n### GERONIMO directory structure\r\n\r\nThe GERONIMO directory structure is designed to produce files in a highly structured manner, ensuring clear insight and facilitating the analysis of results. During a successful run, GERONIMO produces the following folders:\r\n* `/database` - which contains genome assemblies that were downloaded from the *NCBI database* and grouped in subfolders\r\n* `/taxonomy` - where taxonomy information is gathered and stored in the form of tables\r\n* `/results` - the main folder containing all produced results:\r\n * `/infernal_raw` - contains the raw results produced by *Infernal*\r\n * `/infernal` - contains restructured results of *Infernal* in table format\r\n * `/cmdBLAST` - contains results of *cmdblast*, which extracts the extended genomic region\r\n * `/summary` - contains summary files that join results from *Infernal*, *cmdblast*, and attach taxonomy context\r\n * `/plots` - contains two types of summary plots\r\n* `/temp` - folder contains the information necessary to download genome assemblies from *NCBI database*\r\n\r\n* `/env` - stores instructions for dependency installation\r\n* `/models` - where calibrated covariance models can be pasted, *for example, from the Rfam database*\r\n* `/modes_to_built` - where multiple alignments in *.stk* format can be pasted\r\n* `/scripts` - contains developed scripts that perform results structurization\r\n\r\n#### The example GERONIMO directory structure:\r\n\r\n```shell\r\nGERONIMO\r\n├── database\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   └── ...\r\n├── env\r\n├── models\r\n├── model_to_build\r\n├── results\r\n│   ├── cmdBLAST\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── ...\r\n│   ├── infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   ├── plots\r\n│   ├── raw_infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   └── summary\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   ├── ...\r\n├── scripts\r\n├── taxonomy\r\n└── temp\r\n```\r\n\r\n## GERONIMO applicability\r\n\r\n### Expanding the evolutionary context\r\nTo add new genomes or database queries to an existing analysis, please follow the instructions:\r\n1) Rename the `list_of_genomes.txt` file to `previous_list_of_genomes.txt` or any other preferred name.\r\n2) Modify the `config.yaml` file by replacing the previous database query with the new one.\r\n3) Delete:\r\n - `summary_table.xlsx`, `part_summary_table.csv`, `summary_table_models.xlsx` files located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n5) Run GERONIMO to calculate new results using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n ```\r\n7) Once the new results are generated, reviewing them before merging them with the original results is recommended.\r\n8) Copy the contents of the `previous_list_of_genomes.txt` file and paste them into the current `list_of_genomes.txt`.\r\n9) Delete:\r\n - `summary_table.xlsx` located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n10) Run GERONIMO to merge the results from both analyses using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n ```\r\n\r\n### Incorporating new covariance models into existing analysis\r\n1) Copy the new covariance model to `GERONIMO/models`\r\n2) Modify the `config.yaml` file by adding the name of the new model to the line `models: [...]`\r\n3) Run GERONIMO to see the updated analysis outcome\r\n\r\n### Building a new covariance model\r\nWith GERONIMO, building a new covariance model from multiple sequence alignment in the `.stk` format is possible. \r\n\r\nTo do so, simply paste `.stk` file to `GERONIMO/models_to_build` and paste the name of the new covariance model to `config.yaml` file to the line `models: [\"\"]`\r\n\r\nand run GERONIMO.\r\n\r\n\r\n## Questions & Answers\r\n\r\n### How to specify the database query?\r\n- Visit the [NCBI Assemblies] website. \r\n- Follow the instruction on the graphic below:\r\n\r\n[NCBI Assemblies]: https://www.ncbi.nlm.nih.gov/assembly/?term=\r\n\r\n### WSL: problem with creating `snakemake_env`\r\nIn the case of an error similar to the one below:\r\n> CondaError: Unable to create prefix directory '/mnt/c/Windows/system32/env_snakemake'.\r\n> Check that you have sufficient permissions. \r\n \r\nYou might try to delete the cache with: `rm -r ~/.cache/` and try again.\r\n\r\n### When `snakemake` does not seem to be installed properly\r\nIn the case of the following error:\r\n> Command 'snakemake' not found ...\r\n\r\nCheck whether the `env_snakemake` is activated.\r\n> It should result in a change from (base) to (env_snakemake) before your login name in the command line window.\r\n\r\nIf you still see `(base)` before your login name, please try to activate the environment with conda:\r\n`conda activate env_snakemake`\r\n\r\n\r\nPlease note that you might need to specify the full path to the `env_snakemake`, like /home/your user name/env_snakemake\r\n\r\n### How to browse GERONIMO results obtained in WSL?\r\nYou can easily access the results obtained on WSL from your Windows environment by opening `File Explorer` and pasting the following line into the search bar: `\\\\wsl.localhost\\Ubuntu\\home\\`. This will reveal a folder with your username, as specified during the configuration of your Ubuntu system. To locate the GERONIMO results, simply navigate to the folder with your username and then to the `home` folder. (`\\\\wsl.localhost\\Ubuntu\\home\\\\home\\GERONIMO`)\r\n\r\n### GERONIMO occupies a lot of storage space\r\nThrough genome downloads, GERONIMO can potentially consume storage space, rapidly leading to a shortage. Currently, downloading genomes is an essential step for optimal GERONIMO performance.\r\n\r\nRegrettably, if the analysis is rerun without the `/database` folder, it will result in the need to redownload genomes, which is a highly time-consuming process.\r\n\r\nNevertheless, if you do not intend to repeat the analysis and have no requirement for additional genomes or models, you are welcome to retain your results tables and plots while removing the remaining files.\r\n\r\nIt is strongly advised against using local machines for extensive analyses. If you lack access to external storage space, it is recommended to divide the analysis into smaller segments, which can be later merged, as explained in the section titled `Expanding the evolutionary context`.\r\n\r\nConsidering this limitation, I am currently working on implementing a solution that will help circumvent the need for redundant genome downloads without compromising GERONIMO performance in the future.\r\n\r\nYou might consider deleting the `.snakemake` folder to free up storage space. However, please note that deleting this folder will require the reinstallation of GERONIMO dependencies when the analysis is rerun.\r\n\r\n## License\r\nCopyright (c) 2023 Agata M. Kilar\r\n\r\nPermission is hereby granted, free of charge, to any person obtaining a copy\r\nof this software and associated documentation files (the \"Software\"), to deal\r\nin the Software without restriction, including without limitation the rights\r\nto use, copy, modify, merge, publish, distribute, sublicense, and/or sell\r\ncopies of the Software, and to permit persons to whom the Software is\r\nfurnished to do so, subject to the following conditions:\r\n\r\nThe above copyright notice and this permission notice shall be included in all\r\ncopies or substantial portions of the Software.\r\n\r\nTHE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR\r\nIMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,\r\nFITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE\r\nAUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER\r\nLIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,\r\nOUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE\r\nSOFTWARE.\r\n\r\n## Contact\r\nmgr inż. Agata Magdalena Kilar, PhD (agata.kilar@gmail.com)\r\n\r\n" . "# GERONIMO\r\n\r\n## Introduction\r\nGERONIMO is a bioinformatics pipeline designed to conduct high-throughput homology searches of structural genes using covariance models. These models are based on the alignment of sequences and the consensus of secondary structures. The pipeline is built using Snakemake, a workflow management tool that allows for the reproducible execution of analyses on various computational platforms. \r\n\r\nThe idea for developing GERONIMO emerged from a comprehensive search for [telomerase RNA in lower plants] and was subsequently refined through an [expanded search of telomerase RNA across Insecta]. GERONIMO can test hundreds of genomes and ensures the stability and reproducibility of the analyses performed.\r\n\r\n\r\n[telomerase RNA in lower plants]: https://doi.org/10.1093/nar/gkab545\r\n[expanded search of telomerase RNA across Insecta]: https://doi.org/10.1093/nar/gkac1202\r\n\r\n## Scope\r\nThe GERONIMO tool utilises covariance models (CMs) to conduct homology searches of RNA sequences across a wide range of gene families in a broad evolutionary context. Specifically, it can be utilised to:\r\n\r\n* Detect RNA sequences that share a common evolutionary ancestor\r\n* Identify and align orthologous RNA sequences among closely related species, as well as paralogous sequences within a single species\r\n* Identify conserved non-coding RNAs in a genome, and extract upstream genomic regions to characterise potential promoter regions. \r\nIt is important to note that GERONIMO is a computational tool, and as such, it is intended to be run on a computer with a small amount of data. Appropriate computational infrastructure is necessary for analysing hundreds of genomes.\r\n\r\nAlthough GERONIMO was primarily designed for Telomerase RNA identification, its functionality extends to include the detection and alignment of other RNA gene families, including **rRNA**, **tRNA**, **snRNA**, **miRNA**, and **lncRNA**. This can aid in identifying paralogs and orthologs across different species that may carry specific functions, making it useful for phylogenetic analyses. \r\n\r\nIt is crucial to remember that some gene families may exhibit similar characteristics but different functions. Therefore, analysing the data and functional annotation after conducting the search is essential to characterise the sequences properly.\r\n\r\n## Pipeline overview\r\n\r\n\r\nBy default, the GERONIMO pipeline conducts high-throughput searches of homology sequences in downloaded genomes utilizing covariance models. If a significant similarity is detected between the model and genome sequence, the pipeline extracts the upstream region, making it convenient to identify the promoter of the discovered gene. In brief, the pipeline:\r\n- Compiles a list of genomes using the NCBI's [Entrez] database based on a specified query, *e.g. \"Rhodophyta\"[Organism]*\r\n- Downloads and decompresses the requested genomes using *rsync* and *gunzip*, respectively\r\n- *Optionally*, generates a covariance model based on a provided alignment using [Infernal]\r\n- Conducts searches among the genomes using the covariance model [Infernal]\r\n- Supplements genome information with taxonomy data using [rentrez]\r\n- Expands the significant hits sequence by extracting upstream genomic regions using [*blastcmd*]\r\n- Compiles the results, organizes them into a tabular format, and generates a visual summary of the performed analysis.\r\n\r\n[Entrez]: https://www.ncbi.nlm.nih.gov/books/NBK179288/\r\n[Infernal]: http://eddylab.org/infernal/\r\n[rentrez]: https://github.com/ropensci/rentrez\r\n[*blastcmd*]: https://www.ncbi.nlm.nih.gov/books/NBK569853/\r\n\r\n## Quick start\r\nThe GERONIMO is available as a `snakemake pipeline` running on Linux and Windows operating systems.\r\n\r\n### Windows 10\r\nInstal Linux on Windows 10 (WSL) according to [instructions], which bottling down to opening PowerShell or Windows Command Prompt in *administrator mode* and pasting the following:\r\n```shell\r\nwsl --install\r\nwsl.exe --install UBUNTU\r\n```\r\nThen restart the machine and follow the instructions for setting up the Linux environment.\r\n\r\n[instructions]: https://learn.microsoft.com/en-us/windows/wsl/install\r\n\r\n### Linux:\r\n#### Check whether the conda is installed:\r\n```shell\r\nconda -V\r\n```\r\n> GERONIMO was tested on conda 23.3.1\r\n#### 1) If you do not have installed `conda`, please install `miniconda`\r\nPlease follow the instructions for installing [miniconda]\r\n\r\n[miniconda]: https://conda.io/projects/conda/en/stable/user-guide/install/linux.html\r\n\r\n#### 2) Continue with installing `mamba` (recommended but optional)\r\n```shell\r\nconda install -n base -c conda-forge mamba\r\n```\r\n#### 3) Install `snakemake`\r\n```shell\r\nconda activate base\r\nmamba create -p env_snakemake -c conda-forge -c bioconda snakemake\r\nmamba activate env_snakemake\r\nsnakemake --help\r\n```\r\nIn case of complications, please check the section `Questions & Answers` below or follow the [official documentation] for troubleshooting.\r\n\r\n[official documentation]: https://snakemake.readthedocs.io/en/stable/getting_started/installation.html\r\n\r\n### Clone the GERONIMO repository\r\nGo to the path in which you want to run the analysis and clone the repository:\r\n```shell\r\ncd \r\ngit clone https://github.com/amkilar/GERONIMO.git\r\n```\r\n\r\n### Run sample analysis to ensure GERONIMO installation was successful\r\nAll files are prepared for the sample analysis as a default. Please execute the line below:\r\n```shell\r\nsnakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n```\r\n\r\nThis will prompt GERONIMO to quickly scan all modules, verifying the correct setup of the pipeline without executing any analysis.\r\nYou should see the message `Building DAG of jobs...`, followed by `Nothing to be done (all requested files are present and up to date).`, when successfully completed.\r\n\r\nIf you want to run the sample analysis fully, please remove the folder `results` from the GERONIMO directory and execute GERONIMO again with:\r\n\r\n`snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx`\r\n\r\n> You might consider allowing more cores to speed up the analysis, which might take up to several hours.\r\n\r\n#### You might want to clean `GERONIMO/` directory from the files produced by the example analysis. You can safely remove the following:\r\n- `GERONIMO/results`\r\n- `GERONIMO/database`\r\n- `GERONIMO/taxonomy`\r\n- `GERONIMO/temp`\r\n- `.create_genome_list.touch`\r\n- `list_of_genomes.txt`\r\n\r\n## Setup the inputs\r\n\r\n### 1) Prepare the `covariance models`:\r\n\r\n#### Browse the collection of available `covariance models` at [Rfam] (*You can find the covariance model in the tab `Curation`.*) \r\nPaste the covariance model to the folder `GERONIMO/models` and ensure its name follows the convention: `cov_model_`\r\n\r\n[Rfam]: https://rfam.org/\r\n\r\n#### **OR**\r\n\r\n#### Prepare your own `covariance model` using [LocARNA]\r\n1. Paste or upload your sequences to the web server and download the `.stk` file with the alignment result. \r\n \r\n > *Please note that the `.stk` file format is crucial for the analysis, containing sequence alignment and secondary structure consensus.*\r\n \r\n > The LocARNA web service allows you to align 30 sequences at once - if you need to align more sequences, please use the standalone version available [here] \r\n > After installation run: \r\n ```shell\r\n mlocarna my_fasta_sequences.fasta\r\n ```\r\n \r\n2. Paste the `.stk` alignment file to the folder `GERONIMO/model_to_build` and ensure its name follows the convention: `.stk`\r\n\r\n > Please check the example `heterotrichea.stk` format in `GERONIMO/models_to_built` for reference\r\n \r\n\r\n[LocARNA]: http://rna.informatik.uni-freiburg.de/LocARNA/Input.jsp\r\n[here]: http://www.bioinf.uni-freiburg.de/Software/LocARNA/\r\n\r\n\r\n### 2) Adjust the `config.yaml` file\r\nPlease adjust the analysis specifications, as in the following example:\r\n\r\n> - database: ' [Organism]' (in case of difficulties with defining the database query, please follow the instructions below)\r\n> - extract_genomic_region-length: (here you can determine how long the upstream genomic region should be extracted; tested for 200)\r\n> - models: [\"\", \"\"] (here specify the names of models that should be used to perform analysis)\r\n> \r\n> *Here you can also insert the name of the covariance model you want to build with GERONIMO - just be sure you placed `.stk` file in `GERONIMO/models_to_build` before starting analysis*\r\n> - CPU_for_model_building: (specify the number of available CPUs devoted to the process of building model (cannot exceed the CPU number allowed to snakemake with `--cores`)\r\n>\r\n> *You might ignore this parameter when you do not need to create a new covariance model*\r\n\r\n\r\nKeep in mind that the covariance models and alignments must be present in the respective GERONIMO folders.\r\n \r\n### 3) Remove folder `results`, which contains example analysis output\r\n### 4) **Please ensure you have enough storage capacity to download all the requested genomes (in the `GERONIMO/` directory)**\r\n\r\n## Run GERONIMO\r\n```shell\r\nmamba activate env_snakemake\r\ncd ~/GERONIMO\r\nsnakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n```\r\n \r\n## Example results\r\n\r\n### Outputs characterisation\r\n\r\n#### A) Summary table\r\nThe Excel table contains the results arranged by taxonomy information and hit significance. The specific columns include:\r\n* family, organism_name, class, order, phylum (taxonomy context)\r\n* GCA_id - corresponds to the genome assembly in the *NCBI database*\r\n* model - describes which covariance model identified the result\r\n* label - follows the *Infernal* convention of categorizing hits\r\n* number - the counter of the result\r\n* e_value - indicates the significance level of the hit\r\n* HIT_sequence - the exact HIT sequence found by *Infernal*, which corresponds to the covariance model\r\n* HIT_ID - describes in which part of the genome assembly the hit was found, which may help publish novel sequences\r\n* extended_genomic_region - upstream sequence, which may contain a possible promoter sequence\r\n* secondary_structure - the secondary structure consensus of the covariance model\r\n\r\n\r\n#### B) Significant Hits Distribution Across Taxonomy Families\r\nThe plot provides an overview of the number of genomes in which at least one significant hit was identified, grouped by family. The bold black line corresponds to the number of genomes present in each family, helping to minimize bias regarding unequal data representation across the taxonomy.\r\n\r\n\r\n#### C) Hits Distribution in Genomes Across Families\r\nThe heatmap provides information about the most significant hits from the genome, identified by a specific covariance model. Genomes are grouped by families (on the right). Hits are classified into three categories based on their e-values. Generally, these categories correspond to hit classifications (\"HIT,\" \"MAYBE,\" \"NO HIT\"). The \"HIT\" category is further divided to distinguish between highly significant hits and moderately significant ones.\r\n\r\n\r\n\r\n### GERONIMO directory structure\r\n\r\nThe GERONIMO directory structure is designed to produce files in a highly structured manner, ensuring clear insight and facilitating the analysis of results. During a successful run, GERONIMO produces the following folders:\r\n* `/database` - which contains genome assemblies that were downloaded from the *NCBI database* and grouped in subfolders\r\n* `/taxonomy` - where taxonomy information is gathered and stored in the form of tables\r\n* `/results` - the main folder containing all produced results:\r\n * `/infernal_raw` - contains the raw results produced by *Infernal*\r\n * `/infernal` - contains restructured results of *Infernal* in table format\r\n * `/cmdBLAST` - contains results of *cmdblast*, which extracts the extended genomic region\r\n * `/summary` - contains summary files that join results from *Infernal*, *cmdblast*, and attach taxonomy context\r\n * `/plots` - contains two types of summary plots\r\n* `/temp` - folder contains the information necessary to download genome assemblies from *NCBI database*\r\n\r\n* `/env` - stores instructions for dependency installation\r\n* `/models` - where calibrated covariance models can be pasted, *for example, from the Rfam database*\r\n* `/modes_to_built` - where multiple alignments in *.stk* format can be pasted\r\n* `/scripts` - contains developed scripts that perform results structurization\r\n\r\n#### The example GERONIMO directory structure:\r\n\r\n```shell\r\nGERONIMO\r\n├── database\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   └── ...\r\n├── env\r\n├── models\r\n├── model_to_build\r\n├── results\r\n│   ├── cmdBLAST\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── ...\r\n│   ├── infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   ├── plots\r\n│   ├── raw_infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   └── summary\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   ├── ...\r\n├── scripts\r\n├── taxonomy\r\n└── temp\r\n```\r\n\r\n## GERONIMO applicability\r\n\r\n### Expanding the evolutionary context\r\nTo add new genomes or database queries to an existing analysis, please follow the instructions:\r\n1) Rename the `list_of_genomes.txt` file to `previous_list_of_genomes.txt` or any other preferred name.\r\n2) Modify the `config.yaml` file by replacing the previous database query with the new one.\r\n3) Delete:\r\n - `summary_table.xlsx`, `part_summary_table.csv`, `summary_table_models.xlsx` files located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n5) Run GERONIMO to calculate new results using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n ```\r\n7) Once the new results are generated, reviewing them before merging them with the original results is recommended.\r\n8) Copy the contents of the `previous_list_of_genomes.txt` file and paste them into the current `list_of_genomes.txt`.\r\n9) Delete:\r\n - `summary_table.xlsx` located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n10) Run GERONIMO to merge the results from both analyses using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n ```\r\n\r\n### Incorporating new covariance models into existing analysis\r\n1) Copy the new covariance model to `GERONIMO/models`\r\n2) Modify the `config.yaml` file by adding the name of the new model to the line `models: [...]`\r\n3) Run GERONIMO to see the updated analysis outcome\r\n\r\n### Building a new covariance model\r\nWith GERONIMO, building a new covariance model from multiple sequence alignment in the `.stk` format is possible. \r\n\r\nTo do so, simply paste `.stk` file to `GERONIMO/models_to_build` and paste the name of the new covariance model to `config.yaml` file to the line `models: [\"\"]`\r\n\r\nand run GERONIMO.\r\n\r\n\r\n## Questions & Answers\r\n\r\n### How to specify the database query?\r\n- Visit the [NCBI Assemblies] website. \r\n- Follow the instruction on the graphic below:\r\n\r\n[NCBI Assemblies]: https://www.ncbi.nlm.nih.gov/assembly/?term=\r\n\r\n### WSL: problem with creating `snakemake_env`\r\nIn the case of an error similar to the one below:\r\n> CondaError: Unable to create prefix directory '/mnt/c/Windows/system32/env_snakemake'.\r\n> Check that you have sufficient permissions. \r\n \r\nYou might try to delete the cache with: `rm -r ~/.cache/` and try again.\r\n\r\n### When `snakemake` does not seem to be installed properly\r\nIn the case of the following error:\r\n> Command 'snakemake' not found ...\r\n\r\nCheck whether the `env_snakemake` is activated.\r\n> It should result in a change from (base) to (env_snakemake) before your login name in the command line window.\r\n\r\nIf you still see `(base)` before your login name, please try to activate the environment with conda:\r\n`conda activate env_snakemake`\r\n\r\n\r\nPlease note that you might need to specify the full path to the `env_snakemake`, like /home/your user name/env_snakemake\r\n\r\n### How to browse GERONIMO results obtained in WSL?\r\nYou can easily access the results obtained on WSL from your Windows environment by opening `File Explorer` and pasting the following line into the search bar: `\\\\wsl.localhost\\Ubuntu\\home\\`. This will reveal a folder with your username, as specified during the configuration of your Ubuntu system. To locate the GERONIMO results, simply navigate to the folder with your username and then to the `home` folder. (`\\\\wsl.localhost\\Ubuntu\\home\\\\home\\GERONIMO`)\r\n\r\n### GERONIMO occupies a lot of storage space\r\nThrough genome downloads, GERONIMO can potentially consume storage space, rapidly leading to a shortage. Currently, downloading genomes is an essential step for optimal GERONIMO performance.\r\n\r\nRegrettably, if the analysis is rerun without the `/database` folder, it will result in the need to redownload genomes, which is a highly time-consuming process.\r\n\r\nNevertheless, if you do not intend to repeat the analysis and have no requirement for additional genomes or models, you are welcome to retain your results tables and plots while removing the remaining files.\r\n\r\nIt is strongly advised against using local machines for extensive analyses. If you lack access to external storage space, it is recommended to divide the analysis into smaller segments, which can be later merged, as explained in the section titled `Expanding the evolutionary context`.\r\n\r\nConsidering this limitation, I am currently working on implementing a solution that will help circumvent the need for redundant genome downloads without compromising GERONIMO performance in the future.\r\n\r\nYou might consider deleting the `.snakemake` folder to free up storage space. However, please note that deleting this folder will require the reinstallation of GERONIMO dependencies when the analysis is rerun.\r\n\r\n## License\r\nCopyright (c) 2023 Agata M. Kilar\r\n\r\nPermission is hereby granted, free of charge, to any person obtaining a copy\r\nof this software and associated documentation files (the \"Software\"), to deal\r\nin the Software without restriction, including without limitation the rights\r\nto use, copy, modify, merge, publish, distribute, sublicense, and/or sell\r\ncopies of the Software, and to permit persons to whom the Software is\r\nfurnished to do so, subject to the following conditions:\r\n\r\nThe above copyright notice and this permission notice shall be included in all\r\ncopies or substantial portions of the Software.\r\n\r\nTHE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR\r\nIMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,\r\nFITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE\r\nAUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER\r\nLIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,\r\nOUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE\r\nSOFTWARE.\r\n\r\n## Contact\r\nmgr inż. Agata Magdalena Kilar, PhD (agata.kilar@gmail.com)" . "application/ld+json" . . . . . . . . . . . "https://w3id.org/ro-id/892d05bf-da84-4be9-9ab7-60eba5f73c70" . "https://github.com/amkilar/GERONIMO.git" . . "GERONIMO.sm" . "Research Object Crate for GERONIMO" . "https://workflowhub.eu/workflows/547/ro_crate?version=1" . . . . "https://w3id.org/ro-id/1c684234-851e-4998-a0d2-41c66b7a20d5" . "https://w3id.org/ro-id/a15f9eff-f256-4663-87fd-6ab855cd5941" . "https://w3id.org/ro-id/ad9304c3-ccf5-4d0f-8ca0-9f4da0cab28e" . "https://w3id.org/ro-id/005c8422-6fce-4a35-862d-7f2ee4d9c4f5" . "https://w3id.org/ro-id/0f678463-a81a-47a0-a669-079e20fa0125" . "https://w3id.org/ro-id/54c71f37-c543-46a8-84b3-aaa8ae32fb6f" . "https://w3id.org/ro-id/58d0583c-0a1a-4ce8-83ee-3620189c1ddd" . "https://w3id.org/ro-id/7bffe30b-5a6e-4262-80b0-701944edc9f6" . "https://w3id.org/ro-id/7f74e75a-8601-49da-a852-353b1a3994f2" . "https://w3id.org/ro-id/82d5b0b6-0398-46cd-9cc2-8af27f478fc0" . "https://w3id.org/ro-id/83fc731e-3aa1-43c7-9cf5-ab6740c3633d" . "https://w3id.org/ro-id/8ac41430-786a-430f-abe4-37d599cdae43" . "https://w3id.org/ro-id/935b0973-3d62-49ac-9c5a-84cc5b147c8b" . "https://w3id.org/ro-id/9e31b1dd-d3ba-4d69-ad80-5549179acdc5" . "https://w3id.org/ro-id/a028a5da-e9aa-4ae3-b2d4-0ffaa42903f5" . "https://w3id.org/ro-id/acf35671-dadd-4cf7-9241-ed095c2c7de3" . "https://w3id.org/ro-id/af381567-5178-4371-a371-df7607665fc5" . "https://w3id.org/ro-id/b0e0a997-66e3-4631-ae36-b0b54cdee3d9" . "https://w3id.org/ro-id/ba2815fb-b1e0-4bd6-8934-8040df63fe9f" . "https://w3id.org/ro-id/c45ec5d3-c9ab-43ac-857d-f1acd9382c63" . "https://w3id.org/ro-id/c88a7dc6-d3d4-43a4-9de6-b0d62237c95f" . "https://w3id.org/ro-id/d548f2ff-c42c-4cba-b8fe-feee2ff0fb11" . "https://w3id.org/ro-id/3f7be379-7d08-4ee4-b1b9-aa17032ec971" . "https://w3id.org/ro-id/d4da8bc8-008c-4385-9640-9bed35f86494" . "https://w3id.org/ro-id/4fae7e6b-9f40-4a7c-b2aa-e85c1a97e8ec" . "https://w3id.org/ro-id/612078aa-2703-4c7b-843e-0b683ce23f9c" . "https://w3id.org/ro-id/73687a3c-6aaa-429a-9d13-7d4b854f803e" . "https://w3id.org/ro-id/15a63617-f9a8-4753-8830-af3a1bac7b4d" . "https://w3id.org/ro-id/17b944ac-f8bb-403d-bddf-98d90bd13b83" . "https://w3id.org/ro-id/28bfade8-878b-4c72-8135-4078fbd09e42" . "https://w3id.org/ro-id/359d5480-f589-47f7-aaae-ff4c8aeb9a7a" . "https://w3id.org/ro-id/3f7102cd-996b-4e5b-b498-525210a18b3a" . "https://w3id.org/ro-id/4b31f4de-1524-4cef-b8b3-72a240d63dd0" . "https://w3id.org/ro-id/5d9ce283-c305-46ae-99e6-e81b9acf9018" . "https://w3id.org/ro-id/9d5f1beb-ded6-4482-9c66-179de5feff0b" . "https://w3id.org/ro-id/a48c1ce5-f60b-461b-b1a9-249075c973e8" . "https://w3id.org/ro-id/bb5ac70a-e6aa-4213-95d8-a7cf4b58f3f6" . "https://w3id.org/ro-id/c32b7cb6-cdf9-45f8-9acd-64e8609d64f0" . "https://w3id.org/ro-id/d9478a50-7977-42f2-afc6-f9f5016caa29" . "https://w3id.org/ro-id/e4e53726-b83b-4025-ac62-990609bc13db" . "https://w3id.org/ro-id/5b68f8b8-d94d-44ed-b47c-9b30492dd501" . "https://w3id.org/ro-id/cceabc2e-5232-4060-8697-5fcaa87b6f18" . "https://w3id.org/ro-id/05648ad5-fcc8-4f24-882a-0bf77f727e52" . "https://w3id.org/ro-id/0735a772-38d8-45e7-a3d2-e20f62056146" . "https://w3id.org/ro-id/285fc200-65ce-4175-afc6-da0c325375f0" . "https://w3id.org/ro-id/2c76f1de-2578-4660-8926-16dbd0f88985" . "https://w3id.org/ro-id/3a8d21f9-3481-47d7-925e-2c129cf76f8d" . "https://w3id.org/ro-id/5f76202b-8605-4866-9750-2ae44e74cf74" . "https://w3id.org/ro-id/630cd87c-2fd2-4a70-bbad-e85002d3cc15" . "https://w3id.org/ro-id/7589da61-b04f-425f-83c1-6dbbedd2d27a" . "https://w3id.org/ro-id/76e5ce8b-2407-4d11-adba-17e295f4b1bd" . "https://w3id.org/ro-id/869e48d4-c589-4db1-af1e-0a2a00da300b" . "https://w3id.org/ro-id/8cff22b7-ba04-4331-9c66-9839b47eaa32" . "https://w3id.org/ro-id/a20d1313-4ed3-40e3-b3a2-410652e82081" . "https://w3id.org/ro-id/a5ffdb9a-54c2-4a38-946c-77190ed34431" . "https://w3id.org/ro-id/b0a6cec6-4527-4450-8012-5fa257d2ab77" . "https://w3id.org/ro-id/b1d6ebd1-ab77-4337-b253-7b08f406c27e" . "https://w3id.org/ro-id/b591df3e-7def-4e90-a20b-28fcbbff4fc2" . "https://w3id.org/ro-id/c86349a6-01f9-4c71-a872-ce3527630004" . "https://w3id.org/ro-id/cde37df3-f892-4c6e-89a0-f8b532fde95d" . "https://w3id.org/ro-id/d2685e00-63e1-4365-b375-92539c6893ec" . "https://w3id.org/ro-id/14d92214-4a88-46fb-b596-4f2e99bfcc37" . "https://w3id.org/ro-id/87d45696-5e29-461e-b6bf-bc8fcc2a033f" . "https://w3id.org/ro-id/8865ad4d-7c4a-46a5-a12a-89213e7e14bd" . "https://w3id.org/ro-id/c6632496-4497-4e92-869f-1072a707de02" . "Agata Kilar. \"Research Object Crate for GERONIMO.\" ROHub. Sep 08 ,2023. https://w3id.org/ro-id/892d05bf-da84-4be9-9ab7-60eba5f73c70." . . . . . . . . "env" . . . . "models_to_build" . . . . . . . . . . . "models" . . . . . "plots" . . . . . . . . . . "scripts" . . . . . . "results" . . . . "92192"^^ . "https://api.rohub.org/api/resources/04fd81d6-ef3e-4e2f-8e50-ae6a40144ea2/download/" . . "2023-09-08 12:04:05.466271+00:00" . "2023-09-08 12:04:10.301548+00:00" . . "cov_model_U4" . "2023-09-08 12:04:05.466271+00:00" . . . . "111328"^^ . "https://api.rohub.org/api/resources/07ffb5a3-347e-4f5d-883a-78f016ee907a/download/" . . "2023-09-08 12:04:05.460885+00:00" . "2023-09-08 12:04:09.341000+00:00" . . "cov_model_U1" . "2023-09-08 12:04:05.460885+00:00" . . . "https://ror.org/https://doi.org/10.48546/workflowhub.workflow.547.1" . . "72017"^^ . "https://api.rohub.org/api/resources/083e698f-333d-4c5d-9610-cc48cdad9e4d/download/" . . "2023-09-08 12:04:05.481914+00:00" . "2023-09-08 12:05:20.681373+00:00" . "text/html" . . "ro-crate-preview.html" . "2023-09-08 12:04:05.481914+00:00" . . . . . "657"^^ . "https://api.rohub.org/api/resources/1037d046-1068-4888-a38b-ba2ef9ce42eb/download/" . . "2023-09-08 12:04:05.472015+00:00" . "2023-09-08 12:05:12.815668+00:00" . . "infernal_env.yaml" . "2023-09-08 12:04:05.472015+00:00" . . . . "1071"^^ . "https://api.rohub.org/api/resources/17bc9250-5ec3-4696-a1e2-2d8fbfab6736/download/" . . "2023-09-08 12:04:05.440532+00:00" . "2023-09-08 12:04:06.435874+00:00" . . "LICENSE" . "2023-09-08 12:04:05.440532+00:00" . . . . "21805" . "https://api.rohub.org/api/resources/1c28c455-7dac-4035-8898-511615255bcb/download/" . . "2023-09-08 12:04:05.472901+00:00" . "2023-09-08 12:05:20.104078+00:00" . "text/markdown" . . "README.md" . "#cwl" . "2023-09-08 12:04:05.472901+00:00" . . . . . . "159869"^^ . "https://api.rohub.org/api/resources/200efa59-3ae7-4970-aea3-e22645d28d96/download/" . . "2023-09-08 12:04:05.462474+00:00" . "2023-09-08 12:04:09.540442+00:00" . . "cov_model_SRP" . "2023-09-08 12:04:05.462474+00:00" . . . . "6681"^^ . "https://api.rohub.org/api/resources/2062fe32-a75a-4957-a90b-8dafc5ad0367/download/" . . "2023-09-08 12:04:05.470567+00:00" . "2023-09-08 12:05:12.326066+00:00" . . "make_summary_R.yaml" . "2023-09-08 12:04:05.470567+00:00" . . . . "4883"^^ . "https://api.rohub.org/api/resources/21369f0b-77af-4840-8160-0aaab2e6b15f/download/" . . "2023-09-08 12:04:05.477109+00:00" . "2023-09-08 12:05:14.349759+00:00" . . "search_taxonomy.r" . "2023-09-08 12:04:05.477109+00:00" . . . . "2262"^^ . "https://api.rohub.org/api/resources/215f4518-6bf1-41f1-b1f3-72f45c0f387d/download/" . . "2023-09-08 12:04:05.477887+00:00" . "2023-09-08 12:05:14.755680+00:00" . . "transform_cmdBLAST_output.R" . "2023-09-08 12:04:05.477887+00:00" . . . . "213819"^^ . "https://api.rohub.org/api/resources/229330c4-3395-4872-8fe5-94eae1b2b2ed/download/" . . "2023-09-08 12:04:05.454424+00:00" . "2023-09-08 12:04:07.339998+00:00" . "application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" . . "summary_table.xlsx" . "2023-09-08 12:04:05.454424+00:00" . . . . "68255"^^ . "https://api.rohub.org/api/resources/280bae50-1fa1-45c8-ae1b-17fe40b9f300/download/" . . "2023-09-08 12:04:05.458093+00:00" . "2023-09-08 12:04:08.198257+00:00" . . "cov_model_U6" . "2023-09-08 12:04:05.458093+00:00" . . . . "146892"^^ . "https://api.rohub.org/api/resources/3008db58-7827-4bf7-b378-3257e455e597/download/" . . "2023-09-08 12:04:05.459620+00:00" . "2023-09-08 12:04:08.394076+00:00" . . "cov_model_U3" . "2023-09-08 12:04:05.459620+00:00" . . . "https://bioschemas.org/profiles/ComputationalWorkflow/1.0-RELEASE/" . . "6740" . "https://api.rohub.org/api/resources/34d6227c-e493-4c4d-a841-67685c68be6d/download/" . . "2023-08-01 01:34:42+00:00" . "2023-09-08 12:05:20.296455+00:00" . "# GERONIMO\r\n\r\n## Introduction\r\nGERONIMO is a bioinformatics pipeline designed to conduct high-throughput homology searches of structural genes using covariance models. These models are based on the alignment of sequences and the consensus of secondary structures. The pipeline is built using Snakemake, a workflow management tool that allows for the reproducible execution of analyses on various computational platforms. \r\n\r\nThe idea for developing GERONIMO emerged from a comprehensive search for [telomerase RNA in lower plants] and was subsequently refined through an [expanded search of telomerase RNA across Insecta]. GERONIMO can test hundreds of genomes and ensures the stability and reproducibility of the analyses performed.\r\n\r\n\r\n[telomerase RNA in lower plants]: https://doi.org/10.1093/nar/gkab545\r\n[expanded search of telomerase RNA across Insecta]: https://doi.org/10.1093/nar/gkac1202\r\n\r\n## Scope\r\nThe GERONIMO tool utilises covariance models (CMs) to conduct homology searches of RNA sequences across a wide range of gene families in a broad evolutionary context. Specifically, it can be utilised to:\r\n\r\n* Detect RNA sequences that share a common evolutionary ancestor\r\n* Identify and align orthologous RNA sequences among closely related species, as well as paralogous sequences within a single species\r\n* Identify conserved non-coding RNAs in a genome, and extract upstream genomic regions to characterise potential promoter regions. \r\nIt is important to note that GERONIMO is a computational tool, and as such, it is intended to be run on a computer with a small amount of data. Appropriate computational infrastructure is necessary for analysing hundreds of genomes.\r\n\r\nAlthough GERONIMO was primarily designed for Telomerase RNA identification, its functionality extends to include the detection and alignment of other RNA gene families, including **rRNA**, **tRNA**, **snRNA**, **miRNA**, and **lncRNA**. This can aid in identifying paralogs and orthologs across different species that may carry specific functions, making it useful for phylogenetic analyses. \r\n\r\nIt is crucial to remember that some gene families may exhibit similar characteristics but different functions. Therefore, analysing the data and functional annotation after conducting the search is essential to characterise the sequences properly.\r\n\r\n## Pipeline overview\r\n\r\n\r\nBy default, the GERONIMO pipeline conducts high-throughput searches of homology sequences in downloaded genomes utilizing covariance models. If a significant similarity is detected between the model and genome sequence, the pipeline extracts the upstream region, making it convenient to identify the promoter of the discovered gene. In brief, the pipeline:\r\n- Compiles a list of genomes using the NCBI's [Entrez] database based on a specified query, *e.g. \"Rhodophyta\"[Organism]*\r\n- Downloads and decompresses the requested genomes using *rsync* and *gunzip*, respectively\r\n- *Optionally*, generates a covariance model based on a provided alignment using [Infernal]\r\n- Conducts searches among the genomes using the covariance model [Infernal]\r\n- Supplements genome information with taxonomy data using [rentrez]\r\n- Expands the significant hits sequence by extracting upstream genomic regions using [*blastcmd*]\r\n- Compiles the results, organizes them into a tabular format, and generates a visual summary of the performed analysis.\r\n\r\n[Entrez]: https://www.ncbi.nlm.nih.gov/books/NBK179288/\r\n[Infernal]: http://eddylab.org/infernal/\r\n[rentrez]: https://github.com/ropensci/rentrez\r\n[*blastcmd*]: https://www.ncbi.nlm.nih.gov/books/NBK569853/\r\n\r\n## Quick start\r\nThe GERONIMO is available as a `snakemake pipeline` running on Linux and Windows operating systems.\r\n\r\n### Windows 10\r\nInstal Linux on Windows 10 (WSL) according to [instructions], which bottling down to opening PowerShell or Windows Command Prompt in *administrator mode* and pasting the following:\r\n```shell\r\nwsl --install\r\nwsl.exe --install UBUNTU\r\n```\r\nThen restart the machine and follow the instructions for setting up the Linux environment.\r\n\r\n[instructions]: https://learn.microsoft.com/en-us/windows/wsl/install\r\n\r\n### Linux:\r\n#### Check whether the conda is installed:\r\n```shell\r\nconda -V\r\n```\r\n> GERONIMO was tested on conda 23.3.1\r\n#### 1) If you do not have installed `conda`, please install `miniconda`\r\nPlease follow the instructions for installing [miniconda]\r\n\r\n[miniconda]: https://conda.io/projects/conda/en/stable/user-guide/install/linux.html\r\n\r\n#### 2) Continue with installing `mamba` (recommended but optional)\r\n```shell\r\nconda install -n base -c conda-forge mamba\r\n```\r\n#### 3) Install `snakemake`\r\n```shell\r\nconda activate base\r\nmamba create -p env_snakemake -c conda-forge -c bioconda snakemake\r\nmamba activate env_snakemake\r\nsnakemake --help\r\n```\r\nIn case of complications, please check the section `Questions & Answers` below or follow the [official documentation] for troubleshooting.\r\n\r\n[official documentation]: https://snakemake.readthedocs.io/en/stable/getting_started/installation.html\r\n\r\n### Clone the GERONIMO repository\r\nGo to the path in which you want to run the analysis and clone the repository:\r\n```shell\r\ncd \r\ngit clone https://github.com/amkilar/GERONIMO.git\r\n```\r\n\r\n### Run sample analysis to ensure GERONIMO installation was successful\r\nAll files are prepared for the sample analysis as a default. Please execute the line below:\r\n```shell\r\nsnakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n```\r\n\r\nThis will prompt GERONIMO to quickly scan all modules, verifying the correct setup of the pipeline without executing any analysis.\r\nYou should see the message `Building DAG of jobs...`, followed by `Nothing to be done (all requested files are present and up to date).`, when successfully completed.\r\n\r\nIf you want to run the sample analysis fully, please remove the folder `results` from the GERONIMO directory and execute GERONIMO again with:\r\n\r\n`snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx`\r\n\r\n> You might consider allowing more cores to speed up the analysis, which might take up to several hours.\r\n\r\n#### You might want to clean `GERONIMO/` directory from the files produced by the example analysis. You can safely remove the following:\r\n- `GERONIMO/results`\r\n- `GERONIMO/database`\r\n- `GERONIMO/taxonomy`\r\n- `GERONIMO/temp`\r\n- `.create_genome_list.touch`\r\n- `list_of_genomes.txt`\r\n\r\n## Setup the inputs\r\n\r\n### 1) Prepare the `covariance models`:\r\n\r\n#### Browse the collection of available `covariance models` at [Rfam] (*You can find the covariance model in the tab `Curation`.*) \r\nPaste the covariance model to the folder `GERONIMO/models` and ensure its name follows the convention: `cov_model_`\r\n\r\n[Rfam]: https://rfam.org/\r\n\r\n#### **OR**\r\n\r\n#### Prepare your own `covariance model` using [LocARNA]\r\n1. Paste or upload your sequences to the web server and download the `.stk` file with the alignment result. \r\n \r\n > *Please note that the `.stk` file format is crucial for the analysis, containing sequence alignment and secondary structure consensus.*\r\n \r\n > The LocARNA web service allows you to align 30 sequences at once - if you need to align more sequences, please use the standalone version available [here] \r\n > After installation run: \r\n ```shell\r\n mlocarna my_fasta_sequences.fasta\r\n ```\r\n \r\n2. Paste the `.stk` alignment file to the folder `GERONIMO/model_to_build` and ensure its name follows the convention: `.stk`\r\n\r\n > Please check the example `heterotrichea.stk` format in `GERONIMO/models_to_built` for reference\r\n \r\n\r\n[LocARNA]: http://rna.informatik.uni-freiburg.de/LocARNA/Input.jsp\r\n[here]: http://www.bioinf.uni-freiburg.de/Software/LocARNA/\r\n\r\n\r\n### 2) Adjust the `config.yaml` file\r\nPlease adjust the analysis specifications, as in the following example:\r\n\r\n> - database: ' [Organism]' (in case of difficulties with defining the database query, please follow the instructions below)\r\n> - extract_genomic_region-length: (here you can determine how long the upstream genomic region should be extracted; tested for 200)\r\n> - models: [\"\", \"\"] (here specify the names of models that should be used to perform analysis)\r\n> \r\n> *Here you can also insert the name of the covariance model you want to build with GERONIMO - just be sure you placed `.stk` file in `GERONIMO/models_to_build` before starting analysis*\r\n> - CPU_for_model_building: (specify the number of available CPUs devoted to the process of building model (cannot exceed the CPU number allowed to snakemake with `--cores`)\r\n>\r\n> *You might ignore this parameter when you do not need to create a new covariance model*\r\n\r\n\r\nKeep in mind that the covariance models and alignments must be present in the respective GERONIMO folders.\r\n \r\n### 3) Remove folder `results`, which contains example analysis output\r\n### 4) **Please ensure you have enough storage capacity to download all the requested genomes (in the `GERONIMO/` directory)**\r\n\r\n## Run GERONIMO\r\n```shell\r\nmamba activate env_snakemake\r\ncd ~/GERONIMO\r\nsnakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n```\r\n \r\n## Example results\r\n\r\n### Outputs characterisation\r\n\r\n#### A) Summary table\r\nThe Excel table contains the results arranged by taxonomy information and hit significance. The specific columns include:\r\n* family, organism_name, class, order, phylum (taxonomy context)\r\n* GCA_id - corresponds to the genome assembly in the *NCBI database*\r\n* model - describes which covariance model identified the result\r\n* label - follows the *Infernal* convention of categorizing hits\r\n* number - the counter of the result\r\n* e_value - indicates the significance level of the hit\r\n* HIT_sequence - the exact HIT sequence found by *Infernal*, which corresponds to the covariance model\r\n* HIT_ID - describes in which part of the genome assembly the hit was found, which may help publish novel sequences\r\n* extended_genomic_region - upstream sequence, which may contain a possible promoter sequence\r\n* secondary_structure - the secondary structure consensus of the covariance model\r\n\r\n\r\n#### B) Significant Hits Distribution Across Taxonomy Families\r\nThe plot provides an overview of the number of genomes in which at least one significant hit was identified, grouped by family. The bold black line corresponds to the number of genomes present in each family, helping to minimize bias regarding unequal data representation across the taxonomy.\r\n\r\n\r\n#### C) Hits Distribution in Genomes Across Families\r\nThe heatmap provides information about the most significant hits from the genome, identified by a specific covariance model. Genomes are grouped by families (on the right). Hits are classified into three categories based on their e-values. Generally, these categories correspond to hit classifications (\"HIT,\" \"MAYBE,\" \"NO HIT\"). The \"HIT\" category is further divided to distinguish between highly significant hits and moderately significant ones.\r\n\r\n\r\n\r\n### GERONIMO directory structure\r\n\r\nThe GERONIMO directory structure is designed to produce files in a highly structured manner, ensuring clear insight and facilitating the analysis of results. During a successful run, GERONIMO produces the following folders:\r\n* `/database` - which contains genome assemblies that were downloaded from the *NCBI database* and grouped in subfolders\r\n* `/taxonomy` - where taxonomy information is gathered and stored in the form of tables\r\n* `/results` - the main folder containing all produced results:\r\n * `/infernal_raw` - contains the raw results produced by *Infernal*\r\n * `/infernal` - contains restructured results of *Infernal* in table format\r\n * `/cmdBLAST` - contains results of *cmdblast*, which extracts the extended genomic region\r\n * `/summary` - contains summary files that join results from *Infernal*, *cmdblast*, and attach taxonomy context\r\n * `/plots` - contains two types of summary plots\r\n* `/temp` - folder contains the information necessary to download genome assemblies from *NCBI database*\r\n\r\n* `/env` - stores instructions for dependency installation\r\n* `/models` - where calibrated covariance models can be pasted, *for example, from the Rfam database*\r\n* `/modes_to_built` - where multiple alignments in *.stk* format can be pasted\r\n* `/scripts` - contains developed scripts that perform results structurization\r\n\r\n#### The example GERONIMO directory structure:\r\n\r\n```shell\r\nGERONIMO\r\n├── database\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   └── ...\r\n├── env\r\n├── models\r\n├── model_to_build\r\n├── results\r\n│   ├── cmdBLAST\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   │   ├── extended\r\n│   │   │   │   └── filtered\r\n│   │   │   └── ...\r\n│   │   ├── ...\r\n│   ├── infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   ├── plots\r\n│   ├── raw_infernal\r\n│   │   ├── MRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   │   ├── SRP\r\n│   │   │   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   │   │   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   │   │   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   │   │   ├── ...\r\n│   └── summary\r\n│   ├── GCA_000091205.1_ASM9120v1_genomic\r\n│   ├── GCA_000341285.1_ASM34128v1_genomic\r\n│   ├── GCA_000350225.2_ASM35022v2_genomic\r\n│   ├── ...\r\n├── scripts\r\n├── taxonomy\r\n└── temp\r\n```\r\n\r\n## GERONIMO applicability\r\n\r\n### Expanding the evolutionary context\r\nTo add new genomes or database queries to an existing analysis, please follow the instructions:\r\n1) Rename the `list_of_genomes.txt` file to `previous_list_of_genomes.txt` or any other preferred name.\r\n2) Modify the `config.yaml` file by replacing the previous database query with the new one.\r\n3) Delete:\r\n - `summary_table.xlsx`, `part_summary_table.csv`, `summary_table_models.xlsx` files located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n5) Run GERONIMO to calculate new results using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores --use-conda results/summary_table.xlsx\r\n ```\r\n7) Once the new results are generated, reviewing them before merging them with the original results is recommended.\r\n8) Copy the contents of the `previous_list_of_genomes.txt` file and paste them into the current `list_of_genomes.txt`.\r\n9) Delete:\r\n - `summary_table.xlsx` located in the `GERONIMO\\results` directory\r\n - `.create_genome_list.touch` file\r\n10) Run GERONIMO to merge the results from both analyses using the command:\r\n ```shell\r\n snakemake -s GERONIMO.sm --cores 1 --use-conda results/summary_table.xlsx\r\n ```\r\n\r\n### Incorporating new covariance models into existing analysis\r\n1) Copy the new covariance model to `GERONIMO/models`\r\n2) Modify the `config.yaml` file by adding the name of the new model to the line `models: [...]`\r\n3) Run GERONIMO to see the updated analysis outcome\r\n\r\n### Building a new covariance model\r\nWith GERONIMO, building a new covariance model from multiple sequence alignment in the `.stk` format is possible. \r\n\r\nTo do so, simply paste `.stk` file to `GERONIMO/models_to_build` and paste the name of the new covariance model to `config.yaml` file to the line `models: [\"\"]`\r\n\r\nand run GERONIMO.\r\n\r\n\r\n## Questions & Answers\r\n\r\n### How to specify the database query?\r\n- Visit the [NCBI Assemblies] website. \r\n- Follow the instruction on the graphic below:\r\n\r\n[NCBI Assemblies]: https://www.ncbi.nlm.nih.gov/assembly/?term=\r\n\r\n### WSL: problem with creating `snakemake_env`\r\nIn the case of an error similar to the one below:\r\n> CondaError: Unable to create prefix directory '/mnt/c/Windows/system32/env_snakemake'.\r\n> Check that you have sufficient permissions. \r\n \r\nYou might try to delete the cache with: `rm -r ~/.cache/` and try again.\r\n\r\n### When `snakemake` does not seem to be installed properly\r\nIn the case of the following error:\r\n> Command 'snakemake' not found ...\r\n\r\nCheck whether the `env_snakemake` is activated.\r\n> It should result in a change from (base) to (env_snakemake) before your login name in the command line window.\r\n\r\nIf you still see `(base)` before your login name, please try to activate the environment with conda:\r\n`conda activate env_snakemake`\r\n\r\n\r\nPlease note that you might need to specify the full path to the `env_snakemake`, like /home/your user name/env_snakemake\r\n\r\n### How to browse GERONIMO results obtained in WSL?\r\nYou can easily access the results obtained on WSL from your Windows environment by opening `File Explorer` and pasting the following line into the search bar: `\\\\wsl.localhost\\Ubuntu\\home\\`. This will reveal a folder with your username, as specified during the configuration of your Ubuntu system. To locate the GERONIMO results, simply navigate to the folder with your username and then to the `home` folder. (`\\\\wsl.localhost\\Ubuntu\\home\\\\home\\GERONIMO`)\r\n\r\n### GERONIMO occupies a lot of storage space\r\nThrough genome downloads, GERONIMO can potentially consume storage space, rapidly leading to a shortage. Currently, downloading genomes is an essential step for optimal GERONIMO performance.\r\n\r\nRegrettably, if the analysis is rerun without the `/database` folder, it will result in the need to redownload genomes, which is a highly time-consuming process.\r\n\r\nNevertheless, if you do not intend to repeat the analysis and have no requirement for additional genomes or models, you are welcome to retain your results tables and plots while removing the remaining files.\r\n\r\nIt is strongly advised against using local machines for extensive analyses. If you lack access to external storage space, it is recommended to divide the analysis into smaller segments, which can be later merged, as explained in the section titled `Expanding the evolutionary context`.\r\n\r\nConsidering this limitation, I am currently working on implementing a solution that will help circumvent the need for redundant genome downloads without compromising GERONIMO performance in the future.\r\n\r\nYou might consider deleting the `.snakemake` folder to free up storage space. However, please note that deleting this folder will require the reinstallation of GERONIMO dependencies when the analysis is rerun.\r\n\r\n## License\r\nCopyright (c) 2023 Agata M. Kilar\r\n\r\nPermission is hereby granted, free of charge, to any person obtaining a copy\r\nof this software and associated documentation files (the \"Software\"), to deal\r\nin the Software without restriction, including without limitation the rights\r\nto use, copy, modify, merge, publish, distribute, sublicense, and/or sell\r\ncopies of the Software, and to permit persons to whom the Software is\r\nfurnished to do so, subject to the following conditions:\r\n\r\nThe above copyright notice and this permission notice shall be included in all\r\ncopies or substantial portions of the Software.\r\n\r\nTHE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR\r\nIMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,\r\nFITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE\r\nAUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER\r\nLIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,\r\nOUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE\r\nSOFTWARE.\r\n\r\n## Contact\r\nmgr inż. Agata Magdalena Kilar, PhD (agata.kilar@gmail.com)" . "README.md" . "Bioinformatics, Snakemake, rna" . . "GERONIMO" . "https://workflowhub.eu/projects/178" . "#snakemake" . "2023-08-01 01:34:42+00:00" . "https://about.workflowhub.eu/" . "README.md" . "https://workflowhub.eu/workflows/547?version=1" . "1" . . . . . . "354983"^^ . "https://api.rohub.org/api/resources/39ab3c54-396c-4bdc-a0da-fdb563e1a32f/download/" . . "2023-09-08 12:04:05.450498+00:00" . "2023-09-08 12:04:06.925272+00:00" . "image/png" . . "Hits_distribution_across_families.png" . "2023-09-08 12:04:05.450498+00:00" . . . . "129276"^^ . "https://api.rohub.org/api/resources/4b7473cf-f431-4b64-b5c2-199d9fd3f328/download/" . . "2023-09-08 12:04:05.463907+00:00" . "2023-09-08 12:04:09.908441+00:00" . . "cov_model_U2" . "2023-09-08 12:04:05.463907+00:00" . . . . "77501"^^ . "https://api.rohub.org/api/resources/6664e5f4-08cb-4e0b-9436-e59a6719b6e6/download/" . . "2023-09-08 12:04:05.465058+00:00" . "2023-09-08 12:04:10.101003+00:00" . . "cov_model_U5" . "2023-09-08 12:04:05.465058+00:00" . . . . "8159"^^ . "https://api.rohub.org/api/resources/6c59aa5d-1ba7-4bd1-8be4-ff050e3ead30/download/" . . "2023-09-08 12:04:05.469780+00:00" . "2023-09-08 12:05:11.784941+00:00" . . "search_taxonomy_r_env.yaml" . "2023-09-08 12:04:05.469780+00:00" . . . . "3019"^^ . "https://api.rohub.org/api/resources/73838afc-95e9-4778-9397-b1177d41d5a9/download/" . . "2023-09-08 12:04:05.478657+00:00" . "2023-09-08 12:05:14.941074+00:00" . . "create_input_for_cmdBLAST.R" . "2023-09-08 12:04:05.478657+00:00" . . . . "1962"^^ . "https://api.rohub.org/api/resources/76f62c7c-9e00-43e4-8065-99eee7aae198/download/" . . "2023-09-08 12:04:05.480716+00:00" . "2023-09-08 12:05:15.392082+00:00" . . "heterotrichera.stk" . "2023-09-08 12:04:05.480716+00:00" . . . . "16338"^^ . "https://api.rohub.org/api/resources/81147bb0-01a7-4a8d-b2c6-d4aa632ba662/download/" . . "2023-09-08 12:04:05.476362+00:00" . "2023-09-08 12:05:14.172718+00:00" . . "read_results_infernal.R" . "2023-09-08 12:04:05.476362+00:00" . . . . "218"^^ . "https://api.rohub.org/api/resources/8cebbd0d-c2fc-4812-a092-5965ce747c14/download/" . . "2023-09-08 12:04:05.455163+00:00" . "2023-09-08 12:04:07.781299+00:00" . . "config.yaml" . "2023-09-08 12:04:05.455163+00:00" . . . . "176087"^^ . "https://api.rohub.org/api/resources/9cbd68ad-d410-48e1-8bd2-3067f695e62d/download/" . . "2023-09-08 12:04:05.457012+00:00" . "2023-09-08 12:04:08.004058+00:00" . . "cov_model_MRP" . "2023-09-08 12:04:05.457012+00:00" . . . . "1833"^^ . "https://api.rohub.org/api/resources/a2b33867-14fb-4100-b9a5-c8286bcf138e/download/" . . "2023-09-08 12:04:05.467422+00:00" . "2023-09-08 12:04:10.485032+00:00" . . "blast_env.yaml" . "2023-09-08 12:04:05.467422+00:00" . . . . "2567"^^ . "https://api.rohub.org/api/resources/a5812e02-c6e1-4bf0-a13f-9fad72c73f79/download/" . . "2023-09-08 12:04:05.479423+00:00" . "2023-09-08 12:05:15.204459+00:00" . . "join_infernal_and_BLASTcmd_results.R" . "2023-09-08 12:04:05.479423+00:00" . . . . "538"^^ . "https://api.rohub.org/api/resources/af1cd576-0d29-4f18-be29-84e5ec3ccd82/download/" . . "2023-09-08 12:04:05.471287+00:00" . "2023-09-08 12:05:12.619195+00:00" . . "entrez_env.yaml" . "2023-09-08 12:04:05.471287+00:00" . . . . "6627"^^ . "https://api.rohub.org/api/resources/b5c49c93-d1ad-4601-80c4-70277c4fd400/download/" . . "2023-09-08 12:04:05.468996+00:00" . "2023-09-08 12:04:11.117606+00:00" . . "r_tidyverse_env.yaml" . "2023-09-08 12:04:05.468996+00:00" . . . . "6670"^^ . "https://api.rohub.org/api/resources/bf376d44-9a16-4535-92b1-4b66ad2a519a/download/" . . "2023-09-08 12:04:05.468204+00:00" . "2023-09-08 12:04:10.908257+00:00" . . "cmdBLAST_to_R_env.yaml" . "2023-09-08 12:04:05.468204+00:00" . . . . "459"^^ . "https://api.rohub.org/api/resources/c119db00-a9bb-469c-8aef-97d2927ad412/download/" . . "2023-09-08 12:04:05.475571+00:00" . "2023-09-08 12:05:13.878329+00:00" . "text/x-sh" . . "cmdBLAST.sh" . "2023-09-08 12:04:05.475571+00:00" . . . . "851052"^^ . "https://api.rohub.org/api/resources/d642e7aa-836b-4808-a6a7-65043f69a119/download/" . . "2023-09-08 12:04:05.447483+00:00" . "2023-09-08 12:04:06.648816+00:00" . "image/png" . . "Hits_distribution_heatmap.png" . "2023-09-08 12:04:05.447483+00:00" . . . . "227587"^^ . "https://api.rohub.org/api/resources/d7c7b947-a603-404d-bbba-42138827e166/download/" . . "2023-09-08 12:04:05.452494+00:00" . "2023-09-08 12:04:07.138945+00:00" . "application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" . . "summary_table_models.xlsx" . "2023-09-08 12:04:05.452494+00:00" . . . . "7577"^^ . "https://api.rohub.org/api/resources/f985439a-75a1-4c88-9d85-c49b1bf41a9f/download/" . . "2023-09-08 12:04:05.474815+00:00" . "2023-09-08 12:05:13.695278+00:00" . . "make_table_plots.R" . "2023-09-08 12:04:05.474815+00:00" . . . . . . "search" . . "6.270903010033445" . "7.5" . "database' folder" . . "3.6474164133738602" . "2.4" . "ribonucleic acid" . . "7.5250836120401345" . "9.0" . "username" . . "5.5464926590538335" . "3.4" . "username" . . "5.518394648829432" . "6.6" . "covariance" . . "7.859531772575251" . "9.4" . "software" . . "27.474747474747474" . "13.6" . "The software" . . "2.8875379939209727" . "1.9" . "file" . . "7.6672104404567705" . "4.7" . "database query" . . "7.750759878419452" . "5.1" . "directory" . . "3.3444816053511706" . "4.0" . "genetics" . . "13.93939393939394" . "6.9" . "throughput" . . "2.508361204013378" . "3.0" . "homology search" . . "3.4954407294832825" . "2.3" . "folder" . . "8.862876254180602" . "10.6" . "permission notice" . . "2.8875379939209727" . "1.9" . "GERONIMO tool" . . "3.4954407294832825" . "2.3" . "pipeline" . . "8.19397993311037" . "9.8" . "folder" . . "8.97226753670473" . "5.5" . "GERONIMO" . . "12.39804241435563" . "7.6" . "homology" . . "3.5953177257525084" . "4.3" . "However, please note that deleting this folder will require the reinstallation of GERONIMO dependencies when the analysis is rerun." . . "20.88607594936709" . "3.3" . "file explorer" . . "5.471124620060791" . "3.6" . "genome" . . "9.448160535117058" . "11.3" . "computer programming and software" . . "100.0" . "0.41168108582496643" . "Rfam database" . . "2.8875379939209727" . "1.9" . "Detect RNA sequence" . . "3.6474164133738602" . "2.4" . "earth sciences" . . "100.0" . "0.9349703192710876" . "taxonomy" . . "3.3444816053511706" . "4.0" . "search" . . "6.688417618270798" . "4.1" . "name" . . "4.893964110929853" . "3.0" . "Agata Kilar" . . "Mendel Centre for Plant Genomics and Proteomics" . . . . "2025-11-11T16:09:16.035+01:00"^^ . . . "Research Object Crate for GERONIMO" . "RSA" . "MIIBIjANBgkqhkiG9w0BAQEFAAOCAQ8AMIIBCgKCAQEA4pPaESKwmC6l37P86K6TNLq6yeQtc7m9CvcqauLs/1FC0viHvQnFBgxj0a+loPDv/Egwe6OqFpa0iW9Ypnyz9YPoh+pxbRXonbuMOb+8Ry9hXZ+TEKfWjhjVDGEaClwfRwglh2HI/xfV4CD9AgvDOEoZQiyta8a90PYwJ3G6e70oCHTn61+OWTkI9KRYHOYgg3btdy2Z7q/30PTFawb2ZT5aIfIJYobUYv2a7yhtcqWCHZeKv0bxGnRjTFNx1rscBMlLJSzvRtpQc1cCRVEPFZHo1adaXCI9tGvn4cxeNQ96y8dxkN1XhpaJairde+23MDzf42Oe97KG2HYzKiyVnQIDAQAB" . "2wqnBh1ev32QIZtPLmjLfElb3kGQK7i4YbhM2CyEmEh8Bkuib6BYJ9WRECh92ytxjBhveTbeKOAWTGeU9yTzDfvzundH22XYiTkhYm1hJei9tezNxdIE/GwG4f2BWzzegFtaxAkh52+yZJK2HgWfz9Xwo4ZOnZUArbJtxP/S1Nlbl/p6I+Cdwbm3lKRtpINW2X2GE/VVvSpEa3KfJxit11Ot3iFjgOKGfWs3Mrl46kIGi6nxtrMTwllKkyS2FnnV2Ofa7IUl66jR2OB+gi1ivTNs3ztuWolMi+ssMBs9R/dJQLJxDDx4Fzy4DzcPTMIx+gP+CSAz7B1yOoJkUr5eGA==" . . .