A Machine Learning Method for Identifying Critical Interactions Between Gene Pairs in Alzheimer's Disease Prediction

The ROSMAP project: aging and neurodegenerative diseases through omic sciences

The Religious Order Study and Memory and Aging Project (ROSMAP) is an initiative that integrates two longitudinal cohort studies, which have been collecting clinicopathological and molecular data since the early 1990s. This extensive dataset includes a wide array of omic data, revealing the complex interactions between molecular levels in neurodegenerative diseases (ND) and aging. Neurodegenerative diseases (ND) are frequently associated with morbidity and cognitive decline in older adults. Omics research, in conjunction with clinical variables, is crucial for advancing our understanding of the diagnosis and treatment of neurodegenerative diseases. This summary reviews the extensive omics research-encompassing genomics, transcriptomics, proteomics, metabolomics, epigenomics, and multiomics-conducted through the ROSMAP study. It highlights the significant advancements in understanding the mechanisms underlying neurodegenerative diseases, with a particular focus on Alzheimer's disease.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Applying Proteomics and Computational Approaches to Identify Novel Targets in Blast-Associated Post-Traumatic Epilepsy

Traumatic brain injury (TBI) can lead to post-traumatic epilepsy (PTE). Blast TBI (bTBI) found in Veterans presents with several complications, including cognitive and behavioral disturbances and PTE; however, the underlying mechanisms that drive the long-term sequelae are not well understood. Using an unbiased proteomics approach in a mouse model of repeated bTBI (rbTBI), this study addresses this gap in the knowledge. After rbTBI, mice were monitored using continuous, uninterrupted video-EEG for up to four months. Following this period, we collected cortex and hippocampus tissues from three groups of mice: those with post-traumatic epilepsy (PTE+), those without epilepsy (PTE-), and the control group (sham). Hundreds of differentially expressed proteins were identified in the cortex and hippocampus of PTE+ and PTE- relative to sham. Focusing on protein pathways unique to PTE+, pathways related to mitochondrial function, post-translational modifications, and transport were disrupted. Computational metabolic modeling using dysregulated protein expression predicted mitochondrial proton pump dysregulation, suggesting electron transport chain dysregulation in the epileptic tissue relative to PTE-. Finally, data mining enabled the identification of several novel and previously validated TBI and epilepsy biomarkers in our data set, many of which were found to already be targeted by drugs in various phases of clinical testing. These findings highlight novel proteins and protein pathways that may drive the chronic PTE sequelae following rbTBI.

Deep belief network-based approach for detecting Alzheimer's disease using the multi-omics data

Alzheimer's disease (AD) is the most uncertain form of Dementia in terms of finding out the mechanism. AD does not have a vital genetic factor to relate to. There were no reliable techniques and methods to identify the genetic risk factors associated with AD in the past. Most of the data available were from the brain images. However, recently, there have been drastic advancements in the high-throughput techniques in bioinformatics. It has led to focused researches in discovering the AD causing genetic risk factors. Recent analysis has resulted in considerable prefrontal cortex data with which classification and prediction models can be developed for AD. We have developed a Deep Belief Network-based prediction model using the DNA Methylation and Gene Expression Microarray Data, with High Dimension Low Sample Size (HDLSS) issues. To overcome the HDLSS challenge, we performed a two-layer feature selection considering the biological aspects of the features as well. In the two-layered feature selection approach, first the differentially expressed genes and differentially methylated positions are identified, then both the datasets are combined using Jaccard similarity measure. As the second step, an ensemble-based feature selection approach is implemented to further narrow down the gene selection. The results show that the proposed feature selection technique outperforms the existing commonly used feature selection techniques, such as Support Vector Machine Recursive Feature Elimination (SVM-RFE), and Correlation-based Feature Selection (CBS). Furthermore, the Deep Belief Network-based prediction model performs better than the widely used Machine Learning models. Also, the multi-omics dataset shows promising results compared to the single omics.

Using machine learning to identify gene interaction networks associated with breast cancer

BACKGROUND

Breast cancer (BC) is one of the most prevalent cancers worldwide but its etiology remains unclear. Obesity is recognized as a risk factor for BC, and many obesity-related genes may be involved in its occurrence and development. Research assessing the complex genetic mechanisms of BC should not only consider the effect of a single gene on the disease, but also focus on the interaction between genes. This study sought to construct a gene interaction network to identify potential pathogenic BC genes.

