A Personalized Metabolic Modelling Approach through Integrated Analysis of RNA-Seq-Based Genomic Variants and Gene Expression Levels in Alzheimer’s Disease

Motivation Alzheimer’s disease (AD) is known to cause alterations in brain metabolism. Furthermore, genomic variants in enzyme-coding genes may exacerbate AD-linked metabolic changes. Generating condition-specific metabolic models by mapping gene expression data to genome-scale metabolic models is a routine approach to elucidate disease mechanisms from a metabolic perspective. RNAseq data provides both gene expression and genomic variation information. Integrating variants that perturb enzyme functionality from the same RNAseq data may enhance model accuracy, offering insights into genome-wide AD metabolic pathology. Results Our study pioneers the extraction of both transcriptomic and genomic data from the same RNA-seq data to reconstruct personalized metabolic models. We mapped genes with significantly higher load of pathogenic variants in AD onto a human genome-scale metabolic network together with the gene expression data. Comparative analysis of the resulting personalized patient metabolic models with the control models showed enhanced accuracy in detecting AD-associated metabolic pathways compared to the case where only expression data was mapped on the metabolic network. Besides, several otherwise would-be missed pathways were annotated in AD by considering the effect of genomic variants. Implementation The scripts are available at https://github.com/SysBioGTU/GenomicVariantsMetabolicModels . Contact [email protected]