Analysis, Classification and Marker Discovery of Gene Expression Data with Evolving Spiking Neural Networks

Analysis of temporal gene expression data poses a significant challenge due to the combination of high dimensionality and low sample size. The purpose of this paper is to present a methodology for classification, modelling, and analysis of short time-series gene expression data using spiking neural networks (SNN) and to uncover temporal expression patterns for knowledge discovery. The classification is based on the NeuCube SNN model. Time-series gene expression data of mouse primary cortical neurons is examined as a case study. The results of the analysis are promising, indicating that SNN methodologies can be effectively used to model and analyse temporal gene expression data with surpassing performance over traditional machine learning algorithms. Additionally, a gene interaction network is constructed from the temporal gene activity modelled using the NeuCube architecture offering a new way of knowledge discovery. Future work will be directed towards using gene interactions networks to help guide pharmacological research for dementia.

BackgroundSerial Analysis of Gene Expression (SAGE) is a DNA sequencing-based method for large-scale gene expression profiling that provides an alternative to microarray analysis. Most analyses of SAGE data aimed at identifying co-expressed genes have been accomplished using various versions of clustering approaches that often result in a number of false positives.Principal FindingsHere we explore the use of seriation, a statistical approach for ordering sets of objects based on their similarity, for large-scale expression pattern discovery in SAGE data. For this specific task we implement a seriation heuristic we term ‘progressive construction of contigs’ that constructs local chains of related elements by sequentially rearranging margins of the correlation matrix. We apply the heuristic to the analysis of simulated and experimental SAGE data and compare our results to those obtained with a clustering algorithm developed specifically for SAGE data. We show using simulations that the performance of seriation compares favorably to that of the clustering algorithm on noisy SAGE data.ConclusionsWe explore the use of a seriation approach for visualization-based pattern discovery in SAGE data. Using both simulations and experimental data, we demonstrate that seriation is able to identify groups of co-expressed genes more accurately than a clustering algorithm developed specifically for SAGE data. Our results suggest that seriation is a useful method for the analysis of gene expression data whose applicability should be further pursued.

Machine Learning algorithms have been widely used for gene expression data classification, despite the fact that these data have often intrinsic limitations, such as high dimensionality and a small number of examples. Few studies try to characterize to which extent these aspects can influence the performance of the classification models induced. In this paper we compute different measures characterizing the complexity of gene expression data sets for cancer diagnosis. We then investigate how these measures relate to the classification performances achieved by support vector machines, a popular Machine Learning technique usually employed in the analysis of gene expression data. The results obtained indicate that some of the complexity indices utilized are indeed successful in explaining the difficulty involved in the classification of cancer gene expression data.

Analysis of large-scale gene expression data is a research hotspot in the field of bioinformatics, which can be used to study abnormal phenomenon in plant growth process. This paper proposes a biological knowledge integration method based on parallel clustering to select gene subsets effectively. Gene ontology is utilised to obtain the biological functional similarity, and combined with gene expression data. Parallelised affinity propagation algorithm is used to cluster data since it can not only obtain more biologically meaningful subsets, but also avoid the loss of some potential value in genes from simple gene primary selection. The algorithm is verified with four typical plant datasets and compared with other well-known integration methods. Experimental results on plant stress response datasets demonstrate that the proposed method can select genes with stronger classification ability.

