PhenoExam: gene set analyses through integration of different phenotype databases.

BACKGROUND: Gene set enrichment analysis (detecting phenotypic terms that emerge as significant in a set of genes) plays an important role in bioinformatics focused on diseases of genetic basis. To facilitate phenotype-oriented gene set analysis, we developed PhenoExam, a freely available R package for tool developers and a web interface for users, which performs: (1) phenotype and disease enrichment analysis on a gene set; (2) measures statistically significant phenotype similarities between gene sets and (3) detects significant differential phenotypes or disease terms across different databases. RESULTS: PhenoExam generates sensitive and accurate phenotype enrichment analyses. It is also effective in segregating gene sets or Mendelian diseases with very similar phenotypes. We tested the tool with two similar diseases (Parkinson and dystonia), to show phenotype-level similarities but also potentially interesting differences. Moreover, we used PhenoExam to validate computationally predicted new genes potentially associated with epilepsy. CONCLUSIONS: We developed PhenoExam, a freely available R package and Web application, which performs phenotype enrichment and disease enrichment analysis on gene set G, measures statistically significant phenotype similarities between pairs of gene sets G and G' and detects statistically significant exclusive phenotypes or disease terms, across different databases. We proved with simulations and real cases that it is useful to distinguish between gene sets or diseases with very similar phenotypes. Github R package URL is https://github.com/alexcis95/PhenoExam . Shiny App URL is https://alejandrocisterna.shinyapps.io/phenoexamweb/ .

Modeling multifunctionality of genes with secondary gene co-expression networks in human brain provides novel disease insights.

MOTIVATION: Co-expression networks are a powerful gene expression analysis method to study how genes co-express together in clusters with functional coherence that usually resemble specific cell type behavior for the genes involved. They can be applied to bulk-tissue gene expression profiling and assign function, and usually cell type specificity, to a high percentage of the gene pool used to construct the network. One of the limitations of this method is that each gene is predicted to play a role in a specific set of coherent functions in a single cell type (i.e. at most we get a single for each gene). We present here GMSCA (Gene Multifunctionality Secondary Co-expression Analysis), a software tool that exploits the co-expression paradigm to increase the number of functions and cell types ascribed to a gene in bulk-tissue co-expression networks. RESULTS: We applied GMSCA to 27 co-expression networks derived from bulk-tissue gene expression profiling of a variety of brain tissues. Neurons and glial cells (microglia, astrocytes and oligodendrocytes) were considered the main cell types. Applying this approach, we increase the overall number of predicted triplets by 46.73%. Moreover, GMSCA predicts that the SNCA gene, traditionally associated to work mainly in neurons, also plays a relevant function in oligodendrocytes. AVAILABILITYAND IMPLEMENTATION: The tool is available at GitHub, https://github.com/drlaguna/GMSCA as open-source software. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Paired plasma lipidomics and proteomics analysis in the conversion from mild cognitive impairment to Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) is a neurodegenerative condition for which there is currently no available medication that can stop its progression. Previous studies suggest that mild cognitive impairment (MCI) is a phase that precedes the disease. Therefore, a better understanding of the molecular mechanisms behind MCI conversion to AD is needed. METHOD: Here, we propose a machine learning-based approach to detect the key metabolites and proteins involved in MCI progression to AD using data from the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery Study. Proteins and metabolites were evaluated separately in multiclass models (controls, MCI and AD) and together in MCI conversion models (MCI stable vs converter). Only features selected as relevant by 3/4 algorithms proposed were kept for downstream analysis. RESULTS: Multiclass models of metabolites highlighted nine features further validated in an independent cohort (0.726 mean balanced accuracy). Among these features, one metabolite, oleamide, was selected by all the algorithms. Further in-vitro experiments in rodents showed that disease-associated microglia excreted oleamide in vesicles. Multiclass models of proteins stood out with nine features, validated in an independent cohort (0.720 mean balanced accuracy). However, none of the proteins was selected by all the algorithms. Besides, to distinguish between MCI stable and converters, 14 key features were selected (0.872 AUC), including tTau, alpha-synuclein (SNCA), junctophilin-3 (JPH3), properdin (CFP) and peptidase inhibitor 15 (PI15) among others. CONCLUSIONS: This omics integration approach highlighted a set of molecules associated with MCI conversion important in neuronal and glia inflammation pathways.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications.

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.