Metabolomics

Prenatal alcohol exposure (PAE) is a major cause of neurodevelopmental disorders, yet most children are diagnosed late or misdiagnosed. Neuroplacentology suggest that placental factors released into maternal and/or umbilical cord blood contribute to fetal brain development. Consistently, a preclinical inter-organ transcriptomic database revealed that PAE disrupts the expression ratio of angiogenic and inflammatory factors suggesting an angio-inflammatory response. This study aimed i) to assay, by multiplex immunoassay, angiogenic and inflammatory factors in maternal and umbilical cord blood from alcohol-consuming women and ii) to perform a maternofetal analysis according to neonatal sex. Afterwards, dysregulated factors from mothers who gave birth to females or males were submitted to STRING and ShinyGO analyses. Results showed that PAE differently altered the distribution profiles of dysregulated angiogenic and inflammatory factors in maternal and umbilical cord blood. Moreover, sex-specific differences were observed, with 36% of dysregulated proteins specific to males, 48% to females, and 16% common to both. STRING analysis revealed robust functional protein-protein interactions linking together inflammatory and angiogenic clusters while the ShinyGO analysis identified enriched pathways related to vascular shear stress. These findings provide the first maternofetal analysis of combined angiogenic and inflammatory factors from alcohol-consuming mothers.

Background: Early detection of pancreatic ductal adenocarcinoma (PDAC) is critical. While faecal elastase-1 (FE-1) is a standard clinical marker for pancreatic function, its diagnostic accuracy for malignancy is limited. We sought to identify plasma metabolites that enhance FE-1 performance in symptomatic "at-risk" patients. Methods: Using the DEPEND cohort (CRUK C45617/A29908), plasma metabolomics was performed on patients with resectable PDAC (n=23) and healthy volunteers (n=24). Predictive modelling included feature selection and cross-validation, with further validation in an independent external cohort. Results: Citrulline was identified as significantly depleted in PDAC patients across discovery and validation cohorts. In isolation, Citrulline achieved an AUC of 0.86 (internal) and 0.88 (external validation). Standalone FE-1 demonstrated an AUC of 0.67. However, combining Citrulline and FE-1 significantly improved diagnostic performance, achieving a combined AUC of 0.96. Stratification revealed distinct metabolomic signatures associated with poorly differentiated tumours, suggesting a link to histological grade. Conclusions: Integrating Citrulline with FE-1 testing substantially improves PDAC detection in symptomatic patients. This non-invasive panel offers high diagnostic potential, though prospective validation is required to establish clinical cut-offs for routine practice.

BACKGROUND In line with International Federation of Clinical Chemistry and Laboratory Medicine (IFCC) recommendations for high-sensitivity cardiac troponin assays, analytical validation and reference limit assessments are required to confirm that an assay meets performance criteria. This study evaluated the analytical performance and established the 99th percentile upper reference limit (URL) for the SPINCHIP High-Sensitivity Cardiac Troponin I (SPINCHIP hs-cTnI) point-of-care assay. METHODS Analytical performance characteristics, including the limit of blank (LoB), limit of detection (LoD), and limit of quantification (LoQ), were assessed. Additionally, 1,053 plasma samples and 1,055 whole-blood samples were used to determine the URL. Imprecision around the 99th percentile URL was evaluated as part of the analytical validation. High-sensitivity criteria were assessed by confirming measurable cTnI in [&ge;]50% of healthy individuals (n=432 plasma; n=431 whole blood) and achieving imprecision <10% at the 99th percentile (plasma, n=960; whole blood, n=480). RESULTS SPINCHIP hs-cTnI demonstrated a LoB of 0.3 ng/L; LoDs of 0.8 ng/L (plasma) and 0.9 ng/L (whole blood); and LoQs of 1.1 ng/L (plasma) and 1.4 ng/L (whole blood). The analytical measuring range was 1.1-9,000 ng/L. Imprecision at the common 99th percentile URL (14 ng/L) was 5.8%; for men (URL=16 ng/L) 5.6% and for women (URL=10 ng/L) 6.3%. Greater than 85.2% (94.0% and 76.1% in men and women, respectively) of healthy individuals showed measurable cTnI above the LoD. CONCLUSIONS The SPINCHIP hs-cTnI assay meets the IFCC high-sensitivity requirements, demonstrating <10% imprecision at the 99th percentile, reliable low-concentration precision and cTnI detection in more than half of healthy individuals.

