Stratification of Kidney Transplant Recipients Into Five Subgroups Based on Temporal Disease Trajectories.

Background: Kidney transplantation is the treatment of choice for patients with end-stage renal disease. Considerable clinical research has focused on improving graft survival and an increasing number of kidney recipients die with a functioning graft. There is a need to improve patient survival and to better understand the individualized risk of comorbidities and complications. Here, we developed a method to stratify recipients into similar subgroups based on previous comorbidities and subsequently identify complications and for a subpopulation, laboratory test values associated with survival. Methods: First, we identified significant disease patterns based on all hospital diagnoses from the Danish National Patient Registry for 5752 kidney transplant recipients from 1977 to 2018. Using hierarchical clustering, these longitudinal patterns of diseases segregate into 3 main clusters of glomerulonephritis, hypertension, and diabetes. As some recipients are diagnosed with diseases from >1 cluster, recipients are further stratified into 5 more fine-grained trajectory subgroups for which survival, stratified complication patterns as well as laboratory test values are analyzed. Results: The study replicated known associations indicating that diabetes and low levels of albumin are associated with worse survival when investigating all recipients. However, stratification of recipients by trajectory subgroup showed additional associations. For recipients with glomerulonephritis, higher levels of basophils are significantly associated with poor survival, and these patients are more often diagnosed with bacterial infections. Additional associations were also found. Conclusions: This study demonstrates that disease trajectories can confirm known comorbidities and furthermore stratify kidney transplant recipients into clinical subgroups in which we can characterize stratified risk factors. We hope to motivate future studies to stratify recipients into more fine-grained, homogenous subgroups to better discover associations relevant for the individual patient and thereby enable more personalized disease-management and improve long-term outcomes and survival.

Population-wide analysis of hospital laboratory tests to assess seasonal variation and temporal reference interval modification.

We identified mortality-, age-, and sex-associated differences in relation to reference intervals (RIs) for laboratory tests in population-wide data from nearly 2 million hospital patients in Denmark and comprising more than 300 million measurements. A low-parameter mathematical wave-based modification method was developed to adjust for dietary and environment influences during the year. The resulting mathematical fit allowed for improved association rates between re-classified abnormal laboratory tests, patient diagnoses, and mortality. The study highlights the need for seasonally modified RIs and presents an approach that has the potential to reduce over- and underdiagnosis, affecting both physician-patient interactions and electronic health record research as a whole.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models.

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Integrate mechanistic evidence from new approach methodologies (NAMs) into a read-across assessment to characterise trends in shared mode of action.

Read-across approaches often remain inconclusive as they do not provide sufficient evidence on a common mode of action across the category members. This read-across case study on thirteen, structurally similar, branched aliphatic carboxylic acids investigates the concept of using human-based new approach methods, such as in vitro and in silico models, to demonstrate biological similarity. Five out of the thirteen analogues have preclinical in vivo studies. Three out of them induced lipid accumulation or hypertrophy in preclinical studies with repeated exposure, which leads to the read-across hypothesis that the analogues can potentially induce hepatic steatosis. To confirm the selection of analogues, the expression patterns of the induced differentially expressed genes (DEGs) were analysed in a human liver model. With increasing dose, the expression pattern within the tested analogues got more similar, which serves as a first indication of a common mode of action and suggests differences in the potency of the analogues. Hepatic steatosis is a well-known adverse outcome, for which over 55 adverse outcome pathways have been identified. The resulting adverse outcome pathway (AOP) network, comprised a total 43 MIEs/KEs and enabled the design of an in vitro testing battery. From the AOP network, ten MIEs, early and late KEs were tested to systematically investigate a common mode of action among the grouped compounds. The targeted testing of AOP specific MIE/KEs shows that biological activity in the category decreases with side chain length. A similar trend was evident in measuring liver alterations in zebra fish embryos. However, activation of single MIEs or early KEs at in vivo relevant doses did not necessarily progress to the late KE "lipid accumulation". KEs not related to the read-across hypothesis, testing for example general mitochondrial stress responses in liver cells, showed no trend or biological similarity. Testing scope is a key issue in the design of in vitro test batteries. The Dempster-Shafer decision theory predicted those analogues with in vivo reference data correctly using one human liver model or the CALUX reporter assays. The case study shows that the read-across hypothesis is the key element to designing the testing strategy. In the case of a good mechanistic understanding, an AOP facilitates the selection of reliable human in vitro models to demonstrate a common mode of action. Testing DEGs, MIEs and early KEs served to show biological similarity, whereas the late KEs become important for confirmation, as progression from MIEs to AO is not always guaranteed.

