Search

Pangenomics is a relatively new scientific field which investigates the union of all the genomes of a clade. The word pan means everything in ancient Greek; the term pangenomics originally regarded genomes of bacteria and was later intended to refer to human genomes as well. Modern bioinformatics offers several tools to analyze pangenomics data, paving the way to an emerging field that we can call computational pangenomics. Current computational power available for the bioinformatics community has made computational pangenomic analyses easy to perform, but this higher accessibility to pangenomics analysis also increases the chances to make mistakes and to produce misleading or inflated results, especially by beginners. To handle this problem, we present here a few quick tips for efficient and correct computational pangenomic analyses with a focus on bacterial pangenomics, by describing common mistakes to avoid and experienced best practices to follow in this field. We believe our recommendations can help the readers perform more robust and sound pangenomic analyses and to generate more reliable results. [ABSTRACT FROM AUTHOR]

Structural variation (SV) refers to insertions, deletions, inversions, and duplications in human genomes. SVs are present in approximately 1.5% of the human genome. Still, this small subset of genetic variation has been implicated in the pathogenesis of psoriasis, Crohn's disease and other autoimmune disorders, autism spectrum and other neurodevelopmental disorders, and schizophrenia. Since identifying structural variants is an important problem in genetics, several specialized computational techniques have been developed to detect structural variants directly from sequencing data. With advances in whole-genome sequencing (WGS) technologies, a plethora of SV detection methods have been developed. However, dissecting SVs from WGS data remains a challenge, with the majority of SV detection methods prone to a high false-positive rate, and no existing method able to precisely detect a full range of SVs present in a sample. Previous studies have shown that none of the existing SV callers can maintain high accuracy across various SV lengths and genomic coverages. Here, we report an integrated structural variant calling framework, Variant Identification and Structural Variant Analysis (VISTA), that leverages the results of individual callers using a novel and robust filtering and merging algorithm. In contrast to existing consensus-based tools which ignore the length and coverage, VISTA overcomes this limitation by executing various combinations of top-performing callers based on variant length and genomic coverage to generate SV events with high accuracy. We evaluated the performance of VISTA on comprehensive gold-standard datasets across varying organisms and coverage. We benchmarked VISTA using the Genome-in-a-Bottle gold standard SV set, haplotype-resolved de novo assemblies from the Human Pangenome Reference Consortium, along with an in-house polymerase chain reaction (PCR)-validated mouse gold standard set. VISTA maintained the highest F1 score among top consensus-based tools measured using a comprehensive gold standard across both mouse and human genomes. VISTA also has an optimized mode, where the calls can be optimized for precision or recall. VISTA-optimized can attain 100% precision and the highest sensitivity among other variant callers. In conclusion, VISTA represents a significant advancement in structural variant calling, offering a robust and accurate framework that outperforms existing consensus-based tools and sets a new standard for SV detection in genomic research. [ABSTRACT FROM AUTHOR]

Simple Summary: This review provides a comprehensive overview of the current knowledge on mutational signatures in colorectal cancer (CRC), highlighting their potential clinical applications. It discusses the challenges and limitations in translating these analyses into clinical practice and proposes strategies to overcome these obstacles. Additionally, it provides insights into future directions and emerging proof-of-concept studies that highlight the translational potential of mutational signature analysis in improving patient care and outcomes in CRC. A multitude of exogenous and endogenous processes have the potential to result in DNA damage. While the repair mechanisms are typically capable of correcting this damage, errors in the repair process can result in mutations. The findings of research conducted in 2012 indicate that mutations do not occur randomly but rather follow specific patterns that can be attributed to known or inferred mutational processes. The process of mutational signature analysis allows for the inference of the predominant mutational process for a given cancer sample, with significant potential for clinical applications. A deeper comprehension of these mutational signatures in CRC could facilitate enhanced prevention strategies, facilitate the comprehension of genotoxic drug activity, predict responses to personalized treatments, and, in the future, inform the development of targeted therapies in the context of precision oncology. The efforts of numerous researchers have led to the identification of several mutational signatures, which can be categorized into different mutational signature references. In CRC, distinct mutational signatures are identified as correlating with mismatch repair deficiency, polymerase mutations, and chemotherapy treatment. In this context, a mutational signature analysis offers considerable potential for enhancing minimal residual disease (MRD) tests in stage II (high-risk) and stage III CRC post-surgery, stratifying CRC based on the impacts of genetic and epigenetic alterations for precision oncology, identifying potential therapeutic vulnerabilities, and evaluating drug efficacy and guiding therapy, as illustrated in a proof-of-concept clinical trial. [ABSTRACT FROM AUTHOR]

The categorization of human diseases is mainly based on the affected organ system and phenotypic characteristics. This is limiting the view to the pathological manifestations, while it neglects mechanistic relationships that are crucial to develop therapeutic strategies. This work aims to advance the understanding of diseases and their relatedness beyond traditional phenotypic views. Hence, the similarity among 502 diseases is mapped using six different data dimensions encompassing molecular, clinical, and pharmacological information retrieved from public sources. Multiple distance measures and multi‐view clustering are used to assess the patterns of disease relatedness. The integration of all six dimensions into a consensus map of disease relationships reveals a divergent disease view from the International Classification of Diseases (ICD), emphasizing novel insights offered by a multi‐view disease map. Disease features such as genes, pathways, and chemicals that are enriched in distinct disease groups are identified. Finally, an evaluation of the top similar diseases of three candidate diseases common in the Western population shows concordance with known epidemiological associations and reveals rare features shared between Type 2 diabetes (T2D) and Alzheimer's disease. A revision of disease relationships holds promise for facilitating the reconstruction of comorbidity patterns, repurposing drugs, and advancing drug discovery in the future. [ABSTRACT FROM AUTHOR]

Tumor mutational signatures have gained prominence in cancer research, yet the lack of standardized methods hinders reproducibility and robustness. Leveraging colorectal cancer (CRC) as a model, we explored the influence of computational parameters on mutational signature analyses across 230 CRC cell lines and 152 CRC patients. Results were validated in three independent datasets: 483 endometrial cancer patients stratified by mismatch repair (MMR) status, 35 lung cancer patients by smoking status and 12 patient-derived organoids (PDOs) annotated for colibactin exposure. Assessing various bioinformatic tools, reference datasets and input data sizes including whole genome sequencing, whole exome sequencing and a pan-cancer gene panel, we demonstrated significant variability in the results. We report that the use of distinct algorithms and references led to statistically different results, highlighting how arbitrary choices may induce variability in the mutational signature contributions. Furthermore, we found a differential contribution of mutational signatures between coding and intergenic regions and defined the minimum number of somatic variants required for reliable mutational signature assignment. To facilitate the identification of the most suitable workflows, we developed Comparative Mutational Signature analysis on Coding and Extragenic Regions (CoMSCER), a bioinformatic tool which allows researchers to easily perform comparative mutational signature analysis by coupling the results from several tools and public reference datasets and to assess mutational signature contributions in coding and non-coding genomic regions. In conclusion, our study provides a comparative framework to elucidate the impact of distinct computational workflows on mutational signatures. [ABSTRACT FROM AUTHOR]