Recent research has focused on incorporating biological function and pathway information into the analysis of gene expression data, partly as a means of compensating for insufficient experimental replications, low signal to noise, lack of reproducibility and/or multiple testing confounds. A Bayesian approach seems to be ideal for incorporating functional information into gene expression data analysis. In this study, we tested the feasibility of using literature derived gene relationships in a Bayesian model to analyze gene expression data. Prior distributions were constructed based on gene associations derived from the biomedical literature using Latent Semantic Indexing (LSI). The LSI model was built using more than 1 million Medline abstracts corresponding to 22,000 human and mouse genes. A key advantage of LSI is that both explicit and implicit gene relationships can be derived from the literature. Gene neighborhoods were determined using latent Gaussian Markov random fields and logistic transformation of the latent variables. We tested the procedure on a microarray dataset for interferon-stimulated genes in mouse embryonic fibroblasts. By integrating functional information from literature, Bayesian approach identified relevant genes that previously did not meet the 0.05 significance level. In comparison to a standard mixture model, spatial mixture model has more power for identifying direct and indirect interferon regulated genes. The spatial model enhanced the ranks of some genes which are known to be affected by interferon treatment, such as Nmi (NMI N-myc and STAT interactor) and ifi35 (interferon-induced protein 35). It also identified some genes that previously were ignored because of the marginal p-values, such as dpysl2, map2k1, msn, Psck5, and Il6st. Interestingly, these genes appear to be indirectly related to interferon treatment. In summary, we show that our procedure increases statistical power and produces more biologically meaningful gene lists. These results suggest that Bayesian methods which incorporate functional information from the literature may improve analysis of gene expression data.