Optimizing drug selection from a prescription trajectory of one patient.

It is unknown how sequential drug patterns convey information on a patient's health status and treatment guidelines rarely account for this. Drug-agnostic longitudinal analyses of prescription trajectories in a population-wide setting are needed. In this cohort study, we used 24 years of data (1.1 billion prescriptions) from the Danish prescription registry to model the risk of sequentially redeeming a drug after another. Drug pairs were used to build multistep longitudinal prescription trajectories. These were subsequently used to stratify patients and calculate survival hazard ratios between the stratified groups. The similarity between prescription histories was used to determine individuals' best treatment option. Over the course of 122 million person-years of observation, we identified 9 million common prescription trajectories and demonstrated their predictive power using hypertension as a case. Among patients treated with agents acting on the renin-angiotensin system we identified four groups: patients prescribed angiotensin converting enzyme (ACE) inhibitor without change, angiotensin receptor blockers (ARBs) without change, ACE with posterior change to ARB, and ARB posteriorly changed to ACE. In an adjusted time-to-event analysis, individuals treated with ACE compared to those treated with ARB had lower survival probability (hazard ratio, 0.73 [95% CI, 0.64-0.82]; P < 1 × 10-16). Replication in UK Biobank data showed the same trends. Prescription trajectories can provide novel insights into how individuals' drug use change over time, identify suboptimal or futile prescriptions and suggest initial treatments different from first line therapies. Observations of this kind may also be important when updating treatment guidelines.

Impaired Vitamin D Signaling in T Cells From a Family With Hereditary Vitamin D Resistant Rickets.

The active form of vitamin D, 1,25-dihydroxyvitamin D3 (1,25(OH)2D3), mediates its immunomodulatory effects by binding to the vitamin D receptor (VDR). Here, we describe a new point mutation in the DNA-binding domain of the VDR and its consequences for 1,25(OH)2D3 signaling in T cells from heterozygous and homozygous carriers of the mutation. The mutation did not affect the overall structure or the ability of the VDR to bind 1,25(OH)2D3 and the retinoid X receptor. However, the subcellular localization of the VDR was strongly affected and the transcriptional activity was abolished by the mutation. In heterozygous carriers of the mutation, 1,25(OH)2D3-induced gene regulation was reduced by ~ 50% indicating that the expression level of wild-type VDR determines 1,25(OH)2D3 responsiveness in T cells. We show that vitamin D-mediated suppression of vitamin A-induced gene regulation depends on an intact ability of the VDR to bind DNA. Furthermore, we demonstrate that vitamin A inhibits 1,25(OH)2D3-induced translocation of the VDR to the nucleus and 1,25(OH)2D3-induced up-regulation of CYP24A1. Taken together, this study unravels novel aspects of vitamin D signaling and function of the VDR in human T cells.

Disease trajectory browser for exploring temporal, population-wide disease progression patterns in 7.2 million Danish patients.

We present the Danish Disease Trajectory Browser (DTB), a tool for exploring almost 25 years of data from the Danish National Patient Register. In the dataset comprising 7.2 million patients and 122 million admissions, users can identify diagnosis pairs with statistically significant directionality and combine them to linear disease trajectories. Users can search for one or more disease codes (ICD-10 classification) and explore disease progression patterns via an array of functionalities. For example, a set of linear trajectories can be merged into a disease trajectory network displaying the entire multimorbidity spectrum of a disease in a single connected graph. Using data from the Danish Register for Causes of Death mortality is also included. The tool is disease-agnostic across both rare and common diseases and is showcased by exploring multimorbidity in Down syndrome (ICD-10 code Q90) and hypertension (ICD-10 code I10). Finally, we show how search results can be customized and exported from the browser in a format of choice (i.e. JSON, PNG, JPEG and CSV).

Age-stratified longitudinal study of Alzheimer's and vascular dementia patients.