This paper presents a neutrosophic inference model for bioinformatics. The model is used to develop a system for accurate comparisons of human nucleic acids, where the new nucleic acid is compared to a database of old nucleic acids. The comparisons are analyzed in terms of accuracy, certainty, uncertainty, neutrality, and bias. The proposed system achieves good results and provides a reliable standard for future comparisons. It highlights the potential of neutrosophic inference models in bioinformatics applications. Data mining and bioinformatics play a crucial role in computational biology, with applications in scientific research and industrial development. Biological analysts rely on specialized tools and algorithms to collect, store, categorize, and analyze large volumes of unstructured data. Data mining techniques are used to extract valuable information from this data, aiding in the development of new therapies and understanding genetic relationships between organisms. Recent advancements in bioinformatics include gene expression tools, Bio sequencing, and Bio databases, which facilitate the extraction and analysis of vital biological information. These technologies contribute to the analysis of big data, identification of key bioinformatics insights, and generation of new biological knowledge. Data collection, analysis, and interpretation in this field involves the use of modern technologies such as cloud computing, machine learning, and artificial intelligence, enabling more efficient and accurate results. Ultimately, data mining and bioinformatics enhance our understanding of genetic relationships, aid in developing new therapies, and improve healthcare outcomes. [ABSTRACT FROM AUTHOR]

Bioinformatics and computational biology are two fields that have been exploiting GPUs for more than two decades, with being CUDA the most used programming language for them. However, as CUDA is an NVIDIA proprietary language, it implies a strong portability restriction to a wide range of heterogeneous architectures, like AMD or Intel GPUs. To face this issue, the Khronos group has recently proposed the SYCL standard, which is an open, royalty-free, cross-platform abstraction layer that enables the programming of a heterogeneous system to be written using standard, single-source C++ code. Over the past few years, several implementations of this SYCL standard have emerged, being oneAPI the one from Intel. This paper presents the migration process of the SW# suite, a biological sequence alignment tool developed in CUDA, to SYCL using Intel's oneAPI ecosystem. The experimental results show that SW# was completely migrated with a small programmer intervention in terms of hand-coding. In addition, it was possible to port the migrated code between different architectures (considering multiple vendor GPUs and also CPUs), with no noticeable performance degradation on five different NVIDIA GPUs. Moreover, performance remained stable when switching to another SYCL implementation. As a consequence, SYCL and its implementations can offer attractive opportunities for the bioinformatics community, especially considering the vast existence of CUDA-based legacy codes. [ABSTRACT FROM AUTHOR]

Computational methods have significantly transformed biomedical research, offering a comprehensive exploration of disease mechanisms and molecular protein functions. This article reviews a spectrum of computational tools and network analysis databases that play a crucial role in identifying potential interactions and signaling networks contributing to the onset of disease states. The utilization of protein/gene interaction and genetic variation databases, coupled with pathway analysis can facilitate the identification of potential drug targets. By bridging the gap between molecular-level information and disease understanding, this review contributes insights into the impactful utilization of computational methods, paving the way for targeted interventions and therapeutic advancements in biomedical research. [ABSTRACT FROM AUTHOR]

Darwin recognized that 'a grand and almost untrodden field of inquiry will be opened, on the causes and laws of variation.' However, because the Modern Synthesis assumes that the intrinsic probability of any individual mutation is unrelated to that mutation's potential adaptive value, attention has been focused on selection rather than on the intrinsic generation of variation. Yet many examples illustrate that the term 'random' mutation, as widely understood, is inaccurate. The probabilities of distinct classes of variation are neither evenly distributed across a genome nor invariant over time, nor unrelated to their potential adaptive value. Because selection acts upon variation, multiple biochemical mechanisms can and have evolved that increase the relative probability of adaptive mutations. In effect, the generation of heritable variation is in a feedback loop with selection, such that those mechanisms that tend to generate variants that survive recurring challenges in the environment would be captured by this survival and thus inherited and accumulated within lineages of genomes. Moreover, because genome variation is affected by a wide range of biochemical processes, genome variation can be regulated. Biochemical mechanisms that sense stress, from lack of nutrients to DNA damage, can increase the probability of specific classes of variation. A deeper understanding of evolution involves attention to the evolution of, and environmental influences upon, the intrinsic variation generated in gametes, in other words upon the biochemical mechanisms that generate variation across generations. These concepts have profound implications for the types of questions that can and should be asked, as omics databases become more comprehensive, detection methods more sensitive, and computation and experimental analyses even more high throughput and thus capable of revealing the intrinsic generation of variation in individual gametes. These concepts also have profound implications for evolutionary theory, which, upon reflection it will be argued, predicts that selection would increase the probability of generating adaptive mutations, in other words, predicts that the ability to evolve itself evolves. [ABSTRACT FROM AUTHOR]

Objective(s): MicroRNAs, which are micro-coordinators of gene expression, have been recently investigated as a potential treatment for cancer. The study used computational techniques to identify microRNAs that could target a set of genes simultaneously. Due to their multi-target-directed nature, microRNAs have the potential to impact multiple key pathways and their pathogenic cross-talk. Materials and Methods: We identified microRNAs that target a prostate cancer-associated gene set using integrated bioinformatics analyses and experimental validation. The candidate gene set included genes targeted by clinically approved prostate cancer medications. We used STRING, GO, and KEGG web tools to confirm gene-gene interactions and their clinical significance. Then, we employed integrated predicted and validated bioinformatics approaches to retrieve hsa-miR-124-3p, 16-5p, and 27a-3p as the top three relevant microRNAs. KEGG and DIANA-miRPath showed the related pathways for the candidate genes and microRNAs Results: The Real-time PCR results showed that miR-16-5p simultaneously down-regulated all genes significantly except for PIK3CA/CB in LNCaP; miR-27a-3p simultaneously down-regulated all genes significantly, excluding MET in LNCaP and PIK3CA in PC-3; and miR-124-3p could not downregulate significantly PIK3CB, MET, and FGFR4 in LNCaP and FGFR4 in PC-3. Finally, we used a cell cycle assay to show significant G0/G1 arrest by transfecting miR-124-3p in LNCaP and miR-16-5p in both cell lines. Conclusion: Our findings suggest that this novel approach may have therapeutic benefits and these predicted microRNAs could effectively target the candidate genes. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) have emerged as a prominent layer of regulation of gene expression. This article offers the salient and current aspects of machine learning (ML) tools and approaches from genome to phenome in miRNA research. First, we underline that the complexity in the analysis of miRNA function ranges from their modes of biogenesis to the target diversity in diverse biological conditions. Therefore, it is imperative to first ascertain the miRNA coding potential of genomes and understand the regulatory mechanisms of their expression. This knowledge enables the efficient classification of miRNA precursors and the identification of their mature forms and respective target genes. Second, and because one miRNA can target multiple mRNAs and vice versa, another challenge is the assessment of the miRNA-mRNA target interaction network. Furthermore, long-noncoding RNA (lncRNA)and circular RNAs (circRNAs) also contribute to this complexity. ML has been used to tackle these challenges at the high-dimensional data level. The present expert review covers more than 100 tools adopting various ML approaches pertaining to, for example, (1) miRNA promoter prediction, (2) precursor classification, (3) mature miRNA prediction, (4) miRNA target prediction, (5) miRNA- lncRNA and miRNA-circRNA interactions, (6) miRNA-mRNA expression profiling, (7) miRNA regulatory module detection, (8) miRNA-disease association, and (9) miRNA essentiality prediction. Taken together, we unpack, critically examine, and highlight the cutting-edge synergy of ML approaches and miRNA research so as to develop a dynamic and microlevel understanding of human health and diseases. [ABSTRACT FROM AUTHOR]

