Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches

Introduction

The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing.

Methods

Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors.

Results

Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19.

Discussion

The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

Early-to-mid stage idiopathic Parkinson's disease shows enhanced cytotoxicity and differentiation in CD8 T-cells in females

Neuroinflammation in the brain contributes to the pathogenesis of Parkinson's disease (PD), but the potential dysregulation of peripheral immunity has not been systematically investigated for idiopathic PD (iPD). Here we showed an elevated peripheral cytotoxic immune milieu, with more terminally-differentiated effector memory (TEMRA) CD8 T, CD8+ NKT cells and circulating cytotoxic molecules in fresh blood of patients with early-to-mid iPD, especially females, after analyzing > 700 innate and adaptive immune features. This profile, also reflected by fewer CD8+FOXP3+ T cells, was confirmed in another subcohort. Co-expression between cytotoxic molecules was selectively enhanced in CD8 TEMRA and effector memory (TEM) cells. Single-cell RNA-sequencing analysis demonstrated the accelerated differentiation within CD8 compartments, enhanced cytotoxic pathways in CD8 TEMRA and TEM cells, while CD8 central memory (TCM) and naïve cells were already more-active and transcriptionally-reprogrammed. Our work provides a comprehensive map of dysregulated peripheral immunity in iPD, proposing candidates for early diagnosis and treatments.

Age at onset as stratifier in idiopathic Parkinson’s disease – effect of ageing and polygenic risk score on clinical phenotypes

Several phenotypic differences observed in Parkinson’s disease (PD) patients have been linked to age at onset (AAO). We endeavoured to find out whether these differences are due to the ageing process itself by using a combined dataset of idiopathic PD (n = 430) and healthy controls (HC; n = 556) excluding carriers of known PD-linked genetic mutations in both groups. We found several significant effects of AAO on motor and non-motor symptoms in PD, but when comparing the effects of age on these symptoms with HC (using age at assessment, AAA), only positive associations of AAA with burden of motor symptoms and cognitive impairment were significantly different between PD vs HC. Furthermore, we explored a potential effect of polygenic risk score (PRS) on clinical phenotype and identified a significant inverse correlation of AAO and PRS in PD. No significant association between PRS and severity of clinical symptoms was found. We conclude that the observed non-motor phenotypic differences in PD based on AAO are largely driven by the ageing process itself and not by a specific profile of neurodegeneration linked to AAO in the idiopathic PD patients.

PARK7/DJ-1 promotes pyruvate dehydrogenase activity and maintains Treg homeostasis during ageing

