BiPSTP: Sequence feature encoding method for identifying different RNA modifications with bidirectional position-specific trinucleotides propensities.

RNA modification, a posttranscriptional regulatory mechanism, significantly influences RNA biogenesis and function. The accurate identification of modification sites is paramount for investigating their biological implications. Methods for encoding RNA sequence into numerical data play a crucial role in developing robust models for predicting modification sites. However, existing techniques suffer from limitations, including inadequate information representation, challenges in effectively integrating positional and sequential information, and the generation of irrelevant or redundant features when combining multiple approaches. These deficiencies hinder the effectiveness of machine learning models in addressing the performance challenges associated with predicting RNA modification sites. Here, we introduce a novel RNA sequence feature representation method, named BiPSTP, which utilizes bidirectional trinucleotide position-specific propensities. We employ the parameter ξ to denote the interval between the current nucleotide and its adjacent forward or backward dinucleotide, enabling the extraction of positional and sequential information from RNA sequences. Leveraging the BiPSTP method, we have developed the prediction model mRNAPred using support vector machine classifier to identify multiple types of RNA modification sites. We evaluate the performance of our BiPSTP method and mRNAPred model across 12 distinct RNA modification types. Our experimental results demonstrate the superiority of the mRNAPred model compared to state-of-art models in the domain of RNA modification sites identification. Importantly, our BiPSTP method enhances the robustness and generalization performance of prediction models. Notably, it can be applied to feature extraction from DNA sequences to predict other biological modification sites.

Functional support vector machine.

Linear and generalized linear scalar-on-function modeling have been commonly used to understand the relationship between a scalar response variable (e.g. continuous, binary outcomes) and functional predictors. Such techniques are sensitive to model misspecification when the relationship between the response variable and the functional predictors is complex. On the other hand, support vector machines (SVMs) are among the most robust prediction models but do not take account of the high correlations between repeated measurements and cannot be used for irregular data. In this work, we propose a novel method to integrate functional principal component analysis with SVM techniques for classification and regression to account for the continuous nature of functional data and the nonlinear relationship between the scalar response variable and the functional predictors. We demonstrate the performance of our method through extensive simulation experiments and two real data applications: the classification of alcoholics using electroencephalography signals and the prediction of glucobrassicin concentration using near-infrared reflectance spectroscopy. Our methods especially have more advantages when the measurement errors in functional predictors are relatively large.

Identification of species-specific RNA N6-methyladinosine modification sites from RNA sequences.

N6-methyladinosine (m6A) modification is the most abundant co-transcriptional modification in eukaryotic RNA and plays important roles in cellular regulation. Traditional high-throughput sequencing experiments used to explore functional mechanisms are time-consuming and labor-intensive, and most of the proposed methods focused on limited species types. To further understand the relevant biological mechanisms among different species with the same RNA modification, it is necessary to develop a computational scheme that can be applied to different species. To achieve this, we proposed an attention-based deep learning method, adaptive-m6A, which consists of convolutional neural network, bi-directional long short-term memory and an attention mechanism, to identify m6A sites in multiple species. In addition, three conventional machine learning (ML) methods, including support vector machine, random forest and logistic regression classifiers, were considered in this work. In addition to the performance of ML methods for multi-species prediction, the optimal performance of adaptive-m6A yielded an accuracy of 0.9832 and the area under the receiver operating characteristic curve of 0.98. Moreover, the motif analysis and cross-validation among different species were conducted to test the robustness of one model towards multiple species, which helped improve our understanding about the sequence characteristics and biological functions of RNA modifications in different species.

PROSE: phenotype-specific network signatures from individual proteomic samples.

Proteomic studies characterize the protein composition of complex biological samples. Despite recent advancements in mass spectrometry instrumentation and computational tools, low proteome coverage and interpretability remains a challenge. To address this, we developed Proteome Support Vector Enrichment (PROSE), a fast, scalable and lightweight pipeline for scoring proteins based on orthogonal gene co-expression network matrices. PROSE utilizes simple protein lists as input, generating a standard enrichment score for all proteins, including undetected ones. In our benchmark with 7 other candidate prioritization techniques, PROSE shows high accuracy in missing protein prediction, with scores correlating strongly to corresponding gene expression data. As a further proof-of-concept, we applied PROSE to a reanalysis of the Cancer Cell Line Encyclopedia proteomics dataset, where it captures key phenotypic features, including gene dependency. We lastly demonstrated its applicability on a breast cancer clinical dataset, showing clustering by annotated molecular subtype and identification of putative drivers of triple-negative breast cancer. PROSE is available as a user-friendly Python module from https://github.com/bwbio/PROSE.