Bioinformatics has revolutionized biology and medicine by using computational methods to analyze and interpret biological data. Quantum mechanics has recently emerged as a promising tool for the analysis of biological systems, leading to the development of quantum bioinformatics. This new field employs the principles of quantum mechanics, quantum algorithms, and quantum computing to solve complex problems in molecular biology, drug design, and protein folding. However, the intersection of bioinformatics, biology, and quantum mechanics presents unique challenges. One significant challenge is the possibility of confusion among scientists between quantum bioinformatics and quantum biology, which have similar goals and concepts. Additionally, the diverse calculations in each field make it difficult to establish boundaries and identify purely quantum effects from other factors that may affect biological processes. This review provides an overview of the concepts of quantum biology and quantum mechanics and their intersection in quantum bioinformatics. We examine the challenges and unique features of this field and propose a classification of quantum bioinformatics to promote interdisciplinary collaboration and accelerate progress. By unlocking the full potential of quantum bioinformatics, this review aims to contribute to our understanding of quantum mechanics in biological systems. [ABSTRACT FROM AUTHOR]

Bioinformatics resolves biological problems by predicting the structural-functional annotation of genes and substantially contributes to further downstream applications in the field of genomics. Most of the proteins from Colletotrichum gloeosporoides are uncharacterized yet, which could be of novel biological functions and help to understand the biological function. Hence, in this study, we aimed to annotate the zinc-finger dhhc domain containing uncharacterized protein A0A8H4C4Q8_COLGL of C. gloeosporoides. The uncharacterized protein was analyzed for subcellular localization, physicochemical characteristics, secretory nature, antigenicity, allergenicity determination, conserved domain identification, secondary structure prediction, 3D modeling, docking analysis, active site prediction, and protein-protein interaction. The structural-functional annotation predicted it as an effector protein possessing a role in pathogenicity. [ABSTRACT FROM AUTHOR]

The most critical step in the clinical diagnosis workflow is the pathological evaluation of each tumor sample. Deep learning is a powerful approach that is widely used to enhance diagnostic accuracy and streamline the diagnosis process. In our previous study using omics data, we identified 2 distinct subtypes of pure seminoma. Seminoma is the most common histological type of testicular germ cell tumors (TGCTs). Here we developed a deep learning decision making tool for the identification of seminoma subtypes using histopathological slides. We used all available slides for pure seminoma samples from The Cancer Genome Atlas (TCGA). The developed model showed an area under the ROC curve of 0.896. Our model not only confirms the presence of 2 distinct subtypes within pure seminoma but also unveils the presence of morphological differences between them that are imperceptible to the human eye. [ABSTRACT FROM AUTHOR]

Simple Summary: The effects of radiation exposure seem closely related to effects of old age—so much so that the idea of a radiation–age association came about in the 1960s. While not a new idea, modern technology is allowing us to revisit these ideas and explore them with a fresh perspective. Separately, there are gaps in the community's understanding of the effects of radiation and aging, such as with respect to low-level, long-term effects of radiation and estimating someone's biological age. To study their association, a number of tools exist that need to be efficiently integrated to study this complex and interdisciplinary field. This article includes an extensive literature review on the theory of these two topics, providing a detailed foundation for a current understanding. We then present a resource-agnostic approach for researchers in these areas, focusing on studying the association between the two. Primary points of interest are focused on indirect damage of radiation exposure via oxidative stress within a cell, a comprehensive table of functional estimators for biological age, and using modern computational tools and biology to overlap fields of study to develop and exploit a rad–age association. Multi-omics studies have emerged as powerful tools for tailoring individualized responses to various conditions, capitalizing on genome sequencing technologies' increasing affordability and efficiency. This paper delves into the potential of multi-omics in deepening our understanding of biological age, examining the techniques available in light of evolving technology and computational models. The primary objective is to review the relationship between ionizing radiation and biological age, exploring a wide array of functional, physiological, and psychological parameters. This comprehensive review draws upon an extensive range of sources, including peer-reviewed journal articles, government documents, and reputable websites. The literature review spans from fundamental insights into radiation effects to the latest developments in aging research. Ionizing radiation exerts its influence through direct mechanisms, notably single- and double-strand DNA breaks and cross links, along with other critical cellular events. The cumulative impact of DNA damage forms the foundation for the intricate process of natural aging, intersecting with numerous diseases and pivotal biomarkers. Furthermore, there is a resurgence of interest in ionizing radiation research from various organizations and countries, reinvigorating its importance as a key contributor to the study of biological age. Biological age serves as a vital reference point for the monitoring and mitigation of the effects of various stressors, including ionizing radiation. Ionizing radiation emerges as a potent candidate for modeling the separation of biological age from chronological age, offering a promising avenue for tailoring protocols across diverse fields, including the rigorous demands of space exploration. [ABSTRACT FROM AUTHOR]

Enhancers are short DNA segments (50–1500 bp) that effectively activate gene transcription when transcription factors (TFs) are present. There is a correlation between the genetic differences in enhancers and numerous human disorders including cancer and inflammatory bowel disease. In computational biology, the accurate categorization of enhancers can yield important information for drug discovery and development. High-throughput experimental approaches are thought to be vital tools for researching enhancers' key characteristics; however, because these techniques require a lot of labor and time, it might be difficult for researchers to forecast enhancers and their powers. Therefore, computational techniques are considered an alternate strategy for handling this issue. Based on the types of algorithms that have been used to construct predictors, the current methodologies can be divided into three primary categories: ensemble-based methods, deep learning-based approaches, and traditional ML-based techniques. In this study, we developed a novel two-layer deep forest-based predictor for accurate enhancer and strength prediction, namely, NEPERS. Enhancers and non-enhancers are divided at the first level by NEPERS, whereas strong and weak enhancers are divided at the second level. To evaluate the effectiveness of feature fusion, block-wise deep forest and other algorithms were combined with multi-view features such as PSTNPss, PSTNPdss, CKSNAP, and NCP via 10-fold cross-validation and independent testing. Our proposed technique performs better than competing models across all parameters, with an ACC of 0.876, Sen of 0.864, Spe of 0.888, MCC of 0.753, and AUC of 0.940 for layer 1 and an ACC of 0.959, Sen of 0.960, Spe of 0.958, MCC of 0.918, and AUC of 0.990 for layer 2, respectively, for the benchmark dataset. Similarly, for the independent test, the ACC, Sen, Spe, MCC, and AUC were 0.863, 0.865, 0.860, 0.725, and 0.948 for layer 1 and 0.890, 0.940, 0.840, 0.784, and 0.951 for layer 2, respectively. This study provides conclusive insights for the accurate and effective detection and characterization of enhancers and their strengths. [ABSTRACT FROM AUTHOR]

Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA–miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR. [ABSTRACT FROM AUTHOR]