Pyruvate dehydrogenase (PDH) is the gatekeeper enzyme of the tricarboxylic acid (TCA) cycle. Here we show that the deglycase DJ-1 (encoded by PARK7, a key familial Parkinson's disease gene) is a pacemaker regulating PDH activity in CD4⁺ regulatory T cells (T_reg cells). DJ-1 binds to PDHE1-β (PDHB), inhibiting phosphorylation of PDHE1-α (PDHA), thus promoting PDH activity and oxidative phosphorylation (OXPHOS). Park7 (Dj-1) deletion impairs T_reg survival starting in young mice and reduces T_reg homeostatic proliferation and cellularity only in aged mice. This leads to increased severity in aged mice during the remission of experimental autoimmune encephalomyelitis (EAE). Dj-1 deletion also compromises differentiation of inducible T_reg cells especially in aged mice, and the impairment occurs via regulation of PDHB. These findings provide unforeseen insight into the complicated regulatory machinery of the PDH complex. As T_reg homeostasis is dysregulated in many complex diseases, the DJ-1-PDHB axis represents a potential target to maintain or re-establish T_reg homeostasis.

DJ-1 depletion prevents immunoaging in T-cell compartments

Decline in immune function during aging increases susceptibility to different aging‐related diseases. However, the underlying molecular mechanisms, especially the genetic factors contributing to imbalance of naïve/memory T‐cell subpopulations, still remain largely elusive. Here, we show that loss of DJ‐1 encoded by PARK7/DJ‐1, causing early‐onset familial Parkinson’s disease (PD), unexpectedly diminished signs of immunoaging in T‐cell compartments of both human and mice. Compared with two gender‐matched unaffected siblings of similar ages, the index PD patient with DJ‐1 deficiency showed a decline in many critical immunoaging features, including almost doubled non‐senescent T cells. The observation was further consolidated by the results in 45‐week‐old DJ‐1 knockout mice. Our data demonstrated that DJ‐1 regulates several immunoaging features via hematopoietic‐intrinsic and naïve‐CD8‐intrinsic mechanisms. Mechanistically, DJ‐1 depletion reduced oxidative phosphorylation (OXPHOS) and impaired TCR sensitivity in naïve CD8 T cells at a young age, accumulatively leading to a reduced aging process in T‐cell compartments in older mice. Our finding suggests an unrecognized critical role of DJ‐1 in regulating immunoaging, discovering a potent target to interfere with immunoaging‐ and aging‐associated diseases.

Neurodegeneration and neuroinflammation are linked, but independent of alpha‐synuclein inclusions, in a seeding/spreading mouse model of Parkinson's disease

A key pathological process in Parkinson's disease (PD) is the transneuronal spreading of α‐synuclein. Alpha‐synuclein (α‐syn) is a presynaptic protein that, in PD, forms pathological inclusions. Other hallmarks of PD include neurodegeneration and microgliosis in susceptible brain regions. Whether it is primarily transneuronal spreading of α‐syn particles, inclusion formation, or other mechanisms, such as inflammation, that cause neurodegeneration in PD is unclear. We used a model of spreading of α‐syn induced by striatal injection of α‐syn preformed fibrils into the mouse striatum to address this question. We performed quantitative analysis for α‐syn inclusions, neurodegeneration, and microgliosis in different brain regions, and generated gene expression profiles of the ventral midbrain, at two different timepoints after disease induction. We observed significant neurodegeneration and microgliosis in brain regions not only with, but also without α‐syn inclusions. We also observed prominent microgliosis in injured brain regions that did not correlate with neurodegeneration nor with inclusion load. Using longitudinal gene expression profiling, we observed early gene expression changes, linked to neuroinflammation, that preceded neurodegeneration, indicating an active role of microglia in this process. Altered gene pathways overlapped with those typical of PD. Our observations indicate that α‐syn inclusion formation is not the major driver in the early phases of PD‐like neurodegeneration, but that microglia, activated by diffusible, oligomeric α‐syn, may play a key role in this process. Our findings uncover new features of α‐syn induced pathologies, in particular microgliosis, and point to the necessity for a broader view of the process of α‐syn spreading.

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

A look into the future of the COVID-19 pandemic in Europe: an expert consultation

How will the coronavirus disease 2019 (COVID-19) pandemic develop in the coming months and years? Based on an expert survey, we examine key aspects that are likely to influence the COVID-19 pandemic in Europe. The challenges and developments will strongly depend on the progress of national and global vaccination programs, the emergence and spread of variants of concern (VOCs), and public responses to non-pharmaceutical interventions (NPIs). In the short term, many people remain unvaccinated, VOCs continue to emerge and spread, and mobility and population mixing are expected to increase. Therefore, lifting restrictions too much and too early risk another damaging wave. This challenge remains despite the reduced opportunities for transmission given vaccination progress and reduced indoor mixing in summer 2021. In autumn 2021, increased indoor activity might accelerate the spread again, whilst a necessary reintroduction of NPIs might be too slow. The incidence may strongly rise again, possibly filling intensive care units, if vaccination levels are not high enough. A moderate, adaptive level of NPIs will thus remain necessary. These epidemiological aspects combined with economic, social, and health-related consequences provide a more holistic perspective on the future of the COVID-19 pandemic.

Quantitative trait locus mapping identifies a locus linked to striatal dopamine and points to collagen IV alpha-6 chain as a novel regulator of striatal axonal branching in mice

Dopaminergic neurons (DA neurons) are controlled by multiple factors, many involved in neurological disease. Parkinson's disease motor symptoms are caused by the demise of nigral DA neurons, leading to loss of striatal dopamine (DA). Here, we measured DA concentration in the dorsal striatum of 32 members of Collaborative Cross (CC) family and their eight founder strains. Striatal DA varied greatly in founders, and differences were highly heritable in the inbred CC progeny. We identified a locus, containing 164 genes, linked to DA concentration in the dorsal striatum on chromosome X. We used RNAseq profiling of the ventral midbrain of two founders with substantial difference in striatal DA-C56BL/6 J and A/J-to highlight potential protein-coding candidates modulating this trait. Among the five differentially expressed genes within the locus, we found that the gene coding for the collagen IV alpha 6 chain (Col4a6) was expressed nine times less in A/J than in C57BL/6J. Using single cell RNA-seq data from developing human midbrain, we found that COL4A6 is highly expressed in radial glia-like cells and neuronal progenitors, indicating a role in neuronal development. Collagen IV alpha-6 chain (COL4A6) controls axogenesis in simple model organisms. Consistent with these findings, A/J mice had less striatal axonal branching than C57BL/6J mice. We tentatively conclude that DA concentration and axonal branching in dorsal striatum are modulated by COL4A6, possibly during development. Our study shows that genetic mapping based on an easily measured Central Nervous System (CNS) trait, using the CC population, combined with follow-up observations, can parse heritability of such a trait, and nominate novel functions for commonly expressed proteins.

Standard Peripheral Blood Mononuclear Cell Cryopreservation Selectively Decreases Detection of Nine Clinically Relevant T Cell Markers

Biobanking is an operational component of various epidemiological studies and clinical trials. Although peripheral blood is routinely acquired and stored in biobanks, the effects of specimen processing on cell composition and clinically relevant functional markers of T cells still require a systematic evaluation. In this study, we assessed 25 relevant T cell markers in human PBMCs and showed that the detection of nine membrane markers (e.g., PD-1, CTLA4, KLRG1, CD25, CD122, CD127, CCR7, and others reflecting exhaustion, senescence, and other functions) was reduced among at least one T cell subset following standard processing, although the frequency of CD4, CD8, and regulatory T cells was unaffected. Nevertheless, a 6-mo-long cryopreservation did not impair the percentages of cells expressing many other membrane and all the eight tested intracellular lineage or functional T cell markers. Our findings uncover that several clinically relevant markers are particularly affected by processing and the interpretation of those results in clinical trials and translational research should be done with caution.

Identification of VIMP as a gene inhibiting cytokine production in human CD4+ effector T cells

Many players regulating the CD4+ T cell-mediated inflammatory response have already been identified. However, the critical nodes that constitute the regulatory and signaling networks underlying CD4 T cell responses are still missing. Using a correlation-network-guided approach, here we identified VIMP (VCP-interacting membrane protein), one of the 25 genes encoding selenoproteins in humans, as a gene regulating the effector functions of human CD4 T cells, especially production of several cytokines including IL2 and CSF2. We identified VIMP as an endogenous inhibitor of cytokine production in CD4 effector T cells via both the E2F5 transcription regulatory pathway and the Ca2+/NFATC2 signaling pathway. Our work not only indicates that VIMP might be a promising therapeutic target for various inflammation-associated diseases but also shows that our network-guided approach can significantly aid in predicting new functions of the genes of interest.

Evaluating the Use of Circulating MicroRNA Profiles for Lung Cancer Detection in Symptomatic Patients

Importance

The overall low survival rate of patients with lung cancer calls for improved detection tools to enable better treatment options and improved patient outcomes. Multivariable molecular signatures, such as blood-borne microRNA (miRNA) signatures, may have high rates of sensitivity and specificity but require additional studies with large cohorts and standardized measurements to confirm the generalizability of miRNA signatures.

Objective

To investigate the use of blood-borne miRNAs as potential circulating markers for detecting lung cancer in an extended cohort of symptomatic patients and control participants.

Design, Setting, and Participants

This multicenter, cohort study included patients from case-control and cohort studies (TREND and COSYCONET) with 3102 patients being enrolled by convenience sampling between March 3, 2009, and March 19, 2018. For the cohort study TREND, population sampling was performed. Clinical diagnoses were obtained for 3046 patients (606 patients with non-small cell and small cell lung cancer, 593 patients with nontumor lung diseases, 883 patients with diseases not affecting the lung, and 964 unaffected control participants). No samples were removed because of experimental issues. The collected data were analyzed between April 2018 and November 2019.

Main Outcomes and Measures

Sensitivity and specificity of liquid biopsy using miRNA signatures for detection of lung cancer.

Results

A total of 3102 patients with a mean (SD) age of 61.1 (16.2) years were enrolled. Data on the sex of the participants were available for 2856 participants; 1727 (60.5%) were men. Genome-wide miRNA profiles of blood samples from 3046 individuals were evaluated by machine-learning methods. Three classification scenarios were investigated by splitting the samples equally into training and validation sets. First, a 15-miRNA signature from the training set was used to distinguish patients diagnosed with lung cancer from all other individuals in the validation set with an accuracy of 91.4% (95% CI, 91.0%-91.9%), a sensitivity of 82.8% (95% CI, 81.5%-84.1%), and a specificity of 93.5% (95% CI, 93.2%-93.8%). Second, a 14-miRNA signature from the training set was used to distinguish patients with lung cancer from patients with nontumor lung diseases in the validation set with an accuracy of 92.5% (95% CI, 92.1%-92.9%), sensitivity of 96.4% (95% CI, 95.9%-96.9%), and specificity of 88.6% (95% CI, 88.1%-89.2%). Third, a 14-miRNA signature from the training set was used to distinguish patients with early-stage lung cancer from all individuals without lung cancer in the validation set with an accuracy of 95.9% (95% CI, 95.7%-96.2%), sensitivity of 76.3% (95% CI, 74.5%-78.0%), and specificity of 97.5% (95% CI, 97.2%-97.7%).

Conclusions and Relevance

The findings of the study suggest that the identified patterns of miRNAs may be used as a component of a minimally invasive lung cancer test, complementing imaging, sputum cytology, and biopsy tests.

Common diseases alter the physiological age-related blood microRNA profile

Aging is a key risk factor for chronic diseases of the elderly. MicroRNAs regulate post-transcriptional gene silencing through base-pair binding on their target mRNAs. We identified nonlinear changes in age-related microRNAs by analyzing whole blood from 1334 healthy individuals. We observed a larger influence of the age as compared to the sex and provide evidence for a shift to the 5' mature form of miRNAs in healthy aging. The addition of 3059 diseased patients uncovered pan-disease and disease-specific alterations in aging profiles. Disease biomarker sets for all diseases were different between young and old patients. Computational deconvolution of whole-blood miRNAs into blood cell types suggests that cell intrinsic gene expression changes may impart greater significance than cell abundance changes to the whole blood miRNA profile. Altogether, these data provide a foundation for understanding the relationship between healthy aging and disease, and for the development of age-specific disease biomarkers.

ROS networks: designs, aging, Parkinson's disease and precision therapies

How the network around ROS protects against oxidative stress and Parkinson's disease (PD), and how processes at the minutes timescale cause disease and aging after decades, remains enigmatic. Challenging whether the ROS network is as complex as it seems, we built a fairly comprehensive version thereof which we disentangled into a hierarchy of only five simpler subnetworks each delivering one type of robustness. The comprehensive dynamic model described in vitro data sets from two independent laboratories. Notwithstanding its five-fold robustness, it exhibited a relatively sudden breakdown, after some 80 years of virtually steady performance: it predicted aging. PD-related conditions such as lack of DJ-1 protein or increased α-synuclein accelerated the collapse, while antioxidants or caffeine retarded it. Introducing a new concept (aging-time-control coefficient), we found that as many as 25 out of 57 molecular processes controlled aging. We identified new targets for "life-extending interventions": mitochondrial synthesis, KEAP1 degradation, and p62 metabolism.

Pituitary Tumor Transforming Gene 1 Orchestrates Gene Regulatory Variation in Mouse Ventral Midbrain During Aging

Dopaminergic neurons in the midbrain are of particular interest due to their role in diseases such as Parkinson’s disease and schizophrenia. Genetic variation between individuals can affect the integrity and function of dopaminergic neurons but the DNA variants and molecular cascades modulating dopaminergic neurons and other cells types of ventral midbrain remain poorly defined. Three genetically diverse inbred mouse strains – C57BL/6J, A/J, and DBA/2J – differ significantly in their genomes (∼7 million variants), motor and cognitive behavior, and susceptibility to neurotoxins. To further dissect the underlying molecular networks responsible for these variable phenotypes, we generated RNA-seq and ChIP-seq data from ventral midbrains of the 3 mouse strains. We defined 1000–1200 transcripts that are differentially expressed among them. These widespread differences may be due to altered activity or expression of upstream transcription factors. Interestingly, transcription factors were significantly underrepresented among the differentially expressed genes, and only one transcription factor, Pttg1, showed significant differences between all three strains. The changes in Pttg1 expression were accompanied by consistent alterations in histone H3 lysine 4 trimethylation at Pttg1 transcription start site. The ventral midbrain transcriptome of 3-month-old C57BL/6J congenic Pttg1^–/– mutants was only modestly altered, but shifted toward that of A/J and DBA/2J in 9-month-old mice. Principle component analysis (PCA) identified the genes underlying the transcriptome shift and deconvolution of these bulk RNA-seq changes using midbrain single cell RNA-seq data suggested that the changes were occurring in several different cell types, including neurons, oligodendrocytes, and astrocytes. Taken together, our results show that Pttg1 contributes to gene regulatory variation between mouse strains and influences mouse midbrain transcriptome during aging.

COVID-19 Disease Map, building a computational repository of SARS-CoV-2 virus-host interaction mechanisms

Researchers around the world join forces to reconstruct the molecular processes of the virus-host interactions aiming to combat the cause of the ongoing pandemic.

Towards a European health research and innovation cloud (HRIC)

The European Union (EU) initiative on the Digital Transformation of Health and Care (Digicare) aims to provide the conditions necessary for building a secure, flexible, and decentralized digital health infrastructure. Creating a European Health Research and Innovation Cloud (HRIC) within this environment should enable data sharing and analysis for health research across the EU, in compliance with data protection legislation while preserving the full trust of the participants. Such a HRIC should learn from and build on existing data infrastructures, integrate best practices, and focus on the concrete needs of the community in terms of technologies, governance, management, regulation, and ethics requirements. Here, we describe the vision and expected benefits of digital data sharing in health research activities and present a roadmap that fosters the opportunities while answering the challenges of implementing a HRIC. For this, we put forward five specific recommendations and action points to ensure that a European HRIC: i) is built on established standards and guidelines, providing cloud technologies through an open and decentralized infrastructure; ii) is developed and certified to the highest standards of interoperability and data security that can be trusted by all stakeholders; iii) is supported by a robust ethical and legal framework that is compliant with the EU General Data Protection Regulation (GDPR); iv) establishes a proper environment for the training of new generations of data and medical scientists; and v) stimulates research and innovation in transnational collaborations through public and private initiatives and partnerships funded by the EU through Horizon 2020 and Horizon Europe.

