123VCF: an intuitive and efficient tool for filtering VCF files.

BACKGROUND: The advent of Next-Generation Sequencing (NGS) has catalyzed a paradigm shift in medical genetics, enabling the identification of disease-associated variants. However, the vast quantum of data produced by NGS necessitates a robust and dependable mechanism for filtering irrelevant variants. Annotation-based variant filtering, a pivotal step in this process, demands a profound understanding of the case-specific conditions and the relevant annotation instruments. To tackle this complex task, we sought to design an accessible, efficient and more importantly easy to understand variant filtering tool. RESULTS: Our efforts culminated in the creation of 123VCF, a tool capable of processing both compressed and uncompressed Variant Calling Format (VCF) files. Built on a Java framework, the tool employs a disk-streaming real-time filtering algorithm, allowing it to manage sizable variant files on conventional desktop computers. 123VCF filters input variants in accordance with a predefined filter sequence applied to the input variants. Users are provided the flexibility to define various filtering parameters, such as quality, coverage depth, and variant frequency within the populations. Additionally, 123VCF accommodates user-defined filters tailored to specific case requirements, affording users enhanced control over the filtering process. We evaluated the performance of 123VCF by analyzing different types of variant files and comparing its runtimes to the most similar algorithms like BCFtools filter and GATK VariantFiltration. The results indicated that 123VCF performs relatively well. The tool's intuitive interface and potential for reproducibility make it a valuable asset for both researchers and clinicians. CONCLUSION: The 123VCF filtering tool provides an effective, dependable approach for filtering variants in both research and clinical settings. As an open-source tool available at https://project123vcf.sourceforge.io , it is accessible to the global scientific and clinical community, paving the way for the discovery of disease-causing variants and facilitating the advancement of personalized medicine.

AI-Enhanced Detection of Clinically Relevant Structural and Functional Anomalies in MRI: Traversing the Landscape of Conventional to Explainable Approaches.

Anomaly detection in medical imaging, particularly within the realm of magnetic resonance imaging (MRI), stands as a vital area of research with far-reaching implications across various medical fields. This review meticulously examines the integration of artificial intelligence (AI) in anomaly detection for MR images, spotlighting its transformative impact on medical diagnostics. We delve into the forefront of AI applications in MRI, exploring advanced machine learning (ML) and deep learning (DL) methodologies that are pivotal in enhancing the precision of diagnostic processes. The review provides a detailed analysis of preprocessing, feature extraction, classification, and segmentation techniques, alongside a comprehensive evaluation of commonly used metrics. Further, this paper explores the latest developments in ensemble methods and explainable AI, offering insights into future directions and potential breakthroughs. This review synthesizes current insights, offering a valuable guide for researchers, clinicians, and medical imaging experts. It highlights AI's crucial role in improving the precision and speed of detecting key structural and functional irregularities in MRI. Our exploration of innovative techniques and trends furthers MRI technology development, aiming to refine diagnostics, tailor treatments, and elevate patient care outcomes. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

A predictive ensemble classifier for the gene expression diagnosis of ASD at ages 1 to 4 years.