Understanding binding related changes in antibody conformations is important for epitope prediction and antibody refinement. The increase of available data in the PDB allowed a more detailed investigation of the conformational landscape for free and bound antibodies. A dataset containing a total of 835 unique PDB entries of antibodies that were crystallized in complex with their antigen and in a free state was constructed. It was examined for binding related conformation changes. We present further evidence supporting the theory of a pre-existing-equilibrium in experimental data. Multiple sequence alignments did not show binding induced tendencies in the solvent accessibility of residues in any specific position. Evaluating the changes in solvent accessibility per residue revealed a certain binding induced increase for several amino acids. Antibody-antigen interaction statistics were established and quantify a significant directional asymmetry between many interacting antibody and antigen residue pairs, especially a richness in tyrosine in the antibody epitope compared to its paratope. This asymmetry could potentially facilitate an increase in the success rate of computationally guided antibody refinement. [ABSTRACT FROM AUTHOR]

The oral microbiota has an enormous impact on human health, with oral dysbiosis now linked to many oral and systemic diseases. Recent advancements in sequencing, mass spectrometry, bioinformatics, computational biology, and machine learning are revolutionizing oral microbiome research, enabling analysis at an unprecedented scale and level of resolution using omics approaches. This review contains a comprehensive perspective of the current state-of-the-art tools available to perform genomics, metagenomics, phylogenomics, pangenomics, transcriptomics, proteomics, metabolomics, lipidomics, and multi-omics analysis on (all) microbiomes, and then provides examples of how the techniques have been applied to research of the oral microbiome, specifically. Key findings of these studies and remaining challenges for the field are highlighted. Although the methods discussed here are placed in the context of their contributions to oral microbiome research specifically, they are pertinent to the study of any microbiome, and the intended audience of this includes researchers would simply like to get an introduction to microbial omics and/or an update on the latest omics methods. Continued research of the oral microbiota using omics approaches is crucial and will lead to dramatic improvements in human health, longevity, and quality of life. [ABSTRACT FROM AUTHOR]

In the machine learning field, the technique known as ensemble learning aims at combining different base learners in order to increase the quality and the robustness of the predictions. Indeed, this approach has widely been applied to tackle, with success, real world problems from different domains, including computational biology. Nevertheless, despite their potential, ensembles combining results from different base learners have been understudied in the context of gene regulatory network inference. In this paper we applied genetic algorithms and frequent itemset mining, to design small but effective ensembles of gene regulatory network inference methods. These ensembles were evaluated and compared to well-established single and ensemble methods, on both real and synthetic datasets. Results showed that small ensembles, consisting of few but diverse base learners, enhance the exploration of the solution space, and compensate base learners biases, outperforming state-of-the-art methods. Results advocate for the use of such methods as gene regulatory network inference tools. [ABSTRACT FROM AUTHOR]

The bioinformatics discipline seeks to solve problems in biology with computational theories and methods. Formal concept analysis (FCA) is one such theoretical model, based on partial orders. FCA allows the user to examine the structural properties of data based on which subsets of the dataset depend on each other. This article surveys the current literature related to the use of FCA for bioinformatics. The survey begins with a discussion of FCA, its hierarchical advantages, several advanced models of FCA, and lattice management strategies. It then examines how FCA has been used in bioinformatics applications, followed by future prospects of FCA in those areas. The applications addressed include gene data analysis (with next-generation sequencing), biomarkers discovery, protein-protein interaction, disease analysis (including COVID-19, cancer, and others), drug design and development, healthcare informatics, biomedical ontologies, and phylogeny. Some of the most promising prospects of FCA are identifying influential nodes in a network representing protein-protein interactions, determining critical concepts to discover biomarkers, integrating machine learning and deep learning for cancer classification, and pattern matching for next-generation sequencing. [ABSTRACT FROM AUTHOR]

Deep learning, a potent branch of artificial intelligence, is steadily leaving its transformative imprint across multiple disciplines. Within computational biology, it is expediting progress in the understanding of Protein–Protein Interactions (PPIs), key components governing a wide array of biological functionalities. Hence, an in-depth exploration of PPIs is crucial for decoding the intricate biological system dynamics and unveiling potential avenues for therapeutic interventions. As the deployment of deep learning techniques in PPI analysis proliferates at an accelerated pace, there exists an immediate demand for an exhaustive review that encapsulates and critically assesses these novel developments. Addressing this requirement, this review offers a detailed analysis of the literature from 2021 to 2023, highlighting the cutting-edge deep learning methodologies harnessed for PPI analysis. Thus, this review stands as a crucial reference for researchers in the discipline, presenting an overview of the recent studies in the field. This consolidation helps elucidate the dynamic paradigm of PPI analysis, the evolution of deep learning techniques, and their interdependent dynamics. This scrutiny is expected to serve as a vital aid for researchers, both well-established and newcomers, assisting them in maneuvering the rapidly shifting terrain of deep learning applications in PPI analysis. [ABSTRACT FROM AUTHOR]

Simple Summary: Next-generation sequencing (NGS) is a cutting-edge technology that enables rapid, high-throughput sequencing of DNA and RNA. Researchers and clinicians can identify genetic mutations, gene fusions, and other alterations that may drive cancer growth. This is particularly important in precision oncology as it is applied in the context of Molecular Tumor Boards (MTBs). The latter are multidisciplinary teams of experts who use NGS and bioinformatics tools to analyze patients' genetic profiles and develop personalized treatment recommendations for cancer patients. Thus, a crucial process for MTB decision-making is the analysis, compilation, and presentation of high-dimensional sequencing data, which are used for both preparation of and case presentation to all stakeholders. MIRACUM-Pipe precisely addresses these requirements and offers an easy-to-use, one-prompt standardized solution to analyze NGS data, including quality control, variant calling, copy number estimation, annotation, visualization, and report generation. (1) Background: Next-generation sequencing (NGS) of patients with advanced tumors is becoming an established method in Molecular Tumor Boards. However, somatic variant detection, interpretation, and report generation, require in-depth knowledge of both bioinformatics and oncology. (2) Methods: MIRACUM-Pipe combines many individual tools into a seamless workflow for comprehensive analyses and annotation of NGS data including quality control, alignment, variant calling, copy number variation estimation, evaluation of complex biomarkers, and RNA fusion detection. (3) Results: MIRACUM-Pipe offers an easy-to-use, one-prompt standardized solution to analyze NGS data, including quality control, variant calling, copy number estimation, annotation, visualization, and report generation. (4) Conclusions: MIRACUM-Pipe, a versatile pipeline for NGS, can be customized according to bioinformatics and clinical needs and to support clinical decision-making with visual processing and interactive reporting. [ABSTRACT FROM AUTHOR]

