Portal de Periódicos da CAPES

AptaSUITE: A Full-Featured Bioinformatics Framework for the Comprehensive Analysis of Aptamers from HT-SELEX Experiments

2018; Cell Press; Volume: 11; Linguagem: Inglês

10.1016/j.omtn.2018.04.006

2162-2531

Jan Hoinka, Rolf Backofen, Teresa M. Przytycka,

DNA and Nucleic Acid Chemistry

To the editor: The capability of producing and efficiently processing big data has revolutionized virtually every field of science and technology and has enabled the analysis of experimental results at unprecedented resolutions. This trend is also evidenced in the rapid emergence, and subsequent field-wide adoption, of the high-throughput systematic evolution of ligands by exponential enrichment (HT-SELEX) protocol in the study of in vitro selection.1Blind M. Blank M. Aptamer Selection Technology and Recent Advances.Mol. Ther. Nucleic Acids. 2015; 4: e223Abstract Full Text Full Text PDF PubMed Scopus (217) Google Scholar HT-SELEX extends the traditional SELEX protocol, aimed at the generation of high-affinity and specificity oligonucleotides known as aptamers against a molecular target of interest, by coupling this technology with HT sequencing. Selection is typically performed in iterations and consists of incubating an initially random pool of sequences with the target, followed by the partitioning, and subsequent removal, of non-affine species while amplifying the remaining pool to form the input to the next cycle as well as the source material for sequencing. The resulting sequencing data, consisting of a representative sample of the pool composition after each round of selection, are consequently analyzed in silico through dedicated algorithmic approaches, whereby aptamers predicted to possess the desired application-specific properties are typically subjected to further in vitro verification and post-processing. Notably, in order to guarantee an efficient and accurate in silico pipeline, these computational methods must be carefully designed to maximize efficiency while scaling well vertically (guaranteeing a proportional reduction in computation time with a growing number of available processing units and memory) as well as with increasing data volume. Ideally, such tools would additionally require a low learning curve for experimentalists, be platform-independent, and provide integrated means of storing and retrieving, interacting with, and visualizing aptamer-related information. Indeed, over the past decade, typical HT-SELEX datasets have grown 200-fold from 10,000–100,000 reads per selection cycle to current sizes of routinely over 20 million reads per round. This development constitutes an emerging barrier for many well-established algorithmic tools devised before the big data revolution but which are still actively used in RNA bioinformatics analysis pipelines. Prominent examples include clustering sequence species into aptamer families related to each other by sequence similarity2Alam K.K. Chang J.L. Burke D.H. FASTAptamer: A bioinformatic toolkit for high-throughput sequence analysis of combinatorial selections.Mol. Ther. Nucleic Acids. 2015; 4: e230Abstract Full Text Full Text PDF PubMed Scopus (140) Google Scholar, 3Bailey T.L. Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers.Proc Int Conf Intell Syst Mol Biol. 1994; 2: 28-36PubMed Google Scholar, 4Hiller M. Pudimat R. Busch A. Backofen R. Using RNA secondary structures to guide sequence motif finding towards single-stranded regions.Nucleic Acids Res. 2006; 34: e117Crossref PubMed Scopus (117) Google Scholar and the elucidation of shared motifs in primary and/or secondary structure evolving throughout the selection. Both analysis tasks are well established for small-scale datasets but rapidly become computationally intractable with increasing data volume when performed with traditional methods.5Hoinka J. Zotenko E. Friedman A. Sauna Z.E. Przytycka T.M. Identification of sequence-structure RNA binding motifs for SELEX-derived aptamers.Bioinformatics. 2012; 28: i215-i223Crossref PubMed Scopus (66) Google Scholar, 6Caroli J. Taccioli C. De La Fuente A. Serafini P. Bicciato S. APTANI: a computational tool to select aptamers through sequence-structure motif analysis of HT-SELEX data.Bioinformatics. 2016; 32: 161-164PubMed Google Scholar More complex approaches, designed for scalable, HT data processing on multi-core environments, as found in data centers and cloud environments, typically require expert knowledge to set up a sensible pipeline and may depend on numerous, potentially non-portable, third-party software packages, increasing the burden of long term maintainability.7Hoinka J. Berezhnoy A. Sauna Z.E. Gilboa E. Przytycka T.M. AptaCluster - A method to cluster HT-SELEX aptamer pools and lessons from its application.Res Comput Mol Biol. 2014; 8394: 115-128Crossref PubMed Scopus (57) Google Scholar, 8Hoinka J. Berezhnoy A. Dao P. Sauna Z.E. Gilboa E. Przytycka T.M. Large scale analysis of the mutational landscape in HT-SELEX improves aptamer discovery.Nucleic Acids Res. 2015; 43: 5699-5707Crossref PubMed Scopus (75) Google Scholar, 9Hoinka J. Dao P. Przytycka T.M. AptaGUI-A Graphical User Interface for the Efficient Analysis of HT-SELEX Data.Mol. Ther. Nucleic Acids. 2015; 4: e257Abstract Full Text Full Text PDF PubMed Scopus (14) Google Scholar In addition, the resulting processed data are predominantly output in pure text format or stored in relational databases, adding to the stack of challenges in efficiently interpreting the results. Finally, while undoubtedly being of great utility, web-based (and therefore graphical) solutions backed by cloud services such as the Galaxy project10Afgan E. Baker D. van den Beek M. Blankenberg D. Bouvier D. Čech M. Chilton J. Clements D. Coraor N. Eberhard C. et al.The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update.Nucleic Acids Res. 2016; 44: W3-W10Crossref PubMed Scopus (1239) Google Scholar, 11Bair T. Wertz J.S. Miller R.J. Schappet J. Giangrande P.H. Thiel W.H. 148. AptaMetrics: A Web-Based Aptamer Bioinformatics Platform.Mol. Ther. 2018; 24: S59Abstract Full Text Full Text PDF Google Scholar are limited in their flexibility of visualizing and interacting with vast amounts of data as they must adhere to the constraints imposed by current web browsers and technologies. To address these issues, we have developed AptaSuite, a full-featured, open source, and platform-independent software collection for the comprehensive analysis of HT-SELEX experiments. In stark contrast to previous methods, each implementing their individual and frequently rudimentary data workflow, AptaSuite provides a unified and robust framework for managing aptamer-related data and leverages this framework to serve the required data in a standardized manner to any particular algorithm built with the software. In its core, AptaSuite consists of a collection of carefully designed APIs (application programming interfaces) and corresponding reference implementations for facilitating input, output, and manipulation of aptamer data (such as sequences, aptamer counts in individual selection cycles, structure information, and more). On top of this powerful core library, a number of previously published approaches8Hoinka J. Berezhnoy A. Dao P. Sauna Z.E. Gilboa E. Przytycka T.M. Large scale analysis of the mutational landscape in HT-SELEX improves aptamer discovery.Nucleic Acids Res. 2015; 43: 5699-5707Crossref PubMed Scopus (75) Google Scholar, 12Hoinka J. Przytycka T. AptaPLEX - A dedicated, multithreaded demultiplexer for HT-SELEX data.Methods. 2016; 106: 82-85Crossref PubMed Scopus (12) Google Scholar, 13Dao P. Hoinka J. Takahashi M. Zhou J. Ho M. Wang Y. Costa F. Rossi J.J. Backofen R. Burnett J. Przytycka T.M. AptaTRACE Elucidates RNA Sequence-Structure Motifs from Selection Trends in HT-SELEX Experiments.Cell Syst. 2016; 3: 62-70Abstract Full Text Full Text PDF PubMed Scopus (38) Google Scholar, 14Fukunaga T. Ozaki H. Terai G. Asai K. Iwasaki W. Kiryu H. CapR: revealing structural specificities of RNA-binding protein target recognition using CLIP-seq data.Genome Biol. 2014; 15: R16Crossref PubMed Scopus (54) Google Scholar, 15Lorenz R. Bernhart S.H. Höner Zu Siederdissen C. Tafer H. Flamm C. Stadler P.F. Hofacker I.L. ViennaRNA Package 2.0.Algorithms Mol. Biol. 2011; 6: 26Crossref PubMed Scopus (2542) Google Scholar have been implemented from scratch and are now combined into this uniform, easy-to-use framework (see Figure 1). In particular, the selected methods constitute well-established approaches to analyze HT-SELEX data and are specifically designed to leverage particular properties of aptamers and the SELEX process. These algorithms engage with the analysis process at the earliest possible moment by providing importing capabilities and quality control of raw sequencing data through AptaPLEX, a multithreaded demultiplexer and parser for HT-SELEX data.12Hoinka J. Przytycka T. AptaPLEX - A dedicated, multithreaded demultiplexer for HT-SELEX data.Methods. 2016; 106: 82-85Crossref PubMed Scopus (12) Google Scholar Compared to generic demultiplexers, AptaPLEX utilizes the primary structure configuration of aptamer reads to increase the total number of recovered oligonucleotides and effectively partition the raw data into the corresponding selection cycles based on barcoding information contained within the reads. Alternatively, by invoking our SELEX simulator AptaSIM, sensible recreations of the selection process, incorporating features such as error-prone amplification, target affine selection, and aliquotation of the pool after each round, can be created and explored with AptaSuite. This in turn allows for the analysis of higher-order relationships regarding the selection pressures governing a particular experimental setup.8Hoinka J. Berezhnoy A. Dao P. Sauna Z.E. Gilboa E. Przytycka T.M. Large scale analysis of the mutational landscape in HT-SELEX improves aptamer discovery.Nucleic Acids Res. 2015; 43: 5699-5707Crossref PubMed Scopus (75) Google Scholar A common first analysis step, after the initial data import, consists of grouping the sequences into aptamer families related to each other by primary structure, which evolved due to a combination of (hyper)mutagenic amplification followed by subsequent selection. To efficiently perform this clustering, and to trace the evolution of the resulting families throughout the selection, we have ported our previously published approach AptaCLUSTER into this new framework.8Hoinka J. Berezhnoy A. Dao P. Sauna Z.E. Gilboa E. Przytycka T.M. Large scale analysis of the mutational landscape in HT-SELEX improves aptamer discovery.Nucleic Acids Res. 2015; 43: 5699-5707Crossref PubMed Scopus (75) Google Scholar AptaCLUSTER leverages the constant size of the randomized region in conjunction with locality sensitive hashing to outperform traditional clustering algorithms, which rely on expensive all-versus-all comparison techniques for their operations. This, and its multi-threaded architecture therefore enable clustering of datasets stemming from next-generation (and future) sequencing technologies in a computationally tractable manner. Building on the results of AptaCLUSTER, we additionally provide the ability to perform an in-depth analysis of the mutational landscape within aptamer families via AptaMUT,8Hoinka J. Berezhnoy A. Dao P. Sauna Z.E. Gilboa E. Przytycka T.M. Large scale analysis of the mutational landscape in HT-SELEX improves aptamer discovery.Nucleic Acids Res. 2015; 43: 5699-5707Crossref PubMed Scopus (75) Google Scholar based on a theoretic model capable of discriminating favorable mutants from those that decrease the binding affinity to the target. Another challenge paramount to the analysis of in vitro selection data centers on the identification of sequence-structure patterns shared among target-affine aptamers, which are responsible for the binding interaction between the species and target. In AptaSuite, this analysis is made possible through the integration of AptaTRACE,13Dao P. Hoinka J. Takahashi M. Zhou J. Ho M. Wang Y. Costa F. Rossi J.J. Backofen R. Burnett J. Przytycka T.M. AptaTRACE Elucidates RNA Sequence-Structure Motifs from Selection Trends in HT-SELEX Experiments.Cell Syst. 2016; 3: 62-70Abstract Full Text Full Text PDF PubMed Scopus (38) Google Scholar our computational approach that leverages the experimental design of the HT-SELEX protocol, RNA secondary structure information, and the potential presence of many secondary motifs to identify sequence-structure elements that show a signature of selection. Finally, every aspect of the data and analysis results can be exported back to file in well-established formats for further downstream processing using third-party pipelines. To maximize the applicability of AptaSuite to as many use cases as possible, every feature of our approach is made available to the user as both command line interfaces (CLIs) and graphical user interfaces (GUIs). The former can additionally be chained together in a variety of combinations to create seamless data processing pipelines. This allows AptaSuite to be deployed as fully automated applications in high-performance computing (HPC) environments, as point-and-click solutions for small- to mid-sized datasets on modern desktop hardware, or as a combination of the two, in which the heavy computation is outsourced to HPC systems whereas visualization and interpretation of the results can be performed locally using the graphical user interface (Figure 2). This flexibility of AptaSuite is mainly attributed to our solution being implemented in pure Java without requiring any third-party dependencies, making this software truly platform independent and portable across a large array of hardware environments. In addition, AptaSuite is designed to be highly scalable with both, data size and CPU count, while minimizing the memory footprint by providing fast, off-heap data structures and storage solutions. Finally, its modular design and well-documented APIs allow for trivial extension with new algorithmic solutions as research in the field of in vitro selection progresses and novel questions calling for computer-aided problem solving arise. To the best of our knowledge, AptaSuite represents the most comprehensive data engine to store and retrieve, manipulate, and analyze aptamer data to date. Our software is already being used by the aptamer community at large, and we are confident that AptaSuite will continue to establish itself as one of the de facto in silico analysis tools in the field of oligonucleotide in vitro selection. AptaSuite is an open source and the precompiled binaries, source code, and manual are available at https://github.com/drivenbyentropy/aptasuite. This work was supported by the Intramural Research Program of the NIH, National Library of Medicine. Funding for the open access charge was provided by the NIH. This work utilized the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov/).