@article {818, title = {A genomic strategy for precision medicine in rare diseases: integrating customized algorithms into clinical practice.}, journal = {J Transl Med}, volume = {23}, year = {2025}, month = {2025 Jan 20}, pages = {86}, abstract = {

BACKGROUND: Despite the use of Next-Generation Sequencing (NGS) as the gold standard for the diagnosis of rare diseases, its clinical implementation has been challenging, limiting the cost-effectiveness of NGS and the understanding, control and safety essential for decision-making in clinical applications. Here, we describe a personalized NGS-based strategy integrating precision medicine into a public healthcare system and its implementation in the routine diagnosis process during a five-year pilot program.

METHODS: Our approach involved customized probe designs, the generation of virtual panels and the development of a personalized medicine module (PMM) for variant prioritization. This strategy was applied to 6500 individuals including 6267 index patients and 233 NGS-based carrier screenings.

RESULTS: Causative variants were identified in 2061 index patients (average 32.9\%, ranging from 12 to 62\% by condition). Also, 131 autosomal-recessive cases could be partially genetically diagnosed. These results led to over 5000 additional studies including carrier, prenatal and preimplantational tests or pharmacological and gene therapy treatments.

CONCLUSION: This strategy has shown promising improvements in the diagnostic rate, facilitating timely diagnosis and gradually expanding our services portfolio for rare diseases. The steps taken towards the integration of clinical and genomic data are opening new possibilities for conducting both retrospective and prospective healthcare studies. Overall, this study represents a major milestone in the ongoing efforts to improve our understanding and clinical management of rare diseases, a crucial area of medical research and care.

}, keywords = {Algorithms, Female, Genomics, High-Throughput Nucleotide Sequencing, Humans, Precision Medicine, rare diseases}, issn = {1479-5876}, doi = {10.1186/s12967-025-06069-2}, author = {M{\'e}ndez-Vidal, Cristina and Bravo-Gil, Nereida and Perez-Florido, Javier and Marcos-Luque, Irene and Fern{\'a}ndez, Raquel M and Fernandez-Rueda, Jose Luis and Gonz{\'a}lez-del Pozo, Mar{\'\i}a and Mart{\'\i}n-S{\'a}nchez, Marta and Fern{\'a}ndez-Su{\'a}rez, Elena and Mena, Marcela and Carmona, Rosario and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {728, title = {DOME: recommendations for supervised machine learning validation in biology.}, journal = {Nat Methods}, volume = {18}, year = {2021}, month = {2021 10}, pages = {1122-1127}, keywords = {Algorithms, Computational Biology, Guidelines as Topic, Humans, Models, Biological, Research Design, Supervised Machine Learning}, issn = {1548-7105}, doi = {10.1038/s41592-021-01205-4}, author = {Walsh, Ian and Fishman, Dmytro and Garcia-Gasulla, Dario and Titma, Tiina and Pollastri, Gianluca and Harrow, Jennifer and Psomopoulos, Fotis E and Tosatto, Silvio C E} } @article {712, title = {A versatile workflow to integrate RNA-seq genomic and transcriptomic data into mechanistic models of signaling pathways.}, journal = {PLoS Comput Biol}, volume = {17}, year = {2021}, month = {2021 02}, pages = {e1008748}, abstract = {

MIGNON is a workflow for the analysis of RNA-Seq experiments, which not only efficiently manages the estimation of gene expression levels from raw sequencing reads, but also calls genomic variants present in the transcripts analyzed. Moreover, this is the first workflow that provides a framework for the integration of transcriptomic and genomic data based on a mechanistic model of signaling pathway activities that allows a detailed biological interpretation of the results, including a comprehensive functional profiling of cell activity. MIGNON covers the whole process, from reads to signaling circuit activity estimations, using state-of-the-art tools, it is easy to use and it is deployable in different computational environments, allowing an optimized use of the resources available.

}, keywords = {Algorithms, Cell Line, Tumor, Computational Biology, Databases, Factual, Gene Expression Profiling, Genomics, High-Throughput Nucleotide Sequencing, Humans, Models, Theoretical, mutation, RNA-seq, Signal Transduction, Software, Transcriptome, whole exome sequencing, Workflow}, issn = {1553-7358}, doi = {10.1371/journal.pcbi.1008748}, author = {Garrido-Rodriguez, Mart{\'\i}n and L{\'o}pez-L{\'o}pez, Daniel and Ortuno, Francisco M and Pe{\~n}a-Chilet, Maria and Mu{\~n}oz, Eduardo and Calzado, Marco A and Dopazo, Joaquin} } @article {704, title = {Transparency and reproducibility in artificial intelligence.}, journal = {Nature}, volume = {586}, year = {2020}, month = {2020 10}, pages = {E14-E16}, keywords = {Algorithms, Artificial Intelligence, Reproducibility of Results}, issn = {1476-4687}, doi = {10.1038/s41586-020-2766-y}, author = {Haibe-Kains, Benjamin and Adam, George Alexandru and Hosny, Ahmed and Khodakarami, Farnoosh and Waldron, Levi and Wang, Bo and McIntosh, Chris and Goldenberg, Anna and Kundaje, Anshul and Greene, Casey S and Broderick, Tamara and Hoffman, Michael M and Leek, Jeffrey T and Korthauer, Keegan and Huber, Wolfgang and Brazma, Alvis and Pineau, Joelle and Tibshirani, Robert and Hastie, Trevor and Ioannidis, John P A and Quackenbush, John and Aerts, Hugo J W L} } @article {403, title = {A comparison of mechanistic signaling pathway activity analysis methods.}, journal = {Brief Bioinform}, volume = {20}, year = {2019}, month = {2019 Sep 27}, pages = {1655-1668}, abstract = {

Understanding the aspects of cell functionality that account for disease mechanisms or drug modes of action is a main challenge for precision medicine. Classical gene-based approaches ignore the modular nature of most human traits, whereas conventional pathway enrichment approaches produce only illustrative results of limited practical utility. Recently, a family of new methods has emerged that change the focus from the whole pathways to the definition of elementary subpathways within them that have any mechanistic significance and to the study of their activities. Thus, mechanistic pathway activity (MPA) methods constitute a new paradigm that allows recoding poorly informative genomic measurements into cell activity quantitative values and relate them to phenotypes. Here we provide a review on the MPA methods available and explain their contribution to systems medicine approaches for addressing challenges in the diagnostic and treatment of complex diseases.

}, keywords = {Algorithms, Humans, Postmortem Changes, Signal Transduction, Systems biology, Transcriptome}, issn = {1477-4054}, doi = {10.1093/bib/bby040}, author = {Amadoz, Alicia and Hidalgo, Marta R and Cubuk, Cankut and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {472, title = {Fast inexact mapping using advanced tree exploration on backward search methods.}, journal = {BMC Bioinformatics}, volume = {16}, year = {2015}, month = {2015 Jan 28}, pages = {18}, abstract = {

BACKGROUND: Short sequence mapping methods for Next Generation Sequencing consist on a combination of seeding techniques followed by local alignment based on dynamic programming approaches. Most seeding algorithms are based on backward search alignment, using the Burrows Wheeler Transform, the Ferragina and Manzini Index or Suffix Arrays. All these backward search algorithms have excellent performance, but their computational cost highly increases when allowing errors. In this paper, we discuss an inexact mapping algorithm based on pruning strategies for search tree exploration over genomic data.

RESULTS: The proposed algorithm achieves a 13x speed-up over similar algorithms when allowing 6 base errors, including insertions, deletions and mismatches. This algorithm can deal with 400 bps reads with up to 9 errors in a high quality Illumina dataset. In this example, the algorithm works as a preprocessor that reduces by 55\% the number of reads to be aligned. Depending on the aligner the overall execution time is reduced between 20-40\%.

CONCLUSIONS: Although not intended as a complete sequence mapping tool, the proposed algorithm could be used as a preprocessing step to modern sequence mappers. This step significantly reduces the number reads to be aligned, accelerating overall alignment time. Furthermore, this algorithm could be used for accelerating the seeding step of already available sequence mappers. In addition, an out-of-core index has been implemented for working with large genomes on systems without expensive memory configurations.

}, keywords = {Algorithms, Genome, Human, Genomics, High-Throughput Nucleotide Sequencing, Humans, Sequence Alignment, Sequence Analysis, DNA, Software}, issn = {1471-2105}, doi = {10.1186/s12859-014-0438-3}, author = {Salavert, Jos{\'e} and Tom{\'a}s, Andr{\'e}s and T{\'a}rraga, Joaqu{\'\i}n and Medina, Ignacio and Dopazo, Joaquin and Blanquer, Ignacio} } @article {474, title = {Using activation status of signaling pathways as mechanism-based biomarkers to predict drug sensitivity.}, journal = {Sci Rep}, volume = {5}, year = {2015}, month = {2015 Dec 18}, pages = {18494}, abstract = {

Many complex traits, as drug response, are associated with changes in biological pathways rather than being caused by single gene alterations. Here, a predictive framework is presented in which gene expression data are recoded into activity statuses of signal transduction circuits (sub-pathways within signaling pathways that connect receptor proteins to final effector proteins that trigger cell actions). Such activity values are used as features by a prediction algorithm which can efficiently predict a continuous variable such as the IC50 value. The main advantage of this prediction method is that the features selected by the predictor, the signaling circuits, are themselves rich-informative, mechanism-based biomarkers which provide insight into or drug molecular mechanisms of action (MoA).

}, keywords = {Algorithms, Antineoplastic Agents, biomarkers, Cell Line, Tumor, Cell Survival, gene expression, Humans, Lethal Dose 50, Neoplasms, Phosphorylation, Proteins, Signal Transduction}, issn = {2045-2322}, doi = {10.1038/srep18494}, author = {Amadoz, Alicia and Sebasti{\'a}n-Leon, Patricia and Vidal, Enrique and Salavert, Francisco and Dopazo, Joaquin} } @article {522, title = {Using GPUs for the exact alignment of short-read genetic sequences by means of the Burrows-Wheeler transform.}, journal = {IEEE/ACM Trans Comput Biol Bioinform}, volume = {9}, year = {2012}, month = {2012 Jul-Aug}, pages = {1245-56}, abstract = {

General Purpose Graphic Processing Units (GPGPUs) constitute an inexpensive resource for computing-intensive applications that could exploit an intrinsic fine-grain parallelism. This paper presents the design and implementation in GPGPUs of an exact alignment tool for nucleotide sequences based on the Burrows-Wheeler Transform. We compare this algorithm with state-of-the-art implementations of the same algorithm over standard CPUs, and considering the same conditions in terms of I/O. Excluding disk transfers, the implementation of the algorithm in GPUs shows a speedup larger than 12, when compared to CPU execution. This implementation exploits the parallelism by concurrently searching different sequences on the same reference search tree, maximizing memory locality and ensuring a symmetric access to the data. The paper describes the behavior of the algorithm in GPU, showing a good scalability in the performance, only limited by the size of the GPU inner memory.

}, keywords = {Algorithms, Animals, Computational Biology, Computer Graphics, Data Compression, Drosophila melanogaster, Genes, Insect, Image Processing, Computer-Assisted, Models, Genetic, Sequence Alignment, Sequence Analysis, DNA}, issn = {1557-9964}, doi = {10.1109/TCBB.2012.49}, author = {Salavert Torres, Jose and Blanquer Espert, Ignacio and Dom{\'\i}nguez, Andr{\'e}s Tom{\'a}s and Hern{\'a}ndez Garc{\'\i}a, Vicente and Medina Castell{\'o}, Ignacio and T{\'a}rraga Gim{\'e}nez, Joaqu{\'\i}n and Dopazo Bl{\'a}zquez, Joaqu{\'\i}n} } @article {529, title = {Differential expression in RNA-seq: a matter of depth.}, journal = {Genome Res}, volume = {21}, year = {2011}, month = {2011 Dec}, pages = {2213-23}, abstract = {

Next-generation sequencing (NGS) technologies are revolutionizing genome research, and in particular, their application to transcriptomics (RNA-seq) is increasingly being used for gene expression profiling as a replacement for microarrays. However, the properties of RNA-seq data have not been yet fully established, and additional research is needed for understanding how these data respond to differential expression analysis. In this work, we set out to gain insights into the characteristics of RNA-seq data analysis by studying an important parameter of this technology: the sequencing depth. We have analyzed how sequencing depth affects the detection of transcripts and their identification as differentially expressed, looking at aspects such as transcript biotype, length, expression level, and fold-change. We have evaluated different algorithms available for the analysis of RNA-seq and proposed a novel approach--NOISeq--that differs from existing methods in that it is data-adaptive and nonparametric. Our results reveal that most existing methodologies suffer from a strong dependency on sequencing depth for their differential expression calls and that this results in a considerable number of false positives that increases as the number of reads grows. In contrast, our proposed method models the noise distribution from the actual data, can therefore better adapt to the size of the data set, and is more effective in controlling the rate of false discoveries. This work discusses the true potential of RNA-seq for studying regulation at low expression ranges, the noise within RNA-seq data, and the issue of replication.

}, keywords = {Algorithms, Expressed Sequence Tags, Gene Expression Profiling, Gene Expression Regulation, Humans, Models, Genetic, Oligonucleotide Array Sequence Analysis}, issn = {1549-5469}, doi = {10.1101/gr.124321.111}, author = {Tarazona, Sonia and Garc{\'\i}a-Alcalde, Fernando and Dopazo, Joaquin and Ferrer, Alberto and Conesa, Ana} } @article {542, title = {Functional analysis of multiple genomic signatures demonstrates that classification algorithms choose phenotype-related genes.}, journal = {Pharmacogenomics J}, volume = {10}, year = {2010}, month = {2010 Aug}, pages = {310-23}, abstract = {

Gene expression signatures of toxicity and clinical response benefit both safety assessment and clinical practice; however, difficulties in connecting signature genes with the predicted end points have limited their application. The Microarray Quality Control Consortium II (MAQCII) project generated 262 signatures for ten clinical and three toxicological end points from six gene expression data sets, an unprecedented collection of diverse signatures that has permitted a wide-ranging analysis on the nature of such predictive models. A comprehensive analysis of the genes of these signatures and their nonredundant unions using ontology enrichment, biological network building and interactome connectivity analyses demonstrated the link between gene signatures and the biological basis of their predictive power. Different signatures for a given end point were more similar at the level of biological properties and transcriptional control than at the gene level. Signatures tended to be enriched in function and pathway in an end point and model-specific manner, and showed a topological bias for incoming interactions. Importantly, the level of biological similarity between different signatures for a given end point correlated positively with the accuracy of the signature predictions. These findings will aid the understanding, and application of predictive genomic signatures, and support their broader application in predictive medicine.

}, keywords = {Algorithms, Databases, Genetic, Endpoint Determination, Gene Expression Profiling, Genomics, Humans, Neural Networks, Computer, Oligonucleotide Array Sequence Analysis, Phenotype, Predictive Value of Tests, Proteins, Quality Control}, issn = {1473-1150}, doi = {10.1038/tpj.2010.35}, author = {Shi, W and Bessarabova, M and Dosymbekov, D and Dezso, Z and Nikolskaya, T and Dudoladova, M and Serebryiskaya, T and Bugrim, A and Guryanov, A and Brennan, R J and Shah, R and Dopazo, J and Chen, M and Deng, Y and Shi, T and Jurman, G and Furlanello, C and Thomas, R S and Corton, J C and Tong, W and Shi, L and Nikolsky, Y} } @article {576, title = {Serial Expression Analysis: a web tool for the analysis of serial gene expression data.}, journal = {Nucleic Acids Res}, volume = {38}, year = {2010}, month = {2010 Jul}, pages = {W239-45}, abstract = {

Serial transcriptomics experiments investigate the dynamics of gene expression changes associated with a quantitative variable such as time or dosage. The statistical analysis of these data implies the study of global and gene-specific expression trends, the identification of significant serial changes, the comparison of expression profiles and the assessment of transcriptional changes in terms of cellular processes. We have created the SEA (Serial Expression Analysis) suite to provide a complete web-based resource for the analysis of serial transcriptomics data. SEA offers five different algorithms based on univariate, multivariate and functional profiling strategies framed within a user-friendly interface and a project-oriented architecture to facilitate the analysis of serial gene expression data sets from different perspectives. SEA is available at sea.bioinfo.cipf.es.

}, keywords = {Algorithms, Gene Expression Profiling, Internet, Kinetics, Linear Models, Oligonucleotide Array Sequence Analysis, Software}, issn = {1362-4962}, doi = {10.1093/nar/gkq488}, author = {Nueda, Maria Jos{\'e} and Carbonell, Jos{\'e} and Medina, Ignacio and Dopazo, Joaquin and Conesa, Ana} } @article {582, title = {Gene set internal coherence in the context of functional profiling.}, journal = {BMC Genomics}, volume = {10}, year = {2009}, month = {2009 Apr 27}, pages = {197}, abstract = {

BACKGROUND: Functional profiling methods have been extensively used in the context of high-throughput experiments and, in particular, in microarray data analysis. Such methods use available biological information to define different types of functional gene modules (e.g. gene ontology -GO-, KEGG pathways, etc.) whose representation in a pre-defined list of genes is further studied. In the most popular type of microarray experimental designs (e.g. up- or down-regulated genes, clusters of co-expressing genes, etc.) or in other genomic experiments (e.g. Chip-on-chip, epigenomics, etc.) these lists are composed by genes with a high degree of co-expression. Therefore, an implicit assumption in the application of functional profiling methods within this context is that the genes corresponding to the modules tested are effectively defining sets of co-expressing genes. Nevertheless not all the functional modules are biologically coherent entities in terms of co-expression, which will eventually hinder its detection with conventional methods of functional enrichment.

RESULTS: Using a large collection of microarray data we have carried out a detailed survey of internal correlation in GO terms and KEGG pathways, providing a coherence index to be used for measuring functional module co-regulation. An unexpected low level of internal correlation was found among the modules studied. Only around 30\% of the modules defined by GO terms and 57\% of the modules defined by KEGG pathways display an internal correlation higher than the expected by chance.This information on the internal correlation of the genes within the functional modules can be used in the context of a logistic regression model in a simple way to improve their detection in gene expression experiments.

CONCLUSION: For the first time, an exhaustive study on the internal co-expression of the most popular functional categories has been carried out. Interestingly, the real level of coexpression within many of them is lower than expected (or even inexistent), which will preclude its detection by means of most conventional functional profiling methods. If the gene-to-function correlation information is used in functional profiling methods, the results obtained improve the ones obtained by conventional enrichment methods.

}, keywords = {Algorithms, Breast Neoplasms, Carcinoma, Intraductal, Noninfiltrating, Computational Biology, Databases, Nucleic Acid, Female, Gene Expression Profiling, Genomics, Humans, Oligonucleotide Array Sequence Analysis, Papillomavirus Infections, Reproducibility of Results}, issn = {1471-2164}, doi = {10.1186/1471-2164-10-197}, author = {Montaner, David and Minguez, Pablo and Al-Shahrour, F{\'a}tima and Dopazo, Joaquin} } @article {600, title = {Use of estimated evolutionary strength at the codon level improves the prediction of disease-related protein mutations in humans.}, journal = {Hum Mutat}, volume = {29}, year = {2008}, month = {2008 Jan}, pages = {198-204}, abstract = {

Predicting the functional impact of protein variation is one of the most challenging problems in bioinformatics. A rapidly growing number of genome-scale studies provide large amounts of experimental data, allowing the application of rigorous statistical approaches for predicting whether a given single point mutation has an impact on human health. Up until now, existing methods have limited their source data to either protein or gene information. Novel in this work, we take advantage of both and focus on protein evolutionary information by using estimated selective pressures at the codon level. Here we introduce a new method (SeqProfCod) to predict the likelihood that a given protein variant is associated with human disease or not. Our method relies on a support vector machine (SVM) classifier trained using three sources of information: protein sequence, multiple protein sequence alignments, and the estimation of selective pressure at the codon level. SeqProfCod has been benchmarked with a large dataset of 8,987 single point mutations from 1,434 human proteins from SWISS-PROT. It achieves 82\% overall accuracy and a correlation coefficient of 0.59, indicating that the estimation of the selective pressure helps in predicting the functional impact of single-point mutations. Moreover, this study demonstrates the synergic effect of combining two sources of information for predicting the functional effects of protein variants: protein sequence/profile-based information and the evolutionary estimation of the selective pressures at the codon level. The results of large-scale application of SeqProfCod over all annotated point mutations in SWISS-PROT (available for download at http://sgu.bioinfo.cipf.es/services/Omidios/; last accessed: 24 August 2007), could be used to support clinical studies.

}, keywords = {Algorithms, Codon, Computational Biology, Databases, Protein, DNA Mutational Analysis, Evolution, Molecular, Genetic Predisposition to Disease, Genetic Variation, Genome, Human, Humans, Iduronic Acid, Point Mutation, Polymorphism, Single Nucleotide, Proteins, Tumor Suppressor Protein p53}, issn = {1098-1004}, doi = {10.1002/humu.20628}, author = {Capriotti, Emidio and Arbiza, Leonardo and Casadio, Rita and Dopazo, Joaquin and Dopazo, Hern{\'a}n and Marti-Renom, Marc A} } @article {603, title = {DBAli tools: mining the protein structure space.}, journal = {Nucleic Acids Res}, volume = {35}, year = {2007}, month = {2007 Jul}, pages = {W393-7}, abstract = {

The DBAli tools use a comprehensive set of structural alignments in the DBAli database to leverage the structural information deposited in the Protein Data Bank (PDB). These tools include (i) the DBAlit program that allows users to input the 3D coordinates of a protein structure for comparison by MAMMOTH against all chains in the PDB; (ii) the AnnoLite and AnnoLyze programs that annotate a target structure based on its stored relationships to other structures; (iii) the ModClus program that clusters structures by sequence and structure similarities; (iv) the ModDom program that identifies domains as recurrent structural fragments and (v) an implementation of the COMPARER method in the SALIGN command in MODELLER that creates a multiple structure alignment for a set of related protein structures. Thus, the DBAli tools, which are freely accessible via the World Wide Web at http://salilab.org/DBAli/, allow users to mine the protein structure space by establishing relationships between protein structures and their functions.

}, keywords = {Algorithms, Amino Acid Sequence, Computational Biology, Data Interpretation, Statistical, Databases, Protein, Internet, Molecular Sequence Data, Protein Conformation, Proteins, Pseudomonas aeruginosa, Sequence Alignment, Sequence Analysis, Protein, Sequence Homology, Amino Acid, Software, Structure-Activity Relationship}, issn = {1362-4962}, doi = {10.1093/nar/gkm236}, author = {Marti-Renom, Marc A and Pieper, Ursula and Madhusudhan, M S and Rossi, Andrea and Eswar, Narayanan and Davis, Fred P and Al-Shahrour, F{\'a}tima and Dopazo, Joaquin and Sali, Andrej} } @article {606, title = {From genes to functional classes in the study of biological systems.}, journal = {BMC Bioinformatics}, volume = {8}, year = {2007}, month = {2007 Apr 03}, pages = {114}, abstract = {

BACKGROUND: With the popularization of high-throughput techniques, the need for procedures that help in the biological interpretation of results has increased enormously. Recently, new procedures inspired in systems biology criteria have started to be developed.

RESULTS: Here we present FatiScan, a web-based program which implements a threshold-independent test for the functional interpretation of large-scale experiments that does not depend on the pre-selection of genes based on the multiple application of independent tests to each gene. The test implemented aims to directly test the behaviour of blocks of functionally related genes, instead of focusing on single genes. In addition, the test does not depend on the type of the data used for obtaining significance values, and consequently different types of biologically informative terms (gene ontology, pathways, functional motifs, transcription factor binding sites or regulatory sites from CisRed) can be applied to different classes of genome-scale studies. We exemplify its application in microarray gene expression, evolution and interactomics.

CONCLUSION: Methods for gene set enrichment which, in addition, are independent from the original data and experimental design constitute a promising alternative for the functional profiling of genome-scale experiments. A web server that performs the test described and other similar ones can be found at: http://www.babelomics.org.

}, keywords = {Algorithms, Chromosome Mapping, Computer Simulation, Gene Expression Profiling, Models, Biological, Multigene Family, Signal Transduction, Software, Systems biology, User-Computer Interface}, issn = {1471-2105}, doi = {10.1186/1471-2105-8-114}, author = {Al-Shahrour, F{\'a}tima and Arbiza, Leonardo and Dopazo, Hern{\'a}n and Huerta-Cepas, Jaime and Minguez, Pablo and Montaner, David and Dopazo, Joaquin} }