Back

Metabolism

Elsevier BV

Preprints posted in the last 7 days, ranked by how well they match Metabolism's content profile, based on 15 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit.

1
Discordant associations of IGF-binding proteins 1 & 2 with diabetes and cardiovascular disease: insights from UK Biobank

Rolfe-Hammerton, E. R.; Conning-Rowland, M. S.; De Faveri, L. E.; Simmons, K. J.; Meakin, P. J.; Cubbon, R. M.; Wheatcroft, S. B.

2026-07-20 endocrinology 10.64898/2026.07.17.26358347 medRxiv
Top 0.1%
3.1%
Show abstract

The insulin-like growth factor (IGF)/IGF-binding protein (IGFBP) axis has been implicated in diabetes mellitus and the associated burden of cardiovascular complications. Higher circulating levels of IGFBP-1 and IGFBP-2 have been established as markers of protection from incident type 2 diabetes, yet their associations with cardiovascular disease remain unclear. Utilising the UK Biobank (UKB) resource to integrate disease outcomes, plasma proteomics and MRI data, we examined associations of IGFBP-1 and IGFBP-2 with incident diabetes and cardiovascular disease. Approximately 50,000 UKB participants with plasma proteomic measurements for IGFBP-1 and IGFBP-2 were included. Multivariate Cox regression models revealed that participants in the highest quartiles of IGFBP-1 and IGFBP-2 had a substantially lower risk of incident diabetes (hazard ratio (HR) = 0.31 and 0.32 respectively), but, paradoxically, had increased risks of incident macrovascular disease, all-cause and cardiovascular-related mortality (HR = 1.81 and 2.39). Both proteins were negatively associated with HbA1c levels, triglyceride/HDL ratio and abdominal adiposity, yet positively associated with NT-proBNP, troponin I, cardiac chamber size and aortic dimensions. In summary, negative associations of IGFBP-1 and IGFBP-2 with incident diabetes mellitus did not translate to a reduced cardiovascular risk, suggesting potentially complex actions of IGFBP-1 and IGFBP-2 in the pathophysiology of cardiometabolic disease.

2
Hippocampal Volume Predicts Unhealthy Food-Seeking Trajectories in Insulin-Resistant, but Not Insulin-Sensitive, Youth with Obesity and Depression

KHODAYARI, N.; Branchini, J.; Zhao, M.; Valenzuela, R. J. F.; Springs, Z. A.; Khanna, M.; Patron, D.; Singh, M. K.

2026-07-16 endocrinology 10.64898/2026.07.14.26357902 medRxiv
Top 0.3%
0.9%
Show abstract

Insulin resistance, an often-untreated precursor of type 2 diabetes mellitus (T2DM), is implicated in cognitive decline in adults, yet its impact on the developing brain in youth with obesity remains poorly understood. We investigate whether insulin resistance moderates the relation between hippocampal volume and unhealthy food-seeking in overweight and depressed youth ages 9-17 who completed an oral glucose tolerance test and a cognitive task assessing unhealthy food-seeking motivation at baseline, 6-, and 24-months follow-up, and structural MRI at baseline and 6-months follow-up. Insulin sensitivity moderated this relation: smaller baseline hippocampal subfield volumes predicted increased unhealthy food-seeking over 24 months (ps<0.05). Categorical grouping revealed subfield CA2/3 and 4 volumes predicted this relation among insulin-resistant (ps<0.05), but not insulin-sensitive (ps>0.10), youth, suggesting that threshold criteria for insulin resistance are physiologically meaningful. These findings identify a neuro-metabolic risk phenotype that precedes T2DM and may accelerate unhealthy food-seeking severity in youth with obesity.

3
Histidine-rich glycoprotein is not associated with thrombosis in a UK Biobank Mendelian randomisation analysis

Duan, Y.; Aitken-Buck, H. M.; MacCallum, P.; Weitz, J. I.; Kakkar, A. K.; Allen, A. S.

2026-07-21 cardiovascular medicine 10.64898/2026.07.20.26358305 medRxiv
Top 0.5%
0.6%
Show abstract