METHODS

The study included 953 BC patients and 963 control individuals. Chi-square analysis was used to assess the correlation between demographic characteristics and BC. The joint density-based non-parametric differential interaction network analysis and classification (JDINAC) was used to build a BC gene interaction network using single nucleotide polymorphisms (SNP). The odds ratio (OR) and 95% confidence interval (95% CI) of hub gene SNPs were evaluated using a logistic regression model. To assess reliability, the hub genes were quantified by edgeR program using BC RNA-seq data from The Cancer Genome Atlas (TCGA) and identical edges were verified by logistic regression using UK Biobank datasets. Go and KEGG enrichment analysis were used to explore the biological functions of interactive genes.

RESULTS

Body mass index (BMI) and menopause are important risk factors for BC. After adjusting for potential confounding factors, the BC gene interaction network was identified using JDINAC. LEP, LEPR, XRCC6, and RETN were identified as hub genes and both hub genes and edges were verified. LEPR genetic polymorphisms (rs1137101 and rs4655555) were also significantly associated with BC. Enrichment analysis showed that the identified genes were mainly involved in energy regulation and fat-related signaling pathways.

CONCLUSION

We explored the interaction network of genes derived from SNP data in BC progression. Gene interaction networks provide new insight into the underlying mechanisms of BC.

Machine Learning Framework for the Prediction of Alzheimer’s Disease Using Gene Expression Data Based on Efficient Gene Selection

In recent years, much research has focused on using machine learning (ML) for disease prediction based on gene expression (GE) data. However, many diseases have received considerable attention, whereas some, including Alzheimer’s disease (AD), have not, perhaps due to data shortage. The present work is intended to fill this gap by introducing a symmetric framework to predict AD from GE data, with the aim to produce the most accurate prediction using the smallest number of genes. The framework works in four stages after it receives a training dataset: pre-processing, gene selection (GS), classification, and AD prediction. The symmetry of the model is manifested in all of its stages. In the pre-processing stage gene columns in the training dataset are pre-processed identically. In the GS stage, the same user-defined filter metrics are invoked on every gene individually, and so are the same user-defined wrapper metrics. In the classification stage, a number of user-defined ML models are applied identically using the minimal set of genes selected in the preceding stage. The core of the proposed framework is a meticulous GS algorithm which we have designed to nominate eight subsets of the original set of genes provided in the training dataset. Exploring the eight subsets, the algorithm selects the best one to describe AD, and also the best ML model to predict the disease using this subset. For credible results, the framework calculates performance metrics using repeated stratified k-fold cross validation. To evaluate the framework, we used an AD dataset of 1157 cases and 39,280 genes, obtained by combining a number of smaller public datasets. The cases were split in two partitions, 1000 for training/testing, using 10-fold CV repeated 30 times, and 157 for validation. From the testing/training phase, the framework identified only 1058 genes to be the most relevant and the support vector machine (SVM) model to be the most accurate with these genes. In the final validation, we used the 157 cases that were never seen by the SVM classifier. For credible performance evaluation, we evaluated the classifier via six metrics, for which we obtained impressive values. Specifically, we obtained 0.97, 0.97, 0.98, 0.945, 0.972, and 0.975 for the sensitivity (recall), specificity, precision, kappa index, AUC, and accuracy, respectively.

Deep learning-based brain transcriptomic signatures associated with the neuropathological and clinical severity of Alzheimer's disease