Sequence-based prediction model of protein crystallization propensity using machine learning and two-level feature selection.

Protein crystallization is crucial for biology, but the steps involved are complex and demanding in terms of external factors and internal structure. To save on experimental costs and time, the tendency of proteins to crystallize can be initially determined and screened by modeling. As a result, this study created a new pipeline aimed at using protein sequence to predict protein crystallization propensity in the protein material production stage, purification stage and production of crystal stage. The newly created pipeline proposed a new feature selection method, which involves combining Chi-square (${\chi }^{2}$) and recursive feature elimination together with the 12 selected features, followed by a linear discriminant analysisfor dimensionality reduction and finally, a support vector machine algorithm with hyperparameter tuning and 10-fold cross-validation is used to train the model and test the results. This new pipeline has been tested on three different datasets, and the accuracy rates are higher than the existing pipelines. In conclusion, our model provides a new solution to predict multistage protein crystallization propensity which is a big challenge in computational biology.

methylClass: an R package to construct DNA methylation-based classification models.

DNA methylation profiling is a useful tool to increase the accuracy of a cancer diagnosis. However, a comprehensive R package specially for it is lacking. Hence, we developed the R package methylClass for methylation-based classification. Within it, we provide the eSVM (ensemble-based support vector machine) model to achieve much higher accuracy in methylation data classification than the popular random forest model and overcome the time-consuming problem of the traditional SVM. In addition, some novel feature selection methods are included in the package to improve the classification. Furthermore, because methylation data can be converted to other omics, such as copy number variation data, we also provide functions for multi-omics studies. The testing of this package on four datasets shows the accurate performance of our package, especially eSVM, which can be used in both methylation and multi-omics models and outperforms other methods in both cases. methylClass is available at: https://github.com/yuabrahamliu/methylClass.

DBPboost:A method of classification of DNA-binding proteins based on improved differential evolution algorithm and feature extraction.

DNA-binding proteins are a class of proteins that can interact with DNA molecules through physical and chemical interactions. Their main functions include regulating gene expression, maintaining chromosome structure and stability, and more. DNA-binding proteins play a crucial role in cellular and molecular biology, as they are essential for maintaining normal cellular physiological functions and adapting to environmental changes. The prediction of DNA-binding proteins has been a hot topic in the field of bioinformatics. The key to accurately classifying DNA-binding proteins is to find suitable feature sources and explore the information they contain. Although there are already many models for predicting DNA-binding proteins, there is still room for improvement in mining feature source information and calculation methods. In this study, we created a model called DBPboost to better identify DNA-binding proteins. The innovation of this study lies in the use of eight feature extraction methods, the improvement of the feature selection step, which involves selecting some features first and then performing feature selection again after feature fusion, and the optimization of the differential evolution algorithm in feature fusion, which improves the performance of feature fusion. The experimental results show that the prediction accuracy of the model on the UniSwiss dataset is 89.32%, and the sensitivity is 89.01%, which is better than most existing models.

Shared and malignancy-specific functional plasticity of dynamic brain properties for patients with left frontal glioma.

The time-varying brain activity may parallel the disease progression of cerebral glioma. Assessment of brain dynamics would better characterize the pathological profile of glioma and the relevant functional remodeling. This study aims to investigate the dynamic properties of functional networks based on sliding-window approach for patients with left frontal glioma. The generalized functional plasticity due to glioma was characterized by reduced dynamic amplitude of low-frequency fluctuation of somatosensory networks, reduced dynamic functional connectivity between homotopic regions mainly involving dorsal attention network and subcortical nuclei, and enhanced subcortical dynamic functional connectivity. Malignancy-specific functional remodeling featured a chaotic modification of dynamic amplitude of low-frequency fluctuation and dynamic functional connectivity for low-grade gliomas, and attenuated dynamic functional connectivity of the intrahemispheric cortico-subcortical connections and reduced dynamic amplitude of low-frequency fluctuation of the bilateral caudate for high-grade gliomas. Network dynamic activity was clustered into four distinct configuration states. The occurrence and dwell time of the weakly connected state were reduced in patients' brains. Support vector machine model combined with predictive dynamic features achieved an averaged accuracy of 87.9% in distinguishing low- and high-grade gliomas. In conclusion, dynamic network properties are highly predictive of the malignant grade of gliomas, thus could serve as new biomarkers for disease characterization.

