[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Oxidative Stress and Reproductive Function in the Aging Male
Previous Article in Journal
Factors Associated with Increased Morbidity and Mortality of Obese and Overweight COVID-19 Patients
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 9
Issue 9
10.3390/biology9090281
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets

by
Simon Orozco-Arias
1,2,*,†,
Nicolas Tobon-Orozco
1,†,
Johan S. Piña
1,
Cristian Felipe Jiménez-Varón
3,
Reinel Tabares-Soto
4 and
Romain Guyot
4,5,*
1
Department of Computer Science, Universidad Autónoma de Manizales, Manizales 170002, Colombia
2
Department of Systems and Informatics, Universidad de Caldas, Manizales 170002, Colombia
3
Department of Physics and Mathematics, Universidad Autónoma de Manizales, Manizales 170002, Colombia
4
Department of Electronics and Automation, Universidad Autónoma de Manizales, Manizales 170002, Colombia
5
Institut de Recherche pour le Développement (IRD), CIRAD, Université de Montpellier, 34394 Montpellier, France
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally.
Biology 2020, 9(9), 281; https://doi.org/10.3390/biology9090281
Submission received: 27 August 2020 / Revised: 1 September 2020 / Accepted: 7 September 2020 / Published: 9 September 2020
(This article belongs to the Section Bioinformatics)
Browse Figures
Figure 1
<p>TIP_finder methodology and schematic representation of the pipeline. TE: transposable element, TIP: transposon insertion polymorphism.</p> "> Figure 2
<p>Flowchart of the parallel strategy implemented in TIP_finder.</p> "> Figure 3
<p>Total average runtime and speedup of TIP_finder -using NCBI-BLAST and MagicBLAST- and TRACKPOSON using 2, 4, 8, 16, 32, 44, and 56 cores with a randomly selected case dataset (30 million of reads) and executed 10 times (<b>A</b>,<b>B</b>), and a randomly selected control dataset (26.5 million of reads) and executed 10 times (<b>C</b>,<b>D</b>). The times for all executions can be found in <a href="#app1-biology-09-00281" class="html-app">Supplementary Material S2 Table S1</a> (for case dataset), and in <a href="#app1-biology-09-00281" class="html-app">Table S2</a> (for control dataset).</p> "> Figure 4
<p>TIP_finder speedup of each step with 32 cores for one case and one control datasets.</p> "> Figure 5
<p>Total average runtime and speedup of TIP_finder and TRACKPOSON using six randomly selected control datasets and six randomly selected case datasets, each one executed 10 times. (<b>A</b>) Comparison of runtimes between TRACKPOSON and TIP_finder (with B: NCBI-BLAST, and M: Magic-BLAST), and (<b>B</b>) the comparison of the speedup of TIP_finder (NCBI-BLAST and MagicBLAST) compared to TRACKPOSON. Additional information can be consulted in <a href="#app1-biology-09-00281" class="html-app">Table S7</a>.</p> "> Figure 6
<p>Number of TIPs insertions for the cases and controls identified by TIP_finder and grouped into bins of 500. Control and case patients are show in red, and in blue respectively. (<b>A</b>) Distribution of the number of TIPs insertions in cases and controls shown in blue and red, respectively. (<b>B</b>) Distribution of TIPs in controls within a range of 140–220.</p> "> Figure 7
<p>Distribution of TIPs insertions (Human endogenous retroviruses type K (HERV-K) insertions) for the cases and controls. The X axis represents the number of patients in log2 scale and the Y axis is the TIPs insertion number. (<b>A</b>,<b>B</b>) correspond to the cases and controls, respectively. Variations in the peaks’ frequency suggest a change in the insertional activity under a condition of interest, in this case cancer.</p> "> Figure 8
<p>Distribution of the number of TIPs—HERVs—along human chromosomes. The X axis represents the chromosome length (on a scale of 1 × 10<sup>8</sup>), where cases and controls are shown in blue and red, respectively. The Y axis represents the number of TIPs along the chromosome length. The arrows show the end of each chromosome. The graphs for all chromosomes are available in <a href="#app1-biology-09-00281" class="html-app">Supplementary Material S4</a>.</p> "> Figure 9
<p>Number of TIPs statistically associated with breast cancer.</p> ">
Versions Notes

Abstract

:
Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era.

1. Introduction

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another [1,2,3]. These mobile elements can accumulate large copy numbers in their host genomes [4] and have been found in all organisms. The majority of the nuclear DNA content of large genomes is composed of TEs, such as in wheat, barley, and maize [5,6,7] for plants. In humans, these elements (or TE-derived sequences) comprise ~50–70% of the sequenced genome [8]. Several studies have indicated that TEs play crucial genomic roles involved in chromosome structuring, structural variation, the alteration of gene expression [5,7], evolution, the variation of genomic size, and environmental adaptation [9,10,11,12,13]. Nevertheless, these elements can also be associated with human diseases, such as different types of cancer [14,15,16]. TEs are classified into two major classes depending on their replication modes [17]. Accordingly, Class I or retrotransposons use an RNA molecule as an intermediate, while Class II or DNA transposons utilize a DNA intermediate. Each class can be hierarchically sub-classified into orders, super-families, lineages, and families [9,18,19]. Among retrotransposons, an order called LTR-retrotransposons bears structural and functional similarities to retroviruses, including the presence of long terminal repeats (LTR) at both ends that flank central coding domains and a similar replication cycle [20,21].
A particular kind of retrovirus present in the human genome is known as the human endogenous retrovirus (HERVs) and makes up approximatively 8% of the DNA sequenced [22]. HERVs carry genes (gag, pro, pol, and env) encoding essential proteins at the functional and structural levels [12,23]. These genes show overexpression in Mendelian diseases [24] and the etiology of cancer, including breast cancer [25,26], testicular cancer [27], melanoma [28], and germ cell tumors [29]. This overexpression could support the hypothesis that HERVs can be used in comparative studies to find polymorphisms caused by the insertion of these retroelements. Other TEs have been used in comparative studies to find transposon insertion polymorphisms (TIPs) associated with certain biological phenomena, such as the adaptation processes of rice [30], paramutation and gene silencing [31], and gene expression regulation in several organisms [32,33,34].
There are several computational tools to detect TIPs and most of them are based on the discordant location of read pairs and split reads [35,36], such as iTIS_PseTNC [37], Jitterbug [38], Transposeq [39], Metasv [40], DD_Detection [41], and TRACKPOSON [30]. The latter was the most recently developed pipeline for the detection of TIPs in gene pools and it was designed to unravel the transpositional activity of TEs in genomic datasets by applying a faster and validated methodology. Although the bioinformatic programs used by TRACKPOSON can be executed in multiple processors, this pipeline does not use any parallel strategies and as a consequence BLAST is executed in a non-efficient way and the BLAST output is processed serially (i.e., on one processor). In addition, this program cannot be run in more than one computational node (or computer) and has issues with the amount of disk space required. Thus, this pipeline is not scalable to current supercomputers (with hundreds of servers) and remains suboptimal for analyzing large genomes with a substantial amount of resequencing data, such as Arabica coffee with 1.3 Gb [42,43], maize (Zea mays) with 2.3 Gb [44], pine tree (Pinus taeda) with 22 billion base pairs [45], and the human genome with 3.3 Gb [46], as well as large datasets (e.g., 10K plant genomes [47] and the Earth BioGenome [48]) and huge databases such as the NCBI genome collection with sequence assemblies from more than 13,000 different organisms [49].
Therefore, it is important to take advantage of supercomputing, parallel programming [50], and high-performance computing (HPC) approaches to speed up the bioinformatics analyses of large genomes and huge datasets [51,52,53]. These techniques have been applied in several tasks, such as the hierarchical clustering of nucleotide sequences [54], scaled BLAST using CPU [52] or Xeon Phi [53,55] architectures, and a pipeline to analyze LTR-retrotransposons [2]. Message passing interface (MPI) [56] is a well-known library for parallel programming [57], which is used to run many sub-problems that are previously divided given different focuses [58].
Here, we present TIP_finder, a pipeline that can be applied to analyze large genomes and large resequencing datasets to discover TIPs using the classical methodology of inference of discordant read pairs as proposed by several pipelines. TIP_finder, written in python, proposes different BLAST engines (NCBI or MagicBLAST) and works under HPC techniques and parallel programming. It can be scaled on many computational nodes (or servers) and multi-core architectures, which makes it especially functional for applications in massive sequencing projects in the current (post) genomic era [24,34,59]. To test the performance of TIP_finder in a challenging way, we use the very large human genome and available breast cancer datasets, for detecting the insertion of HERVs in relation with this disease. In addition, we developed different utilities that automatically perform analyses, such as graphs over TIPs frequencies, peak charts, and statistical associations analyses related to a certain condition. TIP_finder version 1.0 is publicly available at https://github.com/simonorozcoarias/TIP_finder.

2. Materials and Methods

2.1. Implementation of TIP_finder

TIP_finder follows the strategy of the analysis of discordant read pair as proposed by several algorithms such as by [30] to detect TIPs using (i) short pair reads such as those obtained by Illumina, (ii) a reference genome sequence assembled at the chromosome level, and a database of consensus TE sequences (or other mobile elements like HERVs). With available programs, we identified several bottlenecks dramatically reducing the speed of the analysis and increasing the size of the output files such as: the execution time of the NCBI-BLAST program, the writing time and format of the BLAST results, the absence of parallelization or the non-use of a robust programming language.
The first step consisted of mapping paired reads to a mobile elements database and filtering only the reads showing one of the paired reads mapped. Then, the “unmapped reads” of read pairs are aligned to the reference genome to detect the insertion point of the putative TIP. The second step of TIP_finder processes the alignment output (saved as tabular format and only keeping the following columns: subject sequence id, start and end positions of the alignment in the subject sequence, and the ID of the query sequence) to search for reads with only one hit. This step produces a tabular file (in bed format) that will be used in the third step. Finally, TIP_finder counts how many reads were aligned to the reference genome, previously divided by a 10-Kb window, and generates a tabular-delimited file that is processed using the TIP_finder_utils.py. A general scheme of TIP_finder can be consulted in Figure 1.
TIP_finder_utils.py has four utilities to analyze and process the data generated, as follows: (1) cluster the results in a presence/absence (coded by 1/0) matrix of TIPs for each paired-end sequence file; (2) generate frequency graphs from the number of TIPs for each analyzed sample; (3) make peak charts using the frequency for TIPs present along the entire chromosome length; and (4) perform statistical association tests to identify TIPs displaying the highest association rate with a specific condition.
TIP_finder used the following bioinformatic software: Bowtie2 (v. 2.3.4.1) to map the paired reads of genomic data to the indexed consensus sequence of each TE/HERV family [30,60], Samtools (v. 1.9) to process the bowtie2 output and keep only the unmapped reads [61], bedtools (v. 2.26.0) to split the reference genome into 10 Kb windows and count reads in these windows [62], and NCBI-BLASTn (v. 2.10.0) [63] or Magic-BLAST [64] to align the unmapped reads. Although Magic-BLAST is significantly faster than NCBI-BLASTn, it requires more disk space. TIP_finder was developed using Python3 (3.8) and the following libraries: sys, time, os, subprocess, argparse, and MPI4py [65]. Furthermore, TIP_finder_utils.py requires the additional libraries: math, Pandas [66], matplotlib [67], Seaborn [68], and SciPy [69].

2.2. Parallel Strategy Implemented

To take advantage of multi-core architectures, TIP_finder is parallelized using the MPI standard and executes the same process several times with different data [51]. First, the read files (in FASTQ format) are split into n files (where n is the number of processes); then, each working process maps reads against the TE consensus databases. The mapping output is used to extract the “unmapped reads”, which are aligned against the reference genome. Each process generates a BLAST output file and the master process finally joins all of these into one file. The next step includes the detection of reads with a unique hit in the genome. For this, the number of lines in the BLAST output is divided by the processes number, and a start and end line number is assigned to each work process; the process analyzes its corresponding region and creates a Python dictionary, which is unified by the master process. Finally, each working process extracts reads with a unique hit and writes them to a file in bed format. The master process joins all the partial results, filters the final file, and deletes all the temporary files. An overall flowchart of this parallel strategy can be consulted in Figure 2.

2.3. Statistical Association Analysis

The statistical association tests performed using TIP_finder_utils.py consist of establishing two categorical random variables: X which determines if a patient i (i: 1,...,n) has a disease (cancer) and Y, which determines if a patient i has factor (TIP) j (j: 1,...,m). Given this, n corresponds to the number of patients (cases and controls) and m is the number of TIPs. The possible values associated with the aforementioned random variables are X = 0, 1 where X = 0 means that patient i does not have the disease and X = 1 means that the patient i has the disease, and Y = 0, 1 where Y = 0 means that patient i does not have the factor j and Y = 1 means that the patient i has factor j [70].
Based on a sample of n patients and m factors, the statistical association analysis is performed based on the chi-square distribution. Contingency tables are generated by crossing each of the m factors individually with the probability of having the disease or not. A model of the contingency table used is presented in Table 1.
The Oij values refer to the frequency observed in the data for patient i and the presence or absence of factor j.
Based on the values found in the contingency table, the totals by rows and columns are calculated to determine the expected frequency of each of the Oij observed frequencies. In order to determine whether factor j has a statistically significant influence on having the disease, we decided to propose a non-parametric independence hypothesis test based on the chi-square distribution. This test has the structure presented in (1) and (2) [71].
H 0 :   V a r i a b l e s   X   a n d   Y   a r e   i n d e p e n d e n t .
H 1 :   V a r i a b l e s   X   a n d   Y   a r e   a s s o c i a t e d .
The chi-square statistics for these tests was calculated as follows:
χ 2 S q u a r e d = ∑ i = 1 r ∑ j = 1 s ( O i j − e i j ) 2 e i j ,
where Oij is the observed frequency, eij is the expected frequency, and r, s are the row and column sizes, respectively, which in this case are both worth two. The expected frequency eij for each observed frequency Oij is calculated as follows:
e i j = n i . n . j n ,
where n i .   refers to the total of row i, n . j refers to the total of column j, and n refers to the total of the joint data in the contingency table. Once the contrast statistic is determined χ 2 c a l c u l a t e d , a rejection criterion is proposed based on a significance level α and the degrees of freedom of the test given by F D = ( r − 1 ) ( s − 1 ) . The criterion proposed is:
χ 2 c a l c u l a t e d ≥ χ 2 α , ( r − 1 ) ( s − 1 ) .
The above test requires each of the expected frequencies to have a minimum value of five, in which case, the contrast statistic is recalculated by Yates’ correction for continuity. The new statistic is as follows:
χ 2 c a l c u l a t e d = ∑ i = 1 r ∑ j = 1 s ( | O i j − e i j | − 0.5 ) 2 e i j .
Finally, once the factors associated with having the disease are identified, we proceed to determine the value of the strength of the association, which ranges from 0 to 1 and is a measure of the degree of association found between variables X and Y. To do this, we calculated the contingency coefficient [72], which is estimated based on the chi-square statistic, as follows:
C = χ 2 n + χ 2 .

2.4. Genomic Data Used by TIP_finder

We first use the coffee genome (C. canephora) [73] to test available tools. The pseudomolecule chr2 (with 53 Mb), was used as a reference genome, the LTR-retrotransposons (SIRE lineage) as the repeat library, and one resequencing Illumina data set for C. canephora (SRX5013724) [74].
In order to test the behavior and performance of TIP_finder, we used human endogenous retroviruses (HERVs) as a proof of concept since HERVs have been correlated with many tumors. For example, there are several studies reporting an overexpression of HERVs type K in several cancers, such as melanoma [28], prostate [27], pancreas [75], breast [26], and ovary [76]. Breast cancer was chosen here as a condition of interest [26] since it is the leading cause of death by cancer among women worldwide [77].
We used publicly available datasets from the NCBI Whole Genome Sequencing (WGS) database, corresponding to people with breast cancer (i.e., case patients) and people without the disease (i.e., control patients). Case patients were obtained from the SRA (Sequence Read Archive) database through search Equation (8). We retrieved 512 datasets and discarded those obtained by sequencing technologies other than Illumina and that lacked public permission for download. As a result, ~300 datasets of paired reads were pre-selected and only 10% were used (30 patients) due to the huge disk space required for storage. Control patients were consulted from Chen 2016 [78], obtaining a list of 30 patients from the BioProject PRJNA551447 on a study of schizophrenia in a Chinese population (Supplementary Material S1):
((((((((((breast cancer) AND Homo sapiens[Organism]) AND PAIRED[Layout])) AND GENOMIC[Source])) AND WGS[Strategy]))) NOT exome).
For each patient (i.e., case and control), a maximum of 30 million reads was downloaded using the fastq-dump tool from the SRA toolkit [79]. Overall, the datasets comprised approximately 1.1 Terabytes of disk space (642 Gb for controls and 458 Gb for cases). HERV-K consensus sequences were obtained from Repbase (free version 2017) [80]. Finally, the human reference genome (3.3 Gb) in FASTA format was downloaded from the NCBI (assembly GRCh38.p13).

2.5. Computational Resources

TIP_finder was executed on a server with 56 cores E5-2695 v. 3, 252 GB RAM, and CentOS 7 managed by Slurm. Python 3 and the libraries used were installed using Anaconda 3 environments and pre-required software were loaded using environmental modules [81].

3. Results

3.1. Problems Encountered with Large Genomes and Testing TIP_finder

Several tools have been developed to search for transposable element insertions using the approach of discordant read pair analysis when mapped against a reference genome. Recently, we were interested in the study of the diversity of cultivated coffee trees (C. canephora and C. arabica) induced by the insertion of LTR-retrotransposons, the latter elements representing about 50% of the size of these genomes. The fastest and most accurate tool available (TRACKPOSON) was used with one family of LTR-retrotransposons (SIRE lineage) with a single set of resequencing data and a single pseudomolecule of C. canephora (representing 1/11 and 1/22 of the C. canephora and C. arabica pseudomolecules respectively). The result was obtained in 9474 seconds, and generated a BLAST output file of 7577.6 Mb. Several factors explain this result, such as bottlenecks in the execution time of the BLAST program, the writing time and format of the results, and the lack of overall parallelization of the tools. Based on this observation, we decided to develop a new tool that could overcome these different bottlenecks in data analysis time to obtain a faster program that could be more easily deployed on supercomputers.
Thus, TIP_finder was developed for use in current supercomputers to analyze TIPs in huge genomic datasets. We first tested TIP_finder on coffee data and we obtained a computational time of 594 seconds and a BLAST result file size of 448 Mb, showing a decrease in the execution time of 15.9X and a decrease in file size of 16.9X. To test TIP_finder on very large genomes, we decided to analyze the insertional polymorphisms caused by HERVs of type K in the human genome since several studies have demonstrated that these elements correlate with breast cancer [26], a disease that causes high mortality rates in women [77].
TIP_finder implements the MPI standard to reduce execution times; thus, allowing researchers to use more datasets. We tested TIP_finder (using NCBI-BLAST and MagicBLAST as aligners) and TRACKPOSON using a maximum of 30 million of reads (~1X of the human genome) from one randomly selected case dataset (Figure 3A,B) and one control dataset (Figure 3C,D) using 2, 4, 8, 16, 32, 44, and 56 cores, executed 10 times, to determine the behaviors of runtime (Figure 3A,C) and speedup (Figure 3B,D). To calculate the speedup, we used the execution with two cores as serial time since our MPI implementation used one processor as master and the others as workers.
We analyzed the behavior of TIP_finder for each of the four steps: Step 1: mapping and alignment, Step 2: creation of the dictionary with reads, Step 3: filter reads with one hit, and Step 4: the post-processing of TIPs. Here, we were particularly interested in monitoring the execution time and the speedup for each step with two different numbers of cores: two and 32. We used the same case and control datasets used for the speedup analyses. Figure 4 shows that Step 1 in both data sets benefited the most (up to 11.5X) from the parallel strategy used, because this Step is the most computationally time consuming (Figure S1).
Finally, we compared the execution times between TIP_finder and TRACKPOSON by randomly selecting five case and control datasets, in addition to the two patients used in the speedup analyses. This analysis aimed to measure the runtimes based on different numbers of TIPs in the datasets (Figure 5). All executions of TRACKPOSON and TIP_finder were performed using 32 cores and a maximum of 30 million reads.

3.2. Correlation of HERV-K TIPs with Breast Cancer as a Proof of Concept

After implementing TIP_finder and testing it in several configurations (e.g., different numbers of cores and datasets), we were interested in applying the software to a specific problem, such as the association of TIPs (in this case HERV-K) with breast cancer. We used publicly available Illumina read datasets found in the SRA database (NCBI) [82]. For this proof of concept, we selected 30 datasets for cases and 30 for controls. After obtaining the datasets, we ran TIP_finder to generate a tabular-separated file (in bed format) per individual, which contained the information of TIPs located throughout the chromosomes of each genome. Then, we ran TIP_finder_utils.py to automatically generate the analyses based on these data.
TIP_finder_utils.py (“finalMatrix” utility) was run using the tabular-delimited file produced by TIP_finder.py to generate a presence/absence matrix of TIPs for each patient (using the bed file from TIP_finder produced for each individual). From this matrix, the number of TIPs for cases and controls were identified and grouped into bins. Figure 6A shows only the controls and displays the number of TIPs in a reduced bin with a specific range from 0 to 220. The cumulative number of TIPs in case patients, as seen in Figure 6B, is distributed over all bins (0–4000), while TIPs in control patients are only clustered into the first bin (0–500). These findings suggest a much higher insertional activity (HERV-K) for cases compared to controls (Figure 6).
We executed the “histograms” utility of TIP_finder_utils.py to analyze the insertional activity of cases and controls. Figure 7A shows a high insertional activity of HERVs for case patients while control patients (Figure 7B) show a low insertion activity. Based on this insertional activity, the frequency of TIPs for both categories (cases and controls) were grouped according to their estimated position along each chromosome using bins of 300. TIP_finder split the reference genome into windows of 10 Kb, so the positions of each TIP were given within intervals of the window length. This clustering approach allowed detecting chromosome sections where the insertion of HERVs is markedly more frequent using the “peaks” utility of TIP_finder_utils.py. Figure 8 shows the distribution of TIPs for four of the 23 chromosomes (e.g., chromosomes 6, 8, 12, and 17), in which the peaks of polymorphic presence were most significant.
Finally, we performed a statistical association analysis through the “association” utility of TIP_finder_utils.py to investigate which TIPs may be associated with the condition of interest (i.e., breast cancer). We found a total of 3860 associated TIPs (significance level of 5%) along all the chromosomes. Chromosomes 18, X, 17, 7, 12, and 8 showed the highest associated TIPS with 498, 459, 427, 390, 233, and 299 TIPs, respectively (Figure 9 and Supplementary Material S3).

4. Discussion

The analysis of TIPs offers an efficient way to elucidate the dynamics of genome evolution in all organisms, from archaea to humans, via the activity of mobile elements. TIPs may help to answer multiple questions that arise in the field of genomics and bioinformatics, such as the process of the domestication of plants [30,82,83], how mobile elements have shaped the evolution of mammals [84], and inheritance features between organisms [85]. The study of genome diversity based on the insertional activity of TEs also provides an opportunity in the future to understand genetic diseases [24] caused by TIPs, such as breast cancer and its association with HERVs [26,86].
Although previous algorithms significantly improve the discovery of TIPs [30,87], their programming remains suboptimal considering that a parallel strategy and the use of multi-core architectures were not considered in their implementation. TRACKPOSON displays a speedup of ~1X (Figure 3) on 56 processors compared to two processors, due to the suboptimal execution of NCBI-BLAST and BLAST output processing using a single-processor Perl script. This demonstrates that a parallel strategy is required for better execution and this limitation is a disadvantage for massive sequencing projects that generate huge genomic datasets. To address this deficiency, the HPC and MPI techniques used by TIP_finder allow researchers to use many processors simultaneously in an efficient manner, taking advantage of all the benefits of supercomputers. TIP_finder can accelerate up to 55× the runtime compared to TRACKPOSON (Figure 3B,D). This improvement is important considering that TRACKPOSON is the fastest TIP detector in a recent benchmarking analysis [36]. Furthermore, TIP_finder can process larger datasets, improving its speedup using a high number of processors. Particularly, Step 1, which involves mapping (Bowtie2) and alignment (NCBI-BLAST or Magic-BLAST), requires the highest execution time and computational resources (Figure S1); thus, greatly benefiting from a parallel programming approach to speed up the analysis compared to TRACKPOSON (Figure 4). Furthermore, TIP_finder is a scalable software since it can run single-node applications (i.e., software than can run over multiple processors but are limited to one server), such as Bowtie2, NCBI-BLAST, and MagicBLAST, on multiple servers by splitting data into different inputs (i.e., data is split to run on different servers, see Figure 2). This feature allows taking advantage of current clusters that follow a distributed memory structure and makes the TIP_finder an especially useful tool for the massive sequencing projects required in the post-genomic era. Since TIP_finder applies a validated strategy to find TIPs [30], the results have high quality in terms of precision and sensitivity [36].
In a proof of concept using the sequencing data of 60 patients, TIP_finder demonstrated suitable performance for TIPs identification (represented here as hte consensus sequences of HERV-K) in the large human genome. Despite this not being the objective of the present study, the association between HERV-K and breast cancer in humans was confirmed. HERV-K in human breast cancer are involved in the processes of replication and tumorigenesis [88], derived directly from the protein and gene activities of mobile elements in the host [12]. This association is supported by the literature, indicating that this disease could be caused by the presence of HERVs in chromosomes 6, 7, 8, and 19 [26]. A larger sample of patients is necessary to further study this association and particularly, identify the coding regions involved. Moreover, TIP_finder should be tested on other human genetic diseases involving the mobility of endogenous retroviruses.
Overall, TIP_finder is a software that, besides its computational versatility (scalability and high performance), can run analyses on any group of organisms (plants, bacteria, animals) and uses information from any type of mobile element. Finally, this software will be used in future work on applications of the deep analysis of retrovirus interference or association with diseases in humans. These applications will help unravel potential applications in the industry of personalized medicine and pharmacology, as well as innovate with biological pathways for the treatment and prevention of the condition analyzed.

5. Conclusions

The development of HPC-based applications in bioinformatics is required on a daily basis. Thanks to the reduced cost of sequencing technologies, the challenge is not in obtaining genomic datasets, but rather in the analysis of this information in record time. Studies on TIPs have been performed in large datasets, demonstrating promising results for solving biological questions. Nevertheless, TIPs analyses show limitations in scalability since the available tools do not implement parallel strategies and follow a single-node approach. TIP_finder proves to be a very useful tool to be implemented on an HPC computing cluster architecture with many servers and with genomic datasets of considerable size. Moreover, the analysis of polymorphic activities of HERV-K could be used to perform a genomic-scale study of certain cancer types, such as breast cancer.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-7737/9/9/281/s1, Supplementary Material S1: SRA (Sequence Read Archive) accessions list of the case and control patients used, Supplementary Material S2: Figure S1: The proportion of execution times per step using a randomly selected case dataset and: (A) two cores (serial) and (B) 32 cores; and utilizing a randomly selected control dataset and: (C) two cores (serial) and (D) 32 cores using NCBI-BLAST as aligner, Table S1: TIP_finder execution times using one randomly selected case dataset, Table S2: TIP_finder execution times using one randomly selected control dataset, Table S3: TIP_finder execution time per step using one randomly selected case dataset and two processors (serial), Table S4: TIP_finder execution time per step using one randomly selected case dataset and 32 processors, Table S5: TIP_finder execution time per step using one randomly selected control dataset and two processors (serial), Table S6: TIP_finder execution time per step using one randomly selected control dataset and 32 processors, Table S7: Comparative executions using TRACKPOSON and TIP_finder on five randomly selected case and control datasets and 32 cores, Supplementary Material S3: TIPs showing statistical association with breast cancer, Supplementary Material S4: Chromosomal distribution of TIPs in the human genome.

Author Contributions

Conceptualization, S.O.-A. and R.G.; methodology, S.O.-A., N.T.-O., J.S.P. and C.F.J.-V.; software, S.O.-A.; supervision, R.T.-S. and R.G.; writing—original draft, S.O.-A., N.T.-O., J.S.P., C.F.J.-V. and R.T.-S.; writing—review and editing, S.O.-A., N.T.-O., J.S.P., C.F.J.-V., R.T.-S., and R.G. All authors have read and agreed to the published version of the manuscript.

Funding

Simon Orozco-Arias is supported by a Ph.D. grant from the Ministry of Science, Technology and Innovation (Minciencias) of Colombia, Grant Call 785/2017. The authors were supported by Universidad Autónoma de Manizales, Manizales, Colombia under project 589-089, publication costs were supported by IRD, Montpellier, France, and Romain Guyot was supported by the LMI BIO-INCA. The funding institutions had no role in the study design, data collection and analysis, the decision to publish, or preparation of the manuscript.

Acknowledgments

The authors acknowledge the IFB Core Cluster that is part of the National Network of Compute Resources (NNCR) of the Institut Français de Bioinformatique (https://www.france-bioinformatique.fr), the Genotoul Bioinformatics platform (http://bioinfo.genotoul.fr/), and the IRD itrop (https://bioinfo.ird.fr/) at IRD Montpellier for providing HPC resources that have contributed to the research results reported in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McClintock, B. The Origin and Behavior of Mutable Loci in Maize. Proc. Natl. Acad. Sci. USA 1950, 36, 344–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Orozco-Arias, S.; Liu, J.; Id, R.T.; Ceballos, D.; Silva, D.; Id, D.; Ming, R.; Guyot, R. Inpactor, Integrated and Parallel Analyzer and Classifier of LTR Retrotransposons and Its Application for Pineapple LTR Retrotransposons Diversity and Dynamics. Biology 2018, 7, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bourque, G.; Burns, K.H.; Gehring, M.; Gorbunova, V.; Seluanov, A.; Hammell, M.; Imbeault, M.; Izsvák, Z.; Levin, H.L.; Macfarlan, T.S.; et al. Ten things you should know about transposable elements. Genome Biol. 2018, 19, 199. [Google Scholar] [CrossRef] [PubMed]
  4. Rishishwar, L.; Mariño-Ramírez, L.; Jordan, I.K. Benchmarking computational tools for polymorphic transposable element detection. Brief. Bioinform. 2017, 18, 908–918. [Google Scholar] [CrossRef]
  5. Orozco-Arias, S.; Isaza, G.; Guyot, R.; Tabares-soto, R. A systematic review of the application of machine learning in the detection and classi fi cation of transposable elements. PeerJ 2019, 7, 18311. [Google Scholar] [CrossRef] [Green Version]
  6. Choulet, F.; Alberti, A.; Theil, S.; Glover, N.; Barbe, V.; Daron, J.; Pingault, L.; Sourdille, P.; Couloux, A.; Paux, E.; et al. Structural and functional partitioning of bread wheat chromosome 3B. Science 2014, 345, 1249721. [Google Scholar] [CrossRef]
  7. Su, W.; Gu, X.; Peterson, T. TIR-Learner, a New Ensemble Method for TIR Transposable Element Annotation, Provides Evidence for Abundant New Transposable Elements in the Maize Genome. Mol. Plant 2019, 12, 447–460. [Google Scholar] [CrossRef] [Green Version]
  8. De Koning, A.P.J.; Gu, W.; Castoe, T.A.; Batzer, M.A.; Pollock, D.D. Repetitive elements may comprise over Two-Thirds of the human genome. PLoS Genet. 2011, 7. [Google Scholar] [CrossRef] [Green Version]
  9. Orozco-Arias, S.; Isaza, G.; Guyot, R. Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning. Int. J. Mol. Sci. 2019, 20, 3837. [Google Scholar] [CrossRef] [Green Version]
  10. Todorovska, E. Retrotransposons and their role in plant—Genome evolution. Biotechnol. Biotechnol. Equip. 2014, 21, 294–305. [Google Scholar] [CrossRef] [Green Version]
  11. Casacuberta, E.; González, J. The impact of transposable elements in environmental adaptation. Mol. Ecol. 2013, 22, 1503–1517. [Google Scholar] [CrossRef]
  12. Zhang, M.; Liang, J.Q. Expressional activation and functional roles of human endogenous retroviruses in cancers. Rev. Med. Virol. 2019, 1–11. [Google Scholar] [CrossRef] [PubMed]
  13. Lisch, D. How important are transposons for plant evolution? Nat. Rev. Genet. 2013, 14, 49–61. [Google Scholar] [CrossRef] [PubMed]
  14. Deininger, P.L.; Batzer, M.A. Alu repeats and human disease. Mol. Genet. Metab. 1999, 67, 183–193. [Google Scholar] [CrossRef] [Green Version]
  15. Hancks, D.C.; Kazazian, H.H. Active human retrotransposons: Variation and disease. Curr. Opin. Genet. Dev. 2012, 22, 191–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Beck, C.R.; Garcia-Perez, J.L.; Badge, R.M.; Moran, J.V. LINE-1 Elements in Structural Variation and Disease. Annu. Rev. Genom. Hum. Genet. 2011, 12, 187–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Chaparro, C.; Gayraud, T.; de Souza, R.F.; Domingues, D.S.; Akaffou, S.S.; Vanzela, A.L.L.; de Kochko, A.; Rigoreau, M.; Crouzillat, D.; Hamon, S.; et al. Terminal-repeat retrotransposons with GAG domain in plant genomes: A new testimony on the complex world of transposable elements. Genome Biol. Evol. 2015, 7, 493–504. [Google Scholar] [CrossRef] [Green Version]
  18. Wicker, T.; Sabot, F.; Hua-Van, A.; Bennetzen, J.L.; Capy, P.; Chalhoub, B.; Flavell, A.; Leroy, P.; Morgante, M.; Panaud, O.; et al. A unified classification system for eukaryotic transposable elements. Nat. Rev. Genet. 2007, 8, 973–982. [Google Scholar] [CrossRef]
  19. Neumann, P.; Novák, P.; Hoštáková, N.; MacAs, J. Systematic survey of plant LTR-retrotransposons elucidates phylogenetic relationships of their polyprotein domains and provides a reference for element classification. Mob. DNA 2019, 10, 1. [Google Scholar] [CrossRef]
  20. Laten, H.M.; Gaston, G.D. Plant Endogenous Retroviruses? A Case of Mysterious ORFs. In Plant Transposable Elements; Spriger: Berlin/Heidelberg, Germany, 2012; pp. 89–112. [Google Scholar]
  21. Grandbastien, M.-A.A. LTR retrotransposons, handy hitchhikers of plant regulation and stress response. Biochim. Biophys. Acta Gene Regul. Mech. 2015, 1849, 403–416. [Google Scholar] [CrossRef]
  22. Subramanian, R.P.; Wildschutte, J.H.; Russo, C.; Coffin, J.M. Identification, characterization, and comparative genomic distribution of the HERV-K (HML-2) group of human endogenous retroviruses. Retrovirology 2011, 1–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Havecker, E.R.; Gao, X.; Voytas, D.F. The diversity of LTR retrotransposons. Genome Biol. 2004, 5, 225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rishishwar, L.; Wang, L.; Clayton, E.A.; Mariño-Ramírez, L.; McDonald, J.F.; Jordan, I.K. Population and clinical genetics of human transposable elements in the (post) genomic era. Mob. Genet. Elem. 2017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Asch, H.L.; Eliacin, E.; Fanning, T.G.; Connolly, J.L.; Bratthauer, G.; Asch, B.B. Comparative Expression of the LINE-1 p40 Protein in Human Breast Breast Carcinomas and Normal Breast Tissues. Oncol. Res. Featur. Preclin. Clin. Cancer Ther. 1996, 8, 239–247. [Google Scholar]
  26. Johanning, G.L.; Malouf, G.G.; Zheng, X.; Esteva, F.J.; Weinstein, J.N.; Wang-Johanning, F.; Su, X. Expression of human endogenous retrovirus-K is strongly associated with the basal-like breast cancer phenotype. Sci. Rep. 2017, 7, 41960. [Google Scholar] [CrossRef] [Green Version]
  27. Goering, W.; Schmitt, K.; Dostert, M.; Schaal, H.; Mayer, J.; Schulz, W.A. Human Endogenous Retrovirus HERV-K (HML-2) Activity in Prostate Cancer Is Dominated by a Few Loci. Prostate 2015, 1971, 1958–1971. [Google Scholar] [CrossRef]
  28. Roesch, A.; Meese, E.; Mayer, J.; Schmitt, K. Transcriptional Profiling of Human Endogenous Retrovirus Group HERV-K (HML-2) Loci in Melanoma. Genome Biol. Evol. 2013, 5, 307–328. [Google Scholar] [CrossRef] [Green Version]
  29. Bratthauer, G.L.; Fanning, T.G. LINE-1 retrotransposon expression in pediatric germ cell tumors. Cancer 1993, 71, 2383–2386. [Google Scholar] [CrossRef]
  30. Carpentier, M.C.; Manfroi, E.; Wei, F.J.; Wu, H.P.; Lasserre, E.; Llauro, C.; Debladis, E.; Akakpo, R.; Hsing, Y.I.; Panaud, O. Retrotranspositional landscape of Asian rice revealed by 3000 genomes. Nat. Commun. 2019, 10. [Google Scholar] [CrossRef]
  31. Martienssen, R. Epigenetic phenomena: Paramutation and gene silencing in plants. Curr. Biol. 1996, 6, 810–813. [Google Scholar] [CrossRef] [Green Version]
  32. Drongitis, D.; Aniello, F.; Fucci, L.; Donizetti, A. Roles of Transposable Elements in the Different Layers of Gene Expression Regulation. Int. J. Mol. Sci. 2019, 20, 5755. [Google Scholar] [CrossRef] [Green Version]
  33. Barrón, M.G.; Fiston-Lavier, A.-S.; Petrov, D.A.; González, J. Population Genomics of Transposable Elements in Drosophila. Annu. Rev. Genet. 2014, 48, 561–581. [Google Scholar] [CrossRef] [Green Version]
  34. Rigal, M.; Mathieu, O. A “mille-feuille” of silencing: Epigenetic control of transposable elements. Biochim. Biophys. Acta Gene Regul. Mech. 2011, 1809, 452–458. [Google Scholar] [CrossRef]
  35. Ewing, A.D. Transposable element detection from whole genome sequence data. Mob. DNA 2015, 6. [Google Scholar] [CrossRef] [Green Version]
  36. Vendrell-Mir, P.; Barteri, F.; Merenciano, M.; González, J.; Casacuberta, J.M.; Castanera, R. A benchmark of transposon insertion detection tools using real data. Mob. DNA 2019, 10, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Jiang, C.; Chen, C.; Huang, Z.; Liu, R.; Verdier, J. ITIS, a bioinformatics tool for accurate identification of transposon insertion sites using next-generation sequencing data. BMC Bioinform. 2015, 16. [Google Scholar] [CrossRef] [Green Version]
  38. Hénaff, E.; Zapata, L.; Casacuberta, J.M.; Ossowski, S. Jitterbug: Somatic and germline transposon insertion detection at single-nucleotide resolution. BMC Genom. 2015, 16, 768. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Helman, E.; Lawrence, M.S.; Stewart, C.; Sougnez, C.; Getz, G.; Meyerson, M. Somatic retrotransposition in human cancer revealed by whole-genome and exome sequencing. Genome Res. 2014, 24, 1053–1063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Mohiyuddin, M.; Mu, J.C.; Li, J.; Bani Asadi, N.; Gerstein, M.B.; Abyzov, A.; Wong, W.H.; Lam, H.Y.K. MetaSV: An accurate and integrative structural-variant caller for next generation sequencing. Bioinformatics 2015, 31, 2741–2744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Kroon, M.; Lameijer, E.W.; Lakenberg, N.; Hehir-Kwa, J.Y.; Thung, D.T.; Slagboom, P.E.; Kok, J.N.; Ye, K. Detecting dispersed duplications in high-throughput sequencing data using a database-free approach. Bioinformatics 2016, 32, 505–510. [Google Scholar] [CrossRef] [Green Version]
  42. Tran, H.T.M.; Ramaraj, T.; Furtado, A.; Lee, L.S.; Henry, R.J. Use of a draft genome of coffee (Coffea arabica) to identify SNP s associated with caffeine content. Plant Biotechnol. J. 2018, 16, 1756–1766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Mueller, L.; Strickler, S.; Domingues, D.; Pereira, L.; Andrade, A.; Marraccini, P.; Ming, R.; Wai, J.; Albert, V.; Giuliano, G.; et al. Towards a better understanding of the Coffea Arabica genome structure. In Proceedings of the Embrapa Café-Artigo em Anais de Congresso (ALICE), International Conference on Coffee Science, Armenia, Colombia, 8–13 September 2014. [Google Scholar]
  44. Wu, A.R.; Neff, N.F.; Kalisky, T.; Dalerba, P.; Treutlein, B.; Rothenberg, M.E.; Mburu, F.M.; Mantalas, G.L.; Sim, S.; Clarke, M.F.; et al. Quantitative assessment of single-cell RNA-sequencing methods. Nat. Methods 2014, 11, 41–46. [Google Scholar] [CrossRef] [Green Version]
  45. Berlin, K.; Koren, S.; Chin, C.S.; Drake, J.P.; Landolin, J.M.; Phillippy, A.M. Assembling large genomes with single-molecule sequencing and locality-sensitive hashing. Nat. Biotechnol. 2015, 33, 623–630. [Google Scholar] [CrossRef]
  46. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; FitzHugh, W.; et al. Initial sequencing and analysis of the human genome. Nature 2001, 409, 860–921. [Google Scholar] [PubMed] [Green Version]
  47. Cheng, S.; Melkonian, M.; Smith, S.A.; Brockington, S.; Archibald, J.M.; Delaux, P.-M.; Li, F.-W.; Melkonian, B.; Mavrodiev, E.V.; Sun, W.; et al. 10KP: A phylodiverse genome sequencing plan. Gigascience 2018, 7, giy013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Lewin, H.A.; Robinson, G.E.; Kress, W.J.; Baker, W.J.; Coddington, J.; Crandall, K.A.; Durbin, R.; Edwards, S.V.; Forest, F.; Gilbert, M.T.P.; et al. Earth BioGenome Project: Sequencing life for the future of life. Proc. Natl. Acad. Sci. USA 2018, 115, 4325–4333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Tatusova, T. Update on genomic databases and resources at the national center for biotechnology information. In Methods in Molecular Biology; Humana Press Inc.: Totowa, NJ, USA, 2016; Volume 1415, pp. 3–30. [Google Scholar]
  50. Tabares Soto, R. Programación Paralela Sobre Arquitecturas Heterogéneas. Master’s Thesis, Universidad Nacional de Colombia, Manizales, Colombia, November 2016. [Google Scholar]
  51. Orozco-Arias, S.; Tabares-Soto, R.; Ceballos, D.; Guyot, R. Parallel Programming in Biological Sciences, Taking Advantage of Supercomputing in Genomics. In Advances in Computing; Solano, A., Ordoñez, H., Eds.; Springer: Zurich, Switzerland, 2017; Volume 735, pp. 627–643. ISBN 9781457720819. [Google Scholar]
  52. Mikailov, M.; Luo, F.J.; Barkley, S.; Valleru, L.; Whitney, S.; Liu, Z.; Thakkar, S.; Tong, W.; Petrick, N. Scaling bioinformatics applications on HPC. BMC Bioinform. 2017, 18, 501. [Google Scholar] [CrossRef]
  53. Orozco-Arias, S.; Camargo-forero, L.; Correa, J.C.; Guyot, R.; Cristancho, M. BIOS-ParallelBlast: Paralelización optimizada de alineamiento de secuencias sobre Xeon Phi. Ing. Investig. Technol. 2017, 18, 423–432. [Google Scholar]
  54. Rodrigues, F.M.; von Mering, C. Sequence analysis HPC-CLUST: Distributed hierarchical clustering for large sets of nucleotide sequences. Bioinformatics 2014, 30, 287–288. [Google Scholar] [CrossRef]
  55. Sawyer, S.; Horton, M.; Burdyshaw, C.; Brook, G. HPC-BLAST: Distributed BLAST for Modern HPC Clusters. In Proceedings of the 11th International Conference on Bioinformatics and Computational Biology, Honolulu, HI, USA, 18–20 March 2019. [Google Scholar]
  56. Gropp, W.; Lusk, E. Fault Tolerance in Message Passing Interface Programs. Int. J. High Perform. Comput. Appl. 2004, 18, 363–372. [Google Scholar] [CrossRef]
  57. Gropp, W.; Lusk, E.; Doss, N.; Skjellum, A. A high-performance, portable implementation of the MPI message passing interface standard. Parallel Comput. 1996, 22, 789–828. [Google Scholar] [CrossRef]
  58. Aguilar Castro, J.L.; Leiss, E. Introducción a la Computación Paralela; Universidad de los Andes: Mérida, Venezuela, 2004; ISBN 9801207523. [Google Scholar]
  59. Chen, W.; Feng, P.M.; Deng, E.Z.; Lin, H.; Chou, K.C. iTIS-PseTNC: A sequence-based predictor for identifying translation initiation site in human genes using pseudo trinucleotide composition. Anal. Biochem. 2014, 462, 76–83. [Google Scholar] [CrossRef] [PubMed]
  60. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Data, G.P.; et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [Green Version]
  63. Altschup, S.F.; Gish, W.; Pennsylvania, T.; Park, U. Basic Local Alignment Search Tool 2Department of Computer Science. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  64. Boratyn, G.M.; Thierry-Mieg, J.; Thierry-Mieg, D.; Busby, B.; Madden, T.L. Magic-BLAST, an accurate RNA-seq aligner for long and short reads. BMC Bioinform. 2019, 20, 1–19. [Google Scholar] [CrossRef]
  65. Dalcín, L.; Paz, R.; Storti, M.; D’Elía, J. MPI for Python: Performance improvements and MPI-2 extensions. J. Parallel Distrib. Comput. 2008, 68, 655–662. [Google Scholar] [CrossRef]
  66. McKinney, W. Data Structures for Statistical Computing in Python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 28 June–3 July 2010; pp. 51–56. [Google Scholar]
  67. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  68. Waskom, M.; Botvinnik, O.; O’Kane, D.; Hobson, P.; Lukauskas, S.; Gemperline, D.C.; Augspurger, T.; Halchenko, Y.; Cole, J.B.; Warmenhoven, J.; et al. Mwaskom/seaborn: v0.8.1 (September 2017). Zenodo 2017. [Google Scholar] [CrossRef]
  69. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
  70. Maurya, V.N.; Singh, V.V.; Yusuf, M.U. Statistical Analysis on the Rate of Kidney (Renal) Failure. Am. J. Appl. Math. Stat. 2014, 2, 6–12. [Google Scholar] [CrossRef] [Green Version]
  71. Edition, S. Smooth Tests of Goodness of Fit; John Wiley & Sons: Hoboken, NJ, USA, 2009; ISBN 9780470824429. [Google Scholar]
  72. Cochran, W.G. Some methods for strenghthening the commom χ^2 tests. Biometrics 2012, 10, 417–451. [Google Scholar] [CrossRef]
  73. Denoeud, F.; Carretero-Paulet, L.; Dereeper, A.; Droc, G.; Guyot, R.; Pietrella, M.; Zheng, C.; Alberti, A.; Anthony, F.; Aprea, G.; et al. The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. Science 2014, 345, 1181–1184. [Google Scholar] [CrossRef] [Green Version]
  74. Huang, L.; Wang, X.; Dong, Y.; Long, Y.; Hao, C.; Yan, L.; Shi, T. Resequencing 93 accessions of coffee unveils independent and parallel selection during Coffea species divergence. Plant Mol. Biol. 2020, 103, 1–11. [Google Scholar] [CrossRef]
  75. Li, M.; Radvanyi, L.; Yin, B.; Rycaj, K.; Li, J.; Chivukula, R.; Lin, K.; Lu, Y.; Shen, J.; Chang, D.Z.; et al. Downregulation of Human Endogenous Retrovirus Type K (HERV-K) Viral env RNA in Pancreatic Cancer Cells Decreases Cell Proliferation and Tumor Growth. Clin. Cancer Res. 2017, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Cegolon, L.; Salata, C.; Weiderpass, E.; Vineis, P.; Palù, G.; Mastrangelo, G. Human endogenous retroviruses and cancer prevention: Evidence and prospects. BMC Cancer 2013, 13, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Desantis, C.E.; Ma, J.; Goding Sauer, A.; Newman, L.A.; Jemal, A. Breast Cancer Statistics, 2017, Racial Disparity in Mortality by State. CA Cancer J. Clin. 2017, 67, 439–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Chen, X. Understanding the Genetic Architecture of Schizophrenia in Chinese Population; University of Nevada Las Vegas: Las Vegas, NV, USA, 2016. [Google Scholar]
  79. Sherry, S.; Xiao, C.; Durbrow, K.; Kimelman, M.; Rodarmer, K.; Shumway, M.; Yaschenko, E. NCBI SRA Toolkit Technology for Next Generation Sequence Data. In Proceedings of the Plant and Animal Genome XX Conference, San Diego, CA, USA, 14–18 January 2012. [Google Scholar]
  80. Jurka, J.; Kapitonov, V.V.; Pavlicek, A.; Klonowski, P.; Kohany, O.; Walichiewicz, J. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 2005, 110, 462–467. [Google Scholar] [CrossRef]
  81. Furlani, J.L.; Osel, P.W. Abstract Yourself With Modules. In Proceedings of the 10th USENIX Conference on System Administration, San Jose, CA, USA, 7–12 November1996; USENIX Association: Berkley, CA, USA, 1996; pp. 193–204. [Google Scholar]
  82. Leinonen, R.; Sugawara, H.; Shumway, M.; International Nucleotide Sequence Database Collaboration. The sequence read archive. Nucleic Acids Res. 2010, 39, D19–D21. [Google Scholar] [CrossRef] [Green Version]
  83. Lynch, V.J.; Leclerc, R.D.; May, G.; Wagner, G.P. Transposon-mediated rewiring of gene regulatory networks contributed to the evolution of pregnancy in mammals. Nat. Genet. 2011, 43, 1154–1159. [Google Scholar] [CrossRef]
  84. Chuong, E.B. Retroviruses facilitate the rapid evolution of the mammalian placenta. Bioessays 2013, 35, 853–861. [Google Scholar] [CrossRef] [PubMed]
  85. Hermann, D.; Egue, F.; Tastard, E.; Nguyen, D.-H.; Casse, N.; Caruso, A.; Hiard, S.; Marchand, J.; Chenais, B.; Morant-Manceau, A.; et al. An introduction to the vast world of transposable elements—what about the diatoms? DIATOM Res. 2014, 29, 91–104. [Google Scholar] [CrossRef]
  86. Pericay, C.; Díez, O.; Campos, B.; Balmaña, J.; Domènech, M.; Lerma, E.; Baena, M.; Sabaté, J.M.; Gómez, A.; López, J.J.; et al. Características clinicopatológicas y evolución clínica de pacientes con cáncer de mama y mutaciones en los genes BRCA1 o BRCA2. Med. Clin. 2001, 117, 161–166. [Google Scholar] [CrossRef]
  87. The 3000 Rice Genomes Project. The 3000 rice genomes project. Gigascience 2014, 3, 1–6. [Google Scholar] [CrossRef]
  88. McDowell, J.M.; Meyers, B.C. A transposable element is domesticated for service in the plant immune system. Proc. Natl. Acad. Sci. USA 2013, 110, 14821–14822. [Google Scholar] [CrossRef] [Green Version]
Biology 09 00281 g001 550
Figure 1. TIP_finder methodology and schematic representation of the pipeline. TE: transposable element, TIP: transposon insertion polymorphism.
Figure 1. TIP_finder methodology and schematic representation of the pipeline. TE: transposable element, TIP: transposon insertion polymorphism.
Biology 09 00281 g001
Biology 09 00281 g002 550
Figure 2. Flowchart of the parallel strategy implemented in TIP_finder.
Figure 2. Flowchart of the parallel strategy implemented in TIP_finder.
Biology 09 00281 g002
Biology 09 00281 g003 550
Figure 3. Total average runtime and speedup of TIP_finder -using NCBI-BLAST and MagicBLAST- and TRACKPOSON using 2, 4, 8, 16, 32, 44, and 56 cores with a randomly selected case dataset (30 million of reads) and executed 10 times (A,B), and a randomly selected control dataset (26.5 million of reads) and executed 10 times (C,D). The times for all executions can be found in Supplementary Material S2 Table S1 (for case dataset), and in Table S2 (for control dataset).
Figure 3. Total average runtime and speedup of TIP_finder -using NCBI-BLAST and MagicBLAST- and TRACKPOSON using 2, 4, 8, 16, 32, 44, and 56 cores with a randomly selected case dataset (30 million of reads) and executed 10 times (A,B), and a randomly selected control dataset (26.5 million of reads) and executed 10 times (C,D). The times for all executions can be found in Supplementary Material S2 Table S1 (for case dataset), and in Table S2 (for control dataset).
Biology 09 00281 g003
Biology 09 00281 g004 550
Figure 4. TIP_finder speedup of each step with 32 cores for one case and one control datasets.
Figure 4. TIP_finder speedup of each step with 32 cores for one case and one control datasets.
Biology 09 00281 g004
Biology 09 00281 g005 550
Figure 5. Total average runtime and speedup of TIP_finder and TRACKPOSON using six randomly selected control datasets and six randomly selected case datasets, each one executed 10 times. (A) Comparison of runtimes between TRACKPOSON and TIP_finder (with B: NCBI-BLAST, and M: Magic-BLAST), and (B) the comparison of the speedup of TIP_finder (NCBI-BLAST and MagicBLAST) compared to TRACKPOSON. Additional information can be consulted in Table S7.
Figure 5. Total average runtime and speedup of TIP_finder and TRACKPOSON using six randomly selected control datasets and six randomly selected case datasets, each one executed 10 times. (A) Comparison of runtimes between TRACKPOSON and TIP_finder (with B: NCBI-BLAST, and M: Magic-BLAST), and (B) the comparison of the speedup of TIP_finder (NCBI-BLAST and MagicBLAST) compared to TRACKPOSON. Additional information can be consulted in Table S7.
Biology 09 00281 g005
Biology 09 00281 g006 550
Figure 6. Number of TIPs insertions for the cases and controls identified by TIP_finder and grouped into bins of 500. Control and case patients are show in red, and in blue respectively. (A) Distribution of the number of TIPs insertions in cases and controls shown in blue and red, respectively. (B) Distribution of TIPs in controls within a range of 140–220.
Figure 6. Number of TIPs insertions for the cases and controls identified by TIP_finder and grouped into bins of 500. Control and case patients are show in red, and in blue respectively. (A) Distribution of the number of TIPs insertions in cases and controls shown in blue and red, respectively. (B) Distribution of TIPs in controls within a range of 140–220.
Biology 09 00281 g006
Biology 09 00281 g007 550
Figure 7. Distribution of TIPs insertions (Human endogenous retroviruses type K (HERV-K) insertions) for the cases and controls. The X axis represents the number of patients in log2 scale and the Y axis is the TIPs insertion number. (A,B) correspond to the cases and controls, respectively. Variations in the peaks’ frequency suggest a change in the insertional activity under a condition of interest, in this case cancer.
Figure 7. Distribution of TIPs insertions (Human endogenous retroviruses type K (HERV-K) insertions) for the cases and controls. The X axis represents the number of patients in log2 scale and the Y axis is the TIPs insertion number. (A,B) correspond to the cases and controls, respectively. Variations in the peaks’ frequency suggest a change in the insertional activity under a condition of interest, in this case cancer.
Biology 09 00281 g007
Biology 09 00281 g008 550
Figure 8. Distribution of the number of TIPs—HERVs—along human chromosomes. The X axis represents the chromosome length (on a scale of 1 × 108), where cases and controls are shown in blue and red, respectively. The Y axis represents the number of TIPs along the chromosome length. The arrows show the end of each chromosome. The graphs for all chromosomes are available in Supplementary Material S4.
Figure 8. Distribution of the number of TIPs—HERVs—along human chromosomes. The X axis represents the chromosome length (on a scale of 1 × 108), where cases and controls are shown in blue and red, respectively. The Y axis represents the number of TIPs along the chromosome length. The arrows show the end of each chromosome. The graphs for all chromosomes are available in Supplementary Material S4.
Biology 09 00281 g008
Biology 09 00281 g009 550
Figure 9. Number of TIPs statistically associated with breast cancer.
Figure 9. Number of TIPs statistically associated with breast cancer.
Biology 09 00281 g009
Table
Table 1. Contingency table used for the statistical association analysis.
Table 1. Contingency table used for the statistical association analysis.
Y/XX = 0X = 1
Y = 0O00O01
Y = 1O10O11

Share and Cite

MDPI and ACS Style

Orozco-Arias, S.; Tobon-Orozco, N.; Piña, J.S.; Jiménez-Varón, C.F.; Tabares-Soto, R.; Guyot, R. TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets. Biology 2020, 9, 281. https://doi.org/10.3390/biology9090281

AMA Style

Orozco-Arias S, Tobon-Orozco N, Piña JS, Jiménez-Varón CF, Tabares-Soto R, Guyot R. TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets. Biology. 2020; 9(9):281. https://doi.org/10.3390/biology9090281

Chicago/Turabian Style

Orozco-Arias, Simon, Nicolas Tobon-Orozco, Johan S. Piña, Cristian Felipe Jiménez-Varón, Reinel Tabares-Soto, and Romain Guyot. 2020. "TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets" Biology 9, no. 9: 281. https://doi.org/10.3390/biology9090281

APA Style

Orozco-Arias, S., Tobon-Orozco, N., Piña, J. S., Jiménez-Varón, C. F., Tabares-Soto, R., & Guyot, R. (2020). TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets. Biology, 9(9), 281. https://doi.org/10.3390/biology9090281

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop