scMuffin: an R package to disentangle solid tumor heterogeneity by single-cell gene expression analysis.

INTRODUCTION: Single-cell (SC) gene expression analysis is crucial to dissect the complex cellular heterogeneity of solid tumors, which is one of the main obstacles for the development of effective cancer treatments. Such tumors typically contain a mixture of cells with aberrant genomic and transcriptomic profiles affecting specific sub-populations that might have a pivotal role in cancer progression, whose identification eludes bulk RNA-sequencing approaches. We present scMuffin, an R package that enables the characterization of cell identity in solid tumors on the basis of a various and complementary analyses on SC gene expression data. RESULTS: scMuffin provides a series of functions to calculate qualitative and quantitative scores, such as: expression of marker sets for normal and tumor conditions, pathway activity, cell state trajectories, Copy Number Variations, transcriptional complexity and proliferation state. Thus, scMuffin facilitates the combination of various evidences that can be used to distinguish normal and tumoral cells, define cell identities, cluster cells in different ways, link genomic aberrations to phenotypes and identify subtle differences between cell subtypes or cell states. We analysed public SC expression datasets of human high-grade gliomas as a proof-of-concept to show the value of scMuffin and illustrate its user interface. Nevertheless, these analyses lead to interesting findings, which suggest that some chromosomal amplifications might underlie the invasive tumor phenotype and the presence of cells that possess tumor initiating cells characteristics. CONCLUSIONS: The analyses offered by scMuffin and the results achieved in the case study show that our tool helps addressing the main challenges in the bioinformatics analysis of SC expression data from solid tumors.

Characterization and comparison of gene-centered human interactomes.

The complex web of macromolecular interactions occurring within cells-the interactome-is the backbone of an increasing number of studies, but a clear consensus on the exact structure of this network is still lacking. Different genome-scale maps of human interactome have been obtained through several experimental techniques and functional analyses. Moreover, these maps can be enriched through literature-mining approaches, and different combinations of various 'source' databases have been used in the literature. It is therefore unclear to which extent the various interactomes yield similar results when used in the context of interactome-based approaches in network biology. We compared a comprehensive list of human interactomes on the basis of topology, protein complexes, molecular pathways, pathway cross-talk and disease gene prediction. In a general context of relevant heterogeneity, our study provides a series of qualitative and quantitative parameters that describe the state of the art of human interactomes and guidelines for selecting interactomes in future applications.

Gene relevance based on multiple evidences in complex networks.

MOTIVATION: Multi-omics approaches offer the opportunity to reconstruct a more complete picture of the molecular events associated with human diseases, but pose challenges in data analysis. Network-based methods for the analysis of multi-omics leverage the complex web of macromolecular interactions occurring within cells to extract significant patterns of molecular alterations. Existing network-based approaches typically address specific combinations of omics and are limited in terms of the number of layers that can be jointly analysed. In this study, we investigate the application of network diffusion to quantify gene relevance on the basis of multiple evidences (layers). RESULTS: We introduce a gene score (mND) that quantifies the relevance of a gene in a biological process taking into account the network proximity of the gene and its first neighbours to other altered genes. We show that mND has a better performance over existing methods in finding altered genes in network proximity in one or more layers. We also report good performances in recovering known cancer genes. The pipeline described in this article is broadly applicable, because it can handle different types of inputs: in addition to multi-omics datasets, datasets that are stratified in many classes (e.g., cell clusters emerging from single cell analyses) or a combination of the two scenarios. AVAILABILITY AND IMPLEMENTATION: The R package 'mND' is available at URL: https://www.itb.cnr.it/mnd. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

γ-TRIS: a graph-algorithm for comprehensive identification of vector genomic insertion sites.

SUMMARY: Retroviruses and their vector derivatives integrate semi-randomly in the genome of host cells and are inherited by their progeny as stable genetic marks. The retrieval and mapping of the sequences flanking the virus-host DNA junctions allows the identification of insertion sites in gene therapy or virally infected patients, essential for monitoring the evolution of genetically modified cells in vivo. However, since ∼30% of insertions land in low complexity or repetitive regions of the host cell genome, they cannot be correctly assigned and are currently discarded, limiting the accuracy and predictive power of clonal tracking studies. Here, we present γ-TRIS, a new graph-based genome-free alignment tool for identifying insertion sites even if embedded in low complexity regions. By using γ-TRIS to reanalyze clinical studies, we observed improvements in clonal quantification and tracking. AVAILABILITY AND IMPLEMENTATION: Source code at https://bitbucket.org/bereste/g-tris. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

isma: an R package for the integrative analysis of mutations detected by multiple pipelines.

BACKGROUND: Recent comparative studies have brought to our attention how somatic mutation detection from next-generation sequencing data is still an open issue in bioinformatics, because different pipelines result in a low consensus. In this context, it is suggested to integrate results from multiple calling tools, but this operation is not trivial and the burden of merging, comparing, filtering and explaining the results demands appropriate software. RESULTS: We developed isma (integrative somatic mutation analysis), an R package for the integrative analysis of somatic mutations detected by multiple pipelines for matched tumor-normal samples. The package provides a series of functions to quantify the consensus, estimate the variability, underline outliers, integrate evidences from publicly available mutation catalogues and filter sites. We illustrate the capabilities of isma analysing breast cancer somatic mutations generated by The Cancer Genome Atlas (TCGA) using four pipelines. CONCLUSIONS: Comparing different "points of view" on the same data, isma generates a unique mutation catalogue and a series of reports that underline common patterns, variability, as well as sites already catalogued by other studies (e.g. TCGA), so as to design and apply filtering strategies to screen more reliable sites. The package is available for non-commercial users at the URL https://www.itb.cnr.it/isma .

The 2017 Network Tools and Applications in Biology (NETTAB) workshop: aims, topics and outcomes.

The 17th International NETTAB workshop was held in Palermo, Italy, on October 16-18, 2017. The special topic for the meeting was "Methods, tools and platforms for Personalised Medicine in the Big Data Era", but the traditional topics of the meeting series were also included in the event. About 40 scientific contributions were presented, including four keynote lectures, five guest lectures, and many oral communications and posters. Also, three tutorials were organised before and after the workshop. Full papers from some of the best works presented in Palermo were submitted for this Supplement of BMC Bioinformatics. Here, we provide an overview of meeting aims and scope. We also shortly introduce selected papers that have been accepted for publication in this Supplement, for a complete presentation of the outcomes of the meeting.

Network-Based Integrative Analysis of Genomics, Epigenomics and Transcriptomics in Autism Spectrum Disorders.

Current studies suggest that autism spectrum disorders (ASDs) may be caused by many genetic factors. In fact, collectively considering multiple studies aimed at characterizing the basic pathophysiology of ASDs, a large number of genes has been proposed. Addressing the problem of molecular data interpretation using gene networks helps to explain genetic heterogeneity in terms of shared pathways. Besides, the integrative analysis of multiple omics has emerged as an approach to provide a more comprehensive view of a disease. In this work, we carry out a network-based meta-analysis of the genes reported as associated with ASDs by studies that involved genomics, epigenomics, and transcriptomics. Collectively, our analysis provides a prioritization of the large number of genes proposed to be associated with ASDs, based on genes' relevance within the intracellular circuits, the strength of the supporting evidence of association with ASDs, and the number of different molecular alterations affecting genes. We discuss the presence of the prioritized genes in the SFARI (Simons Foundation Autism Research Initiative) database and in gene networks associated with ASDs by other investigations. Lastly, we provide the full results of our analyses to encourage further studies on common targets amenable to therapy.

An infrastructure for precision medicine through analysis of big data.

BACKGROUND: Nowadays, the increasing availability of omics data, due to both the advancements in the acquisition of molecular biology results and in systems biology simulation technologies, provides the bases for precision medicine. Success in precision medicine depends on the access to healthcare and biomedical data. To this end, the digitization of all clinical exams and medical records is becoming a standard in hospitals. The digitization is essential to collect, share, and aggregate large volumes of heterogeneous data to support the discovery of hidden patterns with the aim to define predictive models for biomedical purposes. Patients' data sharing is a critical process. In fact, it raises ethical, social, legal, and technological issues that must be properly addressed. RESULTS: In this work, we present an infrastructure devised to deal with the integration of large volumes of heterogeneous biological data. The infrastructure was applied to the data collected between 2010-2016 in one of the major diagnostic analysis laboratories in Italy. Data from three different platforms were collected (i.e., laboratory exams, pathological anatomy exams, biopsy exams). The infrastructure has been designed to allow the extraction and aggregation of both unstructured and semi-structured data. Data are properly treated to ensure data security and privacy. Specialized algorithms have also been implemented to process the aggregated information with the aim to obtain a precise historical analysis of the clinical activities of one or more patients. Moreover, three Bayesian classifiers have been developed to analyze examinations reported as free text. Experimental results show that the classifiers exhibit a good accuracy when used to analyze sentences related to the sample location, diseases presence and status of the illnesses. CONCLUSIONS: The infrastructure allows the integration of multiple and heterogeneous sources of anonymized data from the different clinical platforms. Both unstructured and semi-structured data are processed to obtain a precise historical analysis of the clinical activities of one or more patients. Data aggregation allows to perform a series of statistical assessments required to answer complex questions that can be used in a variety of fields, such as predictive and precision medicine. In particular, studying the clinical history of patients that have developed similar pathologies can help to predict or individuate markers able to allow an early diagnosis of possible illnesses.

Genome-based analysis for the identification of genes involved in o-xylene degradation in Rhodococcus opacus R7.

BACKGROUND: Bacteria belonging to the Rhodococcus genus play an important role in the degradation of many contaminants, including methylbenzenes. These bacteria, widely distributed in the environment, are known to be a powerhouse of numerous degradation functions, due to their ability to metabolize a wide range of organic molecules including aliphatic, aromatic, polycyclic aromatic compounds (PAHs), phenols, and nitriles. In accordance with their immense catabolic diversity, Rhodococcus spp. possess large and complex genomes, which contain a multiplicity of catabolic genes, a high genetic redundancy of biosynthetic pathways and a sophisticated regulatory network. The present study aimed to identify genes involved in the o-xylene degradation in R. opacus strain R7 through a genome-based approach. RESULTS: Using genome-based analysis we identified all the sequences in the R7 genome annotated as dioxygenases or monooxygenases/hydroxylases and clustered them into two different trees. The akb, phe and prm sequences were selected as genes encoding respectively for dioxygenases, phenol hydroxylases and monooxygenases and their putative involvement in o-xylene oxidation was evaluated. The involvement of the akb genes in o-xylene oxidation was demonstrated by RT-PCR/qPCR experiments after growth on o-xylene and by the selection of the R7-50 leaky mutant. Although the akb genes are specifically activated for o-xylene degradation, metabolic intermediates of the pathway suggested potential alternative oxidation steps, possibly through monooxygenation. This led us to further investigate the role of the prm and the phe genes. Results showed that these genes were transcribed in a constitutive manner, and that the activity of the Prm monooxygenase was able to transform o-xylene slowly in intermediates as 3,4-dimethylphenol and 2-methylbenzylalcohol. Moreover, the expression level of phe genes, homologous to the phe genes of Rhodococcus spp. 1CP and UPV-1 with a 90% identity, could explain their role in the further oxidation of o-xylene and R7 growth on dimethylphenols. CONCLUSIONS: These results suggest that R7 strain is able to degrade o-xylene by the Akb dioxygenase system leading to the production of the corresponding dihydrodiol. Likewise, the redundancy of sequences encoding for several monooxygenases/phenol hydroxylases, supports the involvement of other oxygenases converging in the o-xylene degradation pathway in R7 strain.

How computer science can help in understanding the 3D genome architecture.

Chromosome conformation capture techniques are producing a huge amount of data about the architecture of our genome. These data can provide us with a better understanding of the events that induce critical regulations of the cellular function from small changes in the three-dimensional genome architecture. Generating a unified view of spatial, temporal, genetic and epigenetic properties poses various challenges of data analysis, visualization, integration and mining, as well as of high performance computing and big data management. Here, we describe the critical issues of this new branch of bioinformatics, oriented at the comprehension of the three-dimensional genome architecture, which we call 'Nucleome Bioinformatics', looking beyond the currently available tools and methods, and highlight yet unaddressed challenges and the potential approaches that could be applied for tackling them. Our review provides a map for researchers interested in using computer science for studying 'Nucleome Bioinformatics', to achieve a better understanding of the biological processes that occur inside the nucleus.

Systems medicine of inflammaging.