INTRODUCTION: Similar symptoms, comorbidities and suboptimal diagnostic tests make the distinction between different types of dementia difficult, although this is essential for improved work-up and treatment optimization. METHODS: We calculated temporal disease trajectories of earlier multi-morbidities in Alzheimer's disease (AD) dementia and vascular dementia (VaD) patients using the Danish National Patient Registry covering all hospital encounters in Denmark (1994 to 2016). Subsequently, we reduced the comorbidity space dimensionality using a non-linear technique, uniform manifold approximation and projection. RESULTS: We found 49,112 and 24,101 patients that were diagnosed with AD or VaD, respectively. Temporal disease trajectories showed very similar disease patterns before the dementia diagnosis. Stratifying patients by age and reducing the comorbidity space to two dimensions, showed better discrimination between AD and VaD patients in early-onset dementia. DISCUSSION: Similar age-associated comorbidities, the phenomenon of mixed dementia, and misdiagnosis create great challenges in discriminating between classical subtypes of dementia.

A Comparative Assessment Study of Known Small-Molecule Keap1-Nrf2 Protein-Protein Interaction Inhibitors: Chemical Synthesis, Binding Properties, and Cellular Activity.

Inhibiting the protein-protein interaction (PPI) between the transcription factor Nrf2 and its repressor protein Keap1 has emerged as a promising strategy to target oxidative stress in diseases, including central nervous system (CNS) disorders. Numerous non-covalent small-molecule Keap1-Nrf2 PPI inhibitors have been reported to date, but many feature suboptimal physicochemical properties for permeating the blood-brain barrier, while others contain problematic structural moieties. Here, we present the first side-by-side assessment of all reported Keap1-Nrf2 PPI inhibitor classes using fluorescence polarization, thermal shift assay, and surface plasmon resonance-and further evaluate the compounds in an NQO1 induction cell assay and in counter tests for nonspecific activities. Surprisingly, half of the compounds were inactive or deviated substantially from reported activities, while we confirm the cross-assay activities for others. Through this study, we have identified the most promising Keap1-Nrf2 inhibitors that can serve as pharmacological probes or starting points for developing CNS-active Keap1 inhibitors.

sAOP: linking chemical stressors to adverse outcomes pathway networks.

MOTIVATION: Adverse outcome pathway (AOP) is a toxicological concept proposed to provide a mechanistic representation of biological perturbation over different layers of biological organization. Although AOPs are by definition chemical-agnostic, many chemical stressors can putatively interfere with one or several AOPs and such information would be relevant for regulatory decision-making. RESULTS: With the recent development of AOPs networks aiming to facilitate the identification of interactions among AOPs, we developed a stressor-AOP network (sAOP). Using the 'cytotoxitiy burst' (CTB) approach, we mapped bioactive compounds from the ToxCast data to a list of AOPs reported in AOP-Wiki database. With this analysis, a variety of relevant connections between chemicals and AOP components can be identified suggesting multiple effects not observed in the simplified 'one-biological perturbation to one-adverse outcome' model. The results may assist in the prioritization of chemicals to assess risk-based evaluations in the context of human health. AVAILABILITY AND IMPLEMENTATION: sAOP is available at http://saop.cpr.ku.dk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prediction of human drug-induced liver injury (DILI) in relation to oral doses and blood concentrations.

Drug-induced liver injury (DILI) cannot be accurately predicted by animal models. In addition, currently available in vitro methods do not allow for the estimation of hepatotoxic doses or the determination of an acceptable daily intake (ADI). To overcome this limitation, an in vitro/in silico method was established that predicts the risk of human DILI in relation to oral doses and blood concentrations. This method can be used to estimate DILI risk if the maximal blood concentration (Cmax) of the test compound is known. Moreover, an ADI can be estimated even for compounds without information on blood concentrations. To systematically optimize the in vitro system, two novel test performance metrics were introduced, the toxicity separation index (TSI) which quantifies how well a test differentiates between hepatotoxic and non-hepatotoxic compounds, and the toxicity estimation index (TEI) which measures how well hepatotoxic blood concentrations in vivo can be estimated. In vitro test performance was optimized for a training set of 28 compounds, based on TSI and TEI, demonstrating that (1) concentrations where cytotoxicity first becomes evident in vitro (EC10) yielded better metrics than higher toxicity thresholds (EC50); (2) compound incubation for 48 h was better than 24 h, with no further improvement of TSI after 7 days incubation; (3) metrics were moderately improved by adding gene expression to the test battery; (4) evaluation of pharmacokinetic parameters demonstrated that total blood compound concentrations and the 95%-population-based percentile of Cmax were best suited to estimate human toxicity. With a support vector machine-based classifier, using EC10 and Cmax as variables, the cross-validated sensitivity, specificity and accuracy for hepatotoxicity prediction were 100, 88 and 93%, respectively. Concentrations in the culture medium allowed extrapolation to blood concentrations in vivo that are associated with a specific probability of hepatotoxicity and the corresponding oral doses were obtained by reverse modeling. Application of this in vitro/in silico method to the rat hepatotoxicant pulegone resulted in an ADI that was similar to values previously established based on animal experiments. In conclusion, the proposed method links oral doses and blood concentrations of test compounds to the probability of hepatotoxicity.