Background: Histidine-rich glycoprotein (HRG) is a regulator of coagulation that has been linked to experimental thrombosis, but its causal role in human thrombotic disease remains unclear. Objectives: To determine whether HRG has a role in thrombosis, we evaluated the association between genetically determined HRG variation and thrombosis risk using Mendelian randomisation (MR). Methods: We performed a two sample MR analysis in UK Biobank using non-overlapping samples. Separate genome wide association studies (GWAS) were conducted to identify single nucleotide polymorphisms (SNPs) associated with circulating HRG protein levels and to estimate SNP associations with thrombosis outcomes. Genetic instruments were derived from the HRG GWAS measurements and applied to assess associations with overall, venous, and arterial thrombosis. Sensitivity analyses using multiple MR methods were undertaken, alongside adjusted logistic regression models in participants with measured HRG levels. Results: Among 30,680 participants with HRG measurements, GWAS identified multiple loci associated with HRG levels, with the strongest signal at the rs9898 SNP ({beta} =0.52; P=1.1x10-306). Single instrument MR found no association between HRG protein levels predicted by the rs9898 SNP and risk of overall thrombosis ({beta} =-0.00660; P=0.622), venous thrombosis ({beta} =0.0101; P=0.650) and arterial thrombosis ({beta} =-0.0137; P=0.384). Similar null findings were observed using multi-instrument MR approaches. In complementary analyses, measured HRG levels were not associated with thrombosis after adjustment for age, sex, and C reactive protein, and results were unchanged after stratification by rs9898 genotype. Conclusion: Genetically determined variation in HRG is not associated with thrombotic risk, indicating that HRG related coagulation phenotypes do not translate into clinically meaningful thrombosis.

4
Shared genetic and molecular architecture between insulin resistance and cognitive performance

Martone, A.; Roth Mota, N.; Sakic, B.; Klein, M.; Franke, B.; Fanelli, G.; Bralten, J.

2026-07-16 psychiatry and clinical psychology 10.64898/2026.07.15.26358124 medRxiv
Top 0.9%
0.3%
Show abstract

Insulin signalling contributes to neurodevelopment and brain function, and insulin resistance (IR)-related traits are associated with cognitive performance. However, the genetic architecture shared across specific cognitive domains and IR-related phenotypes remains insufficiently defined. We analysed large-scale genome-wide association study summary statistics for 11 IR-related traits (N=53,334-933,970) and 10 cognitive measures (N=28,156-436,853) to quantify global and local genetic correlations, fine-map shared association signals, and annotate implicated genes and drug-gene interactions. Pairwise global and local genetic correlations were estimated, and shared high-confidence variants were prioritised using the multivariate Sum of Single Effects model. Positional and expression quantitative trait locus mapping was performed, and implicated genes were examined through functional annotation, tissue enrichment, and drug-gene interaction analyses. Low-to-moderate genetic correlations were observed between six IR-related traits and seven cognitive measures (|rg|=0.08-0.34), with predominantly opposite directions, except for correlations involving visual declarative short-term memory. Local genetic correlations showed mixed effect directions across most trait pairs, and multivariate fine-mapping prioritised 696 shared likely causal variants with high posterior support. Gene annotation indicated enrichment in several pathways, including immune-related, signal transduction, neurogenesis, neurotransmitter metabolism, receptor regulation, and lipid and cholesterol metabolism regulation. Implicated genes were expressed across various brain regions and showed prior associations with neuropsychiatric and cardiometabolic conditions. Several drug-gene interactions were identified, involving immunomodulatory and anti-inflammatory compounds. These findings indicate widespread heterogeneous genetic overlap between IR-related traits, particularly body mass index and waist-to-hip ratio, and cognitive measures of general intelligence, processing speed, and short-term visual declarative memory. The findings prioritise apolipoprotein-related lipid transport and inflammatory and oxidative stress pathways as candidate mechanisms linking cognitive, cardiometabolic, and neuropsychiatric phenotypes.

5
Real-Time Glycaemic and Metabolic Adaptation During Unsupplemented Spiritual Fasting up to 30 Days: A Self-Controlled Observational Study