A novel sand cat swarm optimization algorithm-based SVM for diagnosis imaging genomics in Alzheimer's disease.

In recent years, brain imaging genomics has advanced significantly in revealing underlying pathological mechanisms of Alzheimer's disease (AD) and providing early diagnosis. In this paper, we present a framework for diagnosing AD that integrates magnetic resonance imaging (fMRI) genetic preprocessing, feature selection, and a support vector machine (SVM) model. In particular, a novel sand cat swarm optimization (SCSO) algorithm, named SS-SCSO, which integrates the spiral search strategy and alert mechanism from the sparrow search algorithm, is proposed to optimize the SVM parameters. The optimization efficacy of the SS-SCSO algorithm is evaluated using CEC2017 benchmark functions, with results compared with other metaheuristic algorithms (MAs). The proposed SS-SCSO-SVM framework has been effectively employed to classify different stages of cognitive impairment in Alzheimer's Disease using imaging genetic datasets from the Alzheimer's Disease Neuroimaging Initiative. It has demonstrated excellent classification accuracies for four typical cases, including AD, early mild cognitive impairment, late mild cognitive impairment, and healthy control. Furthermore, experiment results indicate that the SS-SCSO-SVM algorithm has a stronger exploration capability for diagnosing AD compared to other well-established MAs and machine learning techniques.

Integration of graph network with kernel SVM and logistic regression for identification of biomarkers in SCA12 and its diagnosis.

Spinocerebellar ataxia type 12 is a hereditary and neurodegenerative illness commonly found in India. However, there is no established noninvasive automatic diagnostic system for its diagnosis and identification of imaging biomarkers. This work proposes a novel four-phase machine learning-based diagnostic framework to find spinocerebellar ataxia type 12 disease-specific atrophic-brain regions and distinguish spinocerebellar ataxia type 12 from healthy using a real structural magnetic resonance imaging dataset. Firstly, each brain region is represented in terms of statistics of coefficients obtained using 3D-discrete wavelet transform. Secondly, a set of relevant regions are selected using a graph network-based method. Thirdly, a kernel support vector machine is used to capture nonlinear relationships among the voxels of a brain region. Finally, the linear relationship among the brain regions is captured to build a decision model to distinguish spinocerebellar ataxia type 12 from healthy by using the regularized logistic regression method. A classification accuracy of 95% and a harmonic mean of precision and recall, i.e. F1-score of 94.92%, is achieved. The proposed framework provides relevant regions responsible for the atrophy. The importance of each region is captured using Shapley Additive exPlanations values. We also performed a statistical analysis to find volumetric changes in spinocerebellar ataxia type 12 group compared to healthy. The promising result of the proposed framework shows that clinicians can use it for early and timely diagnosis of spinocerebellar ataxia type 12.

Spatial functional reorganizations can serve as potential biomarkers of post facial palsy synkinesis.

Facial palsy can result in a serious complication known as facial synkinesis, causing both physical and psychological harm to the patients. There is growing evidence that patients with facial synkinesis have brain abnormalities, but the brain mechanisms and underlying imaging biomarkers remain unclear. Here, we employed functional magnetic resonance imaging (fMRI) to investigate brain function in 31 unilateral post facial palsy synkinesis patients and 25 healthy controls during different facial expression movements and at rest. Combining surface-based mass-univariate analysis and multivariate pattern analysis, we identified diffused activation and intrinsic connection patterns in the primary motor cortex and the somatosensory cortex on the patient's affected side. Further, we classified post facial palsy synkinesis patients from healthy subjects with favorable accuracy using the support vector machine based on both task-related and resting-state functional magnetic resonance imaging data. Together, these findings indicate the potential of the identified functional reorganizations to serve as neuroimaging biomarkers for facial synkinesis diagnosis.