Background: Despite being in the 21st century, the world has still not been able to vanquish the global AIDS epidemic, and the only foreseeable solution seems to be a safe and effective vaccine. Unfortunately, vaccine trials so far have returned unfruitful results, possibly due to their inability to induce effective cellular, humoral and innate immune responses. The current study aims to tackle these limitations and propose the desired vaccine utilizing immunoinformatic approaches that have returned promising results in designing vaccines against various rapidly mutating organisms. For this, all polyprotein and protein sequences of HIV-1 were retrieved from the LANL (Los Alamos National Laboratory) database. The consensus sequence was generated after alignment and used to predict epitopes. Conserved, antigenic, non-allergenic, T-cell inducing, B-cell inducing, IFN-ɣ inducing, non-human homologous epitopes were selected and combined to propose two vaccine constructs i.e., HIV-1a (without adjuvant) and HIV-1b (with adjuvant). Results: HIV-1a and HIV-1b were subjected to antigenicity, allergenicity, structural quality analysis, immune simulations, and MD (molecular dynamics) simulations. Both proposed multi-epitope vaccines were found to be antigenic, non-allergenic, stable, and induce cellular, humoral, and innate immune responses. TLR-3 docking and in-silico cloning of both constructs were also performed. Conclusion: Our results indicate HIV-1b to be more promising than HIV-1a; experimental validations can confirm the efficacy and safety of both constructs and in-vivo efficacy in animal models. [ABSTRACT FROM AUTHOR]

A comprehensive overview of various omics techniques to understand the biological responsive to engineered biomaterials is provided. Specifically, the role of omics-based approaches in designing next-generation of bioresponsive materials for various regenerative medicine approaches, as well as future outlook on integrated multi-omics platform are discussed. [Display omitted] Recent developments in biomedical engineering have focused on designing smart and bio-responsive materials. However, a critical step in designing this next generation of biomaterials is to evaluate their holistic cellular and molecular interactions. Recent technological advances in biology, specifically in omics techniques, can provide high-throughput functional readouts of the cell genome, epigenome, transcriptome, proteome and metabolome. The ability of omics-based approaches to predict the biological response to engineered biomaterials have the potential to revolutionize biomedical research and replace the traditional "trial and error" approach as the first step in designing the next-generation of biomaterials. In addition, these omics-based approaches can be used to gain insight into in vivo responses of biomaterials, which currently relies heavily on semi-quantitative imaging-based methodologies rather than rigorous computational approaches. In this review, we have outlined various omics techniques that have been utilized to quantify and understand the biological performance of engineered biomaterials. Additionally, we have critically discussed the role of omics-based approaches in designing new biomaterials for regenerative medicine, immune engineering, and drug delivery. Finally, the potential of integrated multi-omics approaches to build a holistic understanding of the biological responses to newly developed biomaterials have been discussed. [ABSTRACT FROM AUTHOR]

Omics data are multidimensional, heterogeneous, and high throughput. Robust computational methods and machine learning (ML)-based models offer new prospects to accelerate the data-to-knowledge trajectory. Deep learning (DL) is a powerful subset of ML inspired by brain structure and has created unprecedented momentum in bioinformatics and computational biology research. This article provides an overview of the current DL models applied to multi-omics data for both the beginner and the expert user. Additionally, COVID-19 will continue to impact planetary health as a pandemic and an endemic disease, with genomic and multi-omic pathophysiology. DL offers, therefore, new ways of harnessing systems biology research on COVID-19 diagnostics and therapeutics. Herein, we discuss, first, the statistical ML algorithms and essential deep architectures. Then, we review DL applications in multi-omics data analysis and their intersection with COVID-19. Finally, challenges and several promising directions are highlighted going forward in the current era of COVID-19. [ABSTRACT FROM AUTHOR]

Signal transducer and activator of transcription 6 (STAT6) is a multifunctional protein that plays critical functions in cell proliferation, apoptosis, differentiation, and angiogenesis. Mutations in STAT6 may contribute to the development of certain complex diseases such as cancer. This study examined single amino acid substitutions in STAT6 to pinpoint deleterious variants and their related structural and functional impairments. Data on STAT6 mutations were obtained from the Ensembl database and analyzed to evaluate the selected mutations for their pathogenicity and destabilizing or harmful effects. Specifically, we analyzed aggregation propensity, nonpacking density, and accessible surface area on the chosen mutations. The results suggest that seven out of eight mutations are less soluble, which might lead to aggregation, disrupt ordered helices, and alter strand propensity. Four mutations lay in the conserved regions of the protein, as revealed by the Consurf analysis. We found that three mutations, E318G, L365F, and R562H, change hydrophobic contacts and lead to frustration of STAT6, which can alter its stability, contributing to disease progression in cancer. In conclusion, these findings inform how single amino acid changes can destabilize STAT6. This has implications for cancer progression which warrants further experimental research. [ABSTRACT FROM AUTHOR]