Background: Plasma phosphorylated tau-217 (pTau-217) measurements for use in Alzheimer disease (AD) identification require thresholds to define positivity and there exist different approaches to operationally defining the boundary. We compared amyloid {beta} (AB) PET-anchored and distribution-based positivity cut-off values and explored how these mapped onto latent biomarker states. Methods: We analysed plasma pTau-217 measured in the Bio-Hermes-001 cohort (N = 990) using an immunoassay (Lilly) and mass spectrometry assay (University of Gothenburg). Gaussian mixture models were used to identify latent classes and thresholds were derived in two ways: achieving 90% specificity for AB PET positivity and exceeding the mean + 2SDs of the lowest latent class. We explore classes in reference to AB PET status and clinical diagnosis, as well as agreement between approaches using Cohen kappa for both assays. Results: In both assays, three latent biomarker classes were identified with monotonic increases in AD clinical diagnosis and AB PET positivity. PET-anchored thresholds showed lower specificity but higher sensitivity to amyloid positivity than distribution-based thresholds. Overall agreement between the approaches was acceptable (k = 0.678 for Lilly and 0.575 for University of Gothenburg), with disagreement concentrated in the intermediate latent class. Classes with the lowest and highest pTau-217 concentrations were classified consistently using both thresholds Discussion: The two thresholding approaches yielded similar classifications at both the negative and positive tail of the observed biomarker distribution, but classify intermediate concentrations differently. The boundary definition influenced pTau-217 positivity more than the analytical platform itself. Thresholding approaches may capture different pTau-217 biomarker states, therefore such methodological decisions should be grounded in the context of the intended application.

Background: Childhood obesity is a prevalent public health concern, particularly among children from marginalized backgrounds. Epigenomic mechanisms, including DNA methylation (DNAm), offer insight into the biological embedding of metabolic risk, yet protective family-level factors remain understudied. This study examined whether participation in a positive parenting intervention was associated with reduced epigenetic cardiometabolic risk among preschool-aged children. Methods: Participants were 74 children (n = 35 intervention; n = 39 services-as-usual; mean age = 36 months). Using a secondary analysis of a randomized controlled trial, DNAm-derived body mass index (BMI) and anthropometric BMI were assessed at baseline and 12-month follow-up, and parenting practices were observed post-treatment. Results: Children in the intervention group demonstrated significantly lower DNAm BMI at 12 months relative to controls, adjusting for child sex, race and/or ethnicity, and anthropometric BMI. No treatment effect was observed for anthropometric BMI, and DNAm BMI was not associated with anthropometric BMI at 12 months. Although parenting practices improved, they did not mediate intervention effects on DNAm-derived BMI. Conclusions: Parenting interventions may influence biological pathways related to cardiometabolic risk, even before changes in anthropometric outcomes, underscoring the potential of family-centered approaches to promote early metabolic health.

