https://github.com/chrisnelsonlab/crispr-tn5
Code for analysis of unbiased sequencing of CRISPR-mediated genome editing
https://github.com/chrisnelsonlab/crispr-tn5
crispr matlab tns
Last synced: about 2 months ago
JSON representation
Code for analysis of unbiased sequencing of CRISPR-mediated genome editing
- Host: GitHub
- URL: https://github.com/chrisnelsonlab/crispr-tn5
- Owner: chrisnelsonlab
- Created: 2018-12-21T19:20:51.000Z (over 6 years ago)
- Default Branch: master
- Last Pushed: 2019-01-13T18:52:22.000Z (over 6 years ago)
- Last Synced: 2025-04-06T15:17:30.305Z (about 2 months ago)
- Topics: crispr, matlab, tns
- Language: MATLAB
- Homepage:
- Size: 379 KB
- Stars: 1
- Watchers: 0
- Forks: 1
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
# CRISPR-TN5
Code for characterizing complex genomic alterations caused by CRISPR gene editing using Illumina Nextera. Code is provided in MATLAB and should run in multiple operating systems but only tested on Windows. We used Illumina Nextera for preparing libraries but purified Tn5 could also work. See Extended data Figure (not published yet) to see details for optimizing the protocol.
# Purpose
The code is provided for reproducibility of the published report (insert link - not published yet). This is not intended as a research tool but if you use it and have questions, please don't hesitate to ask and I'll do my best to respond. The code analyzes genomic editing results of dual nucleases (e.g. CRISPR) for targeted deletion of a region of DNA. The analysis also includes room for two AAV vector genomes to queary for targeted integration, either intentional or unintentional (see parameter file).
## Figure 1 - Possible outcomes of multiplexed genome editing
# How Tn5-based enrichment works
## Figure 2 - Overview of Nextera enrichment for genome editing analysis
# How the code works
## Figure 3 - Outline of how the code works
- Read in fastq file in either .fastq or .fastq.gz
- Remove misprimnig events
- Bin reads as one of the following:
1. Intact unedited allele
2. Indel
3. Targeted deletion
4. AAV vector genome integration
5. Targeted inversion
6. Other (mystery) - these were then analyzed manually or with other alignment software. These could include chromosome translocations, large unexpected deletions, and sequencing/polymerase artificats.
# Experimental design
1. Nuclease (or CRISPR gRNA) design
2. Primer design
3. Optimization and troubleshooting
4. Analysis
# Parameter file
- Line 1 - Amplicon sequence if no editing takes place
- Line 2 - Amplicon sequence if a targeted deletion occurs
- Line 3 - Amplicon sequence if a targeted inversion occurs
- Line 4 - DSB location in bp
- Line 5 - BP window to check for indels on each side of DSB
- Line 6 - Minimum score for aligning
- Line 7 - Minimum score for aligning short portion
- Line 8 - Minimum alignment score for the AAV genome
- Line 9 - Large portion of target gene with inversion for unexpected modifications (e.g. chewback)
- Line 10 - Large portion of target gene for unexpected modifications (e.g. chewback)
- Line 11 - minimum amplicon size. Reject shorter amplicons
- Line 12 - AAV Genome 1
- Line 13 - AAV genome 2
- Line 14 - Value for gap opening when aligning to the dystrophin gene
- Line 15 - Value for gap opening when aligning to the AAV genome
- Line 16 - Value for gap extending for AAV genome
## Example Parameter file
> See linked files here
Parameter file 1
Parameter file 2
## Example data set
> See linked .fastq files here. These paired-end reads provide 1000 reads to ensure the algorithm works. The complete data set is linked through the NCBI (see below).
## Data to recreate figures from (PAPER) can be found here
> https://www.ncbi.nlm.nih.gov/sra/?term=SRP157083