Fibroblast mitochondria in idiopathic Parkinson's disease display morphological changes and enhanced resistance to depolarization

Mitochondrial dysfunction is a hallmark in idiopathic Parkinson's disease (IPD). Here, we established screenable phenotypes of mitochondrial morphology and function in primary fibroblasts derived from patients with IPD. Upper arm punch skin biopsy was performed in 41 patients with mid-stage IPD and 21 age-matched healthy controls. At the single-cell level, the basal mitochondrial membrane potential (Ψm) was higher in patients with IPD than in controls. Similarly, under carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP) stress, the remaining Ψm was increased in patients with IPD. Analysis of mitochondrial morphometric parameters revealed significantly decreased mitochondrial connectivity in patients with IPD, with 9 of 14 morphometric mitochondrial parameters differing from those in controls. Significant morphometric mitochondrial changes included the node degree, mean volume, skeleton size, perimeter, form factor, node count, erosion body count, endpoints, and mitochondria count (all P-values < 0.05). These functional data reveal that resistance to depolarization was increased by treatment with the protonophore FCCP in patients with IPD, whereas morphometric data revealed decreased mitochondrial connectivity and increased mitochondrial fragmentation.

Connecting environmental exposure and neurodegeneration using cheminformatics and high resolution mass spectrometry: potential and challenges