Background. Maximal oxygen uptake (VO2max) is a leading marker of cardiorespiratory fitness and a strong predictor of all-cause mortality. Cardiopulmonary exercise testing (CPET) is the reference method but is resource-intensive, so consumer wearables estimate VO2max from passively collected signals; these estimates compress the fitness range, returning near-correct group averages while ranking individuals poorly. No peer-reviewed validation of a smart-ring VO2max estimate against CPET has been reported, and none in a South Asian cohort. Objective. To validate the Ultrahuman Ring AIR VO2max estimate against laboratory CPET, benchmark it against published prediction equations, and assess its generalization and construct validity. Methods. In a single-site paired ring-CPET cohort (N = 101; mean CPET peak VO2 43.3 mL{middle dot}kg-{superscript 1}{middle dot}min-{superscript 1}, SD 9.9), peak oxygen uptake was measured by treadmill or cycle-ergometer CPET, and the Ultrahuman Ring AIR estimate was computed from passively collected signals using a transparent ensemble based on published equations. Ensemble weights and calibration were selected on an 85-subject development set by an automated search minimizing a composite 5-fold cross-validated error criterion; the locked estimate was evaluated on a 16-subject held-out test set. The calibrated coefficients are proprietary. Agreement was quantified with mean absolute error (MAE), bias, Pearson r, regression slope and Lin's concordance correlation coefficient (CCC; bootstrap 95% CIs), and Bland-Altman limits of agreement. Separately, in 181,133 de-identified Ring users (no CPET reference), construct validity was assessed against ring-measured sleep, continuous glucose monitoring (n = 2,597), and a venous blood panel (n up to 15,203), adjusted for age, sex, and BMI, with lipoprotein(a) as a pre-specified negative control. Reporting followed TRIPOD and STARD. Results. With a self-reported fitness level provided, the estimate agreed with CPET peak VO2 at MAE 4.68 mL{middle dot}kg-{superscript 1}{middle dot}min-{superscript 1} (95% CI 3.93 to 5.49), Pearson r 0.79, CCC 0.79, and slope 0.71. The five published equations were worse on every metric (MAE 6.2 to 10.6, CCC 0.28 to 0.56, slope 0.32 to 0.42), each compressing the fitness range. On the held-out test set (n = 16), agreement held (r 0.84, slope 0.81, MAE essentially unchanged). Without the fitness input, full-cohort MAE was 5.16, still ahead of every published equation. At population scale, higher estimated fitness tracked a healthier profile on measurements the estimate does not use: better ring-measured sleep; higher continuous-glucose time in target range (79.6% versus 61.5%, top versus bottom decile; n = 222 and 399 of 2,597 users); and lower triglycerides, fasting glucose, and HOMA-IR (n up to 15,203 assayed per marker). These associations held after adjustment for age, sex, and BMI, whereas the pre-specified negative control lipoprotein(a) did not separate the deciles. Conclusions. The Ultrahuman Ring AIR VO2max estimate agreed with laboratory CPET substantially better than published prediction equations, held its agreement on held-out subjects, and ordered a large population along independent cardiometabolic gradients consistent with true fitness.

Importance Sleep-disordered breathing (SDB) is common in childhood and is associated with attentional and behavioral impairments despite largely preserved sleep macrostructure and minimal abnormalities in conventional electroencephalographic measures. This discrepancy has contributed to the perception that sleep is relatively preserved in pediatric SDB and has limited understanding of the physiological mechanisms underlying morbidity. Objective To determine whether pediatric SDB is associated with disruption of the regional organization and homeostatic dynamics of slow-wave activity (SWA), a key physiological marker of sleep-dependent neural recovery and development. Design, Setting, and Participants Cross-sectional study of 62 children aged 4 to 12 years who underwent overnight polysomnography with high-density electroencephalography in a laboratory setting. Participants were recruited from clinical referrals and the community, spanning the full spectrum of SDB severity. Exposures SDB severity indexed by hypopnea index (HI), apnea-hypopnea index (AHI), and obstructive apnea index (OAI). Main Outcomes and Measures Regional electroencephalogram-derived SWA (0.5 to 4 Hz) topography and exponential decay parameters derived from frontal and posterior cortical regions. The frontal-to-posterior decay-rate ratio was evaluated as a summary measure of regional sleep homeostasis. Results In children with lower hypopnea index, SWA demonstrated the expected developmental pattern, with posterior predominance in younger children and a progressive shift toward a more balanced anterior-posterior distribution with age. Increasing HI was associated with attenuation or reversal of this spatial organization. Global SWA showed no meaningful association with SDB severity. In contrast, regional frontal and posterior decay parameters were strongly associated with HI (adjusted R2 = 0.53; p < 1e-6) but not OAI (adjusted R2 = 0.05; p = .95). The frontal-to-posterior decay-rate ratio showed the strongest association with HI {beta} = 4.15; 95% CI, 3.17-5.13; p < 1e-10; adjusted R2 = 0.55. Conclusions and Relevance Pediatric SDB was associated with regional disruption of slow-wave sleep homeostasis rather than global loss of deep sleep. These alterations affected both the spatial organization and temporal dynamics of SWA during a period of active cortical maturation and were not captured by conventional sleep metrics. Regional SWA dynamics may provide a developmentally sensitive marker of physiological disease burden in children with SDB.

