Search

Abstract Quantifying the contribution of genetics and environmental effects on disease initiation and progression, as well as the shared genetics of different diseases, is vital for the understanding of the disease etiology of multimorbidities. In this study, we leverage nationwide Danish registries to provide a granular atlas of the genetic origin of disease phenotypes for a cohort of all Danes 1978–2018 with partially known pedigree (n = 6.3 million). We estimate the heritability and genetic correlation between thousands of disease phenotypes using a novel approach that can be scaled to nationwide data. Our findings confirm the importance of genetics for a number of known associations and increase the resolution of heritability by adding numerous associations, some of which point to shared biologically origin of different phenotypes. We also establish the heritability of disease trajectories and the importance of sex-specific genetic contributions. Results can be accessed at https://h2.cpr.ku.dk/ .

Abstract Headache disorders are the most common disorders of the nervous system. The lifetime prevalence of headache disorders show that some individuals never experience headache. The etiology of complete freedom from headache is not known. To assess genetic variants associated with complete freedom from headache, we performed a genome-wide association study of individuals who have never experienced a headache. We included 63,992 individuals (2,998 individuals with complete freedom from headache and 60,994 controls) from the Danish Blood Donor Study Genomic Cohort. Participants were included in two rounds, from 2015 to 2018 and in 2020. We discovered a genome-wide significant association, with the lead variant rs7904615[G] in ADARB2 (EAF = 27%, OR = 1.20 [1.13–1.27], p = 3.92 × 10−9). The genomic locus was replicated in a non-overlapping cohort of 13,032 individuals (539 individuals with complete freedom from headache and 12,493 controls) from the Danish Blood Donor Study Genomic Cohort (p

Abstract Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson’s disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Abstract Background The emerging use of biomarkers in research and tailored care introduces a need for information about the association between biomarkers and basic demographics and lifestyle factors revealing expectable concentrations in healthy individuals while considering general demographic differences. Methods A selection of 47 biomarkers, including markers of inflammation and vascular stress, were measured in plasma samples from 9876 Danish Blood Donor Study participants. Using regression models, we examined the association between biomarkers and sex, age, Body Mass Index (BMI), and smoking. Results Here we show that concentrations of inflammation and vascular stress biomarkers generally increase with higher age, BMI, and smoking. Sex-specific effects are observed for multiple biomarkers. Conclusion This study provides comprehensive information on concentrations of 47 plasma biomarkers in healthy individuals. The study emphasizes that knowledge about biomarker concentrations in healthy individuals is critical for improved understanding of disease pathology and for tailored care and decision support tools.

Abstract Background Variation in laboratory healthcare data due to seasonal changes is a widely accepted phenomenon. Seasonal variation is generally not systematically accounted for in healthcare settings. This study applies a newly developed adjustment method for seasonal variation to analyze the effect seasonality has on machine learning model classification of diagnoses. Methods Machine learning methods were trained and tested on ~ 22 million unique records from ~ 575,000 unique patients admitted to Danish hospitals. Four machine learning models (adaBoost, decision tree, neural net, and random forest) classifying 35 diseases of the circulatory system (ICD-10 diagnosis codes, chapter IX) were run before and after seasonal adjustment of 23 laboratory reference intervals (RIs). The effect of the adjustment was benchmarked via its contribution to machine learning models trained using hyperparameter optimization and assessed quantitatively using performance metrics (AUROC and AUPRC). Results Seasonally adjusted RIs significantly improved cardiovascular disease classification in 24 of the 35 tested cases when using neural net models. Features with the highest average feature importance (via SHAP explainability) across all disease models were sex, C- reactive protein, and estimated glomerular filtration. Classification of diseases of the vessels, such as thrombotic diseases and other atherosclerotic diseases consistently improved after seasonal adjustment. Conclusions As data volumes increase and data-driven methods are becoming more advanced, it is essential to improve data quality at the pre-processing level. This study presents a method that makes it feasible to introduce seasonally adjusted RIs into the clinical research space in any disease domain. Seasonally adjusted RIs generally improve diagnoses classification and thus, ought to be considered and adjusted for in clinical decision support methods.

Abstract Social trust is a heritable trait that has been linked with physical health and longevity. In this study, we performed genome-wide association studies of self-reported social trust in n = 33,882 Danish blood donors. We observed genome-wide and local evidence of genetic similarity with other brain-related phenotypes and estimated the single nucleotide polymorphism-based heritability of trust to be 6% (95% confidence interval = (2.1, 9.9)). In our discovery cohort (n = 25,819), we identified one significantly associated locus (lead variant: rs12776883) in an intronic enhancer region of PLPP4, a gene highly expressed in brain, kidneys, and testes. However, we could not replicate the signal in an independent set of donors who were phenotyped a year later (n = 8063). In the subsequent meta-analysis, we found a second significantly associated variant (rs71543507) in an intergenic enhancer region. Overall, our work confirms that social trust is heritable, and provides an initial look into the genetic factors that influence it.

Artificial intelligence (AI) and algorithms are heralded as significant solutions to the widening gap between the rising healthcare needs of ageing and multi-morbid populations and the scarcity of resources to provide such care. Objective This article investigates how the PMHnet algorithm – an AI prognostication tool developed in Denmark to predict the one-year all-cause mortality risk for patients hospitalized with ischemic heart disease – was presented to cardiologists working in the hospital setting, and how they responded to this novel decision-support tool. Methods Empirically, we draw upon ethnographic fieldwork in the Danish-led international research project, PM Heart, which since 2019 has developed the PMHnet algorithm and implemented the software into the electronic health record system in hospitals in Eastern Denmark (the Capital Region and Region Zealand). Results Paying careful attention to the hopes and concerns of cardiologists who will have to embrace and adapt to algorithmic tools in their everyday work of diagnosing and treating patients, we identify three analytical themes meriting attention when AI is implemented in healthcare: 1) the re-negotiation of agency and autonomy in human-algorithm relations, 2) accountability in algorithmic prognostication and 3) the complex relationship between association and causation actualized by predictive algorithms. From these analytical themes, we elicit methodological questions to guide future ethnographic explorations of how AI and advanced algorithms are put to use in the healthcare system, with what implications, and for whom. Conclusion We conclude that local, qualitative investigations of how algorithms are used, embraced and contested in everyday clinical practice are needed in order to understand their implications – good and bad, intended and unintended – for clinicians, patients and healthcare provision.

