Integrating biomedical research and electronic health records to create knowledge-based biologically meaningful machine-readable embeddings.

In order to advance precision medicine, detailed clinical features ought to be described in a way that leverages current knowledge. Although data collected from biomedical research is expanding at an almost exponential rate, our ability to transform that information into patient care has not kept at pace. A major barrier preventing this transformation is that multi-dimensional data collection and analysis is usually carried out without much understanding of the underlying knowledge structure. Here, in an effort to bridge this gap, Electronic Health Records (EHRs) of individual patients are connected to a heterogeneous knowledge network called Scalable Precision Medicine Oriented Knowledge Engine (SPOKE). Then an unsupervised machine-learning algorithm creates Propagated SPOKE Entry Vectors (PSEVs) that encode the importance of each SPOKE node for any code in the EHRs. We argue that these results, alongside the natural integration of PSEVs into any EHR machine-learning platform, provide a key step toward precision medicine

Peroxisome proliferator-activated receptor (PPAR)alpha expression in T cells mediates gender differences in development of T cell-mediated autoimmunity.

Peroxisome proliferator-activated receptor (PPAR)alpha is a nuclear receptor that mediates gender differences in lipid metabolism. PPARalpha also functions to control inflammatory responses by repressing the activity of nuclear factor kappaB (NF-kappaB) and c-jun in immune cells. Because PPARalpha is situated at the crossroads of gender and immune regulation, we hypothesized that this gene may mediate sex differences in the development of T cell-mediated autoimmune disease. We show that PPARalpha is more abundant in male as compared with female CD4(+) cells and that its expression is sensitive to androgen levels. Genetic ablation of this gene selectively removed the brake on NF-kappaB and c-jun activity in male T lymphocytes, resulting in higher production of interferon gamma and tumor necrosis factor (but not interleukin 17), and lower production of T helper (Th)2 cytokines. Upon induction of experimental autoimmune encephalomyelitis, male but not female PPARalpha(-/-) mice developed more severe clinical signs that were restricted to the acute phase of disease. These results suggest that males are less prone to develop Th1-mediated autoimmunity because they have higher T cell expression of PPARalpha

Silent progression in disease activity-free relapsing multiple sclerosis.

ObjectiveRates of worsening and evolution to secondary progressive multiple sclerosis (MS) may be substantially lower in actively treated patients compared to natural history studies from the pretreatment era. Nonetheless, in our recently reported prospective cohort, more than half of patients with relapsing MS accumulated significant new disability by the 10th year of follow-up. Notably, "no evidence of disease activity" at 2 years did not predict long-term stability. Here, we determined to what extent clinical relapses and radiographic evidence of disease activity contribute to long-term disability accumulation.MethodsDisability progression was defined as an increase in Expanded Disability Status Scale (EDSS) of 1.5, 1.0, or 0.5 (or greater) from baseline EDSS = 0, 1.0-5.0, and 5.5 or higher, respectively, assessed from baseline to year 5 (±1 year) and sustained to year 10 (±1 year). Longitudinal analysis of relative brain volume loss used a linear mixed model with sex, age, disease duration, and HLA-DRB1*15:01 as covariates.ResultsRelapses were associated with a transient increase in disability over 1-year intervals (p = 0.012) but not with confirmed disability progression (p = 0.551). Relative brain volume declined at a greater rate among individuals with disability progression compared to those who remained stable (p < 0.05).InterpretationLong-term worsening is common in relapsing MS patients, is largely independent of relapse activity, and is associated with accelerated brain atrophy. We propose the term silent progression to describe the insidious disability that accrues in many patients who satisfy traditional criteria for relapsing-remitting MS. Ann Neurol 2019;85:653-666

Heterogeneity at the HLA-DRB1 locus and risk for multiple sclerosis.