Prakash, S.; Shekhawat, N.; Bardiya, O.; Garg, R.; Tripathi, V.; Fialoke, S.

2026-07-21 endocrinology 10.64898/2026.07.20.26358472 medRxiv
Top 1%
0.2%
Show abstract

Background: Prolonged unsupplemented spiritual fasting (USF), complete caloric abstinence motivated by spiritual practice, is undertaken by many communities worldwide, yet its physiological consequences remain poorly characterized. No prior study has documented continuous real-time monitoring during free-living fasting beyond 10 days. Methods: We conducted a self-controlled observational study of 23 experienced Jain practitioners undergoing USF. Seventeen completed at least 8 days of fasting (11 completed eight days, six continued to 30 days); six discontinued early. Continuous glucose monitoring (CGM) and blood biomarkers at baseline (Day-0), post-fasting (Day-9), and 60-day follow-up (Day-69) were obtained. Mood was assessed daily using PANAS. Within-participant comparisons used both parametric and non-parametric t-tests. Results: CGM revealed near-complete suppression of glycaemic variability within 24-48 hours of fasting onset, sustained throughout with no clinical hypoglycaemia in either cohort. Blood biomarkers showed transient perturbation during fasting -- including rises in hepatic enzymes, bilirubin, uric acid, creatinine, and lipid fractions -- broadly reversible by follow-up (Day-69). hsCRP rose during fasting then fell below baseline at follow-up (5.52{+/-}13.67 to 3.68{+/-}13.18 mg/L, p=0.034). HDL significantly rose above baseline (45.71{+/-}112.32 to 48.97{+/-}111.19, p=0.045) and LDL similarly declined below baseline (122.33{+/-}132.10 to 106.96{+/-}131.62 mg/dL, p=0.035); and then both significantly improved by follow-up. Thyroid axis suppression fully normalized by follow-up. Psychological wellbeing was maintained throughout. Conclusions: Extended USF produces a safely reversible pattern of acute physiological adaptation with net cardiometabolic benefit. Absence of clinical hypoglycaemia during 30-day water-only USF, documented here for the first time with CGM, provides empirical grounding for future controlled trials.

6
Epidemiology, temporal trends, and fibrosis risk stratification in metabolic dysfunction-associated steatotic liver disease in UK primary care: a population-based cohort and nested case-control study

Huang, H.-T.; Hewitt, M.; Li, W.; Temperley, L.; Sattar, N.; Alazawi, W.

2026-07-16 epidemiology 10.64898/2026.07.15.26358136 medRxiv
Top 2%
0.1%
Show abstract

Background: We sought to determine the changing prevalence, incidence, and temporal trends in real-world, recorded diagnoses of metabolic dysfunction-associated steatotic liver disease (MASLD) and assess availability of fibrosis risk stratification following awareness campaigns and guideline updates over the last decade. Methods: This population-based cohort study identified MASLD diagnoses made between 2003 to 2022 in the UK primary-care Clinical Practice Research Datalink (CPRD) to estimate prevalence and incidence. A nested case-control analysis, utilising 1:4 age-, sex-, and general practice-matched controls, assessed clinical characteristics, availability of Fibrosis-4 (Fib-4) components, and its temporal trend pre- and post-2015. Findings: 11.7 million individuals were active in CPRD in 2022. 365,797 comprised the study cohort of people with a MASLD diagnosis (matched to 1,460,288 controls). From 2012-2022, recorded MASLD prevalence rose from 0.52% [N=51,028] to 2.42% [N=283,762] (p<0.001); recorded incidence doubled from 1.60 to 3.31 per 1000 person-years (p<0.001). People with MASLD diagnosis had a higher prevalence of type 2 diabetes (21.0% [N=76,640] vs 7.7% [N=112,812]) and hypertension (35.3% [N=129,156] vs 18.7% [N=273,502]). People of South Asian ethnicity were overrepresented in MASLD cohort but had the lowest availability of Fib-4 components (14.6% [N=4,467]; adjusted odds ratio 0.67, 95% CI: 0.65-0.70, vs White). Overall, Fib-4 availability increased pre- to post-2015 (4.3% [N=4,945] to 22.8% [N=56,634]). Among those with a calculable score, fewer South Asian individuals had indeterminate/high risk (18.6% [N=833] vs 35.3% [N=15,115] in White individuals, p<0.001). Interpretation: Recorded MASLD prevalence has increased 5-fold in a decade, yet a diagnostic gap persists. Fibrosis risk stratification has improved, but remains low and is potentially inequitable for people of South Asian ethnicity. Funding: Merck Sharp & Dohme LLC, a subsidiary of Merck & Co., Inc., Rahway, NJ, USA; Barts Charity. Keywords: Epidemiology, Real-world data, Real-world evidence, MASLD, Primary Care