Cortical-cerebellar circuits changes in preschool ASD children by multimodal MRI.

OBJECTIVE: To investigate the alterations in cortical-cerebellar circuits and assess their diagnostic potential in preschool children with autism spectrum disorder using multimodal magnetic resonance imaging. METHODS: We utilized diffusion basis spectrum imaging approaches, namely DBSI_20 and DBSI_combine, alongside 3D structural imaging to examine 31 autism spectrum disorder diagnosed patients and 30 healthy controls. The participants' brains were segmented into 120 anatomical regions for this analysis, and a multimodal strategy was adopted to assess the brain networks using a multi-kernel support vector machine for classification. RESULTS: The results revealed consensus connections in the cortical-cerebellar and subcortical-cerebellar circuits, notably in the thalamus and basal ganglia. These connections were predominantly positive in the frontoparietal and subcortical pathways, whereas negative consensus connections were mainly observed in frontotemporal and subcortical pathways. Among the models tested, DBSI_20 showed the highest accuracy rate of 86.88%. In addition, further analysis indicated that combining the 3 models resulted in the most effective performance. CONCLUSION: The connectivity network analysis of the multimodal brain data identified significant abnormalities in the cortical-cerebellar circuits in autism spectrum disorder patients. The DBSI_20 model not only provided the highest accuracy but also demonstrated efficiency, suggesting its potential for clinical application in autism spectrum disorder diagnosis.

m6A-TCPred: a web server to predict tissue-conserved human m6A sites using machine learning approach.

BACKGROUND: N6-methyladenosine (m6A) is the most prevalent post-transcriptional modification in eukaryotic cells that plays a crucial role in regulating various biological processes, and dysregulation of m6A status is involved in multiple human diseases including cancer contexts. A number of prediction frameworks have been proposed for high-accuracy identification of putative m6A sites, however, none have targeted for direct prediction of tissue-conserved m6A modified residues from non-conserved ones at base-resolution level. RESULTS: We report here m6A-TCPred, a computational tool for predicting tissue-conserved m6A residues using m6A profiling data from 23 human tissues. By taking advantage of the traditional sequence-based characteristics and additional genome-derived information, m6A-TCPred successfully captured distinct patterns between potentially tissue-conserved m6A modifications and non-conserved ones, with an average AUROC of 0.871 and 0.879 tested on cross-validation and independent datasets, respectively. CONCLUSION: Our results have been integrated into an online platform: a database holding 268,115 high confidence m6A sites with their conserved information across 23 human tissues; and a web server to predict the conserved status of user-provided m6A collections. The web interface of m6A-TCPred is freely accessible at: www.rnamd.org/m6ATCPred .

Reconstruction of Eriocheir sinensis Protein-Protein Interaction Network Based on DGO-SVM Method.

Eriocheir sinensis is an economically important aquatic animal. Its regulatory mechanisms underlying many biological processes are still vague due to the lack of systematic analysis tools. The protein-protein interaction network (PIN) is an important tool for the systematic analysis of regulatory mechanisms. In this work, a novel machine learning method, DGO-SVM, was applied to predict the protein-protein interaction (PPI) in E. sinensis, and its PIN was reconstructed. With the domain, biological process, molecular functions and subcellular locations of proteins as the features, DGO-SVM showed excellent performance in Bombyx mori, humans and five aquatic crustaceans, with 92-96% accuracy. With DGO-SVM, the PIN of E. sinensis was reconstructed, containing 14,703 proteins and 7,243,597 interactions, in which 35,604 interactions were associated with 566 novel proteins mainly involved in the response to exogenous stimuli, cellular macromolecular metabolism and regulation. The DGO-SVM demonstrated that the biological process, molecular functions and subcellular locations of proteins are significant factors for the precise prediction of PPIs. We reconstructed the largest PIN for E. sinensis, which provides a systematic tool for the regulatory mechanism analysis. Furthermore, the novel-protein-related PPIs in the PIN may provide important clues for the mechanism analysis of the underlying specific physiological processes in E. sinensis.

Structural, static, and dynamic functional MRI predictors for conversion from mild cognitive impairment to Alzheimer's disease: Inter-cohort validation of Shanghai Memory Study and ADNI.