Background - In the NeoTRIPaPDL1 phase III trial, triple-negative breast cancer (TNBC) patients were randomized to receive nab-paclitaxel/carboplatin for 8 cycles (CT arm) with or without atezolizumab (CT/A arm). We previously reported that the majority of patients (65%) with PD-L1- at baseline converted to PD-L1+ (the majority with IC2/IC3) after the first cycle of treatment in the atezolizumab arm. Here we studied treatment associated changes in the molecular tumour features and immune microenvironment by PD-L1 expression. Methods – A total of 158 (56.4%) of patients enrolled in the NeoTRIPaPDL1 trial (76 and 82 from CT/A and CT arm, respectively) were included in this analysis by satisfying the following criteria: i) availability of a paired core biopsy at baseline and at day 1 of the second treatment cycle (D1C2); ii) evaluation of stromal TILs (sTILs) and staining for PD-L1 (Ventana SP142); iii) successful RNA-seq gene expression profiling. PD-L1 groups were defined as IC0/IC1 (<5% PD-L1low) vs IC2/IC3 (>=5% PD-L1high). Presence of different immune cell populations was quantified by gene expression profile deconvolution using ConsensusTME R package. Cancer hallmark gene set collection and custom signature activation status were estimated in each sample using singscore R package. Differences in score distribution were evaluated by 2-sided t-test. Results At baseline, both sTILs and immune cell signatures were upregulated in PD-L1high compared to PD-L1low tumours (p<0.001). No significant differences were found between the two treatment arms in the PD-L1high subpopulation. In the PD-L1low cases, both sTILs and immune related signatures were slightly downregulated in the CT/A arm (p<0.05), suggesting a modest unbalance among treatment arms with CT arm being slightly more inflamed than CT/A arm.At D1C2, PD-L1high tumours in the CT arm systematically had high sTILs (median=70%, range=30-90%), while PD-L1low tumours receiving CT/A had low sTILs in a significant proportion of cases (median=30%, range=0-90%; p<0.001). Similarly, at D1C2 several gene expression-estimated immune cell populations and immune-related signatures were upregulated in the CT arm compared to CT/A arm in PD-L1high tumours, with the weakest association observed for M2 macrophages (p=0.058). No tumour-related signatures were differentially expressed among the two treatment arms within groups with PD-L1high, suggesting that different treatment. modulate PD-L1 by engaging a different immune mileau instead of modulating tumor related features.Considering PD-L1low groups at D1C2, in the CT arm 21 tumors had >30% sTILs (n=21/69, 30.4%), while in the CT/A arm only 2 had sTILs >30% (2/27, 7.4%) (p<0.001). Analysis of gene expression data identified IFN-related signatures as the most upregulated in the CT compared to CT/A arm in PD-L1low cases at D1C2 (p<0.05). Conclusions Integrated dynamic analysis of PD-L1 expression and gene expression data highlighted significant treatment-specific changes of the immune landscape according to PD-L1 expression, when this biomarker is assessed during treatment. This indicates that the immune milieu associated with PD-L1 status is strongly dependent on the time of assessment in relationship to treatment received. Such observation may explain why in the NeoTRIP trial, baseline PD-L1 but not on-treatment PD-L1 was predictive of pCR in CT/A arm. In addition, our findings could have implications related to the use of PD-L1 as a predictive biomarker in pre-treated patients, especially when assessed early on during treatment. Citation Format: Maurizio Callari, Chiun-Sheng Huang, Daniel Egle, Begoña Bermejo, Claudio Zamagni, Matteo Dugo, Marc Thill, Antonio Anton, Marco Barreca, Stefania Russo, Eva Maria Ciruelos, Richard Greil, Stefania Zambelli, Balázs Gyorffy, Chanel Smart, Olivia Biasi, Pinuccia Valagussa, Giuseppe Viale, Luca Gianni, Giampaolo Bianchini. Immune milieu associated with PD-L1 status in TNBC is dependent on time of biomarker assessment and treatment received: A secondary analysis of the NeoTRIPaPDL1 trial [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-04-02.

DNA microarray data sets have been widely explored and used to analyze data without any previous biological background. However, analyzing them becomes challenging if data are missing. Thus, machine learning techniques are applied because microarray technology is promising in genomics, especially in the analysis of gene expression data. Furthermore, gene expression data can describe the transcription and translation processes of each genetic information in detail. In this study, a new system was proposed to impute more realizable values for missing data in a microarray dataset. This system was validated and evaluated on 42 samples of rectal cancer. Several evaluation tests were also conducted to confirm the effectiveness of the new system and compare it with highly known imputing algorithms. The proposed clustering column-mean quantile median technique could predict highly informative missing genes, thereby reducing the difference between the original and imputed datasets and demonstrating its efficiency.

ABSTRACT: The field of machine learning (ML) has made significant contributions to the medical domain, notably in the automatic classification of cancer types such as colon cancer. This type of cancer, while generally more common among the elderly, can affect people at any age. The incorporation of ML techniques has minimized the reliance on manual examination and enhanced the capacity to process large datasets in medical research. Utilizing advanced gene technology and analysis of gene expression data, this study addresses the challenges posed by the high dimensionality and limited size of such datasets. The findings reveal that most of the ML models employed here have surpassed expected accuracy levels, highlighting their effectiveness in the precise classification and diagnosis of colon cancer. Consequently, these developments hold promise for improving patient care and the overall efficiency of healthcare services in dealing with colon cancer.

The purpose of this paper is to introduce transductive inference with support vector machines (TSVM) as a powerful methodology for classification of gene expression data, using training and prediction data sets. The following classification problems will be considered: determination of cancer diagnosis categories and classification of genes from the budding yeast Saccharomyces cerevisiae in functional groups. In the case of training samples, experts have already classified the samples in their respective classes. So, given each prediction sample, the purpose is to determine its corresponding class. The main aspect of TSVM is that the classification task will be implemented in just one step, improving the generalization capability of the classifier. The TSVM will be compared with the traditional inductive method (SVM) in a series of experiments concerning the two classification problems, with promising results.

The application of mathematics in the field of bioinformatics has been widely developed. For example Support Vector Machines (SVM) and Random Forest (RF) are state of the art for classification of cancer in many applications. One of them is Chronic Kidney Disease (CKD). CKD is one of the kidney diseases that sufferers are increasing and have symptoms that are difficult to detect at first. Later, microarrays in gene expression are important tools for this approach. Microarrays gene expression provides an overview of all transcription activities in biological samples. The purpose of this research is a hybrid model combining Random Forest (RF) and Support Vector Machine (SVM) can be used to classify gene expression data. RF can highly accurate, generelize better and are interpretable and SVM (called RF-SVM) to effectively predict gene expression data with very high dimensions. In addition, from the simulation results on data from the Gene Expression Omnibus (GEO) database, it is shown that the proposed RF-SVM is a more accurate algorithm on CKD data than RFE-SVM.

Although there have been several analyses for identifying cancer-associated pathways, based on gene expression data, most of these are based on single pathway analyses, and thus do not consider correlations between pathways. In this paper, we propose a hierarchical structural component model for pathway analysis of gene expression data (HisCoM-PAGE), which accounts for the hierarchical structure of genes and pathways, as well as the correlations among pathways. Specifically, HisCoM-PAGE focuses on the survival phenotype and identifies its associated pathways. Moreover, its application to real biological data analysis of pancreatic cancer data demonstrated that HisCoM-PAGE could successfully identify pathways associated with pancreatic cancer prognosis. Simulation studies comparing the performance of HisCoM-PAGE with other competing methods such as Gene Set Enrichment Analysis (GSEA), Global Test, and Wald-type Test showed HisCoM-PAGE to have the highest power to detect causal pathways in most simulation scenarios.

To identify pathways associated with survival phenotypes using gene expression data, we recently proposed the hierarchical structural component model for pathway analysis of gene expression data (HisCoM-PAGE) method. The HisCoM-PAGE software can consider hierarchical structural relationships between genes and pathways and analyze multiple pathways simultaneously. It can be applied to various types of gene expression data, such as microarray data or RNA sequencing data. We expect that the HisCoM-PAGE software will make our method more easily accessible to researchers who want to perform pathway analysis for survival times.

Evolving spiking neural networks for personalised modelling, classification and prediction of spatio-temporal patterns with a case study on stroke

The emergence of large-scale multi-omics data warrants method development for data integration. Genomic studies from cancer patients have identified epigenetic and genetic regulators – such as methylation marks, somatic mutations, and somatic copy number alterations (SCNAs), among others – as predictive features of cancer outcome. However, identification of “driver genes” associated with a given alteration remains a challenge. To this end, we developed a computational tool, iEDGE, to model cis and trans effects of (epi-)DNA alterations and identify potential cis driver genes, where cis and trans genes denote those genes falling within and outside the genomic boundaries of a given (epi-)genetic alteration, respectively. iEDGE first identifies the cis and trans gene expression signatures associated with the presence/absence of a particular epi-DNA alteration across samples. It then applies tests of statistical mediation to determine the cis genes predictive of the trans gene expression. Finally, cis and trans effects are annotated by pathway enrichment analysis to gain insights into the underlying regulatory networks. We used iEDGE to perform integrative analysis of SCNAs and gene expression data from breast cancer and 18 additional cancer types included in The Cancer Genome Atlas (TCGA). Notably, cis gene drivers identified by iEDGE were found to be significantly enriched for known driver genes from multiple compendia of validated oncogenes and tumor suppressors, suggesting that the remainder are of equal importance. Furthermore, predicted drivers were enriched for functionally relevant cancer genes with amplification-driven dependencies, which are of potential prognostic and therapeutic value. All the analyses results are accessible at https://montilab.bu.edu/iEDGE. In summary, integrative analysis of SCNAs and gene expression using iEDGE successfully identified known cancer driver genes and putative cancer therapeutic targets across 19 cancer types in the TCGA. The proposed method can easily be applied to the integration of gene expression profiles with other epi-DNA assays in a variety of disease contexts.

BackgroundMicroarrays have become extremely useful for analysing genetic phenomena, but establishing a relation between microarray analysis results (typically a list of genes) and their biological significance is often difficult. Currently, the standard approach is to map a posteriori the results onto gene networks in order to elucidate the functions perturbed at the level of pathways. However, integrating a priori knowledge of the gene networks could help in the statistical analysis of gene expression data and in their biological interpretation.ResultsWe propose a method to integrate a priori the knowledge of a gene network in the analysis of gene expression data. The approach is based on the spectral decomposition of gene expression profiles with respect to the eigenfunctions of the graph, resulting in an attenuation of the high-frequency components of the expression profiles with respect to the topology of the graph. We show how to derive unsupervised and supervised classification algorithms of expression profiles, resulting in classifiers with biological relevance. We illustrate the method with the analysis of a set of expression profiles from irradiated and non-irradiated yeast strains.ConclusionIncluding a priori knowledge of a gene network for the analysis of gene expression data leads to good classification performance and improved interpretability of the results.