10-year multimorbidity patterns among people with and without rheumatic and musculoskeletal diseases: an observational cohort study using linked electronic health records from Wales, UK.

OBJECTIVES: To compare the patterns of multimorbidity between people with and without rheumatic and musculoskeletal diseases (RMDs) and to describe how these patterns change by age and sex over time, between 2010 and 2019. PARTICIPANTS: 103 426 people with RMDs and 2.9 million comparators registered in 395 Wales general practices (GPs). Each patient with an RMD aged 0-100 years between January 2010 and December 2019 registered in Clinical Practice Research Welsh practices was matched with up to five comparators without an RMD, based on age, gender and GP code. PRIMARY OUTCOME MEASURES: The prevalence of 29 Elixhauser-defined comorbidities in people with RMDs and comparators categorised by age, gender and GP practices. Conditional logistic regression models were fitted to calculate differences (OR, 95% CI) in associations with comorbidities between cohorts. RESULTS: The most prevalent comorbidities were cardiovascular risk factors, hypertension and diabetes. Having an RMD diagnosis was associated with a significantly higher odds for many conditions including deficiency anaemia (OR 1.39, 95% CI (1.32 to 1.46)), hypothyroidism (OR 1.34, 95% CI (1.19 to 1.50)), pulmonary circulation disorders (OR 1.39, 95% CI 1.12 to 1.73) diabetes (OR 1.17, 95% CI (1.11 to 1.23)) and fluid and electrolyte disorders (OR 1.27, 95% CI (1.17 to 1.38)). RMDs have a higher proportion of multimorbidity (two or more conditions in addition to the RMD) compared with non-RMD group (81% and 73%, respectively in 2019) and the mean number of comorbidities was higher in women from the age of 25 and 50 in men than in non-RMDs group. CONCLUSION: People with RMDs are approximately 1.5 times as likely to have multimorbidity as the general population and provide a high-risk group for targeted intervention studies. The individuals with RMDs experience a greater load of coexisting health conditions, which tend to manifest at earlier ages. This phenomenon is particularly pronounced among women. Additionally, there is an under-reporting of comorbidities in individuals with RMDs.

A scoping review finds a growing trend in studies validating multimorbidity patterns and identifies five broad types of validation methods.

OBJECTIVES: Multimorbidity, the presence of two or more long-term conditions, is a growing public health concern. Many studies use analytical methods to discover multimorbidity patterns from data. We aimed to review approaches used in published literature to validate these patterns. STUDY DESIGN AND SETTING: We systematically searched PubMed and Web of Science for studies published between July 2017 and July 2023 that used analytical methods to discover multimorbidity patterns. RESULTS: Out of 31,617 studies returned by the searches, 172 were included. Of these, 111 studies (64%) conducted validation, the number of studies with validation increased from 53.13% (17 out of 32 studies) to 71.25% (57 out of 80 studies) in 2017-2019 to 2022-2023, respectively. Five types of validation were identified: assessing the association of multimorbidity patterns with clinical outcomes (n = 79), stability across subsamples (n = 26), clinical plausibility (n = 22), stability across methods (n = 7) and exploring common determinants (n = 2). Some studies used multiple types of validation. CONCLUSION: The number of studies conducting a validation of multimorbidity patterns is clearly increasing. The most popular validation approach is assessing the association of multimorbidity patterns with clinical outcomes. Methodological guidance on the validation of multimorbidity patterns is needed.

Using hypergraphs to quantify importance of sets of diseases by healthcare resource utilisation: A retrospective cohort study.

Rates of Multimorbidity (also called Multiple Long Term Conditions, MLTC) are increasing in many developed nations. People with multimorbidity experience poorer outcomes and require more healthcare intervention. Grouping of conditions by health service utilisation is poorly researched. The study population consisted of a cohort of people living in Wales, UK aged 20 years or older in 2000 who were followed up until the end of 2017. Multimorbidity clusters by prevalence and healthcare resource use (HRU) were modelled using hypergraphs, mathematical objects relating diseases via links which can connect any number of diseases, thus capturing information about sets of diseases of any size. The cohort included 2,178,938 people. The most prevalent diseases were hypertension (13.3%), diabetes (6.9%), depression (6.7%) and chronic obstructive pulmonary disease (5.9%). The most important sets of diseases when considering prevalence generally contained a small number of diseases, while the most important sets of diseases when considering HRU were sets containing many diseases. The most important set of diseases taking prevalence and HRU into account was diabetes & hypertension and this combined measure of importance featured hypertension most often in the most important sets of diseases. We have used a single approach to find the most important sets of diseases based on co-occurrence and HRU measures, demonstrating the flexibility of the hypergraph approach. Hypertension, the most important single disease, is silent, underdiagnosed and increases the risk of life threatening co-morbidities. Co-occurrence of endocrine and cardiovascular diseases was common in the most important sets. Combining measures of prevalence with HRU provides insights which would be helpful for those planning and delivering services.

Trajectories in chronic disease accrual and mortality across the lifespan in Wales, UK (2005-2019), by area deprivation profile: linked electronic health records cohort study on 965,905 individuals.