7
Early identification of suboptimal responders to metformin in type 2 diabetes using long-term real-world HbA1c trajectories

Yang, E.; Riselli, A.; Xu, F.; Sridhar, S. B.; Kvale, M.; Giacomini, K. M.; Hedderson, M. M.; Yee, S. W.; Savic, R. M.

2026-07-20 endocrinology 10.64898/2026.07.17.26357984 medRxiv
Top 2%
0.1%
Show abstract

Aims Metformin remains the primary treatment for type 2 diabetes, yet over 40% of patients fail to maintain glycaemic control. We aimed to identify patients unlikely to respond to metformin prior to treatment initiation and to evaluate whether on-treatment management can improve glycaemic outcomes in suboptimal responders, informing early treatment decisions. Materials and Methods We analyzed 59,881 longitudinal HbA1c measurements from 7,105 patients with type 2 diabetes receiving metformin monotherapy using real-world electronic health records from Kaiser Permanente Northern California with up to six years of follow-up. We integrated demographic, clinical, genetic, and pharmacological factors to characterize metformin responder phenotypes and quantify the impact of adherence and weight control on time to glycaemic failure. Results Three distinct trajectory-based phenotypes were identified: good (63.6%), poor (8.9%), and non-responders (27.5%). Poor responders initially achieved glycaemic targets but lost control within 2.5 years, while non-responders showed minimal HbA1c reduction and failed within 1 year. Five baseline factors-HbA1c, age at diagnosis, body mass index, sex, and estimated glomerular filtration rate-classified phenotypes with good discrimination (area under the receiver operating characteristic curve = 0.84). Incorporating on-treatment HbA1c further enhanced identification of non-responders. Among suboptimal responders, weight control and improved adherence delayed glycaemic failure by approximately 7 months; however, eventual glycaemic failure remained likely. Conclusions We characterized three clinically relevant metformin responder phenotypes and showed that suboptimal responders can be identified early using baseline features. Poor and non-responders are unlikely to achieve durable glycaemic control with metformin alone and may require alternative treatment strategies.

8
Longitudinal multiomic network rewiring at the complement coagulation interface in post-acute sequelae of COVID 19 (PASC)

Ward, B.; Belkhir, L.; Balligand, J.-L.; Cani, P. D.; De Greef, J.; Dewulf, J. P.; Gatto, L.; Haufroid, V.; Kabamba, B.; Vertommen, D.; Yombi, J. C.; Elens, L.; Bommer, G.; Bamps, L.

2026-07-16 infectious diseases 10.64898/2026.07.14.26358048 medRxiv
Top 2%
0.1%
Show abstract