Systems Medicine (SM) can be defined as an extension of Systems Biology (SB) to Clinical-Epidemiological disciplines through a shifting paradigm, starting from a cellular, toward a patient centered framework. According to this vision, the three pillars of SM are Biomedical hypotheses, experimental data, mainly achieved by Omics technologies and tailored computational, statistical and modeling tools. The three SM pillars are highly interconnected, and their balancing is crucial. Despite the great technological progresses producing huge amount of data (Big Data) and impressive computational facilities, the Bio-Medical hypotheses are still of primary importance. A paradigmatic example of unifying Bio-Medical theory is the concept of Inflammaging. This complex phenotype is involved in a large number of pathologies and patho-physiological processes such as aging, age-related diseases and cancer, all sharing a common inflammatory pathogenesis. This Biomedical hypothesis can be mapped into an ecological perspective capable to describe by quantitative and predictive models some experimentally observed features, such as microenvironment, niche partitioning and phenotype propagation. In this article we show how this idea can be supported by computational methods useful to successfully integrate, analyze and model large data sets, combining cross-sectional and longitudinal information on clinical, environmental and omics data of healthy subjects and patients to provide new multidimensional biomarkers capable of distinguishing between different pathological conditions, e.g. healthy versus unhealthy state, physiological versus pathological aging.

iSmaRT: a toolkit for a comprehensive analysis of small RNA-Seq data.

Summary: The interest in investigating the biological roles of small non-coding RNAs (sncRNAs) is increasing, due to the pleiotropic effects of these molecules exert in many biological contexts. While several methods and tools are available to study microRNAs (miRNAs), only few focus on novel classes of sncRNAs, in particular PIWI-interacting RNAs (piRNAs). To overcome these limitations, we implemented iSmaRT ( i ntegrative Sm all R NA T ool-kit), an automated pipeline to analyze smallRNA-Seq data. Availability and Implementation: iSmaRT is a collection of bioinformatics tools and own algorithms, interconnected through a Graphical User Interface (GUI). In addition to performing comprehensive analyses on miRNAs, it implements specific computational modules to analyze piRNAs, predicting novel ones and identifying their RNA targets. A smallRNA-Seq dataset generated from brain samples of Huntington's Disease patients was used here to illustrate iSmaRT performances, demonstrating how the pipeline can provide, in a rapid and user friendly way, a comprehensive analysis of different classes of sncRNAs. iSmaRT is freely available on the web at ftp://labmedmolge-1.unisa.it (User: iSmart - Password: password). Contact: aweisz@unisa.it or ggiurato@unisa.it. Supplementary information: Supplementary data are available at Bioinformatics online.

Speeding Up the Identification of Cystic Fibrosis Transmembrane Conductance Regulator-Targeted Drugs: An Approach Based on Bioinformatics Strategies and Surface Plasmon Resonance.

Cystic fibrosis (CF) is mainly caused by the deletion of Phe 508 (ΔF508) in the cystic fibrosis transmembrane conductance regulator (CFTR) protein that is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. New drugs able to rescue ΔF508-CFTR trafficking are eagerly awaited. An integrated bioinformatics and surface plasmon resonance (SPR) approach was here applied to investigate the rescue mechanism(s) of a series of CFTR-ligands including VX809, VX770 and some aminoarylthiazole derivatives (AAT). Computational studies tentatively identified a large binding pocket in the ΔF508-CFTR nucleotide binding domain-1 (NBD1) and predicted all the tested compounds to bind to three sub-regions of this main pocket. Noticeably, the known CFTR chaperone keratin-8 (K8) seems to interact with some residues located in one of these sub-pockets, potentially interfering with the binding of some ligands. SPR results corroborated all these computational findings. Moreover, for all the considered ligands, a statistically significant correlation was determined between their binding capability to ΔF508-NBD1 measured by SPR and the pockets availability measured by computational studies. Taken together, these results demonstrate a strong agreement between the in silico prediction and the SPR-generated binding data, suggesting a path to speed up the identification of new drugs for the treatment of cystic fibrosis.

VISPA2: a scalable pipeline for high-throughput identification and annotation of vector integration sites.