Brain tissue gene expression from donors with and without Alzheimer's disease has been used to help inform the molecular changes associated with the development and potential treatment of this disorder. Here, we use a deep learning method to analyse RNA-seq data from 1114 brain donors from the Accelerating Medicines Project for Alzheimer's Disease consortium to characterize post-mortem brain transcriptome signatures associated with amyloid-β plaque, tau neurofibrillary tangles and clinical severity in multiple Alzheimer's disease dementia populations. Starting from the cross-sectional data in the Religious Orders Study and Memory and Aging Project cohort (n = 634), a deep learning framework was built to obtain a trajectory that mirrors Alzheimer's disease progression. A severity index was defined to quantitatively measure the progression based on the trajectory. Network analysis was then carried out to identify key gene (index gene) modules present in the model underlying the progression. Within this data set, severity indexes were found to be very closely correlated with all Alzheimer's disease neuropathology biomarkers (R ∼ 0.5, P < 1e-11) and global cognitive function (R = -0.68, P < 2.2e-16). We then applied the model to additional transcriptomic data sets from different brain regions (MAYO, n = 266; Mount Sinai Brain Bank, n = 214), and observed that the model remained significantly predictive (P < 1e-3) of neuropathology and clinical severity. The index genes that significantly contributed to the model were integrated with Alzheimer's disease co-expression regulatory networks, resolving four discrete gene modules that are implicated in vascular and metabolic dysfunction in different cell types, respectively. Our work demonstrates the generalizability of this signature to frontal and temporal cortex measurements and additional brain donors with Alzheimer's disease, other age-related neurological disorders and controls, and revealed that the transcriptomic network modules contribute to neuropathological and clinical disease severity. This study illustrates the promise of using deep learning methods to analyse heterogeneous omics data and discover potentially targetable molecular networks that can inform the development, treatment and prevention of neurodegenerative diseases like Alzheimer's disease.

Improving the Classification of Alzheimer's Disease Using Hybrid Gene Selection Pipeline and Deep Learning

Alzheimer’s is a progressive, irreversible, neurodegenerative brain disease. Even with prominent symptoms, it takes years to notice, decode, and reveal Alzheimer’s. However, advancements in technologies, such as imaging techniques, help in early diagnosis. Still, sometimes the results are inaccurate, which delays the treatment. Thus, the research in recent times focused on identifying the molecular biomarkers that differentiate the genotype and phenotype characteristics. However, the gene expression dataset’s generated features are huge, 1,000 or even more than 10,000. To overcome such a curse of dimensionality, feature selection techniques are introduced. We designed a gene selection pipeline combining a filter, wrapper, and unsupervised method to select the relevant genes. We combined the minimum Redundancy and maximum Relevance (mRmR), Wrapper-based Particle Swarm Optimization (WPSO), and Auto encoder to select the relevant features. We used the GSE5281 Alzheimer’s dataset from the Gene Expression Omnibus We implemented an Improved Deep Belief Network (IDBN) with simple stopping criteria after choosing the relevant genes. We used a Bayesian Optimization technique to tune the hyperparameters in the Improved Deep Belief Network. The tabulated results show that the proposed pipeline shows promising results.

MACHINE LEARNING FOR DETECTING EPISTASIS INTERACTIONS AND ITS RELEVANCE TO PERSONALIZED MEDICINE IN ALZHEIMER’S DISEASE: SYSTEMATIC REVIEW

Alzheimer’s disease (AD) is a progressive disease that attacks the brain’s neurons and causes problems in memory, thinking, and reasoning skills. Personalized Medicine (PM) needs a better and more accurate understanding of the relationship between human genetic data and complex diseases like AD. The goal of PM is to tailor the treatment of a case person to his individual properties. PM requires the prediction of a person’s disease from genetic data, and its success depends on the accurate detection of genetic biomarkers. Single Nucleotide polymorphisms (SNPs) are considered the most prevalent type of variation in the human genome. Epistasis has a biological relevance to complex diseases and has an important impact on PM. Detection of the most significant epistasis interactions associated with complex diseases is a big challenge. This paper reviews several machine learning techniques and algorithms to detect the most significant epistasis interactions in Alzheimer’s disease. We discuss many machine learning techniques that can be used for detecting SNPs’ combinations like Random Forests, Support Vector Machines, Multifactor Dimensionality Reduction, Neural Network, and Deep Learning. This review paper highlights the pros and cons of these techniques and explains how they can be applied in an efficient framework to apply knowledge discovery and data mining in AD disease.

Potential associations between immune signaling genes, deactivated microglia, and oligodendrocytes and cortical gray matter loss in patients with long-term remitted Cushing's disease

INTRODUCTION

Cushing's disease (CD) is a rare and severe endocrine disease characterized by hypercortisolemia. Previous studies have found structural brain alterations in remitted CD patients compared to healthy controls, specifically in the anterior cingulate cortex (ACC). However, potential mechanisms through which these persistent alterations may have occurred are currently unknown.

METHODS

