Crispr/Cas9 Target Prediction with Deep Learning
Total Page:16
File Type:pdf, Size:1020Kb
Research Article ISSN: 2574 -1241 DOI: 10.26717/BJSTR.2019.21.003652 Crispr/Cas9 Target Prediction with Deep Learning Özlem Aktaş, Elif Doğan* and Tolga Ensari Department of Computer Engineering, Faculty of Engineering, Turkey *Corresponding author: Elif Doğan, Department of Computer Engineering, Faculty of Engineering, Dokuz Eylül University, Turkey ARTICLE INFO Abstract Received: Published: September 24, 2019 The CRISPR/CAS9 system is a powerful tool for regulating damaged genome October 03, 2019 sequences. Nucleases that are damaged in their sequence are called miRNAs (micro RNAs). The miRNAs targeted by multiple promoter sgRNA (single guide RNA) are Citation: cut or regulated from RNA by the CRISPR/CAS9 method. The sgRNAs targeted to the Özlem Aktaş, Elif Doğan, Tolga incorrect miRNAs may provoke undesired genome abnormalities. In this study, in Ensari. Crispr/Cas9 Target Prediction with order to minimize these genome distortions, sgRNA target estimation was performed Deep Learning. Biomed J Sci & Tech Res for CRISPR/CAS9 with deep learning in this study. In this article, convolutional neural 21(5)-2019. BJSTR. MS.ID.003652. networks (Convolutional Neural Networks-CNN) and multilayer perceptron (Multi- Layer Perceptron-MLP) algorithms are used for experimental analysis. We also compare theAbbreviations: performance of CRISPR/CAS9 system for three algorithms. CNN: Convolutional Neural Networks; BLSTM: Bidirectional Long-Short Term Memory; MLP: Multi Layer Perceptron Introduction E. Coli In in silico silico the genomes to be studied in this respect. Various miRNA target In an Esherichia Coli ( ) bacteria that investigated by “ data sets used in “ ” studies will be used [7]. The aim of this ” (simulation) method which is important role nowadays study, it is made mirRNA target estimation with machine learning has been discovered an immune system named CRISPR (clustered E.Coli algorithms. The result of wrong targeted miRNA may be cause of regularly interspaced short palindromic repeats)/Cas9 [1]. undesirable gene mutations [8]. It is aimed to minimize the Type According to this system an bacteria which is infected by 1 and Type 2 errors of miRNAs by applying machine learning and any virus, add the virus DNA its memory and recognize when any deep learning method of Type 1 and Type 2 errors as output. Thus, other virus attack accrues. This defined virus DNA cuts this virus the wrong target estimate will be minimized. In this way, it will be DNA from its own DNA through the CAS9 enzyme. Thus, DNA repair possible to reliably repair gene damage by correctly targeting mis- will occur. According to recent research the leading factor in gene targeted miRNAs. Genetic diseases will be eliminated by repair of regulation has been observed the “microRNA (miRNA)” targeted gene damage. in silico by the “single-guide RNA (sGRNA)” s. MiRNAs are small and non- coding RNA molecules [2]. In diseases, for example, different The purpose of the “ ” technique is to increase the miRNAs involved in all stages of cancer make cancer diagnosis and accuracy of disease prevention by single pointing with the help of treatment available [3]. CRISPR/CAS9 is used to destroy targeted mechanization. The large base readings provided by mechanization miRNAs in cells [4]. It uses “guideRNA (gRNA)” as a guide to target also contribute greatly. Studies for miRNA targeting, support vector the CRISPR/CAS9 nucleus to the DNA sequence and trigger the machines (Support Vector Machine -SVM) [9], deep learning, double-strand break at the desired location. Breaking and repair of constrained logic programming (Constraint Logic Programming) these threads can cause random addition and alteration of DNA [5]. [10] and a class classification (One Class Classification-OCC) [11] methods were used. Besides, in the future big data may be MiRNA activities are different in each cell type [6]. So, the decreased to its most valuable parts with digital data forgetting miRNA activities of a human cell and a bacterial cell will not be concepts. It makes computations faster and the machines will store the same. From this point, It is important to find a dataset about less data in disks [37]. Copyright@ 16131 Elif Doğan | Biomed J Sci & Tech Res | BJSTR. MS.ID.003652. Volume 21- Issue 5 DOI: 10.26717/BJSTR.2019.21.003652 Materials and Methods CRISPR [22], CRISPOR [23], CRISPR Local [24]. In the algorithm studies performed with these data sets, estimation tools such as mirWalk [25], TargetScan (ID2 PPI analysis network) [26], miRanda [27], CRISPR/CAS9 system is a strong genome editing mechanism mirBase [28], mirTarget [29], TarBase [30] have been developed. used in many biotechnology applications [12]. MiRNAs are RNA molecules that play an important role in gene regulation in In this study, CRISPR Local data set has been used for the animals and plants [13]. Damaged miRNAs are targeted by sgRNAs. analysis [31]. The source of CRISPR Local dataset is ENSEMBL However, the effectiveness of sgRNAs have not been defined firmly Plants. There are approximately 854.610 lines of CRISPR data in to target area [14]. In the systems of CRISPR/CAS9 (clustered the original. In this study, we used 34.200 lines of data because of regularly interspaced short palindromic repeats), Cas9 nuclease technical issues. There are 11 column features in this dataset which creates short array (about 20 nucleotide) with RNA guideline two are examples of “Cyanidioschyzon merolae” alga. This features; the stranded in the determined region at DNA. CRISPR/CAS9 has been gene in which the sgRNA, on target estimated chromosome and a unique technology that allows genetic and medical researchers its coordinate, sgRNA sequence with 23’nt., on target prediction to add, subtract or modify DNA sequences in various parts of the score, off target prediction gene which has the greatest CFD score, genome. In contrast to ZFN (Zinc-Finger Nuclease) and TALENs the chromosome on target prediction, its coordinate and beginning (Transcription Activator-Like Effector Nuclease), CRISPR/CAS9 is position, off-target prediction sequence, the number of sgRNA not man-made; currently the system is part of the bacterial immune and mismatch on the off-target sequence, axon name, axon start system, which helps to protect from invasive phages. Originality in position, all off-target and sgRNA having the highest CFD score. The CRISPR is achieved using an RNA molecule that is complementary sequences having 4 channels like Adenin (A), Sitozin (C), Guanin to the gene of interest. After binding, this RNA molecule (also called (G) ve Timin (T) are used as [1,0,0,0], [0,1,0,0], [0,0,1,0], [0,0,0,11]. guide RNA) traps a CAS9 nuclease, a double-strand break that According these channels, the process has been realized with causes a frame shift when repaired by NHEJ (Non-Homologous End converted RNA sequence to binary system. These sequences were Joining). CRISPR is the most effective process to date in gene repair used as binary in the data set. In this example, each base in the and editing tools. Efficiency varies depending on the organism and sequence is considered as a separate column and feature. The 23nt. the target site. Since it is the simplest, versatile and sensitive method sgRNA sequence, on-target prediction score features were used. of genetic manipulation that is currently available, it attracts great Figure 1 shows an illustration of the sample data set. used. Figure 1 attention in the world of science. In summary, CRISPR/CAS9 differs shows an illustration of the sample data set. from ZFN and TALENs in an important aspect that makes them superior to genomic regulation applications: The ZFN and TALENs bind to DNA through a direct protein-DNA interaction that requires DATAredesigning Set the protein for each new target. Applied “in silico” research and review, BLAST (Basic Local Figure 1: Sample Data Set. Alignment Search Tool) was used for CRISPR [15]. BLAST is a search tool that analyzes the amino acids and DNA sequences of proteins and finds similarities between them. Besides BLAST, Experimental Studies the data set resources have been used such as National Human Genome Research Institute (NHGRI) [16], miRBase [17], Genome In this study, CNN, MLP and BLSTM models were constituted Crispr [18], CrisprInc [19], ENSEMBL [20], ENCODE [21], CRISPRz Multilayer Perceptron -MLP and compared. Figure 2: MLP Model. Copyright@ 16132 Elif Doğan | Biomed J Sci & Tech Res | BJSTR. MS.ID.003652. Volume 21- Issue 5 DOI: 10.26717/BJSTR.2019.21.003652 A 4-layer MLP model has been developed according to the MLP of the matrices. In order to prevent over fitting, maxpooling layer algorithm run using the Google Colab GPU (https://colab.research. was used. It selects the elements with the maximum value from the google.com). In the MLP model a fully connected structure of dense matrix pool of the specified size in the maxpooling layer. Accordingly, layers was formed. Information on the model used is shown in the information obtained when the CNN model is generated can be Figure 2. The rates of logistic regression and accuracy according seen in Figure 6. The rate of accuracy according to the CNN model to the MLP model are shown in Figure 3 and Figure 4. Accuracy can be seen from Figure 7. obtained 80.38% according to MLP model. Figure 5: Convolution Process. Figure 3: Model Accuracy. Figure 4: Model loss. Convolutional Neural Network - CNN Figure 6: CNN Model. Convolutional Neural Networks are a model of artificial neutral network which is used successfully in image processing, bioinformatics, robotics, data mining, finance and many other areas. However, except for image analysis, surprisingly high accuracy ratio was obtained in emotion analysis, text classification and question answering applications. According to this model, it is applied to nxn matrix with nxn filtering method (dot product), with acceptation of n>m.