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.

Circulating blood circular RNA in Parkinson's Disease; a systematic study.

We aimed to identify circRNAs associated with Parkinson's disease (PD) by leveraging 1,848 participants and 1,789 circRNA from two of the largest publicly available studies with longitudinal clinical and blood transcriptomic data. To comprehensively understand changes in circRNAs we performed a cross-sectional study utilizing the last visit of each participant, and a longitudinal (mix model) analysis that included 1,166 participants with at least two time points. We identified 192 circRNAs differentially expressed in PD participants compared to healthy controls, with effects that were consistent in the mixed models, mutation carriers, and diverse ancestry. Finally, we included the 149 circRNA in a model with a ROC AUC of 0.825, showing that have the potential to aid the diagnosis of PD. Overall, we demonstrated that circRNAs play an important role in PD and can be leveraged as biomarkers.

Cell-free RNA signatures predict Alzheimer's disease.

There is a need for affordable, scalable, and specific blood-based biomarkers for Alzheimer's disease that can be applied to a population level. We have developed and validated disease-specific cell-free transcriptomic blood-based biomarkers composed by a scalable number of transcripts that capture AD pathobiology even in the presymptomatic stages of the disease. Accuracies are in the range of the current CSF and plasma biomarkers, and specificities are high against other neurodegenerative diseases.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally.

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

Artificial intelligence for dementia genetics and omics.

Genetics and omics studies of Alzheimer's disease and other dementia subtypes enhance our understanding of underlying mechanisms and pathways that can be targeted. We identified key remaining challenges: First, can we enhance genetic studies to address missing heritability? Can we identify reproducible omics signatures that differentiate between dementia subtypes? Can high-dimensional omics data identify improved biomarkers? How can genetics inform our understanding of causal status of dementia risk factors? And which biological processes are altered by dementia-related genetic variation? Artificial intelligence (AI) and machine learning approaches give us powerful new tools in helping us to tackle these challenges, and we review possible solutions and examples of best practice. However, their limitations also need to be considered, as well as the need for coordinated multidisciplinary research and diverse deeply phenotyped cohorts. Ultimately AI approaches improve our ability to interrogate genetics and omics data for precision dementia medicine. HIGHLIGHTS: We have identified five key challenges in dementia genetics and omics studies. AI can enable detection of undiscovered patterns in dementia genetics and omics data. Enhanced and more diverse genetics and omics datasets are still needed. Multidisciplinary collaborative efforts using AI can boost dementia research.

Genome-wide association study stratified by MAPT haplotypes identifies potential novel loci in Parkinson's disease.

Objective: To identify genetic factors that may modify the effects of the MAPT locus in Parkinson's disease (PD). Methods: We used data from the International Parkinson's Disease Genomics Consortium (IPDGC) and the UK biobank (UKBB). We stratified the IPDGC cohort for carriers of the H1/H1 genotype (PD patients n=8,492 and controls n=6,765) and carriers of the H2 haplotype (with either H1/H2 or H2/H2 genotypes, patients n=4,779 and controls n=4,849) to perform genome-wide association studies (GWASs). Then, we performed replication analyses in the UKBB data. To study the association of rare variants in the new nominated genes, we performed burden analyses in two cohorts (Accelerating Medicines Partnership - Parkinson Disease and UKBB) with a total sample size PD patients n=2,943 and controls n=18,486. Results: We identified a novel locus associated with PD among MAPT H1/H1 carriers near EMP1 (rs56312722, OR=0.88, 95%CI= 0.84-0.92, p= 1.80E-08), and a novel locus associated with PD among MAPT H2 carriers near VANGL1 (rs11590278, OR=1.69 95%CI=1.40-2.03, p=2.72E-08). Similar analysis of the UKBB data did not replicate these results and rs11590278 near VANGL1 did have similar effect size and direction in carriers of H2 haplotype, albeit not statistically significant (OR= 1.32, 95%CI= 0.94-1.86, p=0.17). Rare EMP1 variants with high CADD scores were associated with PD in the MAPT H2 stratified analysis (p=9.46E-05), mainly driven by the p.V11G variant. Interpretation: We identified several loci potentially associated with PD stratified by MAPT haplotype and larger replication studies are required to confirm these associations.