Structural 3T MRI's from 25 remitted CD patients were linked with gene expression data from neurotypical donors, derived from the Allen Human Brain Atlas. Differences in gene expression between the ACC and an unaffected control cortical region were examined, followed by a Gene Ontology (GO) enrichment analysis. A cell type enrichment analysis was conducted on the differentially expressed genes, and a disease association enrichment analysis was conducted to determine possible associations between differentially expressed genes and specific diseases. Subsequently, cortisol sensitivity of these genes in existing datasets was examined.

RESULTS

The gene expression analysis identified 300 differentially expressed genes in the ACC compared to the cortical control region. GO analyses found underexpressed genes to represent immune function. The cell type specificity analysis indicated that underexpressed genes were enriched for deactivated microglia and oligodendrocytes. Neither significant associations with diseases, nor evidence of cortisol sensitivity with the differentially expressed genes were found.

DISCUSSION

Underexpressed genes in the ACC, the area vulnerable to permanent changes in remitted CD patients, were often associated with immune functioning. The specific lack of deactivated microglia and oligodendrocytes implicates protective effects of these cell types against the long-term effects of cortisol overexposure.

Integration and interplay of machine learning and bioinformatics approach to identify genetic interaction related to ovarian cancer chemoresistance

Although chemotherapy is the first-line treatment for ovarian cancer (OCa) patients, chemoresistance (CR) decreases their progression-free survival. This paper investigates the genetic interaction (GI) related to OCa-CR. To decrease the complexity of establishing gene networks, individual signature genes related to OCa-CR are identified using a gradient boosting decision tree algorithm. Additionally, the genetic interaction coefficient (GIC) is proposed to measure the correlation of two signature genes quantitatively and explain their joint influence on OCa-CR. Gene pair that possesses high GIC is identified as signature pair. A total of 24 signature gene pairs are selected that include 10 individual signature genes and the influence of signature gene pairs on OCa-CR is explored. Finally, a signature gene pair-based prediction of OCa-CR is identified. The area under curve (AUC) is a widely used performance measure for machine learning prediction. The AUC of signature gene pair reaches 0.9658, whereas the AUC of individual signature gene-based prediction is 0.6823 only. The identified signature gene pairs not only build an efficient GI network of OCa-CR but also provide an interesting way for OCa-CR prediction. This improvement shows that our proposed method is a useful tool to investigate GI related to OCa-CR.

Unearthing of Key Genes Driving the Pathogenesis of Alzheimer's Disease via Bioinformatics

Alzheimer’s disease (AD) is a neurodegenerative disease with unelucidated molecular pathogenesis. Herein, we aimed to identify potential hub genes governing the pathogenesis of AD. The AD datasets of GSE118553 and GSE131617 were collected from the NCBI GEO database. The weighted gene coexpression network analysis (WGCNA), differential gene expression analysis, and functional enrichment analysis were performed to reveal the hub genes and verify their role in AD. Hub genes were validated by machine learning algorithms. We identified modules and their corresponding hub genes from the temporal cortex (TC), frontal cortex (FC), entorhinal cortex (EC), and cerebellum (CE). We obtained 33, 42, 42, and 41 hub genes in modules associated with AD in TC, FC, EC, and CE tissues, respectively. Significant differences were recorded in the expression levels of hub genes between AD and the control group in the TC and EC tissues (P < 0.05). The differences in the expressions of FCGRT, SLC1A3, PTN, PTPRZ1, and PON2 in the FC and CE tissues among the AD and control groups were significant (P < 0.05). The expression levels of PLXNB1, GRAMD3, and GJA1 were statistically significant between the Braak NFT stages of AD. Overall, our study uncovered genes that may be involved in AD pathogenesis and revealed their potential for the development of AD biomarkers and appropriate AD therapeutics targets.

High-dimensional integrative copula discriminant analysis for multiomics data

Multiomics or integrative omics data have been increasingly common in biomedical studies, holding a promise in better understanding human health and disease. In this article, we propose an integrative copula discrimination analysis classifier in the context of two-class classification, which relaxes the common Gaussian assumption and gains power by borrowing information from multiple omics data types in discriminant analysis. Numerical studies are conducted to assess the finite sample performance of the new classifier. We apply our model to the Religious Orders Study and Memory and Aging Project (ROSMAP) Study, integrating gene expression and DNA methylation data for better prediction.