Connecting chemical exposures over a lifetime to complex chronic diseases with multifactorial causes such as neurodegenerative diseases is an immense challenge requiring a long-term, interdisciplinary approach. Rapid developments in analytical and data technologies, such as non-target high resolution mass spectrometry (NT-HR-MS), have opened up new possibilities to accomplish this, inconceivable 20 years ago. While NT-HR-MS is being applied to increasingly complex research questions, there are still many unidentified chemicals and uncertainties in linking exposures to human health outcomes and environmental impacts. In this perspective, we explore the possibilities and challenges involved in using cheminformatics and NT-HR-MS to answer complex questions that cross many scientific disciplines, taking the identification of potential (small molecule) neurotoxicants in environmental or biological matrices as a case study. We explore capturing literature knowledge and patient exposure information in a form amenable to high-throughput data mining, and the related cheminformatic challenges. We then briefly cover which sample matrices are available, which method(s) could potentially be used to detect these chemicals in various matrices and what remains beyond the reach of NT-HR-MS. We touch on the potential for biological validation systems to contribute to mechanistic understanding of observations and explore which sampling and data archiving strategies may be required to form an accurate, sustained picture of small molecule signatures on extensive cohorts of patients with chronic neurodegenerative disorders. Finally, we reflect on how NT-HR-MS can support unravelling the contribution of the environment to complex diseases.