IntroVerse: a comprehensive database of introns across human tissues.

Dysregulation of RNA splicing contributes to both rare and complex diseases. RNA-sequencing data from human tissues has shown that this process can be inaccurate, resulting in the presence of novel introns detected at low frequency across samples and within an individual. To enable the full spectrum of intron use to be explored, we have developed IntroVerse, which offers an extensive catalogue on the splicing of 332,571 annotated introns and a linked set of 4,679,474 novel junctions covering 32,669 different genes. This dataset has been generated through the analysis of 17,510 human control RNA samples from 54 tissues provided by the Genotype-Tissue Expression Consortium. IntroVerse has two unique features: (i) it provides a complete catalogue of novel junctions and (ii) each novel junction has been assigned to a specific annotated intron. This unique, hierarchical structure offers multiple uses, including the identification of novel transcripts from known genes and their tissue-specific usage, and the assessment of background splicing noise for introns thought to be mis-spliced in disease states. IntroVerse provides a user-friendly web interface and is freely available at https://rytenlab.com/browser/app/introverse.

GenoDrawing: An Autoencoder Framework for Image Prediction from SNP Markers.

Advancements in genome sequencing have facilitated whole-genome characterization of numerous plant species, providing an abundance of genotypic data for genomic analysis. Genomic selection and neural networks (NNs), particularly deep learning, have been developed to predict complex traits from dense genotypic data. Autoencoders, an NN model to extract features from images in an unsupervised manner, has proven to be useful for plant phenotyping. This study introduces an autoencoder framework, GenoDrawing, for predicting and retrieving apple images from a low-depth single-nucleotide polymorphism (SNP) array, potentially useful in predicting traits that are difficult to define. GenoDrawing demonstrates proficiency in its task using a small dataset of shape-related SNPs. Results indicate that the use of SNPs associated with visual traits has substantial impact on the generated images, consistent with biological interpretation. While using substantial SNPs is crucial, incorporating additional, unrelated SNPs results in performance degradation for simple NN architectures that cannot easily identify the most important inputs. The proposed GenoDrawing method is a practical framework for exploring genomic prediction in fruit tree phenotyping, particularly beneficial for small to medium breeding companies to predict economically substantial heritable traits. Although GenoDrawing has limitations, it sets the groundwork for future research in image prediction from genomic markers. Future studies should focus on using stronger models for image reproduction, SNP information extraction, and dataset balance in terms of phenotypes for more precise outcomes.

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/ .

A predictive model for hospitalization and survival to COVID-19 in a retrospective population-based study.

The development of tools that provide early triage of COVID-19 patients with minimal use of diagnostic tests, based on easily accessible data, can be of vital importance in reducing COVID-19 mortality rates during high-incidence scenarios. This work proposes a machine learning model to predict mortality and risk of hospitalization using both 2 simple demographic features and 19 comorbidities obtained from 86,867 electronic medical records of COVID-19 patients, and a new method (LR-IPIP) designed to deal with data imbalance problems. The model was able to predict with high accuracy (90-93%, ROC-AUC = 0.94) the patient's final status (deceased or discharged), while its accuracy was medium (71-73%, ROC-AUC = 0.75) with respect to the risk of hospitalization. The most relevant characteristics for these models were age, sex, number of comorbidities, osteoarthritis, obesity, depression, and renal failure. Finally, to facilitate its use by clinicians, a user-friendly website has been developed ( https://alejandrocisterna.shinyapps.io/PROVIA ).

Bioinformatic approach for the discovery of cis-eQTL signals during fruit ripening of a woody species as grape (Vitis vinifera L.).