Autism Spectrum Disorder (ASD) diagnosis remains behavior-based and the median age of diagnosis is ~52 months, nearly 5 years after its first-trimester origin. Accurate and clinically-translatable early-age diagnostics do not exist due to ASD genetic and clinical heterogeneity. Here we collected clinical, diagnostic, and leukocyte RNA data from 240 ASD and typically developing (TD) toddlers (175 toddlers for training and 65 for test). To identify gene expression ASD diagnostic classifiers, we developed 42,840 models composed of 3570 gene expression feature selection sets and 12 classification methods. We found that 742 models had AUC-ROC ≥ 0.8 on both Training and Test sets. Weighted Bayesian model averaging of these 742 models yielded an ensemble classifier model with accurate performance in Training and Test gene expression datasets with ASD diagnostic classification AUC-ROC scores of 85-89% and AUC-PR scores of 84-92%. ASD toddlers with ensemble scores above and below the overall ASD ensemble mean of 0.723 (on a scale of 0 to 1) had similar diagnostic and psychometric scores, but those below this ASD ensemble mean had more prenatal risk events than TD toddlers. Ensemble model feature genes were involved in cell cycle, inflammation/immune response, transcriptional gene regulation, cytokine response, and PI3K-AKT, RAS and Wnt signaling pathways. We additionally collected targeted DNA sequencing smMIPs data on a subset of ASD risk genes from 217 of the 240 ASD and TD toddlers. This DNA sequencing found about the same percentage of SFARI Level 1 and 2 ASD risk gene mutations in TD (12 of 105) as in ASD (13 of 112) toddlers, and classification based only on the presence of mutation in these risk genes performed at a chance level of 49%. By contrast, the leukocyte ensemble gene expression classifier correctly diagnostically classified 88% of TD and ASD toddlers with ASD risk gene mutations. Our ensemble ASD gene expression classifier is diagnostically predictive and replicable across different toddler ages, races, and ethnicities; out-performs a risk gene mutation classifier; and has potential for clinical translation.

Transcriptomic Analysis of Cyclamen persicum to Identify Invovled Genes in Triterpene Secondary Metabolites Pathway.

Saponins are considered as a diverse group of natural active compounds, which are widely found in crops. Mevalonate pathway (MVA) is regarded as the main pathway for synthesis of saponins in crops. This study aims to compare transcriptome of the leaf with tuber of crop including tubers and roots. First, more than 166 million reads were generated. The existence of 36,678 unigenes in the two samples out of 48,936 assembled ones showed a significant difference in expression. Finally, 310 and 290 highly up-regulated genes in leaf and tuber were selected for the next analysis. In addition, the expression profiles of 13 key genes in the MVA pathway were compared in RNA sequencing and reverse transcription-quantitative polymerase chain reaction analysis. The results indicated that cyclamen tuber has a higher level of expression of MVA pathway genes. The topological analysis for gene co-expression network involved in triterpenoid synthesis represented that the genes at the beginning of such pathway play a critical role so that the reduction of their expression challenges triterpenoid synthesis severely. The tuber of the cyclamen appears to be the major site of triterpene synthesis, and transfer of excess Isopentenyl pyrophosphate (IPP) from tuber to leaf activates downstream genes in leaf of crop.

Reconstructing the cell-cell interaction network among mouse immune cells.

Intercellular interactions and cell-cell communication are critical to regulating cell functions, especially in normal immune cells and immunotherapies. Ligand-receptor pairs mediating these cell-cell interactions can be identified using diverse experimental and computational approaches. Here, we reconstructed the intercellular interaction network between Mus musculus immune cells using publicly available receptor-ligand interaction databases and gene expression data from the immunological genome project. This reconstructed network accounts for 50,317 unique interactions between 16 cell types between 731 receptor-ligand pairs. Analysis of this network shows that cells of hematopoietic lineages use fewer communication pathways for interacting with each other, while nonhematopoietic stromal cells use the most network communications. We further observe that the WNT, BMP, and LAMININ pathways are the most significant contributors to the overall number of cell-cell interactions among the various pathways in the reconstructed communication network. This resource will enable the systematic analysis of normal and pathologic immune cell interactions, along with the study of emerging immunotherapies.

Bioinformatic tools for DNA methylation and histone modification: A survey.

Epigenetic inheritance occurs due to different mechanisms such as chromatin and histone modifications, DNA methylation and processes mediated by non-coding RNAs. It leads to changes in gene expressions and the emergence of new traits in different organisms in many diseases such as cancer. Recent advances in experimental methods led to the identification of epigenetic target sites in various organisms. Computational approaches have enabled us to analyze mass data produced by these methods. Next-generation sequencing (NGS) methods have been broadly used to identify these target sites and their patterns. By using these patterns, the emergence of diseases could be prognosticated. In this study, target site prediction tools for two major epigenetic mechanisms comprising histone modification and DNA methylation are reviewed. Publicly accessible databases are reviewed as well. Some suggestions regarding the state-of-the-art methods and databases have been made, including examining patterns of epigenetic changes that are important in epigenotypes detection.