Summary: Seasonality in laboratory healthcare data is associated with possible under- and overdiagnoses of patients in the clinic. Here, we present a protocol to analyze electronic health record data for seasonality patterns and adjust existing reference intervals for these changes using R software. We describe steps for preprocessing population-wide patient laboratory data into a single dataset. We then detail steps for defining strata, normalizing to median, and fitting data to selected functions.For complete details on the use and execution of this protocol, please refer to Muse et al. (2023).1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Abstract We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

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.

Abstract Objectives Allele counts of sequence variants obtained by whole genome sequencing (WGS) often play a central role in interpreting the results of genetic and genomic research. However, such variant counts are not readily available for individuals in the Danish population. Here, we present a dataset with allele counts for sequence variants (single nucleotide variants (SNVs) and indels) identified from WGS of 8,671 (5,418 females) individuals from the Danish population. The data resource is based on WGS data from three independent research projects aimed at assessing genetic risk factors for cardiovascular, psychiatric, and headache disorders. To enable the sharing of information on sequence variation in Danish individuals, we created summarized statistics on allele counts from anonymized data and made them available through the European Genome-phenome Archive (EGA, https://identifiers.org/ega.dataset:EGAD00001009756 ) and in a dedicated browser, DanMAC5 (available at www.danmac5.dk ). The summary level data and the DanMAC5 browser provide insight into the allelic spectrum of sequence variants segregating in the Danish population, which is important in variant interpretation. Data description Three WGS datasets with an average coverage of 30x were processed independently using the same quality control pipeline. Subsequently, we summarized, filtered, and merged allele counts to create a high-quality summary level dataset of sequence variants.

Abstract We identify biomarkers for disease progression in three type 2 diabetes cohorts encompassing 2,973 individuals across three molecular classes, metabolites, lipids and proteins. Homocitrulline, isoleucine and 2-aminoadipic acid, eight triacylglycerol species, and lowered sphingomyelin 42:2;2 levels are predictive of faster progression towards insulin requirement. Of ~1,300 proteins examined in two cohorts, levels of GDF15/MIC-1, IL-18Ra, CRELD1, NogoR, FAS, and ENPP7 are associated with faster progression, whilst SMAC/DIABLO, SPOCK1 and HEMK2 predict lower progression rates. In an external replication, proteins and lipids are associated with diabetes incidence and prevalence. NogoR/RTN4R injection improved glucose tolerance in high fat-fed male mice but impaired it in male db/db mice. High NogoR levels led to islet cell apoptosis, and IL-18R antagonised inflammatory IL-18 signalling towards nuclear factor kappa-B in vitro. This comprehensive, multi-disciplinary approach thus identifies biomarkers with potential prognostic utility, provides evidence for possible disease mechanisms, and identifies potential therapeutic avenues to slow diabetes progression.

Abstract The intestinal lamina propria contains a diverse network of fibroblasts that provide key support functions to cells within their local environment. Despite this, our understanding of the diversity, location and ontogeny of fibroblasts within and along the length of the intestine remains incomplete. Here we show that the small and large intestinal lamina propria contain similar fibroblast subsets that locate in specific anatomical niches. Nevertheless, we find that the transcriptional profile of similar fibroblast subsets differs markedly between the small intestine and colon suggesting region specific functions. We perform in vivo transplantation and lineage-tracing experiments to demonstrate that adult intestinal fibroblast subsets, smooth muscle cells and pericytes derive from Gli1-expressing precursors present in embryonic day 12.5 intestine. Trajectory analysis of single cell RNA-seq datasets of E12.5 and adult mesenchymal cells suggest that adult smooth muscle cells and fibroblasts derive from distinct embryonic intermediates and that adult fibroblast subsets develop in a linear trajectory from CD81+ fibroblasts. Finally, we provide evidence that colonic subepithelial PDGFRαhi fibroblasts comprise several functionally distinct populations that originate from an Fgfr2-expressing fibroblast intermediate. Our results provide insights into intestinal stromal cell diversity, location, function, and ontogeny, with implications for intestinal development and homeostasis.

Pancreatic cancer is one of the deadliest cancer types with poor treatment options. Better detection of early symptoms and relevant disease correlations could improve pancreatic cancer prognosis. In this retrospective study, we used symptom and disease codes (ICD-10) from the Danish National Patient Registry (NPR) encompassing 6.9 million patients from 1994 to 2018,, of whom 23,592 were diagnosed with pancreatic cancer. The Danish cancer registry included 18,523 of these patients. To complement and compare the registry diagnosis codes with deeper clinical data, we used a text mining approach to extract symptoms from free text clinical notes in electronic health records (3078 pancreatic cancer patients and 30,780 controls). We used both data sources to generate and compare symptom disease trajectories to uncover temporal patterns of symptoms prior to pancreatic cancer diagnosis for the same patients. We show that the text mining of the clinical notes was able to complement the registry-based symptoms by capturing more symptoms prior to pancreatic cancer diagnosis. For example, ‘Blood pressure reading without diagnosis’, ‘Abnormalities of heartbeat’, and ‘Intestinal obstruction’ were not found for the registry-based analysis. Chaining symptoms together in trajectories identified two groups of patients with lower median survival (

Digital health technologies have been in use for many years in a wide spectrum of healthcare scenarios. This narrative review outlines the current use and the future strategies and significance of digital health technologies in modern healthcare applications. It covers the current state of the scientific field (delineating major strengths, limitations, and applications) and envisions the future impact of relevant emerging key technologies. Furthermore, we attempt to provide recommendations for innovative approaches that would accelerate and benefit the research, translation and utilization of digital health technologies.

Polypharmacy has generally been assessed by raw counts of different drugs administered concomitantly to the same patients; not with respect to the likelihood of dosage-adjustments. To address this aspect of polypharmacy, the objective of the present study was to identify co-medications associated with more frequent dosage adjustments. The data foundation was electronic health records from 3.2 million inpatient admissions at Danish hospitals (2008-2016). The likelihood of dosage-adjustments when two drugs were administered concomitantly were computed using Bayesian logistic regressions. We identified 3,993 co-medication pairs that associate significantly with dosage changes when administered together. Of these pairs, 2,412 (60%) did associate with readmission, mortality or longer stays, while 308 (8%) associated with reduced kidney function. In comparison to co-medications pairs that were previously classified as drug-drug interactions, pairs not classified as drug-drug interactions had higher odds ratios of dosage modifications than drug pairs with an established interaction. Drug pairs not corresponding to known drug-drug interactions while still being associated significantly with dosage changes were prescribed to fewer patients and mentioned more rarely together in the literature. We hypothesize that some of these pairs could be associated with yet to be discovered interactions as they may be harder to identify in smaller-scale studies.

BackgroundYKL-40, also known as chitinase-3-like protein 1 (CHI3L1), is a secreted glycoprotein produced by various cell types including stromal, immune, and cancer cells. It contributes to cancer progression through tumor-promoting inflammation and has been shown to inhibit the cytotoxicity of T and NK lymphocytes. In vivo studies have demonstrated synergistic anti-cancer effects of blocking YKL-40 in combination with immune checkpoint inhibitors (ICIs). Biomarkers for the prediction of the response to ICIs are highly needed. We investigated the association between plasma YKL-40 and clinical benefit and survival in patients with metastatic pancreatic cancer (mPC) receiving ICIs and stereotactic body radiotherapy (SBRT).MethodsBlood samples were collected from 84 patients with mPC who participated in the randomized phase II CheckPAC study, in which patients received nivolumab with or without ipilimumab combined with a single fraction of SBRT. Plasma YKL-40 was measured using a commercial ELISA kit.ResultsElevated baseline plasma YKL-40 was an independent predictor of shorter overall survival (OS) (HR 2.19, 95% CI 1.21–3.95). A ≥ 40% decrease in plasma YKL-40 during treatment was associated with longer progression-free survival (p = 0.009) and OS (p = 0.0028). There was no correlation between plasma YKL-40 and the tumor burden marker CA19-9 at baseline or during treatment.ConclusionThis study contributes new knowledge regarding YKL-40 as a predictor of clinical benefit from ICIs and radiotherapy. These exploratory results warrant further investigation of YKL-40 as a biomarker for patients treated with immunotherapies.Clinical trial registrationClinicaltrials.gov, identifier NCT02866383.

Abstract Prediction of survival for patients in intensive care units (ICUs) has been subject to intense research. However, no models exist that embrace the multiverse of data in ICUs. It is an open question whether deep learning methods using automated data integration with minimal pre-processing of mixed data domains such as free text, medical history and high-frequency data can provide discrete-time survival estimates for individual ICU patients. We trained a deep learning model on data from patients admitted to ten ICUs in the Capital Region of Denmark and the Region of Southern Denmark between 2011 and 2018. Inspired by natural language processing we mapped the electronic patient record data to an embedded representation and fed the data to a recurrent neural network with a multi-label output layer representing the chance of survival at different follow-up times. We evaluated the performance using the time-dependent concordance index. In addition, we quantified and visualized the drivers of survival predictions using the SHAP methodology. We included 37,355 admissions of 29,417 patients in our study. Our deep learning models outperformed traditional Cox proportional-hazard models with concordance index in the ranges 0.72–0.73, 0.71–0.72, 0.71, and 0.69–0.70, for models applied at baseline 0, 24, 48, and 72 h, respectively. Deep learning models based on a combination of entity embeddings and survival modelling is a feasible approach to obtain individualized survival estimates in data-rich settings such as the ICU. The interpretable nature of the models enables us to understand the impact of the different data domains.

Novel biomarkers are key to addressing the ongoing pandemic of type 2 diabetes mellitus. While new technologies have improved the potential of identifying such biomarkers, at the same time there is an increasing need for informed prioritization to ensure efficient downstream verification. We have built BALDR, an automated pipeline for biomarker comparison and prioritization in the context of diabetes. BALDR includes protein, gene, and disease data from major public repositories, text-mining data, and human and mouse experimental data from the IMI2 RHAPSODY consortium. These data are provided as easy-to-read figures and tables enabling direct comparison of up to 20 biomarker candidates for diabetes through the public website https://baldr.cpr.ku.dk.

Frequent assessment of the severity of illness for hospitalized patients is essential in clinical settings to prevent outcomes such as in-hospital mortality and unplanned admission to the intensive care unit (ICU). Classical severity scores have been developed typically using relatively few patient features. Recently, deep learning-based models demonstrated better individualized risk assessments compared to classic risk scores, thanks to the use of aggregated and more heterogeneous data sources for dynamic risk prediction. We investigated to what extent deep learning methods can capture patterns of longitudinal change in health status using time-stamped data from electronic health records. We developed a deep learning model based on embedded text from multiple data sources and recurrent neural networks to predict the risk of the composite outcome of unplanned ICU transfer and in-hospital death. The risk was assessed at regular intervals during the admission for different prediction windows. Input data included medical history, biochemical measurements, and clinical notes from a total of 852,620 patients admitted to non-intensive care units in 12 hospitals in Denmark's Capital Region and Region Zealand during 2011-2016 (with a total of 2,241,849 admissions). We subsequently explained the model using the Shapley algorithm, which provides the contribution of each feature to the model outcome. The best model used all data modalities with an assessment rate of 6 hours, a prediction window of 14 days and an area under the receiver operating characteristic curve of 0.898. The discrimination and calibration obtained with this model make it a viable clinical support tool to detect patients at higher risk of clinical deterioration, providing clinicians insights into both actionable and non-actionable patient features.

Summary: Background: Type 1 diabetes is a complex heterogenous autoimmune disease without therapeutic interventions available to prevent or reverse the disease. This study aimed to identify transcriptional changes associated with the disease progression in patients with recent-onset type 1 diabetes. Methods: Whole-blood samples were collected as part of the INNODIA study at baseline and 12 months after diagnosis of type 1 diabetes. We used linear mixed-effects modelling on RNA-seq data to identify genes associated with age, sex, or disease progression. Cell-type proportions were estimated from the RNA-seq data using computational deconvolution. Associations to clinical variables were estimated using Pearson's or point-biserial correlation for continuous and dichotomous variables, respectively, using only complete pairs of observations. Findings: We found that genes and pathways related to innate immunity were downregulated during the first year after diagnosis. Significant associations of the gene expression changes were found with ZnT8A autoantibody positivity. Rate of change in the expression of 16 genes between baseline and 12 months was found to predict the decline in C-peptide at 24 months. Interestingly and consistent with earlier reports, increased B cell levels and decreased neutrophil levels were associated with the rapid progression. Interpretation: There is considerable individual variation in the rate of progression from appearance of type 1 diabetes-specific autoantibodies to clinical disease. Patient stratification and prediction of disease progression can help in developing more personalised therapeutic strategies for different disease endotypes. Funding: A full list of funding bodies can be found under Acknowledgments.

Background: Whether natural selection may have attributed to the observed blood group frequency differences between populations remains debatable. The ABO system has been associated with several diseases and recently also with susceptibility to COVID-19 infection. Associative studies of the RhD system and diseases are sparser. A large disease-wide risk analysis may further elucidate the relationship between the ABO/RhD blood groups and disease incidence. Methods: We performed a systematic log-linear quasi-Poisson regression analysis of the ABO/RhD blood groups across 1,312 phecode diagnoses. Unlike prior studies, we determined the incidence rate ratio for each individual ABO blood group relative to all other ABO blood groups as opposed to using blood group O as the reference. Moreover, we used up to 41 years of nationwide Danish follow-up data, and a disease categorization scheme specifically developed for diagnosis-wide analysis. Further, we determined associations between the ABO/RhD blood groups and the age at the first diagnosis. Estimates were adjusted for multiple testing. Results: The retrospective cohort included 482,914 Danish patients (60.4% females). The incidence rate ratios (IRRs) of 101 phecodes were found statistically significant between the ABO blood groups, while the IRRs of 28 phecodes were found statistically significant for the RhD blood group. The associations included cancers and musculoskeletal-, genitourinary-, endocrinal-, infectious-, cardiovascular-, and gastrointestinal diseases. Conclusions: We found associations of disease-wide susceptibility differences between the blood groups of the ABO and RhD systems, including cancer of the tongue, monocytic leukemia, cervical cancer, osteoarthrosis, asthma, and HIV- and hepatitis B infection. We found marginal evidence of associations between the blood groups and the age at first diagnosis. Funding: Novo Nordisk Foundation and the Innovation Fund Denmark

Abstract Background Patients diagnosed with ischemic heart disease (IHD) are becoming increasingly multi-morbid, and studies designed to analyze the full spectrum are few. Methods Disease trajectories, defined as time-ordered series of diagnoses, were used to study the temporality of multi-morbidity. The main data source was The Danish National Patient Register (NPR) comprising 7,179,538 individuals in the period 1994–2018. Patients with a diagnosis code for IHD were included. Relative risks were used to quantify the strength of the association between diagnostic co-occurrences comprised of two diagnoses that were overrepresented in the same patients. Multiple linear regression models were then fitted to test for temporal associations among the diagnostic co-occurrences, termed length two disease trajectories. Length two disease trajectories were then used as basis for constructing disease trajectories of three diagnoses. Results In a cohort of 570,157 IHD disease patients, we identified 1447 length two disease trajectories and 4729 significant length three disease trajectories. These included 459 distinct diagnoses. Disease trajectories were dominated by chronic diseases and not by common, acute diseases such as pneumonia. The temporal association of atrial fibrillation (AF) and IHD differed in different IHD subpopulations. We found an association between osteoarthritis (OA) and heart failure (HF) among patients diagnosed with OA, IHD, and then HF only. Conclusions The sequence of diagnoses is important in characterization of multi-morbidity in IHD patients as the disease trajectories. The study provides evidence that the timing of AF in IHD marks distinct IHD subpopulations; and secondly that the association between osteoarthritis and heart failure is dependent on IHD.

The underlying genetics of carpal tunnel syndrome is not well understood. Here, the authors perform a GWAS meta-analysis for carpal tunnel syndrome finding variants at 50 loci with connections to the extracellular matrix discovered through various functional analyses.

This study investigates the dynamics of immunological markers after first SARS-CoV-2 vaccination dose in cohort of healthcare professionals in Denmark. Natural infection was associated with higher antibody responses, and IgG decline varied by age, sex, T-cell response, previous infection, and interval between vaccine doses.

Florian Uellendahl-Werth et al. conduct cross-tissue transcriptome-wide association studies to explore genetic mechanisms shared across immune-related and psychiatric traits. Their results identify several genes (including PPP3CA) that could mediate the interplay between psychiatric and inflammatory disease.

People with kidney disease are disproportionately affected by atherosclerosis for unclear reasons. Soluble urokinase plasminogen activator receptor (suPAR) is an immune-derived mediator of kidney disease, levels of which are strongly associated with cardiovascular outcomes. We assessed suPAR’s pathogenic involvement in atherosclerosis using epidemiologic, genetic, and experimental approaches. We found serum suPAR levels to be predictive of coronary artery calcification and cardiovascular events in 5,406 participants without known coronary disease. In a genome-wide association meta-analysis including over 25,000 individuals, we identified a missense variant in the plasminogen activator, urokinase receptor (PLAUR) gene (rs4760), confirmed experimentally to lead to higher suPAR levels. Mendelian randomization analysis in the UK Biobank using rs4760 indicated a causal association between genetically predicted suPAR levels and atherosclerotic phenotypes. In an experimental model of atherosclerosis, proprotein convertase subtilisin/kexin–9 (Pcsk9) transfection in mice overexpressing suPAR (suPARTg) led to substantially increased atherosclerotic plaques with necrotic cores and macrophage infiltration compared with those in WT mice, despite similar cholesterol levels. Prior to induction of atherosclerosis, aortas of suPARTg mice excreted higher levels of CCL2 and had higher monocyte counts compared with WT aortas. Aortic and circulating suPARTg monocytes exhibited a proinflammatory profile and enhanced chemotaxis. These findings characterize suPAR as a pathogenic factor for atherosclerosis acting at least partially through modulation of monocyte function.

Abstract 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

Tuberculosis (TB) presents a serious health problem with approximately a quarter of the world’s population infected with Mycobacterium tuberculosis (M. tuberculosis) in an asymptomatic latent state of which 5–10% develops active TB at some point in their lives. The antimicrobial protein cathelicidin has broad antimicrobial activity towards viruses and bacteria including M. tuberculosis. Vitamin D increases the expression of cathelicidin in many cell types including macrophages, and it has been suggested that the vitamin D-mediated antimicrobial activity against M. tuberculosis is dependent on the induction of cathelicidin. However, unraveling the immunoregulatory effects of vitamin D in humans is hampered by the lack of suitable experimental models. We have previously described a family in which members suffer from hereditary vitamin D-resistant rickets (HVDRR). The family carry a mutation in the DNA-binding domain of the vitamin D receptor (VDR). This mutation leads to a non-functional VDR, meaning that vitamin D cannot exert its effect in family members homozygous for the mutation. Studies of HVDRR patients open unique possibilities to gain insight in the immunoregulatory roles of vitamin D in humans. Here we describe the impaired ability of macrophages to produce cathelicidin in a HVDRR patient, who in her adolescence suffered from extrapulmonary TB. The present case is a rare experiment of nature, which illustrates the importance of vitamin D in the pathophysiology of combating M. tuberculosis.

Dowsett and colleagues used a genome-wide association approach to investigate the genetic influence on soluble urokinase-type plasminogen activator receptor presence in the plasma of humans. Their findings indicate a 60% heritability factor in British twins, and using a wide sample of Northern European genome samples they identify eleven genetic loci associated with an increase or decrease of this chronic inflammation marker.

Abstract Reports of persistent symptoms after hospitalization with COVID-19 have raised concern of a “long COVID” syndrome. This study aimed at determining the prevalence of and risk factors for acute and persistent symptoms in non-hospitalized patients with polymerase chain reaction (PCR) confirmed COVID-19. We conducted a cohort study of non-hospitalized participants identified via the Danish Civil Registration System with a SARS-CoV-2-positive PCR-test and available biobank samples. Participants received a digital questionnaire on demographics and COVID-19-related symptoms. Persistent symptoms: symptoms > 4 weeks (in sensitivity analyses > 12 weeks). We included 445 participants, of whom 34% were asymptomatic. Most common acute symptoms were fatigue, headache, and sneezing, while fatigue and reduced smell and taste were most severe. Persistent symptoms, most commonly fatigue and memory and concentration difficulties, were reported by 36% of 198 symptomatic participants with follow-up > 4 weeks. Risk factors for persistent symptoms included female sex (women 44% vs. men 24%, odds ratio 2.7, 95% CI 1.4–5.1, p = 0.003) and BMI (odds ratio 1.1, 95% CI 1.0–1.2, p = 0.001). In conclusion, among non-hospitalized PCR-confirmed COVID-19 patients one third were asymptomatic while one third of symptomatic participants had persistent symptoms illustrating the heterogeneity of disease presentation. These findings should be considered in health care planning and policy making related to COVID-19.

Type I interferons (IFNs) are essential for antiviral immunity, appear to represent a key component of mRNA vaccine-adjuvanticity, and correlate with severity of systemic autoimmune disease. Relevant to all, type I IFNs can enhance germinal center (GC) B cell responses but underlying signaling pathways are incompletely understood. Here, we demonstrate that a succinct type I IFN response promotes GC formation and associated IgG subclass distribution primarily through signaling in cDCs and B cells. Type I IFN signaling in cDCs, distinct from cDC1, stimulates development of separable Tfh and Th1 cell subsets. However, Th cell-derived IFN-γ induces T-bet expression and IgG2c isotype switching in B cells prior to this bifurcation and has no evident effects once GCs and bona fide Tfh cells developed. This pathway acts in synergy with early B cell-intrinsic type I IFN signaling, which reinforces T-bet expression in B cells and leads to a selective amplification of the IgG2c+ GC B cell response. Despite the strong Th1 polarizing effect of type I IFNs, the Tfh cell subset develops into IL-4 producing cells that control the overall magnitude of the GCs and promote generation of IgG1+ GC B cells. Thus, type I IFNs act on B cells and cDCs to drive GC formation and to coordinate IgG subclass distribution through divergent Th1 and Tfh cell-dependent pathways.

The focus of this paper is on strategic approaches for establishing population-based prospective cohorts that collect and store biological samples from very large numbers of participants to help identify the determinants of common health outcomes. In particular, it aims to address key issues related to investigation of genetic, as well as social, environmental, and ancestral, diversity; generation of detailed genetic and other types of assay data; collection of detailed lifestyle and environmental exposure information; follow-up and characterization of incident health outcomes; and overcoming obstacles to data sharing and access (including capacity building). It concludes that there is a need for strategic planning at an international level (rather than the current ad hoc approach) toward the development of a carefully selected set of deeply characterized large-scale prospective cohorts that are readily accessible by researchers around the world.

Abstract Diabetic retinopathy (DR) is a leading cause of blindness and affects millions of people throughout the world. Early detection and timely checkups are key to reduce the risk of blindness. Automated grading of DR is a cost-effective way to ensure early detection and timely checkups. Deep learning or more specifically convolutional neural network (CNN)—based methods produce state-of-the-art performance in DR detection. Whilst CNN based methods have been proposed, no comparisons have been done between the extracted image features and their clinical relevance. Here we first adopt a CNN visualization strategy to discover the inherent image features involved in the CNN’s decision-making process. Then, we critically analyze those features with respect to commonly known pathologies namely microaneurysms, hemorrhages and exudates, and other ocular components. We also critically analyze different CNNs by considering what image features they pick up during learning to predict and justify their clinical relevance. The experiments are executed on publicly available fundus datasets (EyePACS and DIARETDB1) achieving an accuracy of 89 ~ 95% with AUC, sensitivity and specificity of respectively 95 ~ 98%, 74 ~ 86%, and 93 ~ 97%, for disease level grading of DR. Whilst different CNNs produce consistent classification results, the rate of picked-up image features disagreement between models could be as high as 70%.

Bell et al. report 46 new loci associated with biomarkers of iron homeostasis, including ferritin levels, iron binding capacity, and iron saturation, in the Icelandic, Danish and UK populations. The associated loci point to new iron-regulating proteins and important genetic differences between men and women.

Abstract Diagnostic errors are common and can lead to harmful treatments. We present a data-driven, generic approach for identifying patients at risk of being mis- or overdiagnosed, here exemplified by chronic obstructive pulmonary disease (COPD). It has been estimated that 5–60% of all COPD cases are misdiagnosed. High-throughput methods are therefore needed in this domain. We have used a national patient registry, which contains hospital diagnoses for 6.9 million patients across the entire Danish population for 21 years and identified statistically significant disease trajectories for COPD patients. Using 284,154 patients diagnosed with COPD, we identified frequent disease trajectories comprising time-ordered comorbidities. Interestingly, as many as 42,459 patients did not present with these time-ordered, common comorbidities. Comparison of the individual disease history for each non-follower to the COPD trajectories, demonstrated that 9597 patients were unusual. Survival analysis showed that this group died significantly earlier than COPD patients following a trajectory. Out of the 9597 patients, we identified one subgroup comprising 2185 patients at risk of misdiagnosed COPD without the typical events of COPD patients. In all, 10% of these patients were diagnosed with lung cancer, and it seems likely that they are underdiagnosed for lung cancer as their laboratory test values and survival pattern are similar to such patients. Furthermore, only 4% had a lung function test to confirm the COPD diagnosis. Another subgroup with 2368 patients were found to be at risk of “classically” overdiagnosed COPD that survive >5.5 years after the COPD diagnosis, but without the typical complications of COPD.

Abstract Background The rising prevalence of type 2 diabetes (T2D) poses a major global challenge. It remains unresolved to what extent transcriptomic signatures of metabolic dysregulation and T2D can be observed in easily accessible tissues such as blood. Additionally, large-scale human studies are required to further our understanding of the putative inflammatory component of insulin resistance and T2D. Here we used transcriptomics data from individuals with (n = 789) and without (n = 2127) T2D from the IMI-DIRECT cohorts to describe the co-expression structure of whole blood that mainly reflects processes and cell types of the immune system, and how it relates to metabolically relevant clinical traits and T2D. Methods Clusters of co-expressed genes were identified in the non-diabetic IMI-DIRECT cohort and evaluated with regard to stability, as well as preservation and rewiring in the cohort of individuals with T2D. We performed functional and immune cell signature enrichment analyses, and a genome-wide association study to describe the genetic regulation of the modules. Phenotypic and trans-omics associations of the transcriptomic modules were investigated across both IMI-DIRECT cohorts. Results We identified 55 whole blood co-expression modules, some of which clustered in larger super-modules. We identified a large number of associations between these transcriptomic modules and measures of insulin action and glucose tolerance. Some of the metabolically linked modules reflect neutrophil-lymphocyte ratio in blood while others are independent of white blood cell estimates, including a module of genes encoding neutrophil granule proteins with antibacterial properties for which the strongest associations with clinical traits and T2D status were observed. Through the integration of genetic and multi-omics data, we provide a holistic view of the regulation and molecular context of whole blood transcriptomic modules. We furthermore identified an overlap between genetic signals for T2D and co-expression modules involved in type II interferon signaling. Conclusions Our results offer a large-scale map of whole blood transcriptomic modules in the context of metabolic disease and point to novel biological candidates for future studies related to T2D.

The Danish health system has been collecting health-related data on the entire Danish population for years. Here the authors present the Danish Disease Trajectory Browser (DTB), which allows users to explore population-wide disease progression patterns from data collected between 1994 and 2018.

Abstract Background Congenital heart disease (CHD) occurs in almost 1% of newborn children and is considered a multifactorial disorder. CHD may segregate in families due to significant contribution of genetic factors in the disease etiology. The aim of the study was to identify pathophysiological mechanisms in families segregating CHD. Methods We used whole exome sequencing to identify rare genetic variants in ninety consenting participants from 32 Danish families with recurrent CHD. We applied a systems biology approach to identify developmental mechanisms influenced by accumulation of rare variants. We used an independent cohort of 714 CHD cases and 4922 controls for replication and performed functional investigations using zebrafish as in vivo model. Results We identified 1785 genes, in which rare alleles were shared between affected individuals within a family. These genes were enriched for known cardiac developmental genes, and 218 of these genes were mutated in more than one family. Our analysis revealed a functional cluster, enriched for proteins with a known participation in calcium signaling. Replication in an independent cohort confirmed increased mutation burden of calcium-signaling genes in CHD patients. Functional investigation of zebrafish orthologues of ITPR1, PLCB2, and ADCY2 verified a role in cardiac development and suggests a combinatorial effect of inactivation of these genes. Conclusions The study identifies abnormal calcium signaling as a novel pathophysiological mechanism in human CHD and confirms the complex genetic architecture underlying CHD.

Abstract Background Studies on tumor-secreted microRNAs point to a functional role of these in cellular communication and reprogramming of the tumor microenvironment. Uptake of tumor-secreted microRNAs by neighboring cells may result in the silencing of mRNA targets and, in turn, modulation of the transcriptome. Studying miRNAs externalized from tumors could improve cancer patient diagnosis and disease monitoring and help to pinpoint which miRNA-gene interactions are central for tumor properties such as invasiveness and metastasis. Methods Using a bioinformatics approach, we analyzed the profiles of secreted tumor and normal interstitial fluid (IF) microRNAs, from women with breast cancer (BC). We carried out differential abundance analysis (DAA), to obtain miRNAs, which were enriched or depleted in IFs, from patients with different clinical traits. Subsequently, miRNA family enrichment analysis was performed to assess whether any families were over-represented in the specific sets. We identified dysregulated genes in tumor tissues from the same cohort of patients and constructed weighted gene co-expression networks, to extract sets of co-expressed genes and co-abundant miRNAs. Lastly, we integrated miRNAs and mRNAs to obtain interaction networks and supported our findings using prediction tools and cancer gene databases. Results Network analysis showed co-expressed genes and miRNA regulators, associated with tumor lymphocyte infiltration. All of the genes were involved in immune system processes, and many had previously been associated with cancer immunity. A subset of these, BTLA, CXCL13, IL7R, LAMP3, and LTB, was linked to the presence of tertiary lymphoid structures and high endothelial venules within tumors. Co-abundant tumor interstitial fluid miRNAs within this network, including miR-146a and miR-494, were annotated as negative regulators of immune-stimulatory responses. One co-expression network encompassed differences between BC subtypes. Genes differentially co-expressed between luminal B and triple-negative breast cancer (TNBC) were connected with sphingolipid metabolism and predicted to be co-regulated by miR-23a. Co-expressed genes and TIF miRNAs associated with tumor grade were BTRC, CHST1, miR-10a/b, miR-107, miR-301a, and miR-454. Conclusion Integration of IF miRNAs and mRNAs unveiled networks associated with patient clinicopathological traits, and underlined molecular mechanisms, specific to BC sub-groups. Our results highlight the benefits of an integrative approach to biomarker discovery, placing secreted miRNAs within a biological context.

Abstract Background Medication errors have been identified as the most common preventable cause of adverse events. The lack of granularity in medication error terminology has led pharmacovigilance experts to rely on information in individual case safety reports’ (ICSRs) codes and narratives for signal detection, which is both time consuming and labour intensive. Thus, there is a need for complementary methods for the detection of medication errors from ICSRs. The aim of this study is to evaluate the utility of two natural language processing text mining methods as complementary tools to the traditional approach followed by pharmacovigilance experts for medication error signal detection. Methods The safety surveillance advisor (SSA) method, I2E text mining and University of Copenhagen Center for Protein Research (CPR) text mining, were evaluated for their ability to extract cases containing a type of medication error where patients extracted insulin from a prefilled pen or cartridge by a syringe. A total of 154,209 ICSRs were retrieved from Novo Nordisk’s safety database from January 1987 to February 2018. Each method was evaluated by recall (sensitivity) and precision (positive predictive value). Results We manually annotated 2533 ICSRs to investigate whether these contained the sought medication error. All these ICSRs were then analysed using the three methods. The recall was 90.4, 88.1 and 78.5% for the CPR text mining, the SSA method and the I2E text mining, respectively. Precision was low for all three methods ranging from 3.4% for the SSA method to 1.9 and 1.6% for the CPR and I2E text mining methods, respectively. Conclusions Text mining methods can, with advantage, be used for the detection of complex signals relying on information found in unstructured text (e.g., ICSR narratives) as standardised and both less labour-intensive and time-consuming methods compared to traditional pharmacovigilance methods. The employment of text mining in pharmacovigilance need not be limited to the surveillance of potential medication errors but can be used for the ongoing regulatory requests, e.g., obligations in risk management plans and may thus be utilised broadly for signal detection and ongoing surveillance activities.

Abstract Background Most structured clinical data, such as diagnosis codes, are not sufficient to obtain precise phenotypes and assess disease burden. Text mining of clinical notes could provide a basis for detailed profiles of phenotypic traits. The objective of the current study was to determine whether drug dose, regardless of polypharmacy, is associated with the length of clinical notes, and to determine the frequency of adverse events per word in clinical notes. Methods In this observational study, we utilized restricted-access data from an electronic patient record system. Using three methods (defined daily dose, olanzapine equivalents, and chlorpromazine equivalents) we calculated antipsychotic dose equivalents and compared these with the number of words recorded per treatment day. For each normalization method, the frequencies of adverse events per word in manually curated samples were compared to dose intervals. Results The length of clinical notes per treatment day was positively associated with the prescribed dose for all normalization methods. The number of adverse events per word was stable over the analyzed dose spectrum. Conclusions Assuming that drug dose increases with the severity of disease, the length of clinical notes can serve as a proxy for disease severity. Due to the near-linear relationship, correction of daily word count is unnecessary when text mining for potential adverse drug reactions.

Understanding deregulation of biological pathways in cancer can provide insight into disease etiology and potential therapies. Here, as part of the PanCancer Analysis of Whole Genomes (PCAWG) consortium, the authors present pathway and network analysis of 2583 whole cancer genomes from 27 tumour types.

Disease heritability and genetic correlations between traits depend on genetics, the environment and their interaction. Here, Jia et al. compute disease prevalence curves and disease embeddings from electronic health records and impute heritability for hundreds of diseases and genetic correlations for thousands of disease pairs.

Onychomycosis is a common disease with a significant negative impact on quality of life. While the disease is usually manageable in general practice, a proportion of patients need specialist treatment in academic hospital clinics. However, it is an unknown question whether the incidence in those needing specialist treatments is changing. Furthermore, the comorbidity burden in this patient population severely affected by onychomycosis has never been characterized. We conducted a retrospective study on patients treated for onychomycosis in Danish hospitals from 1994 to 2018. The cohort was observed for 24 years, and the data comprise 7.2 million Danes and their hospital diagnoses. A disease trajectory algorithm was used to examine the comorbidity burden in the cohort. A total of 2,271 patients received hospital treatment for onychomycosis during the time period, of which 1358 (59.8%) were men. The data show an increase in the incidence of hospital-treated cases since 2012 and that the most common comorbidities in this patient population include cardiovascular disease, alcohol-related diagnoses, and diabetes. One explanation of the increase in specialist treatment may include a general increase in patients with decreased resilience to fungal disease. This lack of resilience may both include an increasing elderly population with atherosclerosis, diabetes, and immunosuppression but also a potential increase in patients treated with immunosuppressive agents. Another possible explanation may include a shift in patient expectations in the case of treatment failure. Thus, patients may have an increasing demand for specialist treatment. While our data document an increase in the number of patients in need of specialist treatment for onychomycosis, we suggest future research to examine the general incidence of onychomycosis but also whether this increase in an apparently recalcitrant disease may be attributed to increased antifungal resistance, more specialist treatment options, or increased attention to dermatomycoses.