Expression quantitative trait loci (eQTLs) are associations between genetic variants, such as Single Nucleotide Polymorphisms (SNPs), and gene expression. eQTLs are an important tool to understand the genetic variance of gene expression of complex phenotypes. eQTLs analyses are common in biomedical models but are scarce in woody crop species such as fruit trees or grapes. In this study, a comprehensive bioinformatic analysis was conducted leveraging with expression data from two different growth stages, around ripening onset, of 10 genotypes of grape (Vitis vinifera L.). A total of 2170 cis-eQTL were identified in 212 gene modulated at ripening onset. The 48% of these DEGs have a known function. Among the annotated protein-coding genes, terpene synthase, auxin-regulatory factors, GRFS, ANK_REP_REGION domain-containing protein, Kinesin motor domain-containing protein and flavonol synthase were noted. This new inventory of cis-eQTLs influencing gene expression during fruit ripening will be an important resource to examine variation for this trait and will help to elucidate the complex genetic architecture underlying this process in grape.

Multi-modality machine learning predicting Parkinson's disease.

Personalized medicine promises individualized disease prediction and treatment. The convergence of machine learning (ML) and available multimodal data is key moving forward. We build upon previous work to deliver multimodal predictions of Parkinson's disease (PD) risk and systematically develop a model using GenoML, an automated ML package, to make improved multi-omic predictions of PD, validated in an external cohort. We investigated top features, constructed hypothesis-free disease-relevant networks, and investigated drug-gene interactions. We performed automated ML on multimodal data from the Parkinson's progression marker initiative (PPMI). After selecting the best performing algorithm, all PPMI data was used to tune the selected model. The model was validated in the Parkinson's Disease Biomarker Program (PDBP) dataset. Our initial model showed an area under the curve (AUC) of 89.72% for the diagnosis of PD. The tuned model was then tested for validation on external data (PDBP, AUC 85.03%). Optimizing thresholds for classification increased the diagnosis prediction accuracy and other metrics. Finally, networks were built to identify gene communities specific to PD. Combining data modalities outperforms the single biomarker paradigm. UPSIT and PRS contributed most to the predictive power of the model, but the accuracy of these are supplemented by many smaller effect transcripts and risk SNPs. Our model is best suited to identifying large groups of individuals to monitor within a health registry or biobank to prioritize for further testing. This approach allows complex predictive models to be reproducible and accessible to the community, with the package, code, and results publicly available.

Leveraging omic features with F3UTER enables identification of unannotated 3'UTRs for synaptic genes.

There is growing evidence for the importance of 3' untranslated region (3'UTR) dependent regulatory processes. However, our current human 3'UTR catalogue is incomplete. Here, we develop a machine learning-based framework, leveraging both genomic and tissue-specific transcriptomic features to predict previously unannotated 3'UTRs. We identify unannotated 3'UTRs associated with 1,563 genes across 39 human tissues, with the greatest abundance found in the brain. These unannotated 3'UTRs are significantly enriched for RNA binding protein (RBP) motifs and exhibit high human lineage-specificity. We find that brain-specific unannotated 3'UTRs are enriched for the binding motifs of important neuronal RBPs such as TARDBP and RBFOX1, and their associated genes are involved in synaptic function. Our data is shared through an online resource F3UTER ( https://astx.shinyapps.io/F3UTER/ ). Overall, our data improves 3'UTR annotation and provides additional insights into the mRNA-RBP interactome in the human brain, with implications for our understanding of neurological and neurodevelopmental diseases.

A systems-level analysis highlights microglial activation as a modifying factor in common epilepsies.