Background: Understanding and quantifying the differences in disease development in different socioeconomic groups of people across the lifespan is important for planning healthcare and preventive services. The study aimed to measure chronic disease accrual, and examine the differences in time to individual morbidities, multimorbidity, and mortality between socioeconomic groups in Wales, UK. Methods: Population-wide electronic linked cohort study, following Welsh residents for up to 20 years (2000-2019). Chronic disease diagnoses were obtained from general practice and hospitalisation records using the CALIBER disease phenotype register. Multi-state models were used to examine trajectories of accrual of 132 diseases and mortality, adjusted for sex, age and area-level deprivation. Restricted mean survival time was calculated to measure time spent free of chronic disease(s) or mortality between socioeconomic groups. Findings: In total, 965,905 individuals aged 5-104 were included, from a possible 2.9 m individuals following a 5-year clearance period, with an average follow-up of 13.2 years (12.7 million person-years). Some 673,189 (69.7%) individuals developed at least one chronic disease or died within the study period. From ages 10 years upwards, the individuals living in the most deprived areas consistently experienced reduced time between health states, demonstrating accelerated transitions to first and subsequent morbidities and death compared to their demographic equivalent living in the least deprived areas. The largest difference were observed in 10 and 20 year old males developing multimorbidity (-0.45 years (99% CI: -0.45, -0.44)) and in 70 year old males dying after developing multimorbidity (-1.98 years (99% CI: -2.01, -1.95)). Interpretation: This study adds to the existing literature on health inequalities by demonstrating that individuals living in more deprived areas consistently experience accelerated time to diagnosis of chronic disease and death across all ages, accounting for competing risks. Funding: UK Medical Research Council, Health Data Research UK, and Administrative Data Research Wales.

Adapting Random Forests to Predict Obesity-Associated Gene Expression.

Random forests (RFs) are effective at predicting gene expression from genotype data. However, a comparison of RF regressors and classifiers, including feature selection and encoding, has been under-explored in the context of gene expression prediction. Specifically, we examine the role of ordinal or one-hot encoding and of data balancing via oversam-pling in the prediction of obesity-associated gene expression. Our work shows that RFs compete with PrediXcan in the prediction of obesity-associated gene expression in subcutaneous adipose tissue, a highly relevant tissue to obesity. Additionally, RFs generate predictions for obesity-associated genes where PrediXcan fails to do so.

Comparative Genetic Analysis of Psoriatic Arthritis and Psoriasis for the Discovery of Genetic Risk Factors and Risk Prediction Modeling.

OBJECTIVES: Psoriatic arthritis (PsA) has a strong genetic component, and the identification of genetic risk factors could help identify the ~30% of psoriasis patients at high risk of developing PsA. Our objectives were to identify genetic risk factors and pathways that differentiate PsA from cutaneous-only psoriasis (PsC) and to evaluate the performance of PsA risk prediction models. METHODS: Genome-wide meta-analyses were conducted separately for 5,065 patients with PsA and 21,286 healthy controls and separately for 4,340 patients with PsA and 6,431 patients with PsC. The heritability of PsA was calculated as a single-nucleotide polymorphism (SNP)-based heritability estimate (h2 SNP ) and biologic pathways that differentiate PsA from PsC were identified using Priority Index software. The generalizability of previously published PsA risk prediction pipelines was explored, and a risk prediction model was developed with external validation. RESULTS: We identified a novel genome-wide significant susceptibility locus for the development of PsA on chromosome 22q11 (rs5754467; P = 1.61 × 10-9 ), and key pathways that differentiate PsA from PsC, including NF-κB signaling (adjusted P = 1.4 × 10-45 ) and Wnt signaling (adjusted P = 9.5 × 10-58 ). The heritability of PsA in this cohort was found to be moderate (h2 SNP = 0.63), which was similar to the heritability of PsC (h2 SNP = 0.61). We observed modest performance of published classification pipelines (maximum area under the curve 0.61), with similar performance of a risk model derived using the current data. CONCLUSION: Key biologic pathways associated with the development of PsA were identified, but the investigation of risk classification revealed modest utility in the available data sets, possibly because many of the PsC patients included in the present study were receiving treatments that are also effective in PsA. Future predictive models of PsA should be tested in PsC patients recruited from primary care.

Application of information theoretic feature selection and machine learning methods for the development of genetic risk prediction models.