Background. Post acute sequelae of COVID 19 (PASC) is clinically heterogeneous and mechanistically unresolved, and single-analyte studies have struggled to explain it. Methods. We profiled matched plasma proteomics, metabolomics and whole-blood transcriptomics at acute infection and convalescence (mean 86 days later) in a Belgian cohort, using linear mixed models, multiomic gene-set enrichment, and a degree-matched differential-correlation approach to quantify how each node's interactions were rewired between patients who developed PASC and those who recovered; seven axis proteins were additionally quantified by multiplex immunoassay as orthogonal validation. Findings. Single omic testing yielded few FDR significant features, yet multi-omic enrichment showed sustained complement cascade involvement from acute illness to follow-up in PASC. Correlation networks re-organised topologically toward C3 and lost the immunoglobulin V gene coexpression seen in recovery. The most rewired nodes, heparin cofactor II (SERPIND1), alpha 1 antitrypsin (SERPINA1), complement factor H related 5 (CFHR5), prothrombin/thrombin (F2) and immunoglobulin V gene transcripts (notably IGLV3 21), changed in their co-expression structure rather than in abundance. In multiplex validation, acute CRP was elevated in patients who developed PASC (FDR = 0.012), whereas the directly measured abundances of the network-nominated proteins were unchanged. Interpretation. These trajectory aware, cross omic networks nominate a thrombo inflammatory axis in which complement and coagulation regulation remain dysregulated in PASC at the level of wiring rather than abundance, providing a systems framework for validation and for exploring interventions at the complement coagulation platelet interface.

9
Off-Trial: Real-World Weight Loss on Tirzepatide and Semaglutide

Erly, B.; Raja, S.

2026-07-16 pharmacology and therapeutics 10.64898/2026.07.14.26357502 medRxiv
Top 2%
0.1%
Show abstract

Background. GLP-1 receptor agonist trials are tightly controlled: standardized titration, intensive dietary counseling, frequent in-person follow-up, and rigorous exclusion criteria. The real world is none of those things. In a U.S. telehealth GLP-1 program, diet engagement, exercise, medication choice, dose timing, and out-of-pocket cost vary substantially from patient to patient. Whether trial-level efficacy translates into the outcomes a patient and clinician will actually see is an open question, and the answer matters, because telehealth is now where most GLP-1 prescribing happens. Methods. We conducted a retrospective cohort study of 13,507 adults who used a single GLP-1 agent (tirzepatide or semaglutide) through the Mochi Health telehealth obesity program and had a documented six-month weight observation. The primary outcome was achievement of >=10% total body weight loss at six months. To address selection bias in the tirzepatide-semaglutide comparison, we used 1:1 nearest-neighbor propensity-score matching on age, sex, baseline BMI, baseline weight, and comorbid diabetes, hypertension, dyslipidemia, and prior bariatric surgery (recorded at intake), with a 0.25 SD caliper on the propensity logit. We drew a directed acyclic graph (DAG) with a clinical co-author to make the identifying assumptions explicit and to mark where unobserved variables (insurance, socioeconomic status, concomitant medications such as metformin) limit causal interpretation. We report multivariable predictors via logistic regression, compute an E-value for the matched contrast, and benchmark our point estimates against landmark RCT outcomes. Results. Overall, 59.1% of patients achieved >=10% loss at six months, with mean loss of 11.5% (median 11.3%). Threshold attainment was 86.6% at >=5%, 59.1% at >=10%, 27.5% at >=15%, and 9.1% at >=20%. The unadjusted tirzepatide-semaglutide response gap was +16.0 percentage points (68.8% vs 52.8%); after 1:1 propensity-score matching (3,480 pairs, all post-match |SMD| < 0.05) the gap was +18.1 percentage points (69.6% vs 51.6%, 95% CI +15.9 to +20.3). Matching on the measured covariates did not attenuate the advantage, indicating that selection on those characteristics does not explain it; the matched risk ratio was 1.35 (E-value 2.04). The gap was unchanged when a self-reported insurance indicator was added to the matching (+18.4 pp) and remained large (+14.2 pp) within patients who reached a therapeutic dose. Multivariable predictors of response were tirzepatide (OR 2.10, 1.95-2.26), female sex (OR 1.37, 1.20-1.56), and prior bariatric surgery (OR 1.36, 1.18-1.57); response was lower with comorbid diabetes (OR 0.84, 0.77-0.92) and, modestly, with higher baseline BMI per unit (OR 0.98, 0.97-0.99). Response varied by baseline BMI, from 58.0% in overweight patients (BMI <30) and a peak of 63.5% in Obese I to 52.4% in Obese III. Conclusions. Real-world response to GLP-1 therapy in a telehealth setting is meaningfully attenuated from RCT benchmarks but remains clinically substantial: roughly three in five patients reach the 10% threshold. The tirzepatide advantage over semaglutide is large and, notably, does not shrink under propensity-score matching on measured confounders, so it is not an artifact of the observed selection variables; an unmeasured confounder would need a risk-ratio association of about 2.0 with both drug choice and response to explain it away (E-value 2.04). The findings are observational, conditional on the DAG's identifying assumptions, and unmeasured confounders (insurance, socioeconomic status, concomitant medications) remain possible.