Moonlighting protein prediction using physico-chemical and evolutional properties via machine learning methods.

BACKGROUND: Moonlighting proteins (MPs) are a subclass of multifunctional proteins in which more than one independent or usually distinct function occurs in a single polypeptide chain. Identification of unknown cellular processes, understanding novel protein mechanisms, improving the prediction of protein functions, and gaining information about protein evolution are the main reasons to study MPs. They also play an important role in disease pathways and drug-target discovery. Since detecting MPs experimentally is quite a challenge, most of them are detected randomly. Therefore, introducing an appropriate computational approach to predict MPs seems reasonable. RESULTS: In this study, we introduced a competent model for detecting moonlighting and non-MPs through extracted features from protein sequences. We attempted to set up a well-judged scheme for detecting outlier proteins. Consequently, 37 distinct feature vectors were utilized to study each protein's impact on detecting MPs. Furthermore, 8 different classification methods were assessed to find the best performance. To detect outliers, each one of the classifications was executed 100 times by tenfold cross-validation on feature vectors; proteins which misclassified 90 times or more were grouped. This process was applied to every single feature vector and eventually the intersection of these groups was determined as the outlier proteins. The results of tenfold cross-validation on a dataset of 351 samples (containing 215 moonlighting and 136 non-moonlighting proteins) reveal that the SVM method on all feature vectors has the highest performance among all methods in this study and other available methods. Besides, the study of outliers showed that 57 of 351 proteins in the dataset could be an appropriate candidate for the outlier. Among the outlier proteins, there were non-MPs (such as P69797) that have been misclassified in 8 different classification methods with 16 different feature vectors. Because these proteins have been obtained by computational methods, the results of this study could reduce the likelihood of hypothesizing whether these proteins are non-moonlighting at all. CONCLUSIONS: MPs are difficult to be identified through experimentation. Using distinct feature vectors, our method enabled identification of novel moonlighting proteins. The study also pinpointed that a number of non-MPs are likely to be moonlighting.

Drug repositioning based on gene expression data for human HER2-positive breast cancer.

Human epidermal growth factor receptor 2 (HER2)-positive breast cancer represents approximately 15-30% of all invasive breast cancers. Despite the recent advances in therapeutic practices of HER2 subtype, drug resistance and tumor recurrence still have remained as major problems. Drug discovery is a long and difficult process, so the aim of this study is to find potential new application for existing therapeutic agents. Gene expression data for breast invasive carcinoma were retrieved from The Cancer Genome Atlas (TCGA) database. The normal and tumor samples were analyzed using Linear Models for Microarray Data (LIMMA) R package in order to find the differentially expressed genes (DEGs). These genes were used as entry for the library of integrated network-based cellular signatures (LINCS) L1000CDS2 software and suggested 24 repurposed drugs. According to the obtained results, some of these drugs including vorinostat, mocetinostat, alvocidib, CGP-60474, BMS-387032, AT-7519, and curcumin have significant functional similarity and structural correlation with FDA-approved breast cancer drugs. Based on the drug-target network, which consisted of the repurposed drugs and their target genes, the aforementioned drugs had the highest degrees. Moreover, the experimental approach verified curcumin as an effective therapeutic agent for HER2 positive breast cancer. Hence, our work suggested that some repurposed drugs based on gene expression data can be noticed as potential drugs for the treatment of HER2-positive breast cancer.

Gibberellin causes wide transcriptional modifications in the early stage of grape cluster development.