Variation in major histocompatibility complex genes on chromosome 6p21.3, specifically the human leukocyte antigen HLA-DR2 or DRB1*1501-DQB1*0602 extended haplotype, confers risk for multiple sclerosis (MS). Previous studies of DRB1 variation and both MS susceptibility and phenotypic expression have lacked statistical power to detect modest genotypic influences, and have demonstrated conflicting results. Results derived from analyses of 1339 MS families indicate DRB1 variation influences MS susceptibility in a complex manner. DRB1*15 was strongly associated in families (P=7.8x10(-31)), and a dominant DRB1*15 dose effect was confirmed (OR=7.5, 95% CI=4.4-13.0, P<0.0001). A modest dose effect was also detected for DRB1*03; however, in contrast to DRB1*15, this risk was recessive (OR=1.8, 95% CI=1.1-2.9, P=0.03). Strong evidence for under-transmission of DRB1*14 (P=5.7x10(-6)) even after accounting for DRB1*15 (P=0.03) was present, confirming a protective effect. In addition, a high risk DRB1*15 genotype bearing DRB1*08 was identified (OR=7.7, 95% CI=4.1-14.4, P<0.0001), providing additional evidence for trans DRB1 allelic interactions in MS. Further, a significant DRB1*15 association observed in primary progressive MS families (P=0.0004), similar to relapsing-remitting MS families, suggests that DRB1-related mechanisms are contributing to both phenotypes. In contrast, results obtained from 2201 MS cases argue convincingly that DRB1*15 genotypes do not modulate age of onset, or significantly influence disease severity measured using expanded disease disability score and disease duration. These results contribute substantially to our understanding of the DRB1 locus and MS, and underscore the importance of using large sample sizes to detect modest genetic effects, particularly in studies of genotype-phenotype relationships

Infectious Disease Ontology

Technological developments have resulted in tremendous increases in the volume and diversity of the data and information that must be processed in the course of biomedical and clinical research and practice. Researchers are at the same time under ever greater pressure to share data and to take steps to ensure that data resources are interoperable. The use of ontologies to annotate data has proven successful in supporting these goals and in providing new possibilities for the automated processing of data and information. In this chapter, we describe different types of vocabulary resources and emphasize those features of formal ontologies that make them most useful for computational applications. We describe current uses of ontologies and discuss future goals for ontology-based computing, focusing on its use in the field of infectious diseases. We review the largest and most widely used vocabulary resources relevant to the study of infectious diseases and conclude with a description of the Infectious Disease Ontology (IDO) suite of interoperable ontology modules that together cover the entire infectious disease domain

Pathway and network-based analysis of genome-wide association studies in multiple sclerosis

Genome-wide association studies (GWAS) testing several hundred thousand SNPs have been performed in multiple sclerosis (MS) and other complex diseases. Typically, the number of markers in which the evidence for association exceeds the genome-wide significance threshold is very small, and markers that do not exceed this threshold are generally neglected. Classical statistical analysis of these datasets in MS revealed genes with known immunological functions. However, many of the markers showing modest association may represent false negatives. We hypothesize that certain combinations of genes flagged by these markers can be identified if they belong to a common biological pathway. Here we conduct a pathway-oriented analysis of two GWAS in MS that takes into account all SNPs with nominal evidence of association (P < 0.05). Gene-wise P-values were superimposed on a human protein interaction network and searches were conducted to identify sub-networks containing a higher proportion of genes associated with MS than expected by chance. These sub-networks, and others generated at random as a control, were categorized for membership of biological pathways. GWAS from eight other diseases were analyzed to assess the specificity of the pathways identified. In the MS datasets, we identified sub-networks of genes from several immunological pathways including cell adhesion, communication and signaling. Remarkably, neural pathways, namely axon-guidance and synaptic potentiation, were also over-represented in MS. In addition to the immunological pathways previously identified, we report here for the first time the potential involvement of neural pathways in MS susceptibilit

Genotype-Phenotype correlations in multiple sclerosis: HLA genes influence disease severity inferred by 1HMR spectroscopy and MRI measures