Mild cognitive impairment (MCI) is a critical prodromal stage of Alzheimer's disease (AD), and the mechanism underlying the conversion is not fully explored. Construction and inter-cohort validation of imaging biomarkers for predicting MCI conversion is of great challenge at present, due to lack of longitudinal cohorts and poor reproducibility of various study-specific imaging indices. We proposed a novel framework for inter-cohort MCI conversion prediction, involving comparison of structural, static, and dynamic functional brain features from structural magnetic resonance imaging (sMRI) and resting-state functional MRI (fMRI) between MCI converters (MCI_C) and non-converters (MCI_NC), and support vector machine for construction of prediction models. A total of 218 MCI patients with 3-year follow-up outcome were selected from two independent cohorts: Shanghai Memory Study cohort for internal cross-validation, and Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort for external validation. In comparison with MCI_NC, MCI_C were mainly characterized by atrophy, regional hyperactivity and inter-network hypo-connectivity, and dynamic alterations characterized by regional and connectional instability, involving medial temporal lobe (MTL), posterior parietal cortex (PPC), and occipital cortex. All imaging-based prediction models achieved an area under the curve (AUC) > 0.7 in both cohorts, with the multi-modality MRI models as the best with excellent performances of AUC > 0.85. Notably, the combination of static and dynamic fMRI resulted in overall better performance as relative to static or dynamic fMRI solely, supporting the contribution of dynamic features. This inter-cohort validation study provides a new insight into the mechanisms of MCI conversion involving brain dynamics, and paves a way for clinical use of structural and functional MRI biomarkers in future.

Attenuated total reflection Fourier-transform infrared spectroscopy reveals environment specific phenotypes in clonal Japanese knotweed.

BACKGROUND: Japanese knotweed (Reynoutria japonica var. japonica), a problematic invasive species, has a wide geographical distribution. We have previously shown the potential for attenuated total reflection Fourier-transform infrared (ATR-FTIR) spectroscopy and chemometrics to segregate regional differentiation between Japanese knotweed plants. However, the contribution of environment to spectral differences remains unclear. Herein, the response of Japanese knotweed to varied environmental habitats has been studied. Eight unique growth environments were created by manipulation of the red: far-red light ratio (R: FR), water availability, nitrogen, and micronutrients. Their impacts on plant growth, photosynthetic parameters, and ATR-FTIR spectral profiles, were explored using chemometric techniques, including principal component analysis (PCA), linear discriminant analysis, support vector machines (SVM) and partial least squares regression. Key wavenumbers responsible for spectral differences were identified with PCA loadings, and molecular biomarkers were assigned. Partial least squared regression (PLSR) of spectral absorbance and root water potential (RWP) data was used to create a predictive model for RWP. RESULTS: Spectra from plants grown in different environments were differentiated using ATR-FTIR spectroscopy coupled with SVM. Biomarkers highlighted through PCA loadings corresponded to several molecules, most commonly cell wall carbohydrates, suggesting that these wavenumbers could be consistent indicators of plant stress across species. R: FR most affected the ATR-FTIR spectra of intact dried leaf material. PLSR prediction of root water potential achieved an R2 of 0.8, supporting the potential use of ATR-FTIR spectrometers as sensors for prediction of plant physiological parameters. CONCLUSIONS: Japanese knotweed exhibits environmentally induced phenotypes, indicated by measurable differences in their ATR-FTIR spectra. This high environmental plasticity reflected by key biomolecular changes may contribute to its success as an invasive species. Light quality (R: FR) appears critical in defining the growth and spectral response to environment. Cross-species conservation of biomarkers suggest that they could function as indicators of plant-environment interactions including abiotic stress responses and plant health.

Novel Artificial Intelligence Combining Convolutional Neural Network and Support Vector Machine to Predict Colorectal Cancer Prognosis and Mutational Signatures From Hematoxylin and Eosin Images.

Reducing recurrence following radical resection of colon cancer without overtreatment or undertreatment remains a challenge. Postoperative adjuvant chemotherapy (Adj) is currently administered based solely on pathologic TNM stage. However, prognosis can vary significantly among patients with the same disease stage. Therefore, novel classification systems in addition to the TNM are necessary to inform decision-making regarding postoperative treatment strategies, especially stage II and III disease, and minimize overtreatment and undertreatment with Adj. We developed a prognostic prediction system for colorectal cancer using a combined convolutional neural network and support vector machine approach to extract features from hematoxylin and eosin staining images. We combined the TNM and our artificial intelligence (AI)-based classification system into a modified TNM-AI classification system with high discriminative power for recurrence-free survival. Furthermore, the cancer cell population recognized by this system as low risk of recurrence exhibited the mutational signature SBS87 as a genetic phenotype. The novel AI-based classification system developed here is expected to play an important role in prognostic prediction and personalized treatment selection in oncology.