Genes associated with Parkinson's disease respond to increasing polychlorinated biphenyl levels in the blood of healthy females

Polychlorinated biphenyls (PCBs) are a class of widespread environmental pollutants, commonly found in human blood, that have been suggested to be linked to the occurrence of sporadic Parkinson's disease (PD). It has been reported that some non-coplanar PCBs accumulate in the brains of female PD patients. To improve our understanding of the association between PCB exposure and PD risk we have applied whole transcriptome gene expression analysis in blood cells from 594 PCB-exposed subjects (369 female, 225 male). Interestingly, we observe that in females, blood levels of non-coplanar PCBs appear to be associated with expression levels of PD-specific genes. However, no such association was detected in males. Among the 131 PD-specific genes affected, 39 have been shown to display similar changes in expression levels in the substantia nigra of deceased PD patients. Especially among the down-regulated genes, transcripts of genes involved in neurotransmitter vesicle-related functions were predominant.

Presenting and sharing clinical data using the eTRIKS Standards Master Tree for tranSMART

Motivation

Standardization and semantic alignment have been considered one of the major challenges for data integration in clinical research. The inclusion of the CDISC SDTM clinical data standard into the tranSMART i2b2 via a guiding master ontology tree positively impacts and supports the efficacy of data sharing, visualization and exploration across datasets.

Results

We present here a schema for the organization of SDTM variables into the tranSMART i2b2 tree along with a script and test dataset to exemplify the mapping strategy. The eTRIKS master tree concept is demonstrated by making use of fictitious data generated for four patients, including 16 SDTM clinical domains. We describe how the usage of correct visit names and data labels can help to integrate multiple readouts per patient and avoid ETL crashes when running a tranSMART loading routine.

Availability and implementation

The eTRIKS Master Tree package and test datasets are publicly available at https://doi.org/10.5281/zenodo.1009098 and a functional demo installation at https://public.etriks.org/transmart/datasetExplorer/ under eTRIKS—Master Tree branch, where the discussed examples can be visualized.

MIC-MAC: An automated pipeline for high-throughput characterization and classification of three-dimensional microglia morphologies in mouse and human postmortem brain samples

The phenotypic changes of microglia in brain diseases are particularly diverse and their role in disease progression, beneficial, or detrimental, is still elusive. High‐throughput molecular approaches such as single‐cell RNA‐sequencing can now resolve the high heterogeneity in microglia population for a specific physiological condition, however, the relation between the different microglial signatures and their surrounding brain microenvironment is barely understood. Thus, better tools to characterize the phenotypic variations of microglia in situ are needed, particularly for human brain postmortem samples analysis. To address this challenge, we developed MIC‐MAC, a Microglia and Immune Cells Morphologies Analyser and Classifier pipeline that semiautomatically segments, extracts, and classifies all microglia and immune cells labeled in large three‐dimensional (3D) confocal image stacks of mouse and human brain samples. Our imaging‐based approach enables automatic 3D‐morphology characterization and classification of thousands of individual microglia in situ and revealed species‐ and disease‐specific morphological phenotypes in mouse aging, human Alzheimer's disease, and dementia with Lewy Bodie's samples. MIC‐MAC is a precision diagnostic tool that allows a rapid, unbiased, and large‐scale analysis of microglia morphological states in mouse models and patient brain samples.

Community-driven roadmap for integrated disease maps

The Disease Maps Project builds on a network of scientific and clinical groups that exchange best practices, share information and develop systems biomedicine tools. The project aims for an integrated, highly curated and user-friendly platform for disease-related knowledge. The primary focus of disease maps is on interconnected signaling, metabolic and gene regulatory network pathways represented in standard formats. The involvement of domain experts ensures that the key disease hallmarks are covered and relevant, up-to-date knowledge is adequately represented. Expert-curated and computer readable, disease maps may serve as a compendium of knowledge, allow for data-supported hypothesis generation or serve as a scaffold for the generation of predictive mathematical models. This article summarizes the 2nd Disease Maps Community meeting, highlighting its important topics and outcomes. We outline milestones on the roadmap for the future development of disease maps, including creating and maintaining standardized disease maps; sharing parts of maps that encode common human disease mechanisms; providing technical solutions for complexity management of maps; and Web tools for in-depth exploration of such maps. A dedicated discussion was focused on mathematical modeling approaches, as one of the main goals of disease map development is the generation of mathematically interpretable representations to predict disease comorbidity or drug response and to suggest drug repositioning, altogether supporting clinical decisions.

Large-scale validation of miRNAs by disease association, evolutionary conservation and pathway activity