AIMS: The causes of distinct patterns of reduced cortical thickness in the common human epilepsies, detectable on neuroimaging and with important clinical consequences, are unknown. We investigated the underlying mechanisms of cortical thinning using a systems-level analysis. METHODS: Imaging-based cortical structural maps from a large-scale epilepsy neuroimaging study were overlaid with highly spatially resolved human brain gene expression data from the Allen Human Brain Atlas. Cell-type deconvolution, differential expression analysis and cell-type enrichment analyses were used to identify differences in cell-type distribution. These differences were followed up in post-mortem brain tissue from humans with epilepsy using Iba1 immunolabelling. Furthermore, to investigate a causal effect in cortical thinning, cell-type-specific depletion was used in a murine model of acquired epilepsy. RESULTS: We identified elevated fractions of microglia and endothelial cells in regions of reduced cortical thickness. Differentially expressed genes showed enrichment for microglial markers and, in particular, activated microglial states. Analysis of post-mortem brain tissue from humans with epilepsy confirmed excess activated microglia. In the murine model, transient depletion of activated microglia during the early phase of the disease development prevented cortical thinning and neuronal cell loss in the temporal cortex. Although the development of chronic seizures was unaffected, the epileptic mice with early depletion of activated microglia did not develop deficits in a non-spatial memory test seen in epileptic mice not depleted of microglia. CONCLUSIONS: These convergent data strongly implicate activated microglia in cortical thinning, representing a new dimension for concern and disease modification in the epilepsies, potentially distinct from seizure control.

A genetic link between risk for Alzheimer's disease and severe COVID-19 outcomes via the OAS1 gene.

Recently, we reported oligoadenylate synthetase 1 (OAS1) contributed to the risk of Alzheimer's disease, by its enrichment in transcriptional networks expressed by microglia. However, the function of OAS1 within microglia was not known. Using genotyping from 1313 individuals with sporadic Alzheimer's disease and 1234 control individuals, we confirm the OAS1 variant, rs1131454, is associated with increased risk for Alzheimer's disease. The same OAS1 locus has been recently associated with severe coronavirus disease 2019 (COVID-19) outcomes, linking risk for both diseases. The single nucleotide polymorphisms rs1131454(A) and rs4766676(T) are associated with Alzheimer's disease, and rs10735079(A) and rs6489867(T) are associated with severe COVID-19, where the risk alleles are linked with decreased OAS1 expression. Analysing single-cell RNA-sequencing data of myeloid cells from Alzheimer's disease and COVID-19 patients, we identify co-expression networks containing interferon (IFN)-responsive genes, including OAS1, which are significantly upregulated with age and both diseases. In human induced pluripotent stem cell-derived microglia with lowered OAS1 expression, we show exaggerated production of TNF-α with IFN-γ stimulation, indicating OAS1 is required to limit the pro-inflammatory response of myeloid cells. Collectively, our data support a link between genetic risk for Alzheimer's disease and susceptibility to critical illness with COVID-19 centred on OAS1, a finding with potential implications for future treatments of Alzheimer's disease and COVID-19, and development of biomarkers to track disease progression.

ASL expression in ALDH1A1+ neurons in the substantia nigra metabolically contributes to neurodegenerative phenotype.

Argininosuccinate lyase (ASL) is essential for the NO-dependent regulation of tyrosine hydroxylase (TH) and thus for catecholamine production. Using a conditional mouse model with loss of ASL in catecholamine neurons, we demonstrate that ASL is expressed in dopaminergic neurons in the substantia nigra pars compacta, including the ALDH1A1 + subpopulation that is pivotal for the pathogenesis of Parkinson disease (PD). Neuronal loss of ASL results in catecholamine deficiency, in accumulation and formation of tyrosine aggregates, in elevation of α-synuclein, and phenotypically in motor and cognitive deficits. NO supplementation rescues the formation of aggregates as well as the motor deficiencies. Our data point to a potential metabolic link between accumulations of tyrosine and seeding of pathological aggregates in neurons as initiators for the pathological processes involved in neurodegeneration. Hence, interventions in tyrosine metabolism via regulation of NO levels may be therapeutic beneficial for the treatment of catecholamine-related neurodegenerative disorders.

An integrated genomic approach to dissect the genetic landscape regulating the cell-to-cell transfer of α-synuclein.