10
Proteogenomic mapping of multimorbidity identifies C1R linking coronary artery disease and dementia

Li, L.; Tang, Z.; Zhong, Z.; Geng, T.; Guo, Y.; Liao, Y.; Demirkan, A.; Bowden, J.; Bragg, F.; Pan, A.; Sun, X.; Liu, J.; Liu, G.; Liu, J.

2026-07-16 genetic and genomic medicine 10.64898/2026.07.14.26358022 medRxiv
Top 2%
0.1%
Show abstract

Multimorbidity is highly prevalent in ageing populations, yet its shared molecular basis remains poorly defined, limiting the development of therapies that target multiple conditions. We systematically integrated measurements of 1,954 circulating proteins from 54,219 individuals in discovery and 35,559 in replication, focusing on ten common age-related diseases: coronary artery disease, chronic kidney disease, chronic obstructive pulmonary disease, dementia, heart failure, major depressive disorder, osteoarthritis, Parkinson's disease, stroke, and type 2 diabetes. Coronary artery disease emerged as a central condition in the multimorbidity network, sharing circulating protein signatures with seven other diseases. Through genetic causal-inference analyses, we identified 40 circulating proteins with cross-disease relevance, of which four were further supported by colocalization of genetic variant associations. Among these, complement C1r, encoded by C1R, emerged as a key link between coronary artery disease and dementia, supported by independent colocalization evidence (PP.H4 = 0.86). Phenome-wide association analyses of C1R variants suggested that this signal was not driven by widespread unrelated genetic effects, but instead may reflect a more specific contribution to coronary artery disease-dementia pathogenesis. In vitro experiments further suggested that fibroblast-derived C1R promotes endothelial inflammation and neuronal apoptosis, providing mechanistic plausibility. Together, these findings position C1R as a biologically plausible and therapeutically relevant molecular link between coronary artery disease and dementia.

11
An ancestry-matched Mendelian randomisation analysis of kidney function and heart failure subtypes in African ancestry populations

Gaye, N. D.; Diawara, A.

2026-07-17 genetic and genomic medicine 10.64898/2026.07.15.26358145 medRxiv
Top 2%
0.1%
Show abstract

