- Research
- Open access
- Published:
A deep learning method for lincRNA detection using auto-encoder algorithm
BMC Bioinformatics volume 18, Article number: 511 (2017)
Abstract
Background
RNA sequencing technique (RNA-seq) enables scientists to develop novel data-driven methods for discovering more unidentified lincRNAs. Meantime, knowledge-based technologies are experiencing a potential revolution ignited by the new deep learning methods. By scanning the newly found data set from RNA-seq, scientists have found that: (1) the expression of lincRNAs appears to be regulated, that is, the relevance exists along the DNA sequences; (2) lincRNAs contain some conversed patterns/motifs tethered together by non-conserved regions. The two evidences give the reasoning for adopting knowledge-based deep learning methods in lincRNA detection. Similar to coding region transcription, non-coding regions are split at transcriptional sites. However, regulatory RNAs rather than message RNAs are generated. That is, the transcribed RNAs participate the biological process as regulatory units instead of generating proteins. Identifying these transcriptional regions from non-coding regions is the first step towards lincRNA recognition.
Results
The auto-encoder method achieves 100% and 92.4% prediction accuracy on transcription sites over the putative data sets. The experimental results also show the excellent performance of predictive deep neural network on the lincRNA data sets compared with support vector machine and traditional neural network. In addition, it is validated through the newly discovered lincRNA data set and one unreported transcription site is found by feeding the whole annotated sequences through the deep learning machine, which indicates that deep learning method has the extensive ability for lincRNA prediction.
Conclusions
The transcriptional sequences of lincRNAs are collected from the annotated human DNA genome data. Subsequently, a two-layer deep neural network is developed for the lincRNA detection, which adopts the auto-encoder algorithm and utilizes different encoding schemes to obtain the best performance over intergenic DNA sequence data. Driven by those newly annotated lincRNA data, deep learning methods based on auto-encoder algorithm can exert their capability in knowledge learning in order to capture the useful features and the information correlation along DNA genome sequences for lincRNA detection. As our knowledge, this is the first application to adopt the deep learning techniques for identifying lincRNA transcription sequences.
Background
LincRNA refers to long intergenic non-coding RNA with the length greater than 200 nucleotides that are transcribed from non-coding DNA sequences between protein-coding regions. These intergenic regions were referred as junk DNA, however, now it is discovered that intergenic regions can be transcribed and provide functional non-coding RNA genes within intergenic regions [1]. Various classes of transposable elements are embeded in lincRNAs and lincRNAs are viewed as a tool box of elements with some regulatory functions in transcription and translation. For example some lincRNAs attach to messenger RNA to block protein production [2] and families of transposable elements-derived lincRNAs have been implicated in the regulation of pluripotency [3]. In addition, lincRNA is highly tissue-specific, indicating that it might be closely related to epigenetic regulation. Thus, identifying these lincRNAs is the critical step towards understanding complicated regulatory mechanisms.
Non-coding RNA regions are four times longer than coding RNA sequences. However, currently only 21 thousand lincRNAs (about 2M bytes) are computationally discovered [4]. This is also one of the most important findings in lincRNA identification. The latest work are mostly based on RNA-seq data and heavily rely on the RNA-seq assembly technology [4, 5]. As the long intergenic non-coding RNAs are differentially expressed in different tissues and multiple conditions, the RNA-seq data sets allow to detect both rare and tissue-specific transcription events that would be undetectable in other limited studies, such as tiling array studies [6]. Thus, it establishes a philosophy that RNA-seq data can be used for lincRNA detection as the large volume of sequencing data are comprehensive and detailed. A general procedure of the state-of-the-art method to identify lincRNA is composed of the following main steps [5]: (1) Acquiring RNA-seq data set, (2) De novo RNA-seq assembly, (3) filtering and expression analysis.
Acquiring RNA-seq data sets is to collect the RNA sequencing data of different tissues under multiple conditions. Single RNA-seq data set cannot be used for the evidence of lincRNA detection. For example, in [4], more than one hundred previously published RNA-seq data sets covering more than twenty human tissues under multiple conditions and consisting of about four billion uniquely mapped reads. Subsequently, De novo RNA-seq transcriptome assembly [7] is used as the key technology to discover novel lincRNAs in a currently adopted model, which creates a transcriptome without the use of a reference genome. On the contrary, although the reference-based assembly method is a robust way of identifying transcript sequences using genome alignment, it is not able to account for incidents of structural alterations of mRNA transcripts, such as rare splicing sites and alternative splicing [8]. Instead, spliced variants are not actual proteins and they do not align continuously along the genome. An assembled transcript can be represented as introns and exons that are characterized as one of the features of an lincRNA. Thus, finding the alternative splicing transcripts from RNA-seq is regarded as one of the most important factors to the detection of novel lincRNAs. From the assembled transcripts, all known genes, pseudogenes, short ncRNAs, novel protein coding transcripts, novel UTRs, and non-lincRNA non-coding RNAs must be filtered to identify actual lincRNAs. Only intergenic non-coding transcripts with at least 200 nucleotides in length and expressed at least at one copy per cell are kept as ultimately annotated lincRNAs. A set of filters can be designed to achieve this goal.
The aforementioned techniques ensure the quality of annotated lincRNA data and provide the probability to develop a knowledge-based discovery method, although currently knowledge-based discovery methods for identifying the lincRNA remain on the preliminary stage. Driven by the newly found data set, scientists have found some hints that can corroborate their previous speculations: (1) the expression of lincRNAs appears to be regulated, that is, the relevance exists along the DNA sequences; (2) lincRNAs contain some conversed patterns/motifs tethered together by non-conserved regions [9]. The two evidences give the reasoning for developing knowledge-based deep learning methods in lincRNA detection.
The latest findings show that the expression of lincRNAs appears to be specifically regulated, although a widely accepted concept is that the degree to which intergenic transcription is functional remains uncertain and controversial [9]. According to the reasoning that negative transcripts (non-lincRNA) should lack coherent epigenetic patterns, the evaluation of lincRNAs depends on whether lincRNAs contains epigenetic markers. The catalog of lincRNAs shows some patterns of epigenetic modification similar to protein coding genes [10, 11]. For example, activating histone markers including H3K4me3 and H3K36me3 are both significantly contained within highly expressed lincRNAs; similarly, the repressive mark H3K27me3 is significantly enriched within lowly expressed lincRNAs.
The recent studies further reveal that the majority of the lincRNAs identified display a level of conservation consistent with known functional lincRNAs. This studies was performed through a 50 nt window to scan the sequences for the evaluation of conserved patterns [4]. Consistent with prior studies, lincRNAs display detectable but modest conservation [12]. Thus, by taking advantage of these patterns and conservations along DNA sequence, the knowledge-based discovery systems such as deep learning can discover more unidentified lincRNAs as long as the sufficient knowledge can be acquired. Fortunately, those newly found lincRNA data are able to provide such opportunities.
The preliminary concepts of deep learning including deep neural network were proposed in mid-2000s although the ideas of deep neural network had been discussed for long time since 90s [13–15]. After that, deep learning techniques have been applied to life sciences and shown tremendous promise [16–19]. Thus, deep-learning based technologies are regarded as potential tools for computational discovery of lincRNA. Deep neural network uses complicated algorithms, such as convolution, auto-encoder and Boltzmann machine etc., to constrain the error between layers and eliminate the back-propagation problem. Relying on a multiple-layer perceptron architecture, the estimation of input data through the hidden layer can be calculated by iterative encoding-decoding processing so that the minimum difference can be achieved between the input data and the estimation.
Deep learning related methods are barely seen in the methodology of lincRNA annotation. Based on those annotated data, deep learning based methods can exert their capability in knowledge learning in order to improve the aforementioned method and discover novel lincRNAs in DNA genomes.
In this project, three goals are set. The first one is developing a deep learning method for lincRNA transcription splicing sites. Second, validating the annotated lincRNAs transcription sites and testing the performance of deep learning method by comparing with conventional methods such as support vector machine (SVM) and traditional neural network based method. Third, computationally discovering other unidentified splicing sites. For the first goal, auto-encoder method achieves 100% prediction accuracy illustrated in next section. For the second and third goal, one unreported splicing site is found during re-scanning the whole annotated human lincRNA data sets through the deep learning method.
Methods
Auto-encoder
Auto-encoder (AE) is a layer-wise training algorithm we adopt on an artificial neural network that can be used to constitute a multiple-layer percetron architectures for deep learning machine shown in Fig. 1 a. The hidden layer h and the iterative estimation of x ∗ can be expressed as Eq. 1 by calculating the weights as illustrated in Fig. 1 b. The iteration becomes stable when it has the minimum distance between x and x ∗, as shown in Eq. 2. The preliminary ideas of shallow/deep neural network had been discussed for long time since 90s, however, mature concepts of deep learning including deep neural network were proposed in mid-2000s [13–15]. Since then, it has been applied to life sciences and shown tremendous promise [16–19].
Architecture of Deep Neural Network. a An Illustration of Deep Neural Network Architecture. b An Illustration of Auto-encoder
The simplest auto-encoder is based on a feedforward, non-recurrent neural network similar to the multiple-layer perceptron (MLP). The difference is that the output layer of auto-encoder has the same number of nodes as the input layer and an auto-encoder is trained to reconstruct their own inputs instead of being trained to predict the output value. Thus, training the neighboring set of two layers minimizes the errors between layers and eliminates the problem of error propagation that occurs in conventional neural network.
Our auto-encoder method is composed of three main steps as shown in Fig. 2: building, pre-training and validating. In the first step, the basic architecture including input layer, hidden layer and activation functions is built; secondly, the encoder and the decoder are trained layer by layer following the pre-configured iterations; thirdly, fine-grained training/validation is performed through the entire model. In other words, the first step constructs the basic framework of the deep neural network, the second one trains the layer-wise nodes and the last one flows through all layers for validation.
Flow Chart for Auto-encoder Method
As the core of auto-encoder, the pseudo-code of cost update algorithm is shown in Algorithm 1 following the Eqs. 1 and 2.
Transcription Sites
Similar to coding region transcription, non-coding regions are split at transcription sites. However, regulatory RNAs rather than message RNAs are generated. That is, the transcribed RNAs participate the biological process as regulatory units instead of generating proteins. Thus, identifying these transcriptional regions is the first step towards lincRNA recognition. Similar to gene structures, lincRNAs have the complicated exon/intron structures, whereas the difference from gene structures is that many of them have two exons or three exons only.
Benefiting from the increasing annotation data in lincRNAs, lincRNA transcriptional splicing site sequences are collected from the annotated human DNA genome data. However, the annotated data sets of lincRNAs are not so many as that of mRNAs. Thus, all of annotated lincRNAs are used for training, validation and testing.
In the same vein to detection of protein-coding splicing sites, auto-encoder neural network method is used for the lincRNA application. A 2-layer auto-encoder model is used for lincRNA detection and various encoding schemes are used for evaluating the best performance. The similar knowledge-based deep learning methods in lincRNA detection is barely mentioned in literature so far. The experimental results show an excellent predictive performance of deep neural network method on lincRNA data sets.
Encoding Schemes of DNA Sequence
Data representation, particularly the encoding scheme of DNA sequence, is one of important factors that can largely impact on the performance of knowledge-based discovery systems. Different from other data format, the DNA nucleotide sequences are recorded as human readable characters, C, T, A and G. Adopting the improper encoding schemes to feed the learning machine can lead to the failure of prediction task. The encoding schemes we test are shown as Fig. 3, including DAX [20], EIIP [21], Complementary [22], Enthalpy [23], and Galois(4) [24] schemes.
Five Encoding Schemes
Algorithm Implementation and Validation
The auto-encoder algorithm for lincRNA detection is implemented on open source Python libraries, Theano and Keras. The training and validation data sets including the known lincRNA data are collected from UCSC Genome Browser database. The existing methods, including NNSplice [25] and Libsvm [26], are used for validating the proposed deep learning method by the comparisons with traditional Neural Network and Support Vector Machine.
According to the latest findings [4], totally 46,983 lincRNA sequences containing 90 nucleotides and 89,287 lincRNA sequences containing 15 nucleotides are extracted and collected as transcriptional sites, Acceptors and Donors respectively, including 5,000 sequences as validation in each data set. Based on the auto-encoder algorithm, a 2-layer neural network is constructed for the experiments. Five aforementioned encoding schemes are used for comparing and acquiring the best performance.
Results
Tables 1 and 2 respectively show the comparison results for the two data sets. It shows that 100% predictive rate of deep neural network method with complementary encoding scheme on the acceptor data, meaning that complementary scheme has the strong ability on more-feature data sets. Similar performances among all encoding schemes show the similar ability on less-feature data set.
Moreover, we compare the deep learning method with Support Vector Machine (SVM) using the same data sets. SVM software is tested on the latest version of libsvm [26]. Figures 4 and 5 show the comparative results that auto-encoder based deep learning method has an extraordinary ability over conventional SVM method. On the data set with more features in Fig. 4, the deep learning method shows the large superiority over SVM while their performances are very close on the data set with less features in Fig. 5.
Comparison between Support Vector Machine and Deep Learning on lincRNA Acceptor Data Set
Comparison between Support Vector Machine and Deep Learning on lincRNA Donor Data Set
In addition, a comparison between the deep learning method and the traditional neural network (NN) based method [25] is also conducted. Figures 6 and 7 show that DL outperforms the conventional NN based method for detection of transcriptional sites using lincRNA data sets. Similarly, on the data set with more features in Fig. 6, the deep learning method distinguishes itself from the NN based method while their performances are very close on the data set with less features in Fig. 7. It means that various methods have the similar performance on handling the less-feature data set while deep learning can have a large superiority over others on processing the more-feature data set. Such experimental results also manifest that deep learning based method can have better performance than other conventional methods for prediction of lincRNAs on DNA sequence data. The reason that we separate the comparison between SVM-DL group and NN-DL group is that the SVM tool we use for the experiment can accept all encoding schemes as its input while the NN-based web tool accepts only the DNA sequence as its input.
Comparison between Conventional Neural Network Method and Deep Learning Method on lincRNA Acceptor Data Set
Comparison between Conventional Neural Network Method and Deep Learning Method on lincRNA Donor Data Set
Figure 8 shows an unreported splicing site is found by re-scanning the whole human genome through the deep learning method, which is located at 90,763,154 chromosome 12 (hg38) within the annotated lincRNA chr12_90761911_90806776. This result is based on the aforementioned deep learning method that was tested with 100% accuracy on acceptor data set.
An Unidentified lincRNA Acceptor Site
Discussion
Although a deep learning based method has been illustrated for lincRNA detection, distinguishing the coding and non-coding transcription is still an open problem because the transcribed regions have the similar structures of exon and intron in both coding and non-coding regions. Practically, it is hard to find an effective way to differentiate the two types of transcripts. Thus, the intergenic regions have to be selected and the pre-processing is necessary for detection, which is the downside of our method and partially limits the use of the proposed deep learning based method.
In addition, the development of deep learning method for lincRNA detection is still on preliminary stage and the prototype of the auto-encoder based method has more spaces to improve. For example, function modules need to be uniformed and parameters in the work flow have to be optimized.
Conclusion
RNA-seq technologies generate a large volume of transcriptional data that scientists can utilize for lincRNA annotation. Derived from the observations from the newly found lincRNA data set, two evidences can provide the reasoning for adopting knowledge-based deep learning methods in lincRNA detection: (1) the expression of lincRNAs appears to be regulated, indicating that the relevance exists along the DNA sequences; (2) lincRNAs contain some conversed patterns/motifs tethered together by non-conserved regions [9]. In this project, a knowledge-based discovery method using the emerging deep learning technology for lincRNA detection is proposed and developed on DNA genome analysis. It takes advantage of the latest findings of lincRNA data set and aims to utilize the cutting-edge knowledge-based method, namely auto-encoder algorithm, in order to extract the features of lincRNA transcription sites in a more accurate way than conventional methods. The results show its superiority over the support vector machine and the conventional neural network based method.
In the future, developing a generic framework based on deep learning for lincRNA prediction will be focused on, which can provide an uniform platform for user interfaces. Meanwhile, the studies on lincRNA detection will be carried out on other species such as mouse and other mammals.
References
Mercer TR, Gerhardt DJ, Dinger ME, Crawford J, Trapnell C, Jeddeloh JA, Mattick JS, Rinn JL. Targeted rna sequencing reveals the deep complexity of the human transcriptome. Nat Biotechnol. 2012; 30:99–104.
Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, Suzuki H, Carninci P, Hayashizaki Y, Wells C, Frith M, Ravasi T, Pang KC, Hallinan J, Mattick J, Hume DA, Lipovich L, Batalov S, Engström PG, Mizuno Y, Faghihi MA, Sandelin A, Chalk AM, Mottagui-Tabar S, Liang Z, Lenhard B, Wahlestedt C. Antisense transcription in the mammalian transcriptome. Science. 2005; 309(5740):1564–6. doi:10.1126/science.1112009.
Durruthy-Durruthy J, Sebastiano V, Wossidlo M, Cepeda D, Cui J, Grow EJ, Davila J, Mall M, Wong WH, Wysocka J, Au KF, Reijo Pera RA. The primate-specific noncoding rna hpat5 regulates pluripotency during human preimplantation development and nuclear reprogramming. Nat Genet. 2016; 48(1):44–52.
Hangauer MJ, Vaughn IW, McManus MT. Pervasive transcription of the human genome produces thousands of previously unidentified long intergenic noncoding rnas. PLoS Genet. 2013; 9(6):1–13. doi:10.1371/journal.pgen.1003569.
Luo H, Bu D, Sun L, Fang S, Liu Z, Zhao Y. Identification and function annotation of long intervening noncoding rnas. Brief Bioinform. 2016. doi:10.1093/bib/bbw046.
Kapranov P, Cheng J, Dike S, Nix DA, Duttagupta R, Willingham AT, Stadler PF, Hertel J, Hackermüller J, Hofacker IL, Bell I, Cheung E, Drenkow J, Dumais E, Patel S, Helt G, Ganesh M, Ghosh S, Piccolboni A, Sementchenko V, Tammana H, Gingeras TR. Rna maps reveal new rna classes and a possible function for pervasive transcription. Science. 2007; 316(5830):1484–8. doi:10.1126/science.1138341.
Xuan G, Ning Y, Xiaojun D, Jianxin W, Yi P. Dime: A novel framework for de novo metagenomic sequence assembly. J Comput Biol. 2015; 22(2):159–77.
Birol I, Jackman SD, Nielsen CB, Qian JQ, Varhol R, Stazyk G, Morin RD, Zhao Y, Hirst M, Schein JE, Horsman DE, Connors JM, Gascoyne RD, Marra MA, Jones SJM. De novo transcriptome assembly with abyss. Bioinformatics. 2009; 25(21):2872–7. doi:10.1093/bioinformatics/btp367.
Ulitsky I, Shkumatava A, Jan CH, Sive H, Bartel DP. Conserved function of lincrnas in vertebrate embryonic development despite rapid sequence evolution. Cell. 2011; 147(7):1537–50.
Sati S, Ghosh S, Jain V, Scaria V, Sengupta S. Genome-wide analysis reveals distinct patterns of epigenetic features in long non-coding rna loci. Nucleic Acids Res. 2012; 40(20):10018–31. doi:10.1093/nar/gks776.
Ponjavic J, Ponting CP, Lunter G. Functionality or transcriptional noise? evidence for selection within long noncoding rnas. Genome Res. 2007; 17(5):556–65. doi:10.1101/gr.6036807.
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, Lagarde J, Veeravalli L, Ruan X, Ruan Y, Lassmann T, Carninci P, Brown JB, Lipovich L, Gonzalez JM, Thomas M, Davis CA, Shiekhattar R, Gingeras TR, Hubbard TJ, Notredame C, Harrow J, Guigó R. The gencode v7 catalog of human long noncoding rnas: Analysis of their gene structure, evolution, and expression. Genome Res. 2012; 22(9):1775–89. doi:10.1101/gr.132159.111.
Hinton G, Dayan P, Frey B, Neal R. The “wake-sleep” algorithm for unsupervised neural networks. Science. 1995; 268(5214):1158–61.
Hintonemail GE. Learning multiple layers of representation. Trends Cogn Sci. 2007; 11(10):428–34.
Deng L, Hinton G, Kingsbury B. New types of deep neural network learning for speech recognition and related applications: an overview. In: Acoustics, Speech and Signal Processing (ICASSP), 2013 IEEE International Conference On: 2013. p. 8599–603. doi:10.1109/ICASSP.2013.6639344.
Bengio Y, Courville A, Vincent P. Representation learning: A review and new perspectives. IEEE Trans. Pattern Anal Mach Intell. 2013; 35(8):1798–828.
Di Lena P, Nagata K, Baldi P. Deep architectures for protein contact map prediction. Bioinformatics. 2012; 28(19):2449–57. doi:10.1093/bioinformatics/bts475.
Eickholt J, Cheng J. Predicting protein residue-residue contacts using deep networks and boosting. Bioinformatics. 2012; 28(23):3066–72. doi:10.1093/bioinformatics/bts598.
Leung MKK, Xiong HY, Lee LJ, Frey BJ. Deep learning of the tissue-regulated splicing code. Bioinformatics. 2014; 30(12):121–9. doi:10.1093/bioinformatics/btu277.
Yu N, Guo X, Gu F, Pan Y. DNA AS X: An information-coding-based model to improve the sensitivity in comparative gene analysis. In: Bioinformatics Research and Applications: 11th International Symposium, ISBRA 2015 Norfolk, USA, June 7-10, 2015 Proceedings.Cham: Springer International Publishing: 2015. p. 366–377.
Nair AS, Sreenadhan SP. A coding measure scheme employing electron-ion interaction pseudopotential (EIIP). Bioinformation. 2006; 1(6):197–202.
Akhtar M, Epps J, Ambikairajah E. Signal processing in sequence analysis: Advances in Eukaryotic gene prediction. IEEE J Sel Top Signal Process. 2008; 2(3):310–21.
Kauer G, Blöcker H. Applying signal theory to the analysis of biomolecules. Bioinformatics. 2003; 19(16):2016–21. doi:10.1093/bioinformatics/btg273. http://bioinformatics.oxfordjournals.org/content/19/16/2016.full.pdf+html.
Rosen GL. Signal processing for bibiological-inspired gradient source localization and dna sequence analysis. PhD thesis, Georgia Institute of Technology, School of Electrical and Computer Engineering. 2006.
Reese MG, Eeckman FH, Kulp D, Haussler D. Improved splice site detection in genie. J Comput Biol. 1997; 4(3):311–323.
Chang CC, Lin CJ. LIBSVM: A library for support vector machines. ACM Trans Intell Syst Technol. 2011; 2:27–12727.
Acknowledgements
We give thanks to the supports from the Department of Computing Sciences, SUNY Brockport.
Funding
Publication costs were funded by the College at Brockport, State University of New York.
Availability of data and materials
All genome information and annotated data are collected from UCSC Genome Browser database (https://genome.ucsc.edu/). The software and source code can be downloaded from GitHub (https://github.com/ningyu12/lincRNA_predict/).
About this supplement
This article has been published as part of BMC Bioinformatics Volume 18 Supplement 15, 2017: Selected articles from the 6th IEEE International Conference on Computational Advances in Bio and Medical Sciences (ICCABS): bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-18-supplement-15.
Authors’ contributions
NY designs the experiments and performs the implementation; ZY conducts the literature review and the theoretical design; YP coordinates the project and provides the significant advice on the method design. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Yu, N., Yu, Z. & Pan, Y. A deep learning method for lincRNA detection using auto-encoder algorithm. BMC Bioinformatics 18 (Suppl 15), 511 (2017). https://doi.org/10.1186/s12859-017-1922-3
Published:
DOI: https://doi.org/10.1186/s12859-017-1922-3