Summary It is challenging to simulate realistic tracts of genetic ancestry on a scale suitable for simulation-based inference. We present an algorithm that enables this information to be extracted efficiently from tree sequences produced by simulations run with msprime and SLiM. Availability and implementation A C-based implementation of the link-ancestors algorithm is in tskit (https://tskit.dev/tskit/docs/stable/). We also provide a user-friendly wrapper for link-ancestors in tspop, a Python-based utility package. [ABSTRACT FROM AUTHOR]

Motivation There is a rapid growth in the production of omics datasets collected by the diabetes research community. However, such published data are underutilized for knowledge discovery. To make bioinformatics tools and published omics datasets from the diabetes field more accessible to biomedical researchers, we developed the Diabetes Data and Hypothesis Hub (D2H2). Results D2H2 contains hundreds of high-quality curated transcriptomics datasets relevant to diabetes, accessible via a user-friendly web-based portal. The collected and processed datasets are curated from the Gene Expression Omnibus (GEO). Each curated study has a dedicated page that provides data visualization, differential gene expression analysis, and single-gene queries. To enable the investigation of these curated datasets and to provide easy access to bioinformatics tools that serve gene and gene set-related knowledge, we developed the D2H2 chatbot. Utilizing GPT, we prompt users to enter free text about their data analysis needs. Parsing the user prompt, together with specifying information about all D2H2 available tools and workflows, we answer user queries by invoking the most relevant tools via the tools' API. D2H2 also has a hypotheses generation module where gene sets are randomly selected from the bulk RNA-seq precomputed signatures. We then find highly overlapping gene sets extracted from publications listed in PubMed Central with abstract dissimilarity. With the help of GPT, we speculate about a possible explanation of the high overlap between the gene sets. Overall, D2H2 is a platform that provides a suite of bioinformatics tools and curated transcriptomics datasets for hypothesis generation. Availability and implementation D2H2 is available at: https://d2h2.maayanlab.cloud/ and the source code is available from GitHub at https://github.com/MaayanLab/D2H2-site under the CC BY-NC 4.0 license. [ABSTRACT FROM AUTHOR]

Motivation Development of bioinformatics methods is a long, complex and resource-hungry process. Hundreds of these tools were released. While some methods are highly cited and used, many suffer relatively low citation rates. We empirically analyze a large collection of recently released methods in three diverse protein function and disorder prediction areas to identify key factors that contribute to increased citations. Results We show that provision of a working web server significantly boosts citation rates. On average, methods with working web servers generate three times as many citations compared to tools that are available as only source code, have no code and no server, or are no longer available. This observation holds consistently across different research areas and publication years. We also find that differences in predictive performance are unlikely to impact citation rates. Overall, our empirical results suggest that a relatively low-cost investment into the provision and long-term support of web servers would substantially increase the impact of bioinformatics tools. [ABSTRACT FROM AUTHOR]

Motivation In recent years, network models have gained prominence for their ability to capture complex associations. In statistical omics, networks can be used to model and study the functional relationships between genes, proteins, and other types of omics data. If a Gaussian graphical model is assumed, a gene association network can be determined from the non-zero entries of the inverse covariance matrix of the data. Due to the high-dimensional nature of such problems, integrative methods that leverage similarities between multiple graphical structures have become increasingly popular. The joint graphical lasso is a powerful tool for this purpose, however, the current AIC-based selection criterion used to tune the network sparsities and similarities leads to poor performance in high-dimensional settings. Results We propose stabJGL, which equips the joint graphical lasso with a stable and well-performing penalty parameter selection approach that combines the notion of model stability with likelihood-based similarity selection. The resulting method makes the powerful joint graphical lasso available for use in omics settings, and outperforms the standard joint graphical lasso, as well as state-of-the-art joint methods, in terms of all performance measures we consider. Applying stabJGL to proteomic data from a pan-cancer study, we demonstrate the potential for novel discoveries the method brings. Availability and implementation A user-friendly R package for stabJGL with tutorials is available on Github https://github.com/Camiling/stabJGL. [ABSTRACT FROM AUTHOR]

Summary The creation of effective visualizations is a fundamental component of data analysis. In biomedical research, new challenges are emerging to visualize multi-dimensional data in a 2D space, but current data visualization tools have limited capabilities. To address this problem, we leverage Gestalt principles to improve the design and interpretability of multi-dimensional data in 2D data visualizations, layering aesthetics to display multiple variables. The proposed visualization can be applied to spatially-resolved transcriptomics data, but also broadly to data visualized in 2D space, such as embedding visualizations. We provide an open source R package escheR , which is built off of the state-of-the-art ggplot2 visualization framework and can be seamlessly integrated into genomics toolboxes and workflows. Availability and implementation The open source R package escheR is freely available on Bioconductor (https://bioconductor.org/packages/escheR). [ABSTRACT FROM AUTHOR]

Summary PICAFlow is a R-written integrative workflow dedicated to flow/mass cytometry data handling, from pre-processing to deep and comprehensive analysis. It is designed as a powerful all-in-one tool which contains all the necessary functions and packages presented in a user-friendly and ease-to-use fashion. PICAFlow also includes important features that are very frequently lacking in other close software, such as interactive R Shiny applications for real-time data transformation and compensation as well as normalization methods aiming to remove batch effects and unwanted inter- and intra-group heterogeneity. It also allows to perform dimensionality reduction, cell clustering (using different available approaches), as well as complementary statistical analyses and export different support for data interpretation and visualization. Availability PICAFlow is available as a R-written package hosted at the following GitHub repository: https://github.com/PaulRegnier/PICAFlow and is complemented by a fully detailed tutorial available at the following URL: https://paul-regnier.fr/tutoriel-picaflow/. [ABSTRACT FROM AUTHOR]

Motivation Single-cell technologies allow deep characterization of different molecular aspects of cells. Integrating these modalities provides a comprehensive view of cellular identity. Current integration methods rely on overlapping features or cells to link datasets measuring different modalities, limiting their application to experiments where different molecular layers are profiled in different subsets of cells. Results We present scTopoGAN, a method for unsupervised manifold alignment of single-cell datasets with non-overlapping cells or features. We use topological autoencoders (topoAE) to obtain latent representations of each modality separately. A topology-guided Generative Adversarial Network then aligns these latent representations into a common space. We show that scTopoGAN outperforms state-of-the-art manifold alignment methods in complete unsupervised settings. Interestingly, the topoAE for individual modalities also showed better performance in preserving the original structure of the data in the low-dimensional representations when compared to other manifold projection methods. Taken together, we show that the concept of topology preservation might be a powerful tool to align multiple single modality datasets, unleashing the potential of multi-omic interpretations of cells. Availability and implementation Implementation available on GitHub (https://github.com/AkashCiel/scTopoGAN). All datasets used in this study are publicly available. [ABSTRACT FROM AUTHOR]

Summary ExRec (Exclusion of Recombined DNA) is a dependency-free Python pipeline that implements the four-gamete test to automatically filter out recombined DNA blocks from thousands of DNA sequence loci. This procedure helps all loci better meet the "no intralocus recombination" assumption common to many coalescent-based analyses in population genomic, phylogeographic, and shallow-scale phylogenomic studies. The user-friendly pipeline contains five standalone applications—four file conversion scripts and one main script that performs the recombination filtering procedures. The pipeline outputs recombination-filtered data in a variety of common formats and a tab-delimited table that displays descriptive statistics for all loci and the analysis results. A novel feature of this software is that the user can select whether to output the longest nonrecombined sequence blocks from recombined loci (current best practice) or randomly select nonrecombined blocks from loci (a newer approach). We tested ExRec with six published phylogenomic datasets that ranged in size from 27 to 2237 loci and came in a variety of input file formats. In all trials the data could be easily analyzed in only seconds for the smaller datasets and <30 min for the largest using a simple laptop computer. Availability and implementation ExRec was written in Python 3 under the MIT license. The program applications, user manual (including step-by-step tutorials), and sample data are freely available at https://github.com/Sammccarthypotter/ExRec. [ABSTRACT FROM AUTHOR]

Motivation Enzymes are key targets to biosynthesize functional substances in metabolic engineering. Therefore, various machine learning models have been developed to predict Enzyme Commission (EC) numbers, one of the enzyme annotations. However, the previously reported models might predict the sequences with numerous consecutive identical amino acids, which are found within unannotated sequences, as enzymes. Results Here, we propose EnzymeNet for prediction of complete EC numbers using residual neural networks. EnzymeNet can exclude the exceptional sequences described above. Several EnzymeNet models were built and optimized to explore the best conditions for removing such sequences. As a result, the models exhibited higher prediction accuracy with macro F 1 score up to 0.850 than previously reported models. Moreover, even the enzyme sequences with low similarity to training data, which were difficult to predict using the reported models, could be predicted extensively using EnzymeNet models. The robustness of EnzymeNet models will lead to discover novel enzymes for biosynthesis of functional compounds using microorganisms. Availability and implementation The source code of EnzymeNet models is freely available at https://github.com/nwatanbe/enzymenet. [ABSTRACT FROM AUTHOR]

Summary Reconstructing haplotypes of an organism from a set of sequencing reads is a computationally challenging (NP-hard) problem. In reference-guided settings, at the core of haplotype assembly is the task of clustering reads according to their origin, i.e. grouping together reads that sample the same haplotype. Read length limitations and sequencing errors render this problem difficult even for diploids; the complexity of the problem grows with the ploidy of the organism. We present XHap, a novel method for haplotype assembly that aims to learn correlations between pairs of sequencing reads, including those that do not overlap but may be separated by large genomic distances, and utilize the learned correlations to assemble the haplotypes. This is accomplished by leveraging transformers, a powerful deep-learning technique that relies on the attention mechanism to discover dependencies between non-overlapping reads. Experiments on semi-experimental and real data demonstrate that the proposed method significantly outperforms state-of-the-art techniques in diploid and polyploid haplotype assembly tasks on both short and long sequencing reads. Availability and implementation The code for XHap and the included experiments is available at https://github.com/shoryaconsul/XHap. [ABSTRACT FROM AUTHOR]

Motivation Compared to eukaryotes, prokaryote genomes are more diverse through different mechanisms, including a higher mutation rate and horizontal gene transfer. Therefore, using a linear representative reference can cause a reference bias. Graph-based pangenome methods have been developed to tackle this problem. However, comparisons in DNA space are still challenging due to this high diversity. In contrast, amino acid sequences have higher similarity due to evolutionary constraints, whereby a single amino acid may be encoded by several synonymous codons. Coding regions cover the majority of the genome in prokaryotes. Thus, panproteomes present an attractive alternative leveraging the higher sequence similarity while not losing much of the genome in non-coding regions. Results We present PanPA , a method that takes a set of multiple sequence alignments of protein sequences, indexes them, and builds a graph for each multiple sequence alignment. In the querying step, it can align DNA or amino acid sequences back to these graphs. We first showcase that PanPA generates correct alignments on a panproteome from 1350 Escherichia coli. To demonstrate that panproteomes allow comparisons at longer phylogenetic distances, we compare DNA and protein alignments from 1073 Salmonella enterica assemblies against E.coli reference genome, pangenome, and panproteome using BWA , GraphAligner , and PanPA , respectively; with PanPA aligning around 22% more sequences. We also aligned a DNA short-reads whole genome sequencing (WGS) sample from S.enterica against the E.coli reference with BWA and the panproteome with PanPA , where PanPA was able to find alignment for 68% of the reads compared to 5% with BWA. Availalability and implementation PanPA is available at https://github.com/fawaz-dabbaghieh/PanPA. [ABSTRACT FROM AUTHOR]

Summary Quantifying genetic clusters (=populations) from genotypic data is a fundamental, but non-trivial task for population geneticists that is compounded by: hierarchical population structure, diverse analytical methods, and complex software dependencies. AdmixPipe v3 ameliorates many of these issues in a single bioinformatic pipeline that facilitates all facets of population structure analysis by integrating outputs generated by several popular packages (i.e. CLUMPAK, EvalAdmix). The pipeline interfaces disparate software packages to parse Admixture outputs and conduct EvalAdmix analyses in the context of multimodal population structure results identified by CLUMPAK. We further streamline these tasks by packaging AdmixPipe v3 within a Docker container to create a standardized analytical environment that allows for complex analyses to be replicated by different researchers. This also grants operating system flexibility and mitigates complex software dependencies. Availability and implementation Source code, documentation, example files, and usage examples are freely available at https://github.com/stevemussmann/admixturePipeline. Installation is facilitated via Docker container available from https://hub.docker.com/r/mussmann/admixpipe. Usage under Windows operating systems requires the Windows Subsystem for Linux. [ABSTRACT FROM AUTHOR]

Motivation Predictive computational models must be accurate, robust, and interpretable to be considered reliable in important areas such as biology and medicine. A sufficiently robust model should not have its output affected significantly by a slight change in the input. Also, these models should be able to explain how a decision is made to support user trust in the results. Efforts have been made to improve the robustness and interpretability of predictive computational models independently; however, the interaction of robustness and interpretability is poorly understood. Results As an example task, we explore the computational prediction of cell type based on single-cell RNA-seq data and show that it can be made more robust by adversarially training a deep learning model. Surprisingly, we find this also leads to improved model interpretability, as measured by identifying genes important for classification using a range of standard interpretability methods. Our results suggest that adversarial training may be generally useful to improve deep learning robustness and interpretability and that it should be evaluated on a range of tasks. Availability and implementation Our Python implementation of all analysis in this publication can be found at: https://github.com/MehrshadSD/robustness-interpretability. The analysis was conducted using numPy 0.2.5, pandas 2.0.3, scanpy 1.9.3, tensorflow 2.10.0, matplotlib 3.7.1, seaborn 0.12.2, sklearn 1.1.1, shap 0.42.0, lime 0.2.0.1, matplotlib_venn 0.11.9. [ABSTRACT FROM AUTHOR]

Motivation K -mer hashing is a common operation in many foundational bioinformatics problems. However, generic string hashing algorithms are not optimized for this application. Strings in bioinformatics use specific alphabets, a trait leveraged for nucleic acid sequences in earlier work. We note that amino acid sequences, with complexities and context that cannot be captured by generic hashing algorithms, can also benefit from a domain-specific hashing algorithm. Such a hashing algorithm can accelerate and improve the sensitivity of bioinformatics applications developed for protein sequences. Results Here, we present aaHash, a recursive hashing algorithm tailored for amino acid sequences. This algorithm utilizes multiple hash levels to represent biochemical similarities between amino acids. aaHash performs ∼10× faster than generic string hashing algorithms in hashing adjacent k -mers. Availability and implementation aaHash is available online at https://github.com/bcgsc/btllib and is free for academic use. [ABSTRACT FROM AUTHOR]

Summary Gene set scoring (or enrichment) is a common dimension reduction task in bioinformatics that can be focused on the differences between groups or at the single sample level. Gene sets can represent biological functions, molecular pathways, cell identities, and more. Gene set scores are context dependent values that are useful for interpreting biological changes following experiments or perturbations. Single sample scoring produces a set of scores, one for each member of a group, which can be analyzed with statistical models that can include additional clinically important factors such as gender or age. However, the sparsity and technical noise of single-cell expression measures create difficulties for these methods, which were originally designed for bulk expression profiling (microarrays, RNAseq). This can be greatly remedied by first applying a smoothing transformation that shares gene measure information within transcriptomic neighborhoods. In this work, we use the nearest neighbor graph of cells for matrix smoothing to produce high quality gene set scores on a per-cell, per-group, level which is useful for visualization and statistical analysis. Availability and implementation The gssnng software is available using the python package index (PyPI) and works with Scanpy AnnData objects. It can be installed using "pip install gssnng." More information and demo notebooks: see https://github.com/IlyaLab/gssnng. [ABSTRACT FROM AUTHOR]

Motivation The fast and accurate detection of similar geometrical arrangements of protein residues, known as 3D structural motifs, is highly relevant for many applications such as binding region and catalytic site detection, drug discovery and structure conservation analyses. With the recent publication of new protein structure prediction methods, the number of available protein structures is exploding, which makes efficient and easy-to-use tools for identifying 3D structural motifs essential. Results We present an open-source Python package that enables the search for both exact and mutated motifs with position-specific residue substitutions. The tool is efficient, flexible, accurate, and suitable to run both on computer clusters and personal laptops. Two successful applications of pyScoMotif for catalytic site identification are showcased. Availability and implementation The pyScoMotif package can be installed from the PyPI repository and is also available at https://github.com/3BioCompBio/pyScoMotif. It is free to use for non-commercial purposes. [ABSTRACT FROM AUTHOR]

Motivation Extended connectivity interaction features (ECIF) is a method developed to predict protein–ligand binding affinity, allowing for detailed atomic representation. It performed very well in terms of Comparative Assessment of Scoring Functions 2016 (CASF-2016) scoring power. However, ECIF has the limitation of not being able to adequately account for interatomic distances. Results To investigate what kind of distance representation is effective for P-L binding affinity prediction, we have developed two algorithms that improved ECIF's feature extraction method to take distance into account. One is multi-shelled ECIF, which takes into account the distance between atoms by dividing the distance between atoms into multiple layers. The other is weighted ECIF, which weights the importance of interactions according to the distance between atoms. A comparison of these two methods shows that multi-shelled ECIF outperforms weighted ECIF and the original ECIF, achieving a CASF-2016 scoring power Pearson correlation coefficient of 0.877. Availability and implementation All the codes and data are available on GitHub (https://github.com/koji11235/MSECIFv2). [ABSTRACT FROM AUTHOR]

Motivation ProbioMinServer is a platform designed to help researchers access information on probiotics regarding a wide variety of characteristics, such as safety (e.g. antimicrobial resistance, virulence, pathogenic, plasmid, and prophage genes) and functionality (e.g. functional classes, carbohydrate-active enzyme, and metabolite gene cluster profile). Because probiotics are functional foods, their safety and functionality are a crucial part of health care. Genomics has become a crucial methodology for investigating the safety and functionality of probiotics in food and feed. This shift is primarily attributed to the growing affordability of next-generation sequencing technologies. However, no integrated platform is available for simultaneously evaluating probiotic strain safety, investigating probiotic functionality, and identifying known phylogenetically related strains. Results Thus, we constructed a new platform, ProbioMinServer, which incorporates these functions. ProbioMinServer accepts whole-genome sequence files in the FASTA format. If the query genome belongs to the 25 common probiotic species collected in our database, the server performs a database search and analyzes the core-genome multilocus sequence typing. Front-end applications were implemented in JavaScript with a bootstrap framework, and back-end programs were implemented using PHP, Perl, and Python. ProbioMinServer can help researchers quickly and easily retrieve information on the safety and functionality of various probiotics. Availability and implementation The platform is available at https://probiomindb.imst.nsysu.edu.tw. [ABSTRACT FROM AUTHOR]

Motivation Three-dimensional chromatin structure plays an important role in gene regulation by connecting regulatory regions and gene promoters. The ability to detect the formation and loss of these loops in various cell types and conditions provides valuable information on the mechanisms driving these cell states and is critical for understanding long-range gene regulation. Hi-C is a powerful technique for characterizing 3D chromatin structure; however, Hi-C can quickly become costly and labor-intensive, and proper planning is required to ensure efficient use of time and resources while maintaining experimental rigor and well-powered results. Results To facilitate better planning and interpretation of human Hi-C experiments, we conducted a detailed evaluation of statistical power using publicly available Hi-C datasets, paying particular attention to the impact of loop size on Hi-C contacts and fold change compression. In addition, we have developed Hi-C Poweraid, a publicly hosted web application to investigate these findings. For experiments involving well-replicated cell lines, we recommend a total sequencing depth of at least 6 billion contacts per condition, split between at least two replicates to achieve the power to detect differences in the majority of loops. For experiments with higher variation, more replicates and deeper sequencing depths are required. Values for specific cases can be determined by using Hi-C Poweraid. This tool simplifies Hi-C power calculations, allowing for more efficient use of time and resources and more accurate interpretation of experimental results. Availability and implementation Hi-C Poweraid is available as an R Shiny application deployed at http://phanstiel-lab.med.unc.edu/poweraid/ , with code available at https://github.com/sarmapar/poweraid. [ABSTRACT FROM AUTHOR]

Motivation Positive-unlabeled data consists of points with either positive or unknown labels. It is widespread in medical, genetic, and biological settings, creating a high demand for predictive positive-unlabeled models. The performance of such models is usually estimated using validation sets, assumed to be selected completely at random (SCAR) from known positive examples. For certain metrics, this assumption enables unbiased performance estimation when treating positive-unlabeled data as positive/negative. However, the SCAR assumption is often adopted without proper justifications, simply for the sake of convenience. Results We provide an algorithm that under the weak assumptions of a lower bound on the number of positive examples can test for the violation of the SCAR assumption. Applying it to the problem of gene prioritization for complex genetic traits, we illustrate that the SCAR assumption is often violated there, causing the inflation of performance estimates, which we refer to as validation bias. We estimate the potential impact of validation bias on performance estimation. Our analysis reveals that validation bias is widespread in gene prioritization data and can significantly overestimate the performance of models. This finding elucidates the discrepancy between the reported good performance of models and their limited practical applications. Availability and implementation Python code with examples of application of the validation bias detection algorithm is available at github.com/ArtomovLab/ValidationBias. [ABSTRACT FROM AUTHOR]

Summary Maximum likelihood (ML) is a widely used phylogenetic inference method. ML implementations heavily rely on numerical optimization routines that use internal numerical thresholds to determine convergence. We systematically analyze the impact of these threshold settings on the log-likelihood and runtimes for ML tree inferences with RAxML-NG, IQ-TREE, and FastTree on empirical datasets. We provide empirical evidence that we can substantially accelerate tree inferences with RAxML-NG and IQ-TREE by changing the default values of two such numerical thresholds. At the same time, altering these settings does not significantly impact the quality of the inferred trees. We further show that increasing both thresholds accelerates the RAxML-NG bootstrap without influencing the resulting support values. For RAxML-NG, increasing the likelihood thresholds ϵ LnL and ϵ brlen to 10 and 103, respectively, results in an average tree inference speedup of 1.9 ± 0.6 on Data collection 1 , 1.8 ± 1.1 on Data collection 2 , and 1.9 ± 0.8 on Data collection 2 for the RAxML-NG bootstrap compared to the runtime under the current default setting. Increasing the likelihood threshold ϵ LnL to 10 in IQ-TREE results in an average tree inference speedup of 1.3 ± 0.4 on Data collection 1 and 1.3 ± 0.9 on Data collection 2. Availability and implementation All MSAs we used for our analyses, as well as all results, are available for download at https://cme.h-its.org/exelixis/material/freeLunch%5fdata.tar.gz. Our data generation scripts are available at https://github.com/tschuelia/ml-numerical-analysis. [ABSTRACT FROM AUTHOR]

Summary We present bubbleHeatmap, an R plotting package which combines elements of a bubble plot and heatmap to conveniently display two numerical variables for each data point across a categorical two dimensional grid. This has particular advantages for visualizing the 251 metabolomic measures produced by the automated, high-throughput, 1H-NMR-based platform provided by Nightingale Health, which includes 12 measures repeated across each of 14 lipoprotein subclasses. As these metabolomic profiles are currently available for large biobanks, we provide a figure template to aid the use of bubbleHeatmap in displaying results from analyses using these data. Availability and implementation https://cran.r-project.org/web/packages/bubbleHeatmap [ABSTRACT FROM AUTHOR]