Background: Lead exposure is estimated to cause approximately 3.5 million premature deaths a year, yet the key ongoing sources of lead exposure are unclear. Methods: We estimated the contribution of dietary lead intake to global blood lead levels (BLLs) for 7-year-old children and 22-year-old adults by applying the All-Ages Lead Model (AALM) to calculate blood lead levels (BLLs) based on 25 total diet studies (TDS) that quantify dietary lead intake across 46 countries. Results: For children, the population-weighted average dietary lead intake in low- and middle-income countries (LMICs) (32.0 g/day) was found to be more than three times higher than in high-income countries (HICs) (9.3 g/day), and more than 10 times higher than the FDA reference level for children (2.2 g/day). The average impact on BLLs for children is estimated to be near 29 g/L in LMICs and near 12 g/L in HICs. Averaged across the TDS data, vegetables (27%) and cereals (24%) were found to contribute the most to dietary lead. Conclusions: While there are limitations associated with biokinetic modelling and the TDS data from LMICs, these results suggest that the contribution of dietary lead intake to global lead exposure is in the region of 40 to 50%, suggesting, in turn, that dietary lead intake is likely a major global driver of lead poisoning. Lead absorbed from the environment into food crops is expected to be the key driver of dietary lead. Current regulatory levels for maximum lead concentrations in foods (0.05-0.3 mg/kg) are out-of-date and may imply a dietary lead intake of 200 g/day, far higher than the FDA reference level (2.2 g/day). Collecting representative TDS data in high lead burden countries should be a priority. Further research is also recommended on upstream lead sources and pathways of lead uptake in plants, driving global food contamination.

Autonomic dysfunction is a common sequela of mild traumatic brain injury (mTBI). Physical activity progression is an integral component of mTBI rehabilitation, particularly in addressing autonomic dysfunction. However, clinicians often rely on point-in-time evaluation of orthostatic and exercise intolerance to guide activity recommendations. Commercially available wearable devices (e.g., Fitbits) provide an opportunity to evaluate heart rate response to activity in a real-world setting. Previous work has used physiologic (heart rate) and activity (step count) data to identify subgroups of adults with stroke that may be used to guide activity recommendations. This method may be useful to subgroup youth post-mTBI to identify those who have abnormal physiologic responses to activity. We aimed to identify subgroups using heart rate and step count data in adolescents presenting for specialty care after diagnosed mTBI. Eighty participants aged 13-18 within six months of mTBI diagnosis were recruited to wear a Fitbit Sense 2. Data from seven days and two nights collected within fourteen days of enrollment were included. A group-based steps per minute (SPM) threshold (25th percentile; 10 SPM) and individualized heart rate threshold (20% heart rate reserve (HRR)) were used to classify each minute of active daytime data into one of four quadrants: SPM>10 & HRR>20% (QI), SPM<10 & HRR>20% (QII), SPM<10 & HRR<20% (QIII), and SPM>10 & HRR<20% (QIV). We used percentage of minutes in each quadrant, mean steps per day, percentage of minutes with zero steps, mean SPM in QI, and resting heart rate in a k-means clustering algorithm to identify subgroups. We evaluated subgroup differences by clustering variables using Kruskal-Wallis tests. Sixty-one participants were included. Three subgroups emerged: Sedentary (n=12), Active (n=23), and Atypically Elevated Heart Rate (AEHR; n=26). Subgroups varied significantly on all clustering variables (p<0.01). The Active subgroup took a high number of steps per day, had lower sedentary time, and had the highest activity intensity (mean SPM in QI). The Sedentary subgroup took fewer steps per day compared to the Active subgroup, had high sedentary time, and showed the highest resting heart rate. The AEHR subgroup took fewer steps per day compared to the Active subgroup and had high sedentary time. The AEHR subgroup also spent a higher percentage of time with an atypically high heart rate response to low levels of activity compared to the other subgroups. Our findings suggest that data from wearable devices can identify subgroups of adolescents with mTBI with distinct physiologic/physical activity profiles, which may ultimately be used to inform personalized activity prescriptions. Future work should aim to understand how the identified subgroups relate to longitudinal outcomes.

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.