Genetic susceptibility to multiple sclerosis (MS) is associated with the human leukocyte antigen (HLA) DRB1*1501 allele. Here we show a clear association between DRB1*1501 carrier status and four domains of disease severity in an investigation of genotype-phenotype associations in 505 robust, clinically well characterized MS patients evaluated cross-sectionally: (i) a reduction in the N-acetyl-aspartate (NAA) concentration within normal appearing white matter (NAWM) via 1HMR spectroscopy (P = 0.025), (ii) an increase in the volume of white matter (WM) lesions utilizing conventional anatomical MRI techniques (1,127 mm3; P = 0.031), (iii) a reduction in normalized brain parenchymal volume (nBPV) (P = 0.023), and (iv) impairments in cognitive function as measured by the Paced Auditory Serial Addition Test (PASAT-3) performance (Mean Z Score: DRB1*1501+: 0.110 versus DRB1*1501-: 0.048; P = 0.004). In addition, DRB1*1501+ patients had significantly more women (74% versus 63%; P = 0.009) and a younger mean age at disease onset (32.4 years versus 34.3 years; P = 0.025). Our findings suggest that DRB1*1501 increases disease severity in MS by facilitating the development of more T2-foci, thereby increasing the potential for irreversible axonal compromise and subsequent neuronal degeneration, as suggested by the reduction of NAA concentrations in NAWM, ultimately leading to a decline in brain volume. These structural aberrations may explain the significant differences in cognitive performance observed between DRB1*1501 groups. The overall goal of a deep phenotypic approach to MS is to develop an array of meaningful biomarkers to monitor the course of the disease, predict future disease behaviour, determine when treatment is necessary, and perhaps to more effectively recommend an available therapeutic interventio

Modeling the Cumulative Genetic Risk for Multiple Sclerosis from Genome-Wide Association Data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics.Version of Recor

Spatio-Temporal Changes of Extracellular Matrix (ECM) Stiffness in the Development of the Leech Hirudo verbana

The invertebrate leech Hirudo verbana represents a powerful experimental animal model for improving the knowledge about the functional interaction between the extracellular matrix (ECM) and cells within the tissue microenvironment (TME), and the key role played by ECM stiffness during development and growth. Indeed, the medicinal leech is characterized by a simple anatomical organization reproducing many aspects of the basic biological processes of vertebrates and in which a rapid spatiotemporal development is well established and easily assessed. Our results show that ECM structural organization, as well as the amount of fibrillar and non-fibrillar collagen are deeply different from hatching leeches to adult ones. In addition, the changes in ECM remodelling occurring during the different leech developmental stages, leads to a gradient of stiffness regulating both the path of migratory cells and their fates. The ability of cells to perceive and respond to changes in ECM composition and mechanics strictly depend on nuclear or cytoplasmic expression of Yes-Associated Protein 1 (YAP1), a key mediator converting mechanical signals into transcriptional outputs, expression, and activation

Dynamic relationship among extracellular matrix and body wall cells in Hirudo verbana morphogenesis

A great bulk of recent experimental evidence suggests the key role of the complex crosstalk between the extracellular matrix (ECM) and the cellular component of tissues during morphogenesis and embryogenesis. In particular, remodeling of the ECM and of its physical interactions pattern with surrounding cells represent two crucial processes that might be involved in muscle development. However, little information is available on this topic, especially on invertebrate species. To obtain new insights on how tuning the ECM microenvironment might drive cellular fate during embryonic development, we used the invertebrate medicinal leech Hirudo verbana as a valuable experimental model, due to its simple anatomy and the recapitulation of many aspects of the basic biological processes of vertebrates. Our previous studies on leech post-embryonic development have already shown the pivotal role of ECM changes during the growth of the body wall and the role of Yes-associated protein 1 (YAP1) in mechanotransduction. Here, we suggest that the interactions between stromal cell telocytes and ECM might be crucial in driving the organization of muscle layers during embryogenesis. Furthermore, we propose a possible role of the pleiotropic enzyme HvRNASET2 as a possible modulator of collagen deposition and ECM remodeling not only during regenerative processes (as previously demonstrated) but also in embryogenesis