Toxicogenomics directory of rat hepatotoxicants in vivo and in cultivated hepatocytes.

Transcriptomics is developing into an invaluable tool in toxicology. The aim of this study was, using a transcriptomics approach, to identify genes that respond similar to many different chemicals (including drugs and industrial compounds) in both rat liver in vivo and in cultivated hepatocytes. For this purpose, we analyzed Affymetrix microarray expression data from 162 compounds that were previously tested in a concentration-dependent manner in rat livers in vivo and in rat hepatocytes cultivated in sandwich culture. These data were obtained from the Japanese Toxicogenomics Project (TGP) and North Rhine-Westphalian (NRW) data sets, which represent 138 and 29 compounds, respectively, and have only 5 compounds in common between them. The in vitro gene expression data from the NRW data set were generated in the present study, while TGP is publicly available. For each of the data sets, the overlap between up- or down-regulated genes in vitro and in vivo was identified, and named in vitro-in vivo consensus genes. Interestingly, the in vivo-in vitro consensus genes overlapped to a remarkable extent between both data sets, and were 21-times (upregulated genes) or 12-times (down-regulated genes) enriched compared to random expectation. Finally, the genes in the TGP and NRW overlap were used to identify the upregulated genes with the highest compound coverage, resulting in a seven-gene set of Cyp1a1, Ugt2b1, Cdkn1a, Mdm2, Aldh1a1, Cyp4a3, and Ehhadh. This seven-gene set was then successfully tested with structural analogues of valproic acid that are not present in the TGP and NRW data sets. In conclusion, the seven-gene set identified in the present study responds similarly in vitro and in vivo to a wide range of different chemicals. Despite these promising results with the seven-gene set, transcriptomics with cultivated rat hepatocytes remains a challenge, because in general many genes are up- or downregulated by in vitro culture per se, respond differently to test compounds in vitro and in vivo, and/or show higher variability in the in vitro system compared to the corresponding in vivo data.

Analysis of Time-Series Gene Expression Data to Explore Mechanisms of Chemical-Induced Hepatic Steatosis Toxicity.

Non-alcoholic fatty liver disease (NAFLD) represents a wide spectrum of disease, ranging from simple fatty liver through steatosis with inflammation and necrosis to cirrhosis. One of the most challenging problems in biomedical research and within the chemical industry is to understand the underlying mechanisms of complex disease, and complex adverse outcome pathways (AOPs). Based on a set of 28 steatotic chemicals with gene expression data measured on primary hepatocytes at three times (2, 8, and 24 h) and three doses (low, medium, and high), we identified genes and pathways, defined as molecular initiating events (MIEs) and key events (KEs) of steatosis using a combination of a time series and pathway analyses. Among the genes deregulated by these compounds, the study highlighted OSBPL9, ALDH7A1, MYADM, SLC51B, PRDX6, GPAT3, TMEM135, DLGDA5, BCO2, APO10LA, TSPAN6, NEURL1B, and DUSP1. Furthermore, pathway analysis indicated deregulation of pathways related to lipid accumulation, such as fat digestion and absorption, linoleic and linolenic acid metabolism, calcium signaling pathway, fatty acid metabolism, peroxisome, retinol metabolism, and steroid metabolic pathways in a time dependent manner. Such transcription profile analysis can help in the understanding of the steatosis evolution over time generated by chemical exposure.

In Silico Systems Pharmacology to Assess Drug's Therapeutic and Toxic Effects.

For many years, the "one target, one drug" paradigm has been the driving force behind developments in pharmaceutical research. With the recent advances in molecular biology and genomics technologies, the focus is shifting toward "drug-holistic" systems based approaches (i.e. systems pharmacology). The integration of large and diverse amount of data from chemistry and biology coupled with the development and the application of network-based approaches to cope with these data is the next paradigm of drug discovery. Systems pharmacology offers a novel way of approaching drug discovery by developing models that consider the global physiological environment of protein targets and their modification by drugs. Studying drug action across multiple scales of complexity from molecular and cellular to tissue and organism levels may help identify new druggable disease genes and to design new drugs with a better efficacy and clinical safety.