Copper (Cu) is an essential metal involved in neurobiological processes including energy metabolism and neurotransmission, yet dysregulated Cu levels may adversely affect brain health and olfactory performance. Although olfactory dysfunction has primarily been studied in older adults and neurodegenerative disease, adolescence is a critical period of brain maturation during which the olfactory system may be particularly vulnerable. This cross-sectional study examined associations between Cu exposure, olfactory bulb (OB) volume, and olfactory performance in 200 adolescents and young adults (64% female; ages 13 - 25) from the Public Health Impact of Metals Exposure cohort. Cu concentrations in blood, urine, hair, and saliva were measured using inductively coupled plasma mass spectrometry. T2-weighted magnetic resonance imaging scans estimated left, right, and total OB volumes using a three-stage deep learning pipeline. Olfactory performance was assessed using the Sniffin Sticks test. Weighted quantile sum regression evaluated associations between a Cu mixture index and OB outcomes, while standard linear regression models assessed individual Cu biomarkers. Models were adjusted for age and sex. A higher Cu index was associated with reduced left (Beta= -0.72, 95% CI [-1.42, -0.02]), right (Beta = -0.79, 95% CI [-1.43, -0.15]), and total OB volume (Beta= -1.55, 95% CI [-2.85, -0.25]), as well as lower odor threshold scores (Beta = -0.23, 95% CI [-0.42, -0.03]). Individual biomarkers were not independently associated with outcomes. These findings suggest that Cu exposure may adversely affect olfactory neurodevelopment during adolescence and highlight the importance of studying environmental exposures relevant to long-term neurological health.

Antimicrobial resistance (AMR) is an ever-increasing threat to population health. Industrial, environmental and societal factors are increasingly recognized as important contributors to AMR within communities. Here, we investigated the spatial distribution of AMR genes (ARGs) across Alberta, Canada and their association with socio-economic, immigration-related, and agro-industrial characteristics using municipal wastewater-based surveillance. We analyzed monthly wastewater metagenomes collected between March 2022 and March 2023 across eleven municipalities, representing 39% of Alberta's population. Integration with census data enabled multivariate analysis, revealing that municipal resistome profiles were strongly structured along income and immigration-related population gradients. ARGs spanning 14 resistance classes exhibited distinct distributional patterns across income and immigration gradients, including contrasting associations among beta-lactam, aminoglycoside, and macrolide-lincosamide-streptogramin ARGs, consistent with heterogeneous selection pressures across sub-populations. These findings demonstrate the capacity of longitudinal wastewater surveillance to identify persistent population-level resistome patterns and highlight the importance of incorporating sociodemographic context into AMR surveillance and mitigation strategies.

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.

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

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.

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.

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.

Lithium has been the gold standard for the treatment and prevention of relapse in bipolar disorder for over 60 years. Guidance from the National Institute for Health and Clinical Excellence states explicitly to 'offer lithium as a first-line, long-term pharmacological treatment for bipolar disorder'. Yet, in the last two decades its use has been in decline with clinicians favouring anticonvulsants or antipsychotics when treating this condition. In this study, we have used three openly available datasets containing prescribing data from primary and secondary care to explore trends in the use of lithium in England, showing both regional and temporal variance between 2015-2024. We have shown that lithium use declined in primary care by 20.9% in the last ten years (2015-2024) and 10.9% overall in the last five years (2019 to 2025). We have also shown how there is some regional variation in the source of lithium for patients, although the vast majority is prescribed in primary care. Further research into clinical behaviour is needed to understand what is driving the decrease in lithium usage, and what barriers and enablers may influence its use across the country.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

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.