FCM marker importance for MRD assessment in T-cell acute lymphoblastic leukemia: An AIEOP-BFM-ALL-FLOW study group report.

T-lineage acute lymphoblastic leukemia (T-ALL) accounts for about 15% of pediatric and about 25% of adult ALL cases. Minimal/measurable residual disease (MRD) assessed by flow cytometry (FCM) is an important prognostic indicator for risk stratification. In order to assess the MRD a limited number of antibodies directed against the most discriminative antigens must be selected. We propose a pipeline for evaluating the influence of different markers for cell population classification in FCM data. We use linear support vector machine, fitted to each sample individually to avoid issues with patient and laboratory variations. The best separating hyperplane direction as well as the influence of omitting specific markers is considered. Ninety-one bone marrow samples of 43 pediatric T-ALL patients from five reference laboratories were analyzed by FCM regarding marker importance for blast cell identification using combinations of eight different markers. For all laboratories, CD48 and CD99 were among the top three markers with strongest contribution to the optimal hyperplane, measured by median separating hyperplane coefficient size for all samples per center and time point (diagnosis, Day 15, Day 33). Based on the available limited set tested (CD3, CD4, CD5, CD7, CD8, CD45, CD48, CD99), our findings prove that CD48 and CD99 are useful markers for MRD monitoring in T-ALL. The proposed pipeline can be applied for evaluation of other marker combinations in the future.

Label-free cell detection of acute leukemia using ghost cytometry.

Early diagnosis and prompt initiation of appropriate treatment are critical for improving the prognosis of acute leukemia. Acute leukemia is diagnosed by microscopic morphological examination of bone marrow smears and flow cytometric immunophenotyping of bone marrow cells stained with fluorophore-conjugated antibodies. However, these diagnostic processes require trained professionals and are time and resource-intensive. Here, we present a novel diagnostic approach using ghost cytometry, a recently developed high-content flow cytometric approach, which enables machine vision-based, stain-free, high-speed analysis of cells, leveraging their detailed morphological information. We demonstrate that ghost cytometry can detect leukemic cells from the bone marrow cells of patients diagnosed with acute lymphoblastic leukemia and acute myeloid leukemia without relying on biological staining. The approach presented here holds promise as a precise, simple, swift, and cost-effective diagnostic method for acute leukemia in clinical practice.

Image-free recognition of moderate ROP from mild with machine learning algorithm on plasma Raman spectrum.

The retinopathy of prematurity (ROP) can cause serious clinical consequences and, fortunately, it is remediable while the time window for treatment is relatively narrow. Therefore, it is urgent to screen all premature infants and diagnose ROP degree timely, which has become a large workload for pediatric ophthalmologists. We developed a retinal image-free procedure using small amount of blood samples based on the plasma Raman spectrum with the machine learning model to automatically classify ROP cases before medical intervention was performed. Statistical differences in infrared Raman spectra of plasma samples were found among the control, mild (ZIIIS1), moderate (ZIIIS2 & ZIIS1), and advanced (ZIIS2) ROP groups. With the different wave points of Raman spectra as the inputs, the outputs of our support vector machine showed that the area under the curves in the receiver operating characteristic (AUC) were 0.763 for the pair comparisons of the control with the mild groups, 0.821 between moderate and advanced groups (ZIIS2), while more than 90% in comparisons of the other four pairs: control vs. moderate (0.981), control vs. advanced (0.963), mild vs. moderate (0.936), and mild vs. advanced (0.953), respectively. Our study could advance principally the ROP diagnosis in two dimensions: the moderate ROPs have been classified remarkably from the mild ones, which leaves more time for the medical treatments, and the procedure of Raman spectrum with a machine learning model based on blood samples can be conveniently promoted to those hospitals lacking of the pediatric ophthalmologists with experience in reading retinal images.