<div align="center">
<img src="https://github.com/Goosang-Yu/genet/blob/main/docs/en/assets/images/logo.png?raw=true" width="300"/>
**Genome Editing Toolkit** </br>
**Since 2022. 08. 19.** </br>
[](https://badge.fury.io/py/genet)
[](https://badge.fury.io/py/genet)
[](https://genethq.slack.com/archives/C04DP727E4E)
[](https://goosang-yu.github.io/genet/getting_started/)
[](https://img.shields.io/pypi/l/ansicolortags.svg)
<div align="left">
## Welcome to GenET
GenET (Genome Editing Toolkit) is a library of various python functions for the purpose of analyzing and evaluating data from genome editing experiments. GenET is still in its early stages of development and continue to improve and expand. Currently planned functions include guideRNA design, saturation library design, deep sequenced data analysis, and guide RNA activity prediction.
Please see the [documentation](https://goosang-yu.github.io/genet/).
## Installation
```python
# Create virtual env for genet
conda create -n genet python=3.10
conda activate genet
# Install genet
pip install genet
```
## Who should use GenET?
GenET was developed for anyone interested in the field of genome editing. Especially, Genet can provide aid to those with the following objectives.: <br />
- Develop a quick and easy to design an genome editing experiment for a specific gene.
- Perform genome editing analysis based on sequening data
- Predict the activtiy of specific guideRNAs or all guideRNAs designed for editing a specific product.
- Design a saturation library for a specific gene.
## Example 1: Download genomic data from NCBI database
The genomic information required for research is often downloaded from public databases like NCBI. When only basic information is needed, searching on the NCBI website is usually sufficient. However, when a large amount of data is required or when specific reference sequence files are needed for certain analysis pipelines, it may be necessary to find and download files containing the specific information.
The GenET database module provides functions to easily download frequently used data files from NCBI. For example, to download the genomic assembly of `Homo sapiens`, you can use `GetGenome` as follows:
```python
from genet.database import GetGenome
```
To create a GetGenome instance for the desired species:
```python
# Specify the species
species = "Homo sapiens"
# Create a GetGenome instance
genome = GetGenome(species)
```
To check the available files related to the assembly of the specified species:
```python
# Check available files
available_files = genome.contents()
print("Available files:", available_files)
```
> Available files:
['README.txt',
'Annotation_comparison',
'GCF_000001405.40_GRCh38.p14_assembly_structure',
'GCF_000001405.40-RS_2023_10_annotation_report.xml',
'annotation_hashes.txt',
'RefSeq_transcripts_alignments',
'GCF_000001405.40_GRCh38.p14_assembly_regions.txt',
...]
To download the desired file with the specified name to the desired path:
```python
# Specify the desired file name
file_name = "example_genome.fasta"
# Specify the desired download path
download_path = "/desired/download/path/"
# Download the file
genome.download(file_name, download_path)
```
## Example 2: Prediction of prime editing efficiency by DeepPrime
DeepPrime is a prediction model for evaluating prime editing guideRNAs (pegRNAs) that target specific target sites for prime editing ([Yu et al. Cell 2023](https://doi.org/10.1016/j.cell.2023.03.034)). DeepSpCas9 prediction score is calculated simultaneously and requires tensorflow (version >=2.6). DeepPrime was developed on pytorch. For more details, please see the [documentation](https://goosang-yu.github.io/genet/).
```python
from genet.predict import DeepPrime
seq = 'CCGAGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCACGCTCCATTATC(C/T)AGCCCCAAAGCGCAACAAGCCCACTGTCTATGGTGTGTCCCCCAACTACGACAAGTGGGA'
pegrna = DeepPrime(seq)
# check designed pegRNAs
pegrna.features.head()
```
| | ID | Spacer | RT-PBS | PBS_len | RTT_len | RT-PBS_len | Edit_pos | Edit_len | RHA_len | Target | ... | deltaTm_Tm4-Tm2 | GC_count_PBS | GC_count_RTT | GC_count_RT-PBS | GC_contents_PBS | GC_contents_RTT | GC_contents_RT-PBS | MFE_RT-PBS-polyT | MFE_Spacer | DeepSpCas9_score |
| - | ---------- | -------------------- | ------------------------------------------------- | ------- | ------- | ---------- | -------- | -------- | ------- | ------------------------------------------------- | --- | --------------- | ------------ | ------------ | --------------- | --------------- | --------------- | ------------------ | ---------------- | ---------- | ---------------- |
| 0 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATG | 7 | 38 | 45 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \-510.285 | 2 | 23 | 25 | 28.57143 | 60.52632 | 55.55556 | \-12.7 | 0 | 76.43662 |
| 1 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGA | 8 | 38 | 46 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \-510.285 | 2 | 23 | 25 | 25 | 60.52632 | 54.34783 | \-11.4 | 0 | 76.43662 |
| 2 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGAT | 9 | 38 | 47 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \-510.285 | 2 | 23 | 25 | 22.22222 | 60.52632 | 53.19149 | \-11.4 | 0 | 76.43662 |
| 3 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATG | 10 | 38 | 48 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \-510.285 | 3 | 23 | 26 | 30 | 60.52632 | 54.16667 | \-11.2 | 0 | 76.43662 |
| 4 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATGA | 11 | 38 | 49 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \-510.285 | 3 | 23 | 26 | 27.27273 | 60.52632 | 53.06122 | \-11.2 | 0 | 76.43662 |
Next, select model PE system and run DeepPrime
```python
pe2max_output = pegrna.predict(pe_system='PE2max', cell_type='HEK293T')
pe2max_output.head()
```
| | ID | PE2max_score | Spacer | RT-PBS | PBS_len | RTT_len | RT-PBS_len | Edit_pos | Edit_len | RHA_len | Target |
| - | ---------- | ------------ | -------------------- | ------------------------------------------------- | ------- | ------- | ---------- | -------- | -------- | ------- | ------------------------------------------------- |
| 0 | SampleName | 2.143 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATG | 7 | 38 | 45 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |
| 1 | SampleName | 3.140197 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGA | 8 | 38 | 46 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |
| 2 | SampleName | 2.541219 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGAT | 9 | 38 | 47 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |
| 3 | SampleName | 6.538445 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATG | 10 | 38 | 48 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |
| 4 | SampleName | 7.436117 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATGA | 11 | 38 | 49 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |
Please send all comments and questions to gsyu93@gmail.com
Raw data
{
"_id": null,
"home_page": null,
"name": "genet",
"maintainer": null,
"docs_url": null,
"requires_python": ">=3.9",
"maintainer_email": null,
"keywords": "CRISPR, ai, analysis, bioinformatics, deep-learning, genetics, genome-editing, python",
"author": null,
"author_email": "Goosang Yu <gsyu93@gmail.com>",
"download_url": "https://files.pythonhosted.org/packages/8a/fd/1ef1c002cf1bef3979d8b24048f8be0afa9737fb0827945875ebc8b3819c/genet-0.15.0.tar.gz",
"platform": null,
"description": "<div align=\"center\">\n \n <img src=\"https://github.com/Goosang-Yu/genet/blob/main/docs/en/assets/images/logo.png?raw=true\" width=\"300\"/>\n\n**Genome Editing Toolkit** </br>\n**Since 2022. 08. 19.** </br>\n\n[](https://badge.fury.io/py/genet) \n[](https://badge.fury.io/py/genet) \n[](https://genethq.slack.com/archives/C04DP727E4E)\n[](https://goosang-yu.github.io/genet/getting_started/)\n[](https://img.shields.io/pypi/l/ansicolortags.svg) \n\n\n<div align=\"left\">\n\n## Welcome to GenET\nGenET (Genome Editing Toolkit) is a library of various python functions for the purpose of analyzing and evaluating data from genome editing experiments. GenET is still in its early stages of development and continue to improve and expand. Currently planned functions include guideRNA design, saturation library design, deep sequenced data analysis, and guide RNA activity prediction.\n\nPlease see the [documentation](https://goosang-yu.github.io/genet/).\n\n\n## Installation\n```python\n# Create virtual env for genet\nconda create -n genet python=3.10\nconda activate genet\n\n# Install genet\npip install genet\n```\n\n## Who should use GenET?\nGenET was developed for anyone interested in the field of genome editing. Especially, Genet can provide aid to those with the following objectives.: <br />\n\n- Develop a quick and easy to design an genome editing experiment for a specific gene.\n- Perform genome editing analysis based on sequening data\n- Predict the activtiy of specific guideRNAs or all guideRNAs designed for editing a specific product.\n- Design a saturation library for a specific gene.\n\n\n\n## Example 1: Download genomic data from NCBI database\n\nThe genomic information required for research is often downloaded from public databases like NCBI. When only basic information is needed, searching on the NCBI website is usually sufficient. However, when a large amount of data is required or when specific reference sequence files are needed for certain analysis pipelines, it may be necessary to find and download files containing the specific information.\n\nThe GenET database module provides functions to easily download frequently used data files from NCBI. For example, to download the genomic assembly of `Homo sapiens`, you can use `GetGenome` as follows:\n\n```python \nfrom genet.database import GetGenome\n```\n\nTo create a GetGenome instance for the desired species:\n```python\n# Specify the species\nspecies = \"Homo sapiens\"\n\n# Create a GetGenome instance\ngenome = GetGenome(species)\n```\n\nTo check the available files related to the assembly of the specified species:\n```python\n# Check available files\navailable_files = genome.contents()\nprint(\"Available files:\", available_files)\n```\n\n> Available files: \n ['README.txt',\n 'Annotation_comparison',\n 'GCF_000001405.40_GRCh38.p14_assembly_structure',\n 'GCF_000001405.40-RS_2023_10_annotation_report.xml',\n 'annotation_hashes.txt',\n 'RefSeq_transcripts_alignments',\n 'GCF_000001405.40_GRCh38.p14_assembly_regions.txt',\n ...]\n\nTo download the desired file with the specified name to the desired path:\n```python\n# Specify the desired file name\nfile_name = \"example_genome.fasta\"\n\n# Specify the desired download path\ndownload_path = \"/desired/download/path/\"\n\n# Download the file\ngenome.download(file_name, download_path)\n```\n\n\n## Example 2: Prediction of prime editing efficiency by DeepPrime\n\nDeepPrime is a prediction model for evaluating prime editing guideRNAs (pegRNAs) that target specific target sites for prime editing ([Yu et al. Cell 2023](https://doi.org/10.1016/j.cell.2023.03.034)). DeepSpCas9 prediction score is calculated simultaneously and requires tensorflow (version >=2.6). DeepPrime was developed on pytorch. For more details, please see the [documentation](https://goosang-yu.github.io/genet/).\n\n```python \nfrom genet.predict import DeepPrime\n\nseq = 'CCGAGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCACGCTCCATTATC(C/T)AGCCCCAAAGCGCAACAAGCCCACTGTCTATGGTGTGTCCCCCAACTACGACAAGTGGGA'\n\npegrna = DeepPrime(seq)\n\n# check designed pegRNAs\npegrna.features.head()\n```\n\n| | ID | Spacer | RT-PBS | PBS_len | RTT_len | RT-PBS_len | Edit_pos | Edit_len | RHA_len | Target | ... | deltaTm_Tm4-Tm2 | GC_count_PBS | GC_count_RTT | GC_count_RT-PBS | GC_contents_PBS | GC_contents_RTT | GC_contents_RT-PBS | MFE_RT-PBS-polyT | MFE_Spacer | DeepSpCas9_score |\n| - | ---------- | -------------------- | ------------------------------------------------- | ------- | ------- | ---------- | -------- | -------- | ------- | ------------------------------------------------- | --- | --------------- | ------------ | ------------ | --------------- | --------------- | --------------- | ------------------ | ---------------- | ---------- | ---------------- |\n| 0 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATG | 7 | 38 | 45 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \\-510.285 | 2 | 23 | 25 | 28.57143 | 60.52632 | 55.55556 | \\-12.7 | 0 | 76.43662 |\n| 1 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGA | 8 | 38 | 46 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \\-510.285 | 2 | 23 | 25 | 25 | 60.52632 | 54.34783 | \\-11.4 | 0 | 76.43662 |\n| 2 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGAT | 9 | 38 | 47 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \\-510.285 | 2 | 23 | 25 | 22.22222 | 60.52632 | 53.19149 | \\-11.4 | 0 | 76.43662 |\n| 3 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATG | 10 | 38 | 48 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \\-510.285 | 3 | 23 | 26 | 30 | 60.52632 | 54.16667 | \\-11.2 | 0 | 76.43662 |\n| 4 | SampleName | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATGA | 11 | 38 | 49 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... | ... | \\-510.285 | 3 | 23 | 26 | 27.27273 | 60.52632 | 53.06122 | \\-11.2 | 0 | 76.43662 |\n\n\nNext, select model PE system and run DeepPrime\n```python \npe2max_output = pegrna.predict(pe_system='PE2max', cell_type='HEK293T')\n\npe2max_output.head()\n```\n\n| | ID | PE2max_score | Spacer | RT-PBS | PBS_len | RTT_len | RT-PBS_len | Edit_pos | Edit_len | RHA_len | Target |\n| - | ---------- | ------------ | -------------------- | ------------------------------------------------- | ------- | ------- | ---------- | -------- | -------- | ------- | ------------------------------------------------- |\n| 0 | SampleName | 2.143 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATG | 7 | 38 | 45 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |\n| 1 | SampleName | 3.140197 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGA | 8 | 38 | 46 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |\n| 2 | SampleName | 2.541219 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGAT | 9 | 38 | 47 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |\n| 3 | SampleName | 6.538445 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATG | 10 | 38 | 48 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |\n| 4 | SampleName | 7.436117 | GGTTCATCATCATTCAACGG | TAGATAATGGAGCGTGGTGATGAGCCCGTCGGCCACCGTTGAATGATGA | 11 | 38 | 49 | 37 | 1 | 1 | AGTTGGTTCATCATCATTCAACGGTGGCCGACGGGCTCATCACCAC... |\n\n\nPlease send all comments and questions to gsyu93@gmail.com",
"bugtrack_url": null,
"license": null,
"summary": "GenET: Genome Editing Toolkit",
"version": "0.15.0",
"project_urls": {
"Homepage": "https://github.com/Goosang-Yu/genet",
"Repository": "https://github.com/Goosang-Yu/genet",
"Source": "https://github.com/Goosang-Yu/genet",
"Tracker": "https://github.com/Goosang-Yu/genet/issues"
},
"split_keywords": [
"crispr",
" ai",
" analysis",
" bioinformatics",
" deep-learning",
" genetics",
" genome-editing",
" python"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "fbe6b2e3827226ab7e5b135bb78997958ffa4d29ff1a590dd5fcae6cefb5faa5",
"md5": "91f82b6188b1fd5e5906f82ccd194e08",
"sha256": "b9a9baf83d4e247219e10db1feb4d7f43b9230faa82debbf7f5d48561aee1eeb"
},
"downloads": -1,
"filename": "genet-0.15.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "91f82b6188b1fd5e5906f82ccd194e08",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.9",
"size": 104755,
"upload_time": "2024-04-25T09:09:35",
"upload_time_iso_8601": "2024-04-25T09:09:35.097104Z",
"url": "https://files.pythonhosted.org/packages/fb/e6/b2e3827226ab7e5b135bb78997958ffa4d29ff1a590dd5fcae6cefb5faa5/genet-0.15.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "8afd1ef1c002cf1bef3979d8b24048f8be0afa9737fb0827945875ebc8b3819c",
"md5": "c196c85919e10b3bdaafa4515305daaf",
"sha256": "5e591b7260d8e84b7eb1b8631016cde65e7499d71ac0a97c142d6d5682281e9a"
},
"downloads": -1,
"filename": "genet-0.15.0.tar.gz",
"has_sig": false,
"md5_digest": "c196c85919e10b3bdaafa4515305daaf",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.9",
"size": 106647,
"upload_time": "2024-04-25T09:09:37",
"upload_time_iso_8601": "2024-04-25T09:09:37.689491Z",
"url": "https://files.pythonhosted.org/packages/8a/fd/1ef1c002cf1bef3979d8b24048f8be0afa9737fb0827945875ebc8b3819c/genet-0.15.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-04-25 09:09:37",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "Goosang-Yu",
"github_project": "genet",
"travis_ci": false,
"coveralls": false,
"github_actions": true,
"lcname": "genet"
}
JSON