In view of the growth of clinical risk prediction models using genetic data, there is an increasing need for studies that use appropriate methods to select the optimum number of features from a large number of genetic variants with a high degree of redundancy between features due to linkage disequilibrium (LD). Filter feature selection methods based on information theoretic criteria, are well suited to this challenge and will identify a subset of the original variables that should result in more accurate prediction. However, data collected from cohort studies are often high-dimensional genetic data with potential confounders presenting challenges to feature selection and risk prediction machine learning models. Patients with psoriasis are at high risk of developing a chronic arthritis known as psoriatic arthritis (PsA). The prevalence of PsA in this patient group can be up to 30% and the identification of high risk patients represents an important clinical research which would allow early intervention and a reduction of disability. This also provides us with an ideal scenario for the development of clinical risk prediction models and an opportunity to explore the application of information theoretic criteria methods. In this study, we developed the feature selection and psoriatic arthritis (PsA) risk prediction models that were applied to a cross-sectional genetic dataset of 1462 PsA cases and 1132 cutaneous-only psoriasis (PsC) cases using 2-digit HLA alleles imputed using the SNP2HLA algorithm. We also developed stratification method to mitigate the impact of potential confounder features and illustrate that confounding features impact the feature selection. The mitigated dataset was used in training of seven supervised algorithms. 80% of data was randomly used for training of seven supervised machine learning methods using stratified nested cross validation and 20% was selected randomly as a holdout set for internal validation. The risk prediction models were then further validated in UK Biobank dataset containing data on 1187 participants and a set of features overlapping with the training dataset.Performance of these methods has been evaluated using the area under the curve (AUC), accuracy, precision, recall, F1 score and decision curve analysis(net benefit). The best model is selected based on three criteria: the 'lowest number of feature subset' with the 'maximal average AUC over the nested cross validation' and good generalisability to the UK Biobank dataset. In the original dataset, with over 100 different bootstraps and seven feature selection (FS) methods, HLA_C_*06 was selected as the most informative genetic variant. When the dataset is mitigated the single most important genetic features based on rank was identified as HLA_B_*27 by the seven different feature selection methods, consistent with previous analyses of this data using regression based methods. However, the predictive accuracy of these single features in post mitigation was found to be moderate (AUC= 0.54 (internal cross validation), AUC=0.53 (internal hold out set), AUC=0.55(external data set)). Sequentially adding additional HLA features based on rank improved the performance of the Random Forest classification model where 20 2-digit features selected by Interaction Capping (ICAP) demonstrated (AUC= 0.61 (internal cross validation), AUC=0.57 (internal hold out set), AUC=0.58 (external dataset)). The stratification method for mitigation of confounding features and filter information theoretic feature selection can be applied to a high dimensional dataset with the potential confounders.

Ranking sets of morbidities using hypergraph centrality.

Multi-morbidity, the health state of having two or more concurrent chronic conditions, is becoming more common as populations age, but is poorly understood. Identifying and understanding commonly occurring sets of diseases is important to inform clinical decisions to improve patient services and outcomes. Network analysis has been previously used to investigate multi-morbidity, but a classic application only allows for information on binary sets of diseases to contribute to the graph. We propose the use of hypergraphs, which allows for the incorporation of data on people with any number of conditions, and also allows us to obtain a quantitative understanding of the centrality, a measure of how well connected items in the network are to each other, of both single diseases and sets of conditions. Using this framework we illustrate its application with the set of conditions described in the Charlson morbidity index using data extracted from routinely collected population-scale, patient level electronic health records (EHR) for a cohort of adults in Wales, UK. Stroke and diabetes were found to be the most central single conditions. Sets of diseases featuring diabetes; diabetes with Chronic Pulmonary Disease, Renal Disease, Congestive Heart Failure and Cancer were the most central pairs of diseases. We investigated the differences between results obtained from the hypergraph and a classic binary graph and found that the centrality of diseases such as paraplegia, which are connected strongly to a single other disease is exaggerated in binary graphs compared to hypergraphs. The measure of centrality is derived from the weighting metrics calculated for disease sets and further investigation is needed to better understand the effect of the metric used in identifying the clinical significance and ranked centrality of grouped diseases. These initial results indicate that hypergraphs can be used as a valuable tool for analysing previously poorly understood relationships and information available in EHR data.

Cross-phenotype association mapping of the MHC identifies genetic variants that differentiate psoriatic arthritis from psoriasis.

OBJECTIVES: Psoriatic arthritis (PsA) is a chronic inflammatory arthritis, with a strong heritable component, affecting patients with psoriasis. Here we attempt to identify genetic variants within the major histocompatibility complex (MHC) that differentiate patients with PsA from patients with cutaneous psoriasis alone (PsC). METHODS: 2808 patients with PsC, 1945 patients with PsA and 8920 population controls were genotyped. We imputed SNPs, amino acids and classical HLA alleles across the MHC and tested for association with PsA compared to population controls and the PsC patient group. In addition we investigated the impact of the age of disease onset on associations. RESULTS: HLA-C*06:02 was protective of PsA compared to PsC (p=9.57×10-66, OR 0.37). The HLA-C*06:02 risk allele was associated with a younger age of psoriasis onset in all patients (p=1.01×10-59). After controlling for the age of psoriasis onset no association of PsA to HLA-C*06:02 (p=0.07) was observed; instead, the most significant association was to amino acid at position 97 of HLA-B (p=1.54×10-9) where the presence of asparagine or serine residue increased PsA risk. Asparagine at position 97 of HLA-B defines the HLA-B*27 alleles. CONCLUSIONS: By controlling for the age of psoriasis onset, we show, for the first time, that HLA-C*06:02 is not associated with PsA and that amino acid position 97 of HLA-B differentiates PsA from PsC. This amino acid also represents the largest genetic effect for ankylosing spondylitis, thereby refining the genetic overlap of these two spondyloarthropathies. Correcting for bias has important implications for cross-phenotype genetic studies.