Neuropathological and experimental evidence suggests that the cell-to-cell transfer of α-synuclein has an important role in the pathogenesis of Parkinson's disease (PD). However, the mechanism underlying this phenomenon is not fully understood. We undertook a small interfering RNA (siRNA), genome-wide screen to identify genes regulating the cell-to-cell transfer of α-synuclein. A genetically encoded reporter, GFP-2A-αSynuclein-RFP, suitable for separating donor and recipient cells, was transiently transfected into HEK cells stably overexpressing α-synuclein. We find that 38 genes regulate the transfer of α-synuclein-RFP, one of which is ITGA8, a candidate gene identified through a recent PD genome-wide association study (GWAS). Weighted gene co-expression network analysis (WGCNA) and weighted protein-protein network interaction analysis (WPPNIA) show that those hits cluster in networks that include known PD genes more frequently than expected by random chance. The findings expand our understanding of the mechanism of α-synuclein spread.

CoExp: A Web Tool for the Exploitation of Co-expression Networks.

Gene co-expression networks are a powerful type of analysis to construct gene groupings based on transcriptomic profiling. Co-expression networks make it possible to discover modules of genes whose mRNA levels are highly correlated across samples. Subsequent annotation of modules often reveals biological functions and/or evidence of cellular specificity for cell types implicated in the tissue being studied. There are multiple ways to perform such analyses with weighted gene co-expression network analysis (WGCNA) amongst one of the most widely used R packages. While managing a few network models can be done manually, it is often more advantageous to study a wider set of models derived from multiple independently generated transcriptomic data sets (e.g., multiple networks built from many transcriptomic sources). However, there is no software tool available that allows this to be easily achieved. Furthermore, the visual nature of co-expression networks in combination with the coding skills required to explore networks, makes the construction of a web-based platform for their management highly desirable. Here, we present the CoExp Web application, a user-friendly online tool that allows the exploitation of the full collection of 109 co-expression networks provided by the CoExpNets suite of R packages. We describe the usage of CoExp, including its contents and the functionality available through the family of CoExpNets packages. All the tools presented, including the web front- and back-ends are available for the research community so any research group can build its own suite of networks and make them accessible through their own CoExp Web application. Therefore, this paper is of interest to both researchers wishing to annotate their genes of interest across different brain network models and specialists interested in the creation of GCNs looking for a tool to appropriately manage, use, publish, and share their networks in a consistent and productive manner.

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.

Trem2 promotes anti-inflammatory responses in microglia and is suppressed under pro-inflammatory conditions.

Genome-wide association studies have reported that, amongst other microglial genes, variants in TREM2 can profoundly increase the incidence of developing Alzheimer's disease (AD). We have investigated the role of TREM2 in primary microglial cultures from wild type mice by using siRNA to decrease Trem2 expression, and in parallel from knock-in mice heterozygous or homozygous for the Trem2 R47H AD risk variant. The prevailing phenotype of Trem2 R47H knock-in mice was decreased expression levels of Trem2 in microglia, which resulted in decreased density of microglia in the hippocampus. Overall, primary microglia with reduced Trem2 expression, either by siRNA or from the R47H knock-in mice, displayed a similar phenotype. Comparison of the effects of decreased Trem2 expression under conditions of lipopolysaccharide (LPS) pro-inflammatory or IL-4 anti-inflammatory stimulation revealed the importance of Trem2 in driving a number of the genes up-regulated in the anti-inflammatory phenotype. RNA-seq analysis showed that IL-4 induced the expression of a program of genes including Arg1 and Ap1b1 in microglia, which showed an attenuated response to IL-4 when Trem2 expression was decreased. Genes showing a similar expression profile to Arg1 were enriched for STAT6 transcription factor recognition elements in their promoter, and Trem2 knockdown decreased levels of STAT6. LPS-induced pro-inflammatory stimulation suppressed Trem2 expression, thus preventing TREM2's anti-inflammatory drive. Given that anti-inflammatory signaling is associated with tissue repair, understanding the signaling mechanisms downstream of Trem2 in coordinating the pro- and anti-inflammatory balance of microglia, particularly mediating effects of the IL-4-regulated anti-inflammatory pathway, has important implications for fighting neurodegenerative disease.