Yaghooti grape of Sistan is seedless, early ripening but has compact clusters. To study gibberellin effect on cluster compactness of Yaghooti grape, it has been studied transcriptomic changes in three developmental stages (cluster formation, berry formation and final size of cluster). We found out that 5409 of 22,756 genes in cluster tissue showed significant changes under gibberellin. Finally, it was showed that 2855, 2862 and 497 genes have critically important role on above developmental stages, respectively. GO enrichment analysis showed that gibberellin enhances biochemical pathways activity. Moreover, genes involved in ribosomal structure and photosynthesis rate in cluster tissue were up- and down- regulated, respectively. In addition, we observed location of metabolomic activities was transferred from nucleus to cytoplasm and from cytoplasm to cell wall and intercellular spaces during cluster development; but there is not such situation in gibberellin treated samples.

Protein complex prediction: A survey.

Protein complexes are one of the most important functional units for deriving biological processes within the cell. Experimental methods have provided valuable data to infer protein complexes. However, these methods have inherent limitations. Considering these limitations, many computational methods have been proposed to predict protein complexes, in the last decade. Almost all of these in-silico methods predict protein complexes from the ever-increasing protein-protein interaction (PPI) data. These computational approaches usually use the PPI data in the format of a huge protein-protein interaction network (PPIN) as input and output various sub-networks of the given PPIN as the predicted protein complexes. Some of these methods have already reached a promising efficiency in protein complex detection. Nonetheless, there are challenges in prediction of other types of protein complexes, specially sparse and small ones. New methods should further incorporate the knowledge of biological properties of proteins to improve the performance. Additionally, there are several challenges that should be considered more effectively in designing the new complex prediction algorithms in the future. This article not only reviews the history of computational protein complex prediction but also provides new insight for improvement of new methodologies. In this article, most important computational methods for protein complex prediction are evaluated and compared. In addition, some of the challenges in the reconstruction of the protein complexes are discussed. Finally, various tools for protein complex prediction and PPIN analysis as well as the current high-throughput databases are reviewed.

RepCOOL: computational drug repositioning via integrating heterogeneous biological networks.

BACKGROUND: It often takes more than 10 years and costs more than 1 billion dollars to develop a new drug for a particular disease and bring it to the market. Drug repositioning can significantly reduce costs and time in drug development. Recently, computational drug repositioning attracted a considerable amount of attention among researchers, and a plethora of computational drug repositioning methods have been proposed. This methodology has widely been used in order to address various medical challenges, including cancer treatment. The most common cancers are lung and breast cancers. Thus, suggesting FDA-approved drugs via drug repositioning for breast cancer would help us to circumvent the approval process and subsequently save money as well as time. METHODS: In this study, we propose a novel network-based method, named RepCOOL, for drug repositioning. RepCOOL integrates various heterogeneous biological networks to suggest new drug candidates for a given disease. RESULTS: The proposed method showed a promising performance on benchmark datasets via rigorous cross-validation. The final drug repositioning model has been built based on a random forest classifier after examining various machine learning algorithms. Finally, in a case study, four FDA approved drugs were suggested for breast cancer stage II. CONCLUSION: Results show the potency of the proposed method in detecting true drug-disease relationships. RepCOOL suggested four new drugs for breast cancer stage II namely Doxorubicin, Paclitaxel, Trastuzumab, and Tamoxifen.

Performance evaluation measures for protein complex prediction.

Protein complexes play a dominant role in cellular organization and function. Prediction of protein complexes from the network of physical interactions between proteins (PPI networks) has thus become one of the important research areas. Recently, many computational approaches have been developed to identify these complexes. Various performance assessment measures have been proposed for evaluating the efficiency of these methods. However, there are many inconsistencies in the definitions and usage of the measures across the literature. To address this issue, we have gathered and presented the most important performance evaluation measures and developed a tool, named CompEvaluator, to critically assess the protein complex prediction methods. The tool and documentation are publicly available at https://sourceforge.net/projects/compevaluator/files/.

AntAngioCOOL: computational detection of anti-angiogenic peptides.