BACKGROUND: Bioinformatics tools designed to identify lentiviral or retroviral vector insertion sites in the genome of host cells are used to address the safety and long-term efficacy of hematopoietic stem cell gene therapy applications and to study the clonal dynamics of hematopoietic reconstitution. The increasing number of gene therapy clinical trials combined with the increasing amount of Next Generation Sequencing data, aimed at identifying integration sites, require both highly accurate and efficient computational software able to correctly process "big data" in a reasonable computational time. RESULTS: Here we present VISPA2 (Vector Integration Site Parallel Analysis, version 2), the latest optimized computational pipeline for integration site identification and analysis with the following features: (1) the sequence analysis for the integration site processing is fully compliant with paired-end reads and includes a sequence quality filter before and after the alignment on the target genome; (2) an heuristic algorithm to reduce false positive integration sites at nucleotide level to reduce the impact of Polymerase Chain Reaction or trimming/alignment artifacts; (3) a classification and annotation module for integration sites; (4) a user friendly web interface as researcher front-end to perform integration site analyses without computational skills; (5) the time speedup of all steps through parallelization (Hadoop free). CONCLUSIONS: We tested VISPA2 performances using simulated and real datasets of lentiviral vector integration sites, previously obtained from patients enrolled in a hematopoietic stem cell gene therapy clinical trial and compared the results with other preexisting tools for integration site analysis. On the computational side, VISPA2 showed a > 6-fold speedup and improved precision and recall metrics (1 and 0.97 respectively) compared to previously developed computational pipelines. These performances indicate that VISPA2 is a fast, reliable and user-friendly tool for integration site analysis, which allows gene therapy integration data to be handled in a cost and time effective fashion. Moreover, the web access of VISPA2 ( http://openserver.itb.cnr.it/vispa/ ) ensures accessibility and ease of usage to researches of a complex analytical tool. We released the source code of VISPA2 in a public repository ( https://bitbucket.org/andreacalabria/vispa2 ).

RNA-Generated and Gene-Edited Induced Pluripotent Stem Cells for Disease Modeling and Therapy.

Cellular reprogramming by epigenomic remodeling of chromatin holds great promise in the field of human regenerative medicine. As an example, human-induced Pluripotent Stem Cells (iPSCs) obtained by reprograming of patient somatic cells are sufficiently similar to embryonic stem cells (ESCs) and can generate all cell types of the human body. Clinical use of iPSCs is dependent on methods that do not utilize genome altering transgenic technologies that are potentially unsafe and ethically unacceptable. Transient delivery of exogenous RNA into cells provides a safer reprogramming system to transgenic approaches that rely on exogenous DNA or viral vectors. RNA reprogramming may prove to be more suitable for clinical applications and provide stable starting cell lines for gene-editing, isolation, and characterization of patient iPSC lines. The introduction and rapid evolution of CRISPR/Cas9 gene-editing systems has provided a readily accessible research tool to perform functional human genetic experiments. Similar to RNA reprogramming, transient delivery of mRNA encoding Cas9 in combination with guide RNA sequences to target specific points in the genome eliminates the risk of potential integration of Cas9 plasmid constructs. We present optimized RNA-based laboratory procedure for making and editing iPSCs. In the near-term these two powerful technologies are being harnessed to dissect mechanisms of human development and disease in vitro, supporting both basic, and translational research. J. Cell. Physiol. 232: 1262-1269, 2017. © 2016 Wiley Periodicals, Inc.

FGF2 and EGF Are Required for Self-Renewal and Organoid Formation of Canine Normal and Tumor Breast Stem Cells.

Recent studies suggest that human tumors are generated from cancer cells with stem cell (SC) properties. Spontaneously occurring cancers in dogs contain a diversity of cells that like for human tumors suggest that certain canine tumors are also generated from cancer stem cells (CSCs). CSCs, like normal SCs, have the capacity for self-renewal as mammospheres in suspension cultures. To understand how cells with SC properties contribute to canine mammary gland tumor development and progression, comparative analysis between normal SCs and CSCs, obtained from the same individual, is essential. We have utilized the property of sphere formation to develop culture conditions for propagating stem/progenitor cells from canine normal and tumor tissue. We show that cells from dissociated mammospheres retain sphere reformation capacity for several serial passages and have the capacity to generate organoid structures ex situ. Utilizing various culture conditions for passaging SCs and CSCs, fibroblast growth factor 2 (FGF2) and epidermal growth factor (EGF) were found to positively or negatively regulate mammosphere regeneration, organoid formation, and multi-lineage differentiation potential. The response of FGF2 and EGF on SCs and CSCs was different, with increased FGF2 and EGF self-renewal promoted in SCs and repressed in CSCs. Our protocol for propagating SCs from normal and tumor canine breast tissue will provide new opportunities in comparative mammary gland stem cell analysis between species and anticancer treatment and therapies for dogs. J. Cell. Biochem. 118: 570-584, 2017. © 2016 Wiley Periodicals, Inc.

BITS 2015: the annual meeting of the Italian Society of Bioinformatics.

This preface introduces the content of the BioMed Central journal Supplements related to the BITS 2015 meeting, held in Milan, Italy, from the 3th to the 5th of June, 2015.

A novel molecular dynamics approach to evaluate the effect of phosphorylation on multimeric protein interface: the αB-Crystallin case study.

BACKGROUND: Phosphorylation is one of the most important post-translational modifications (PTM) employed by cells to regulate several cellular processes. Studying the effects of phosphorylations on protein structures allows to investigate the modulation mechanisms of several proteins including chaperones, like the small HSPs, which display different multimeric structures according to the phosphorylation of a few serine residues. In this context, the proposed study is aimed at finding a method to correlate different PTM patterns (in particular phosphorylations at the monomers interface of multimeric complexes) with the dynamic behaviour of the complex, using physicochemical parameters derived from molecular dynamics simulations in the timescale of nanoseconds. RESULTS: We have developed a methodology relying on computing nine physicochemical parameters, derived from the analysis of short MD simulations, and combined with N identifiers that characterize the PTMs of the analysed protein. The nine general parameters were validated on three proteins, with known post-translational modified conformation and unmodified conformation. Then, we applied this approach to the case study of αB-Crystallin, a chaperone which multimeric state (up to 40 units) is supposed to be controlled by phosphorylation of Ser45 and Ser59. Phosphorylation of serines at the dimer interface induces the release of hexamers, the active state of αB-Crystallin. 30 ns of MD simulation were obtained for each possible combination of dimer phosphorylation state and average values of structural, dynamic, energetic and functional features were calculated on the equilibrated portion of the trajectories. Principal Component Analysis was applied to the parameters and the first five Principal Components, which summed up to 84 % of the total variance, were finally considered. CONCLUSIONS: The validation of this approach on multimeric proteins, which structures were known both modified and unmodified, allowed us to propose a new approach that can be used to predict the impact of PTM patterns in multi-modified proteins using data collected from short molecular dynamics simulations. Analysis on the αB-Crystallin case study clusters together all-P dimers with all-P hexamers and no-P dimer with no-P hexamer and results suggest a great influence of Ser59 phosphorylation on chain B.

Removing duplicate reads using graphics processing units.

BACKGROUND: During library construction polymerase chain reaction is used to enrich the DNA before sequencing. Typically, this process generates duplicate read sequences. Removal of these artifacts is mandatory, as they can affect the correct interpretation of data in several analyses. Ideally, duplicate reads should be characterized by identical nucleotide sequences. However, due to sequencing errors, duplicates may also be nearly-identical. Removing nearly-identical duplicates can result in a notable computational effort. To deal with this challenge, we recently proposed a GPU method aimed at removing identical and nearly-identical duplicates generated with an Illumina platform. The method implements an approach based on prefix-suffix comparison. Read sequences with identical prefix are considered potential duplicates. Then, their suffixes are compared to identify and remove those that are actually duplicated. Although the method can be efficiently used to remove duplicates, there are some limitations that need to be overcome. In particular, it cannot to detect potential duplicates in the event that prefixes are longer than 27 bases, and it does not provide support for paired-end read libraries. Moreover, large clusters of potential duplicates are split into smaller with the aim to guarantees a reasonable computing time. This heuristic may affect the accuracy of the analysis. RESULTS: In this work we propose GPU-DupRemoval, a new implementation of our method able to (i) cluster reads without constraints on the maximum length of the prefixes, (ii) support both single- and paired-end read libraries, and (iii) analyze large clusters of potential duplicates. CONCLUSIONS: Due to the massive parallelization obtained by exploiting graphics cards, GPU-DupRemoval removes duplicate reads faster than other cutting-edge solutions, while outperforming most of them in terms of amount of duplicates reads.

Methods for the integration of multi-omics data: mathematical aspects.

BACKGROUND: Methods for the integrative analysis of multi-omics data are required to draw a more complete and accurate picture of the dynamics of molecular systems. The complexity of biological systems, the technological limits, the large number of biological variables and the relatively low number of biological samples make the analysis of multi-omics datasets a non-trivial problem. RESULTS AND CONCLUSIONS: We review the most advanced strategies for integrating multi-omics datasets, focusing on mathematical and methodological aspects.