The validation of microRNAs (miRNAs) identified by next generation sequencing involves amplification-free and hybridization-based detection of transcripts as criteria for confirming valid miRNAs. Since respective validation is frequently not performed, miRNA repositories likely still contain a substantial fraction of false positive candidates while true miRNAs are not stored in the repositories yet. Especially if downstream analyses are performed with these candidates (e.g. target or pathway prediction), the results may be misleading. In the present study, we evaluated 558 mature miRNAs from miRBase and 1,709 miRNA candidates from next generation sequencing experiments by amplification-free hybridization and investigated their distributions in patients with various disease conditions. Notably, the most significant miRNAs in diseases are often not contained in the miRBase. However, these candidates are evolutionary highly conserved. From the expression patterns, target gene and pathway analyses and evolutionary conservation analyses, we were able to shed light on the complexity of miRNAs in humans. Our data also highlight that a more thorough validation of miRNAs identified by next generation sequencing is required. The results are available in miRCarta ( https://mircarta.cs.uni-saarland.de ).

The Luxembourg Parkinson's Study: A Comprehensive Approach for Stratification and Early Diagnosis

While genetic advances have successfully defined part of the complexity in Parkinson's disease (PD), the clinical characterization of phenotypes remains challenging. Therapeutic trials and cohort studies typically include patients with earlier disease stages and exclude comorbidities, thus ignoring a substantial part of the real-world PD population. To account for these limitations, we implemented the Luxembourg PD study as a comprehensive clinical, molecular and device-based approach including patients with typical PD and atypical parkinsonism, irrespective of their disease stage, age, comorbidities, or linguistic background. To provide a large, longitudinally followed, and deeply phenotyped set of patients and controls for clinical and fundamental research on PD, we implemented an open-source digital platform that can be harmonized with international PD cohort studies. Our interests also reflect Luxembourg-specific areas of PD research, including vision, gait, and cognition. This effort is flanked by comprehensive biosampling efforts assuring high quality and sustained availability of body liquids and tissue biopsies. We provide evidence for the feasibility of such a cohort program with deep phenotyping and high quality biosampling on parkinsonism in an environment with structural specificities and alert the international research community to our willingness to collaborate with other centers. The combination of advanced clinical phenotyping approaches including device-based assessment will create a comprehensive assessment of the disease and its variants, its interaction with comorbidities and its progression. We envision the Luxembourg Parkinson's study as an important research platform for defining early diagnosis and progression markers that translate into stratified treatment approaches.

Single-cell transcriptomics reveals distinct inflammation-induced microglia signatures

Microglia are specialized parenchymal‐resident phagocytes of the central nervous system (CNS) that actively support, defend and modulate the neural environment. Dysfunctional microglial responses are thought to worsen CNS diseases; nevertheless, their impact during neuroinflammatory processes remains largely obscure. Here, using a combination of single‐cell RNA sequencing and multicolour flow cytometry, we comprehensively profile microglia in the brain of lipopolysaccharide (LPS)‐injected mice. By excluding the contribution of other immune CNS‐resident and peripheral cells, we show that microglia isolated from LPS‐injected mice display a global downregulation of their homeostatic signature together with an upregulation of inflammatory genes. Notably, we identify distinct microglial activated profiles under inflammatory conditions, which greatly differ from neurodegenerative disease‐associated profiles. These results provide insights into microglial heterogeneity and establish a resource for the identification of specific phenotypes in CNS disorders, such as neuroinflammatory and neurodegenerative diseases.

From hype to reality: data science enabling personalized medicine

BACKGROUND

Personalized, precision, P4, or stratified medicine is understood as a medical approach in which patients are stratified based on their disease subtype, risk, prognosis, or treatment response using specialized diagnostic tests. The key idea is to base medical decisions on individual patient characteristics, including molecular and behavioral biomarkers, rather than on population averages. Personalized medicine is deeply connected to and dependent on data science, specifically machine learning (often named Artificial Intelligence in the mainstream media). While during recent years there has been a lot of enthusiasm about the potential of 'big data' and machine learning-based solutions, there exist only few examples that impact current clinical practice. The lack of impact on clinical practice can largely be attributed to insufficient performance of predictive models, difficulties to interpret complex model predictions, and lack of validation via prospective clinical trials that demonstrate a clear benefit compared to the standard of care. In this paper, we review the potential of state-of-the-art data science approaches for personalized medicine, discuss open challenges, and highlight directions that may help to overcome them in the future.

CONCLUSIONS

There is a need for an interdisciplinary effort, including data scientists, physicians, patient advocates, regulatory agencies, and health insurance organizations. Partially unrealistic expectations and concerns about data science-based solutions need to be better managed. In parallel, computational methods must advance more to provide direct benefit to clinical practice.

Systems medicine disease maps: community-driven comprehensive representation of disease mechanisms

The development of computational approaches in systems biology has reached a state of maturity that allows their transition to systems medicine. Despite this progress, intuitive visualisation and context-dependent knowledge representation still present a major bottleneck. In this paper, we describe the Disease Maps Project, an effort towards a community-driven computationally readable comprehensive representation of disease mechanisms. We outline the key principles and the framework required for the success of this initiative, including use of best practices, standards and protocols. We apply a modular approach to ensure efficient sharing and reuse of resources for projects dedicated to specific diseases. Community-wide use of disease maps will accelerate the conduct of biomedical research and lead to new disease ontologies defined from mechanism-based disease endotypes rather than phenotypes.

Rare ABCA7 variants in 2 German families with Alzheimer disease

Objective

The aim of this study was to identify variants associated with familial late-onset Alzheimer disease (AD) using whole-genome sequencing.

Methods

Several families with an autosomal dominant inheritance pattern of AD were analyzed by whole-genome sequencing. Variants were prioritized for rare, likely pathogenic variants in genes already known to be associated with AD and confirmed by Sanger sequencing using standard protocols.

Results

We identified 2 rare ABCA7 variants (rs143718918 and rs538591288) with varying penetrance in 2 independent German AD families, respectively. The single nucleotide variant (SNV) rs143718918 causes a missense mutation, and the deletion rs538591288 causes a frameshift mutation of ABCA7. Both variants have previously been reported in larger cohorts but with incomplete segregation information. ABCA7 is one of more than 20 AD risk loci that have so far been identified by genome-wide association studies, and both common and rare variants of ABCA7 have previously been described in different populations with higher frequencies in AD cases than in controls and varying penetrance. Furthermore, ABCA7 is known to be involved in several AD-relevant pathways.

Conclusions

We conclude that both SNVs might contribute to the development of AD in the examined family members. Together with previous findings, our data confirm ABCA7 as one of the most relevant AD risk genes.

A roadmap towards personalized immunology

Big data generation and computational processing will enable medicine to evolve from a "one-size-fits-all" approach to precise patient stratification and treatment. Significant achievements using "Omics" data have been made especially in personalized oncology. However, immune cells relative to tumor cells show a much higher degree of complexity in heterogeneity, dynamics, memory-capability, plasticity and "social" interactions. There is still a long way ahead on translating our capability to identify potentially targetable personalized biomarkers into effective personalized therapy in immune-centralized diseases. Here, we discuss the recent advances and successful applications in "Omics" data utilization and network analysis on patients' samples of clinical trials and studies, as well as the major challenges and strategies towards personalized stratification and treatment for infectious or non-communicable inflammatory diseases such as autoimmune diseases or allergies. We provide a roadmap and highlight experimental, clinical, computational analysis, data management, ethical and regulatory issues to accelerate the implementation of personalized immunology.

Exome-wide analysis of mutational burden in patients with typical and atypical Rolandic epilepsy

Rolandic epilepsy (RE) is the most common focal epilepsy in childhood. To date no hypothesis-free exome-wide mutational screen has been conducted for RE and atypical RE (ARE). Here we report on whole-exome sequencing of 194 unrelated patients with RE/ARE and 567 ethnically matched population controls. We identified an exome-wide significantly enriched burden for deleterious and loss-of-function variants only for the established RE/ARE gene GRIN2A. The statistical significance of the enrichment disappeared after removing ARE patients. For several disease-related gene-sets, an odds ratio >1 was detected for loss-of-function variants.

Commentaries on “Informatics and Medicine: From Molecules to Populations”

Objective: To discuss interdisciplinary research and education in the context of informatics and medicine by commenting on the paper of Kuhn et al. “Informatics and Medicine: From Molecules to Populations”.

Method: Inviting an international group of experts in biomedical and health informatics and related disciplines to comment on this paper.

Results and Conclusions: The commentaries include a wide range of reasoned arguments and original position statements which, while strongly endorsing the educational needs identified by Kuhn et al., also point out fundamental challenges that are very specific to the unusual combination of scientific, technological, personal and social problems characterizing biomedical informatics. They point to the ultimate objectives of managing difficult human health problems, which are unlikely to yield to technological solutions alone. The psychological, societal, and environmental components of health and disease are emphasized by several of the commentators, setting the stage for further debate and constructive suggestions.

Similar α-Synuclein staining in the colon mucosa in patients with Parkinson's disease and controls

BACKGROUND

The gut is proposed as a starting point of idiopathic IPD, but the presence of α-synuclein in the IPD colon mucosa is debated.

OBJECTIVES

The objective of this study was to evaluate if α-synuclein in the colon mucosa can serve as a biomarker of IPD.

METHODS

Immunohistochemistry was used to locate and quantify in a blinded approach α-synuclein in the mucosa from biopsies of the right and left colon in 19 IPD patients and 8 controls.

RESULTS

Total α-synuclein was present in all but 1 IPD patients and in all controls; phosphorylated α-synuclein was present in all subjects. There was no intensity difference depending on disease status. Staining of total α-synuclein was stronger in the right colon (p = .04).

CONCLUSIONS

Conventional immunohistochemistry α-synuclein staining in colon mucosal biopsies cannot serve as a biomarker of idiopathic PD. These findings do not contradict the assumption of disease starting in the colon, and a colon segment-specific risk for disease initiation can still be hypothesized. © 2016 International Parkinson and Movement Disorder Society.

Rare GABRA3 variants are associated with epileptic seizures, encephalopathy and dysmorphic features

Genetic epilepsies are caused by mutations in a range of different genes, many of them encoding ion channels, receptors or transporters. While the number of detected variants and genes increased dramatically in the recent years, pleiotropic effects have also been recognized, revealing that clinical syndromes with various degrees of severity arise from a single gene, a single mutation, or from different mutations showing similar functional defects. Accordingly, several genes coding for GABAA receptor subunits have been linked to a spectrum of benign to severe epileptic disorders and it was shown that a loss of function presents the major correlated pathomechanism. Here, we identified six variants in GABRA3 encoding the α3-subunit of the GABAA receptor. This gene is located on chromosome Xq28 and has not been previously associated with human disease. Five missense variants and one microduplication were detected in four families and two sporadic cases presenting with a range of epileptic seizure types, a varying degree of intellectual disability and developmental delay, sometimes with dysmorphic features or nystagmus. The variants co-segregated mostly but not completely with the phenotype in the families, indicating in some cases incomplete penetrance, involvement of other genes, or presence of phenocopies. Overall, males were more severely affected and there were three asymptomatic female mutation carriers compared to only one male without a clinical phenotype. X-chromosome inactivation studies could not explain the phenotypic variability in females. Three detected missense variants are localized in the extracellular GABA-binding NH2-terminus, one in the M2-M3 linker and one in the M4 transmembrane segment of the α3-subunit. Functional studies in Xenopus laevis oocytes revealed a variable but significant reduction of GABA-evoked anion currents for all mutants compared to wild-type receptors. The degree of current reduction correlated partially with the phenotype. The microduplication disrupted GABRA3 expression in fibroblasts of the affected patient. In summary, our results reveal that rare loss-of-function variants in GABRA3 increase the risk for a varying combination of epilepsy, intellectual disability/developmental delay and dysmorphic features, presenting in some pedigrees with an X-linked inheritance pattern.

Recessive mutations in SLC35A3 cause early onset epileptic encephalopathy with skeletal defects

We describe the clinical and whole genome sequencing (WGS) study of a non-consanguineous Italian family in which two siblings, a boy and a girl, manifesting a severe epileptic encephalopathy (EE) with skeletal abnormalities, carried novel SLC35A3 compound heterozygous mutations. Both siblings exhibited infantile spasms, associated with focal, and tonic vibratory seizures from early infancy. EEG recordings showed a suppression-burst (SB) pattern and multifocal paroxysmal activity in both. In addition both had quadriplegia, acquired microcephaly, and severe intellectual disability. General examination showed distal arthrogryposis predominant in the hands in both siblings and severe left dorso-lumbar convex scoliosis in one. WGS of the siblings-parents quartet identified novel compound heterozygous mutations in SLC35A3 in both children. SLC35A3 encodes the major Golgi uridine diphosphate N-acetylglucosamine transporter. With this study, we add SLC35A3 to the gene list of epilepsies. Neurological symptoms and skeletal abnormalities might result from impaired glycosylation of proteins involved in normal development and function of the central nervous system and skeletal apparatus.

Embryonic development of selectively vulnerable neurons in Parkinson's disease

A specific set of brainstem nuclei are susceptible to degeneration in Parkinson's disease. We hypothesise that neuronal vulnerability reflects shared phenotypic characteristics that confer selective vulnerability to degeneration. Neuronal phenotypic specification is mainly the cumulative result of a transcriptional regulatory program that is active during the development. By manual curation of the developmental biology literature, we comprehensively reconstructed an anatomically resolved cellular developmental lineage for the adult neurons in five brainstem regions that are selectively vulnerable to degeneration in prodromal or early Parkinson's disease. We synthesised the literature on transcription factors that are required to be active, or required to be inactive, in the development of each of these five brainstem regions, and at least two differentially vulnerable nuclei within each region. Certain transcription factors, e.g., Ascl1 and Lmx1b, seem to be required for specification of many brainstem regions that are susceptible to degeneration in early Parkinson's disease. Some transcription factors can even distinguish between differentially vulnerable nuclei within the same brain region, e.g., Pitx3 is required for specification of the substantia nigra pars compacta, but not the ventral tegmental area. We do not suggest that Parkinson's disease is a developmental disorder. In contrast, we consider identification of shared developmental trajectories as part of a broader effort to identify the molecular mechanisms that underlie the phenotypic features that are shared by selectively vulnerable neurons. Systematic in vivo assessment of fate determining transcription factors should be completed for all neuronal populations vulnerable to degeneration in early Parkinson's disease.

Biallelic Variants in OTUD6B Cause an Intellectual Disability Syndrome Associated with Seizures and Dysmorphic Features

Ubiquitination is a posttranslational modification that regulates many cellular processes including protein degradation, intracellular trafficking, cell signaling, and protein-protein interactions. Deubiquitinating enzymes (DUBs), which reverse the process of ubiquitination, are important regulators of the ubiquitin system. OTUD6B encodes a member of the ovarian tumor domain (OTU)-containing subfamily of deubiquitinating enzymes. Herein, we report biallelic pathogenic variants in OTUD6B in 12 individuals from 6 independent families with an intellectual disability syndrome associated with seizures and dysmorphic features. In subjects with predicted loss-of-function alleles, additional features include global developmental delay, microcephaly, absent speech, hypotonia, growth retardation with prenatal onset, feeding difficulties, structural brain abnormalities, congenital malformations including congenital heart disease, and musculoskeletal features. Homozygous Otud6b knockout mice were subviable, smaller in size, and had congenital heart defects, consistent with the severity of loss-of-function variants in humans. Analysis of peripheral blood mononuclear cells from an affected subject showed reduced incorporation of 19S subunits into 26S proteasomes, decreased chymotrypsin-like activity, and accumulation of ubiquitin-protein conjugates. Our findings suggest a role for OTUD6B in proteasome function, establish that defective OTUD6B function underlies a multisystemic human disorder, and provide additional evidence for the emerging relationship between the ubiquitin system and human disease.