BACKGROUND: Angiogenesis inhibition research is a cutting edge area in angiogenesis-dependent disease therapy, especially in cancer therapy. Recently, studies on anti-angiogenic peptides have provided promising results in the field of cancer treatment. METHODS: A non-redundant dataset of 135 anti-angiogenic peptides (positive instances) and 135 non anti-angiogenic peptides (negative instances) was used in this study. Also, 20% of each class were selected to construct an independent test dataset (see Additional files 1, 2). We proposed an effective machine learning based R package (AntAngioCOOL) to predict anti-angiogenic peptides. We have examined more than 200 different classifiers to build an efficient predictor. Also, more than 17,000 features were extracted to encode the peptides. RESULTS: Finally, more than 2000 informative features were selected to train the classifiers for detecting anti-angiogenic peptides. AntAngioCOOL includes three different models that can be selected by the user for different purposes; it is the most sensitive, most specific and most accurate. According to the obtained results AntAngioCOOL can effectively suggest anti-angiogenic peptides; this tool achieved sensitivity of 88%, specificity of 77% and accuracy of 75% on the independent test set. AntAngioCOOL can be accessed at https://cran.r-project.org/ . CONCLUSIONS: Only 2% of the extracted descriptors were used to build the predictor models. The results revealed that physico-chemical profile is the most important feature type in predicting anti-angiogenic peptides. Also, atomic profile and PseAAC are the other important features.

rpiCOOL: A tool for In Silico RNA-protein interaction detection using random forest.

Understanding the principle of RNA-protein interactions (RPIs) is of critical importance to provide insights into post-transcriptional gene regulation and is useful to guide studies about many complex diseases. The limitations and difficulties associated with experimental determination of RPIs, call an urgent need to computational methods for RPI prediction. In this paper, we proposed a machine learning method to detect RNA-protein interactions based on sequence information. We used motif information and repetitive patterns, which have been extracted from experimentally validated RNA-protein interactions, in combination with sequence composition as descriptors to build a model to RPI prediction via a random forest classifier. About 20% of the "sequence motifs" and "nucleotide composition" features have been selected as the informative features with the feature selection methods. These results suggest that these two feature types contribute effectively in RPI detection. Results of 10-fold cross-validation experiments on three non-redundant benchmark datasets show a better performance of the proposed method in comparison with the current state-of-the-art methods in terms of various performance measures. In addition, the results revealed that the accuracy of the RPI prediction methods could vary considerably across different organisms. We have implemented the proposed method, namely rpiCOOL, as a stand-alone tool with a user friendly graphical user interface (GUI) that enables the researchers to predict RNA-protein interaction. The rpiCOOL is freely available at http://biocool.ir/rpicool.html for non-commercial uses.

LocFuse: human protein-protein interaction prediction via classifier fusion using protein localization information.

Protein-protein interaction (PPI) detection is one of the central goals of functional genomics and systems biology. Knowledge about the nature of PPIs can help fill the widening gap between sequence information and functional annotations. Although experimental methods have produced valuable PPI data, they also suffer from significant limitations. Computational PPI prediction methods have attracted tremendous attentions. Despite considerable efforts, PPI prediction is still in its infancy in complex multicellular organisms such as humans. Here, we propose a novel ensemble learning method, LocFuse, which is useful in human PPI prediction. This method uses eight different genomic and proteomic features along with four types of different classifiers. The prediction performance of this classifier selection method was found to be considerably better than methods employed hitherto. This confirms the complex nature of the PPI prediction problem and also the necessity of using biological information for classifier fusion. The LocFuse is available at: http://lbb.ut.ac.ir/Download/LBBsoft/LocFuse. BIOLOGICAL SIGNIFICANCE: The results revealed that if we divide proteome space according to the cellular localization of proteins, then the utility of some classifiers in PPI prediction can be improved. Therefore, to predict the interaction for any given protein pair, we can select the most accurate classifier with regard to the cellular localization information. Based on the results, we can say that the importance of different features for PPI prediction varies between differently localized proteins; however in general, our novel features, which were extracted from position-specific scoring matrices (PSSMs), are the most important ones and the Random Forest (RF) classifier performs best in most cases. LocFuse was developed with a user-friendly graphic interface and it is freely available for Linux, Mac OSX and MS Windows operating systems.

PPIevo: protein-protein interaction prediction from PSSM based evolutionary information.