Chronic kidney disease and heart failure disproportionately burden populations of African ancestry, yet Mendelian randomisation (MR) studies of the causal relationship between kidney function and heart failure subtypes have been conducted exclusively in European ancestry populations. We performed a forward two-sample MR analysis to evaluate the causal effect of genetically predicted estimated glomerular filtration rate (eGFR) on heart failure with preserved ejection fraction (HFpEF) and heart failure with reduced ejection fraction (HFrEF) in individuals of African ancestry. Genetic instruments were selected from an African ancestry eGFR genome-wide association study (N = 67,943) at genome-wide significance, with linkage disequilibrium clumping using an African ancestry reference panel. Heart failure subtype summary statistics were obtained from the Million Veteran Program (HFpEF: 5,379 cases / 113,041 controls; HFrEF: 9,104 cases / 109,632 controls). Six independent SNPs (F-statistics 30.5 &#8211 107.3; R&#178 = 0.62%) were retained as instruments. The primary inverse-variance weighted analysis provided no evidence of a causal effect of eGFR on HFpEF (OR 0.92, 95% CI 0.80 &#8211 1.06, p = 0.248) or HFrEF (OR 0.98, 95% CI 0.78 &#8211 1.23, p = 0.878). Sensitivity analyses were directionally consistent. There was no evidence of heterogeneity or directional pleiotropy. Minimum detectable effects at 80% power were OR 1.28 for HFpEF and OR 1.22 for HFrEF. These null findings should be interpreted as inconclusive given current power constraints; larger ancestry-matched studies are needed.

12
Neonatal admission as a marker of risk for poor educational attainment and special educational needs in children aged 5-11 years

John, A.; Pike, C.; Olga, L.; Sovio, U.; Wong, H. S.; Smith, G. C.; Aiken, C.

2026-07-17 pediatrics 10.64898/2026.07.15.26358132 medRxiv
Top 3%
0.0%
Show abstract

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.

13
Bridging surveillance gaps in dengue: a hierarchical model integrating mixed data sources for transmission estimation and vaccine targeting

Djaafara, B. A.; Elyazar, I. R.; Yosephine, P.; Surya, A.; Silalahi, F. S.; Handito, A.; Thohir, B.; Aryani, D.; Gunawan, D.; Nisa, A. K.; Prianto, E.; Samad, I.; Cook, A. R.; Huang, A. T.; Clapham, H. E.; Bhatt, S.; Mishra, S.

2026-07-17 epidemiology 10.64898/2026.07.15.26358208 medRxiv
Top 3%
0.0%
Show abstract

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

14
Efficient stochastic epidemic simulation via the Sellke construction

van Boven, M.; Bootsma, M. C.

2026-07-17 epidemiology 10.64898/2026.07.16.26358219 medRxiv
Top 3%
0.0%
Show abstract

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

15
General Practice Perspectives on Post-Infection Conditions: Scoping Review and UK Survey

Aung, K. W.; Scuffell, J.; Podlasek, A.; Engamba, S.; Jones, F.; Edwards, A.; Chew-Graham, C. A.; Sanyaolu, L.; Busse-Morris, M.

2026-07-17 primary care research 10.64898/2026.07.15.26358157 medRxiv
Top 3%
0.0%
Show abstract

Background Post-infection conditions (PICs), such as Long Covid, are associated with heterogeneous, fluctuating symptoms that profoundly affect daily functioning. Despite moderate-certainty evidence from the NIHR-funded LISTEN trial (COV-LT2-0009) that personalised self management support improves outcomes and may reduce societal and economic impacts of Long Covid, many people living with PICs still receive condition-specific services, generic advice, or stand-alone digital tools that do not address their complex needs. Aim To map care approaches in general practice and synthesise UK evidence for PIC management. Design and setting Scoping review and online survey. Method A two-phase study was conducted: (1) a scoping review of UK evidence on PIC management in general practice; and (2) a supplementary online survey of practitioners working in UK general practice to provide contextual insights. Results The scoping review identified 32 studies focused on Long Covid. One study included a comparator group (ME/CFS). Study populations were predominantly white ethnicity and female. Evidence for non-Covid PICs in UK general practice was largely absent. The supplementary survey (n=46) provided preliminary practice-level insights. Healthcare practitioners reported varied PIC presentations, diagnostic uncertainty, limited referral pathways, inequitable access, and low confidence in managing PICs. Conclusion Evidence informing PIC management in UK general practice remains predominantly Long Covid-focused and may not reflect the range of PICs encountered in practice. While survey findings are preliminary and require confirmation in larger samples, they highlight uncertainty around PIC management. Further research is needed to evaluate whether existing Long Covid pathways should be expanded or complemented by broader PIC models. Keywords general practice; Long Covid; self-management; post-viral syndromes

16
Complex intra-host SARS-CoV-2 evolution following monoclonal antibody pre-exposure prophylaxis

Kamelian, K.; Pascall, D. J.; Cheng, M. T. K.; Meng, B.; Altaf, M.; Morse, R. M.; Aggio, J. B.; Egan, D. J. S.; Chen-Xu, M.; Trivioli, G.; Sutton, B.; Richter, A.; Gonzalez-Vazquez, L. D.; Cormie, C.; Kemp, S.; Yeadon, R.; Hyatt, B.; Wong, A.; Thesin Pelamkulangara, N.; Fraser, E.; McCarthy, B.; Novaes, F.; Stott, S.; Galvin, A.; Bellis, K. L.; De Angelis, D.; Harrison, E. M.; Martin, D.; Smith, R. M.; Gupta, R. K.

2026-07-17 infectious diseases 10.64898/2026.07.14.26356329 medRxiv
Top 3%
0.0%
Show abstract

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

17
Nationwide Mpox Genomic Surveillance Reveals Clade Ib Introductions, APOBEC3-Driven Evolution, and Terminal Deletions

Brochu, H. N.; Shi, Q.; Song, K.; Zhang, Q.; Munroe, J.; Harris, N. J.; Britt, N.; Zeng, Q.; Kapuria, K.; Chappell, J.; Norvell, B. M.; Peavy, L.; Williams, J. D.; Harris, A. B.; Chaitram, J.; Hutson, C. L.; Deng, J.; McGrath, D.; Boles, D.; Dale, S. E.; Gigante, C. M.; Iyer, L. K.

2026-07-17 infectious diseases 10.64898/2026.07.15.26357894 medRxiv
Top 3%
0.0%
Show abstract

Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.

18
Temporal relationships between distress and pain in people living with HIV

Arendse, G.; Kamerman, P.; Wadley, A.; Edwards, R. R.; Joska, J.; Parker, R.; Madden, V. J.

2026-07-17 primary care research 10.64898/2026.07.15.26358133 medRxiv
Top 3%
0.0%
Show abstract

Objective: There is a bidirectional relationship between emotional distress and pain. However, this relationship is understudied in people with HIV in low-resource settings. This study sought to describe the temporal relationship between emotional distress and pain in people with HIV. Design: Longitudinal observational study. Methods: Participants with virally suppressed HIV, reporting either no pain or persistent pain at baseline, provided weekly remote ratings of distress, worst pain, and average pain using 0-10 visual analogue scales. Within-individual fluctuations in distress and pain were visualised over time. Group-level correlations were determined using Spearman's correlation tests. Cumulative link mixed models assessed whether distress and pain each predicted the other in the following week. Results: 72 participants provided responses over 49 weeks. The participants had a median (IQR) age of 43 (37-51) years, 63% (n=45) were unemployed and most were females (n=51;71%). Distress and pain fluctuated concurrently within individuals: distress was positively correlated with worst pain ({rho}=0.66, 95% CI= 0.60-0.72, p<0.001) and average pain ({rho}=0.70, 95% CI=0.64-0.75, p<0.001) intensity within the same week. Worst pain (OR=1.42, 95% CI=1.17-1.71, p<0.001) and average pain (OR=1.43, 95% CI=1.20-1.71, p<0.001) intensity both predicted distress in the next week. Distress predicted worst pain intensity (OR=1.25, 95% CI=1.07-1.46, p=0.023) but not average pain intensity (OR=1.19, 95% CI=1.01-1.40, p=0.152) in the next week. Conclusions: The temporal relationship between distress and worst pain intensity was bidirectional, whereas distress did not temporally predict average pain intensity. Both pain and emotional distress should receive attention from HIV research and clinical care in low-resource settings.

19
Comparative Efficacy of Vancomycin and Fidaxomicin Regimens for the Prevention of Recurrent Clostridioides difficile Infection: A Systematic Review and Network Meta-Analysis of Randomized Controlled Trials

Prosty, C.; Butler-Laporte, G.; Brophy, J.; Frenette, C.; Loo, V.; Coburn, B.; Hota, S.; Longtin, Y.; Kong, L.; Muller, M.; Steiner, T.; Valiquette, L.; Daneman, N.; Daley, P.; Nott, C.; MacFadden, D. R.; Kandel, C.; Chen, Y.; Perez- Patrigeon, S.; Lee, T. C.; McDonald, E.

2026-07-17 infectious diseases 10.64898/2026.07.14.26358112 medRxiv
Top 3%
0.0%
Show abstract

Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.

20
Rationale and guidance for implementing the continual reassessment method for dose-finding in controlled human infection model studies

Weerasinghe, C.; Osowicki, J.; Simpson, J. A.; Crocker-Buque, T.; McCarthy, J.; Williams, E.; Price, D. J.

2026-07-17 infectious diseases 10.64898/2026.07.16.26358128 medRxiv
Top 3%
0.0%
Show abstract

Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.