Protein-protein interactions regulate a variety of cellular processes. There is a great need for computational methods as a complement to experimental methods with which to predict protein interactions due to the existence of many limitations involved in experimental techniques. Here, we introduce a novel evolutionary based feature extraction algorithm for protein-protein interaction (PPI) prediction. The algorithm is called PPIevo and extracts the evolutionary feature from Position-Specific Scoring Matrix (PSSM) of protein with known sequence. The algorithm does not depend on the protein annotations, and the features are based on the evolutionary history of the proteins. This enables the algorithm to have more power for predicting protein-protein interaction than many sequence based algorithms. Results on the HPRD database show better performance and robustness of the proposed method. They also reveal that the negative dataset selection could lead to an acute performance overestimation which is the principal drawback of the available methods.

Predictive three-biomarker panel in peripheral blood mononuclear cells for detecting hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) ranks among the most prevalent cancers and accounts for a significant proportion of cancer-associated deaths worldwide. This disease, marked by multifaceted etiology, often poses diagnostic challenges. Finding a reliable and non-invasive diagnostic method seems to be necessary. In this study, we analyzed the gene expression profiles of 20 HCC patients, 12 individuals with chronic hepatitis, and 15 healthy controls. Enrichment analysis revealed that platelet aggregation, secretory granule lumen, and G-protein-coupled purinergic nucleotide receptor activity were common biological processes, cellular components, and molecular function in HCC and chronic hepatitis B (CHB) compared to healthy controls, respectively. Furthermore, pathway analysis demonstrated that "estrogen response" was involved in the pathogenesis of HCC and CHB conditions, while, "apoptosis" and "coagulation" pathways were specific for HCC. Employing computational feature selection and logistic regression classification, we identified candidate genes pivotal for diagnostic panel development and evaluated the performance of these panels. Subsequent machine learning evaluations assessed these panels' performance in an independent cohort. Remarkably, a 3-marker panel, comprising RANSE2, TNF-α, and MAP3K7, demonstrated the best performance in qRT-PCR-validated experimental data, achieving 98.4% accuracy and an area under the curve of 1. Our findings highlight this panel's promising potential as a non-invasive approach not only for detecting HCC but also for distinguishing HCC from CHB patients.

Differences in regional brain structure in toddlers with autism are related to future language outcomes.

Language and social symptoms improve with age in some autistic toddlers, but not in others, and such outcome differences are not clearly predictable from clinical scores alone. Here we aim to identify early-age brain alterations in autism that are prognostic of future language ability. Leveraging 372 longitudinal structural MRI scans from 166 autistic toddlers and 109 typical toddlers and controlling for brain size, we find that, compared to typical toddlers, autistic toddlers show differentially larger or thicker temporal and fusiform regions; smaller or thinner inferior frontal lobe and midline structures; larger callosal subregion volume; and smaller cerebellum. Most differences are replicated in an independent cohort of 75 toddlers. These brain alterations improve accuracy for predicting language outcome at 6-month follow-up beyond intake clinical and demographic variables. Temporal, fusiform, and inferior frontal alterations are related to autism symptom severity and cognitive impairments at early intake ages. Among autistic toddlers, brain alterations in social, language and face processing areas enhance the prediction of the child's future language ability.

Computational Prediction of Protein-Protein Interaction Networks: Algo-rithms and Resources.

Protein interactions play an important role in the discovery of protein functions and pathways in biological processes. This is especially true in case of the diseases caused by the loss of specific protein-protein interactions in the organism. The accuracy of experimental results in finding protein-protein interactions, however, is rather dubious and high throughput experimental results have shown both high false positive beside false negative information for protein interaction. Computational methods have attracted tremendous attention among biologists because of the ability to predict protein-protein interactions and validate the obtained experimental results. In this study, we have reviewed several computational methods for protein-protein interaction prediction as well as describing major databases, which store both predicted and detected protein-protein interactions, and the tools used for analyzing protein interaction networks and improving protein-protein interaction reliability.