Globally resolved surface temperatures since the Last Glacial Maximum.

Climate changes across the past 24,000 years provide key insights into Earth system responses to external forcing. Climate model simulations1,2 and proxy data3-8 have independently allowed for study of this crucial interval; however, they have at times yielded disparate conclusions. Here, we leverage both types of information using paleoclimate data assimilation9,10 to produce the first proxy-constrained, full-field reanalysis of surface temperature change spanning the Last Glacial Maximum to present at 200-year resolution. We demonstrate that temperature variability across the past 24 thousand years was linked to two primary climatic mechanisms: radiative forcing from ice sheets and greenhouse gases; and a superposition of changes in the ocean overturning circulation and seasonal insolation. In contrast with previous proxy-based reconstructions6,7 our results show that global mean temperature has slightly but steadily warmed, by ~0.5 °C, since the early Holocene (around 9 thousand years ago). When compared with recent temperature changes11, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Glacial cooling and climate sensitivity revisited.

The Last Glacial Maximum (LGM), one of the best studied palaeoclimatic intervals, offers an excellent opportunity to investigate how the climate system responds to changes in greenhouse gases and the cryosphere. Previous work has sought to constrain the magnitude and pattern of glacial cooling from palaeothermometers1,2, but the uneven distribution of the proxies, as well as their uncertainties, has challenged the construction of a full-field view of the LGM climate state. Here we combine a large collection of geochemical proxies for sea surface temperature with an isotope-enabled climate model ensemble to produce a field reconstruction of LGM temperatures using data assimilation. The reconstruction is validated with withheld proxies as well as independent ice core and speleothem δ18O measurements. Our assimilated product provides a constraint on global mean LGM cooling of -6.1 degrees Celsius (95 per cent confidence interval: -6.5 to -5.7 degrees Celsius). Given assumptions concerning the radiative forcing of greenhouse gases, ice sheets and mineral dust aerosols, this cooling translates to an equilibrium climate sensitivity of 3.4 degrees Celsius (2.4-4.5 degrees Celsius), a value that is higher than previous LGM-based estimates but consistent with the traditional consensus range of 2-4.5 degrees Celsius3,4.

Predictive Value of Microfilariae-Based Stop-MDA Thresholds After Triple Drug Therapy With IDA Against Lymphatic Filariasis in Treatment-Naive Indian Settings.

Mass drug administration (MDA) of antifilarial drugs is the main strategy for the elimination of lymphatic filariasis (LF). Recent clinical trials indicated that the triple-drug therapy with ivermectin, diethylcarbamazine, and albendazole (IDA) is much more effective against LF than the widely used two-drug combinations (albendazole plus either ivermectin or diethylcarbamazine). For IDA-based MDA, the stop-MDA decision is made based on microfilariae (mf) prevalence in adults. In this study, we assess how the probability of eventually reaching elimination of transmission depends on the critical threshold used in transmission assessment surveys (TAS-es) to define whether transmission was successfully suppressed and triple-drug MDA can be stopped. This analysis focuses on treatment-naive Indian settings. We do this for a range of epidemiological and programmatic contexts, using the established LYMFASIM model for transmission and control of LF. Based on our simulations, a single TAS, one year after the last MDA round, provides limited predictive value of having achieved suppressed transmission, while a higher MDA coverage increases elimination probability, thus leading to a higher predictive value. Every additional TAS, conditional on previous TAS-es being passed with the same threshold, further improves the predictive value for low values of stop-MDA thresholds. An mf prevalence threshold of 0.5% corresponding to TAS-3 results in ≥95% predictive value even when the MDA coverage is relatively low.

Reducing the Antigen Prevalence Target Threshold for Stopping and Restarting Mass Drug Administration for Lymphatic Filariasis Elimination: A Model-Based Cost-effectiveness Simulation in Tanzania, India and Haiti.

BACKGROUND: The Global Programme to Eliminate Lymphatic Filariasis (GPELF) aims to reduce and maintain infection levels through mass drug administration (MDA), but there is evidence of ongoing transmission after MDA in areas where Culex mosquitoes are the main transmission vector, suggesting that a more stringent criterion is required for MDA decision making in these settings. METHODS: We use a transmission model to investigate how a lower prevalence threshold (<1% antigenemia [Ag] prevalence compared with <2% Ag prevalence) for MDA decision making would affect the probability of local elimination, health outcomes, the number of MDA rounds, including restarts, and program costs associated with MDA and surveys across different scenarios. To determine the cost-effectiveness of switching to a lower threshold, we simulated 65% and 80% MDA coverage of the total population for different willingness to pay per disability-adjusted life-year averted for India ($446.07), Tanzania ($389.83), and Haiti ($219.84). RESULTS: Our results suggest that with a lower Ag threshold, there is a small proportion of simulations where extra rounds are required to reach the target, but this also reduces the need to restart MDA later in the program. For 80% coverage, the lower threshold is cost-effective across all baseline prevalences for India, Tanzania, and Haiti. For 65% MDA coverage, the lower threshold is not cost-effective due to additional MDA rounds, although it increases the probability of local elimination. Valuing the benefits of elimination to align with the GPELF goals, we find that a willingness to pay per capita government expenditure of approximately $1000-$4000 for 1% increase in the probability of local elimination would be required to make a lower threshold cost-effective. CONCLUSIONS: Lower Ag thresholds for stopping MDAs generally mean a higher probability of local elimination, reducing long-term costs and health impacts. However, they may also lead to an increased number of MDA rounds required to reach the lower threshold and, therefore, increased short-term costs. Collectively, our analyses highlight that lower target Ag thresholds have the potential to assist programs in achieving lymphatic filariasis goals.

How Does the Proportion of Never Treatment Influence the Success of Mass Drug Administration Programs for the Elimination of Lymphatic Filariasis?

BACKGROUND: Mass drug administration (MDA) is the cornerstone for the elimination of lymphatic filariasis (LF). The proportion of the population that is never treated (NT) is a crucial determinant of whether this goal is achieved within reasonable time frames. METHODS: Using 2 individual-based stochastic LF transmission models, we assess the maximum permissible level of NT for which the 1% microfilaremia (mf) prevalence threshold can be achieved (with 90% probability) within 10 years under different scenarios of annual MDA coverage, drug combination and transmission setting. RESULTS: For Anopheles-transmission settings, we find that treating 80% of the eligible population annually with ivermectin + albendazole (IA) can achieve the 1% mf prevalence threshold within 10 years of annual treatment when baseline mf prevalence is 10%, as long as NT <10%. Higher proportions of NT are acceptable when more efficacious treatment regimens are used. For Culex-transmission settings with a low (5%) baseline mf prevalence and diethylcarbamazine + albendazole (DA) or ivermectin + diethylcarbamazine + albendazole (IDA) treatment, elimination can be reached if treatment coverage among eligibles is 80% or higher. For 10% baseline mf prevalence, the target can be achieved when the annual coverage is 80% and NT ≤15%. Higher infection prevalence or levels of NT would make achieving the target more difficult. CONCLUSIONS: The proportion of people never treated in MDA programmes for LF can strongly influence the achievement of elimination and the impact of NT is greater in high transmission areas. This study provides a starting point for further development of criteria for the evaluation of NT.

Accelerating Progress Towards the 2030 Neglected Tropical Diseases Targets: How Can Quantitative Modeling Support Programmatic Decisions?

Over the past decade, considerable progress has been made in the control, elimination, and eradication of neglected tropical diseases (NTDs). Despite these advances, most NTD programs have recently experienced important setbacks; for example, NTD interventions were some of the most frequently and severely impacted by service disruptions due to the coronavirus disease 2019 (COVID-19) pandemic. Mathematical modeling can help inform selection of interventions to meet the targets set out in the NTD road map 2021-2030, and such studies should prioritize questions that are relevant for decision-makers, especially those designing, implementing, and evaluating national and subnational programs. In September 2022, the World Health Organization hosted a stakeholder meeting to identify such priority modeling questions across a range of NTDs and to consider how modeling could inform local decision making. Here, we summarize the outputs of the meeting, highlight common themes in the questions being asked, and discuss how quantitative modeling can support programmatic decisions that may accelerate progress towards the 2030 targets.

A complete pipeline enables haplotyping and phasing macrohaplotype in long sequencing reads for polyploidy samples and a multi-source DNA mixture.

Macrohaplotype combines multiple types of phased DNA variants, increasing forensic discrimination power. High-quality long-sequencing reads, for example, PacBio HiFi reads, provide data to detect macrohaplotypes in multiploidy and DNA mixtures. However, the bioinformatics tools for detecting macrohaplotypes are lacking. In this study, we developed a bioinformatics software, MacroHapCaller, in which targeted loci (i.e., short TRs [STRs], single nucleotide polymorphisms, and insertion and deletions) are genotyped and combined with novel algorithms to call macrohaplotypes from long reads. MacroHapCaller uses physical phasing (i.e., read-backed phasing) to identify macrohaplotypes, and thus it can detect multi-allelic macrohaplotypes for a given sample. MacroHapCaller was validated with data generated from our designed targeted PacBio HiFi sequencing pipeline, which sequenced ∼8-kb amplicon regions harboring 20 core forensic STR loci in human benchmark samples HG002 and HG003. MacroHapCaller also was validated in whole-genome long-read sequencing data. Robust and accurate genotyping and phased macrohaplotypes were obtained with MacroHapCaller compared with the known ground truth. MacroHapCaller achieved a higher or consistent genotyping accuracy and faster speed than existing tools HipSTR and DeepVar. MacroHapCaller enables efficient macrohaplotype analysis from high-throughput sequencing data and supports applications using discriminating macrohaplotypes.

Neoadjuvant chemotherapy is associated with improved disease-free survival in pancreatic cancer patients undergoing pancreaticoduodenectomy with vascular resection.

BACKGROUND AND OBJECTIVES: Neoadjuvant chemotherapy (NAC) is becoming favored for all pancreatic adenocarcinoma (PDAC). Patients with seemingly resectable disease infrequently still display vascular involvement intraoperatively. Outcomes following NAC versus upfront surgery in patients undergoing pancreaticoduodenectomy (PD) with vascular resection are unknown. METHODS: We performed a retrospective cohort study of PDAC patients who underwent PD with vascular resection between January 1, 2013, to December 31, 2020, within a single academic center. Clinicopathologic characteristics and disease-free survival (DFS) were compared between NAC versus upfront surgery cohorts using the Kaplan-Meier estimate and Cox proportional-hazards regression model. RESULTS: Eighty-one patients who underwent PD with vascular resection for PDAC were included. Forty-six patients (56%) received NAC. The NAC cohort more often had pathologic N0 status (47.8% vs. 8.6%, p < 0.001), had decreased vascular invasion (11% vs. 40%, p = 0.002), and completed chemotherapy (80% vs. 40%, p < 0.01). The NAC cohort demonstrated improved DFS (40.5 vs. 14.3 months, p = 0.007). In multivariable analysis, NAC remained independently associated with increased DFS (HR = 0.48, p = 0.02). CONCLUSIONS: NAC was associated with improved clinicopathologic outcomes and DFS in PD with vascular resection. These findings demonstrate the advantage of NAC in PDAC patients undergoing PD with vascular resection.

A pilot study to evaluate clinical factors associated with iron and ferritin elevations during pediatric extracorporeal membrane oxygenation.

Introduction: Elevations in serum ferritin and serum iron occur during pediatric extracorporeal membrane oxygenation (ECMO). Previous reports attribute the elevation to frequent red blood cell transfusions and/or hemolysis. Chronic transfusion can cause iron deposition in tissues leading to multisystem organ dysfunction. This study aims identify clinical factors associated with elevated ferritin and iron in pediatric ECMO patients, along with post-decannulation magnetic resonance imaging (MRI) assessment of iron deposition in liver and brain.Methods: Prospective, pilot study, using descriptive statistics to investigate potential associations between patient characteristics, serum ferritin and iron levels, and post-decannulation hepatic and basal ganglia iron deposition.Results: In this study, nine patients (100%) had elevated serum ferritin levels during ECMO. High ferritin levels were more common with veno-arterial than with veno-venous cannulation (p = 0.026) and were also associated with high plasma free hemoglobin levels (p < 0.001). Five patients presented with elevated serum iron levels. High serum iron levels were associated with higher daily (p = 0.016) and cumulative transfusion volumes (p = 0.013) as well ECMO duration beyond 7 days. MRI scans were performed on three patients with no evidence of abnormal iron deposition detected in the liver or brain.Conclusions: This pilot study shows that during pediatric ECMO, elevations in serum ferritin and serum iron occur and those elevations may be related to the cannulation modality, ECMO duration, amount of hemolysis, and volume of red blood cell transfusions. Further investigation is warranted to fully understand the implications of elevated serum iron and ferritin in pediatric ECMO.

Effect of meloxicam or robenacoxib administration timing on renal function and postoperative analgesia in cats undergoing ovariohysterectomy: A randomized, blinded, controlled clinical trial.

We evaluated the effect of administration timing of meloxicam and robenacoxib on renal function, platelet cyclo-oxygenase and perioperative analgesia in 60 cats undergoing ovariohysterectomy, in a prospective randomized blinded controlled study. Twelve cats were randomly allocated to one subcutaneous treatment group: meloxicam (0.2 mg/kg) or robenacoxib (2 mg/kg) at admission (MA, RA), at induction (MI, RI) and robenacoxib at the end of surgery (RE). All cats received the same anaesthesia protocol. Plasma renin activity (PRA), plasma creatinine, drug concentrations and serum thromboxane (TxB2) were measured sequentially. Anaesthesia significantly increased PRA, as activity at end of the surgery was higher than 2 h later (mean ± SD: 26.6 ± 2.8 versus 10.0 ± 3.9 ng/mL/h). PRA remained higher at 2 h post-surgery in admission groups compared to induction groups (p = .01). Serum TxB2 was lower with meloxicam than robenacoxib (p = .001), and was lower in the MA than each robenacoxib group at catheter placement. Admission groups (16/24 from RA and MA groups) received earlier rescue analgesia than other groups (p = .033). In conclusion, the renin-angiotensin system was activated during anaesthesia despite cyclo-oxygenase inhibition, possibly due to hypotension or surgical stimulation. There was no effect of drug or timing on the markers of renal function but one cat receiving meloxicam at induction had suspected IRIS grade II acute kidney injury.

Interpreting forces as deep learning gradients improves quality of predicted protein structures.

Protein structure predictions from deep learning models like AlphaFold2, despite their remarkable accuracy, are likely insufficient for direct use in downstream tasks like molecular docking. The functionality of such models could be improved with a combination of increased accuracy and physical intuition. We propose a new method to train deep learning protein structure prediction models using molecular dynamics force fields to work toward these goals. Our custom PyTorch loss function, OpenMM-Loss, represents the potential energy of a predicted structure. OpenMM-Loss can be applied to any all-atom representation of a protein structure capable of mapping into our software package, SidechainNet. We demonstrate our method's efficacy by finetuning OpenFold. We show that subsequently predicted protein structures, both before and after a relaxation procedure, exhibit comparable accuracy while displaying lower potential energy and improved structural quality as assessed by MolProbity metrics.

Using T4 genetics and Laemmli's development of high-resolution SDS gel electrophoresis to reveal structural protein interactions controlling protein folding and phage self-assembly.

One of the most transformative experimental techniques in the rise of modern molecular biology and biochemistry was the development of high-resolution sodium dodecyl sulfate polyacrylamide gel electrophoresis, which allowed separation of proteins-including structural proteins-in complex mixtures according to their molecular weights. Its development was intimately tied to investigations of the control of virus assembly within phage-infected cells. The method was developed by Ulrich K. Laemmli working in the virus structural group led by Aaron Klug at the famed Medical Research Council Laboratory for Molecular Biology at Cambridge, UK. While Laemmli was tackling T4 head assembly, I sat at the next bench working on T4 tail assembly. To date, Laemmli's original paper has been cited almost 300,000 times. His gel procedure and our cooperation allowed us to sort out the sequential protein-protein interactions controlling the viral self-assembly pathways. It is still not fully appreciated that this control involved protein conformational change induced by interaction with an edge of the growing structure. Subsequent efforts of my students and I to understand how temperature-sensitive mutations interfered with assembly were important in revealing the intracellular off-pathway aggregation processes competing with productive protein folding. These misfolding processes slowed the initial productivity of the biotechnology industry. The article below describes the scientific origin, context, and sociology that supported these advances in protein biochemistry, protein expression, and virus assembly. The cooperation and collaboration that was integral to both the Laboratory for Molecular Biology culture and phage genetics fields were key to these endeavors.

Cingulin binds to the ZU5 domain of scaffolding protein ZO-1 to promote its extended conformation, stabilization, and tight junction accumulation.

Zonula occludens-1 (ZO-1), the major scaffolding protein of tight junctions (TJs), recruits the cytoskeleton-associated proteins cingulin (CGN) and paracingulin (CGNL1) to TJs by binding to their N-terminal ZO-1 interaction motif. The conformation of ZO-1 can be either folded or extended, depending on cytoskeletal tension and intramolecular and intermolecular interactions, and only ZO-1 in the extended conformation recruits the transcription factor DbpA to TJs. However, the sequences of ZO-1 that interact with CGN and CGNL1 and the role of TJ proteins in ZO-1 TJ assembly are not known. Here, we used glutathione-S-transferase pulldowns and immunofluorescence microscopy to show that CGN and CGNL1 bind to the C-terminal ZU5 domain of ZO-1 and that this domain is required for CGN and CGNL1 recruitment to TJs and to phase-separated ZO-1 condensates in cells. We show that KO of CGN, but not CGNL1, results in decreased accumulation of ZO-1 at TJs. Furthermore, ZO-1 lacking the ZU5 domain showed decreased accumulation at TJs, was detectable along lateral contacts, had a higher mobile fraction than full-length ZO-1 by fluorescence recovery after photobleaching analysis, and had a folded conformation, as determined by structured illumination microscopy of its N-terminal and C-terminal ends. The CGN-ZU5 interaction promotes the extended conformation of ZO-1, since binding of the CGN-ZO-1 interaction motif region to ZO-1 resulted in its interaction with DbpA in cells and in vitro. Together, these results show that binding of CGN to the ZU5 domain of ZO-1 promotes ZO-1 stabilization and accumulation at TJs by promoting its extended conformation.

Copper Reductase Activity and Free Radical Chemistry by Cataract-Associated Human Lens γ-Crystallins.

Cataracts are caused by high-molecular-weight aggregates of human eye lens proteins that scatter light, causing lens opacity. Metal ions have emerged as important potential players in the etiology of cataract disease, as human lens γ-crystallins are susceptible to metal-induced aggregation. Here, the interaction of Cu2+ ions with γD-, γC-, and γS-crystallins, the three most abundant γ-crystallins in the lens, has been evaluated. Cu2+ ions induced non-amyloid aggregation in all three proteins. Solution turbidimetry, sodium dodecyl sulfate poly(acrylamide) gel electrophoresis (SDS-PAGE), circular dichroism, and differential scanning calorimetry showed that the mechanism for Cu-induced aggregation involves: (i) loss of ß-sheet structure in the N-terminal domain; (ii) decreased thermal and kinetic stability; (iii) formation of metal-bridged species; and (iv) formation of disulfide-bridged dimers. Isothermal titration calorimetry (ITC) revealed distinct Cu2+ binding affinities in the γ-crystallins. Electron paramagnetic resonance (EPR) revealed two distinct Cu2+ binding sites in each protein. Spin quantitation demonstrated the reduction of γ-crystallin-bound Cu2+ ions to Cu+ under aerobic conditions, while X-ray absorption spectroscopy (XAS) confirmed the presence of linear or trigonal Cu+ binding sites in γ-crystallins. Our EPR and XAS studies revealed that γ-crystallins' Cu2+ reductase activity yields a protein-based free radical that is likely a Tyr-based species in human γD-crystallin. This unique free radical chemistry carried out by distinct redox-active Cu sites in human lens γ-crystallins likely contributes to the mechanism of copper-induced aggregation. In the context of an aging human lens, γ-crystallins could act not only as structural proteins but also as key players for metal and redox homeostasis.

Robotic versus open pancreaticoduodenectomy in elderly patients: a propensity score-matched analysis.

BACKGROUND: Pancreaticoduodenectomy (PD) is complex procedure with high morbidity in the elderly. This retrospective study aimed to compare post-operative outcomes in patients ≥75 years of age who underwent robot-assisted (RA)PD and open PD. METHODS: We analyzed 2502 patients ≥75 years of age who underwent PD from 2015 to 2018 in the National Surgical Quality Improvement Program (NSQIP) database. RAPD and open PD patients were propensity score matched 1:5 to assess the 30-day outcomes of interest: postoperative complications, length of stay, discharge destination, and readmissions. RESULTS: Of 725 matched patients, 110 underwent RAPD, 615 OPD, and 12 were converted to an open operation. Post-operative outcomes were largely similar between cohorts. RAPD was associated a shorter length of stay (median 8 days, interquartile range [IQR] 6 to 11) than OPD (median 8 days, IQR 7 to 13) (p = 0.003). However, RAPD was associated with more readmissions (28.1% vs. 17.7%; p = 0.02). CONCLUSIONS: RAPD in patients ≥75 years of age appears to be safe and has a similar complication profile to open PD. Randomized or well-designed prospective matched studies are needed to confirm these findings.

Kilohertz frame-rate two-photon tomography.

Point-scanning two-photon microscopy enables high-resolution imaging within scattering specimens such as the mammalian brain, but sequential acquisition of voxels fundamentally limits its speed. We developed a two-photon imaging technique that scans lines of excitation across a focal plane at multiple angles and computationally recovers high-resolution images, attaining voxel rates of over 1 billion Hz in structured samples. Using a static image as a prior for recording neural activity, we imaged visually evoked and spontaneous glutamate release across hundreds of dendritic spines in mice at depths over 250 µm and frame rates over 1 kHz. Dendritic glutamate transients in anesthetized mice are synchronized within spatially contiguous domains spanning tens of micrometers at frequencies ranging from 1-100 Hz. We demonstrate millisecond-resolved recordings of acetylcholine and voltage indicators, three-dimensional single-particle tracking and imaging in densely labeled cortex. Our method surpasses limits on the speed of raster-scanned imaging imposed by fluorescence lifetime.

Author Correction: Kilohertz frame-rate two-photon tomography.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Correcting thermal-emission-induced detector saturation in infrared spectroscopy.

We found that temperature-dependent infrared spectroscopy measurements (i.e., reflectance or transmittance) using a Fourier-transform spectrometer can have substantial errors, especially for elevated sample temperatures and collection using an objective lens. These errors can arise as a result of partial detector saturation due to thermal emission from the measured sample reaching the detector, resulting in nonphysical apparent reduction of reflectance or transmittance with increasing sample temperature. Here, we demonstrate that these temperature-dependent errors can be corrected by implementing several levels of optical attenuation that enable convergence testing of the measured reflectance or transmittance as the thermal-emission signal is reduced, or by applying correction factors that can be inferred by looking at the spectral regions where the sample is not expected to have a substantial temperature dependence.

Metabotyping the Welsh population of badgers based on thoracic fluid.

INTRODUCTION: The European badger (Meles meles) is a known wildlife reservoir for bovine tuberculosis (bTB) and a better understanding of the epidemiology of bTB in this wildlife species is required for disease control in both wild and farmed animals. Flow infusion electrospray-high-resolution mass spectrometry (FIE-HRMS) may potentially identify novel metabolite biomarkers based on which new, rapid, and sensitive point of care tests for bTB infection could be developed. OBJECTIVES: In this foundational study, we engaged on assessing the baseline metabolomic variation in the non-bTB infected badger population ("metabotyping") across Wales. METHODS: FIE-HRMS was applied on thoracic fluid samples obtained by post-mortem of bTB negative badgers (n = 285) which were part of the Welsh Government 'All Wales Badger Found Dead' study. RESULTS: Using principal component analysis and partial least squares-discriminant analyses, the major sources of variation were linked to sex, and to a much lesser extent age, as indicated by tooth wear. Within the female population, variation was seen between lactating and non-lactating individuals. No significant variation linked to the presence of bite wounds, obvious lymphatic lesions or geographical region of origin was observed. CONCLUSION: Future metabolomic work when making comparisons between bTB infected and non-infected badger samples will only need be sex-matched and could focus on males only, to avoid lactation bias.

Untargeted metabolomic analysis of thoracic blood from badgers indicate changes linked to infection with bovine tuberculosis (Mycobacterium bovis): a pilot study.

INTRODUCTION: Mycobacterium bovis, the causative agent of bovine tuberculosis (bTB) in cattle, represents a major disease burden to UK cattle farming, with considerable costs associated with its control. The European badger (Meles meles) is a known wildlife reservoir for bTB and better knowledge of the epidemiology of bTB through testing wildlife is required for disease control. Current tests available for the diagnosis of bTB in badgers are limited by cost, processing time or sensitivities. MATERIALS AND METHODS: We assessed the ability of flow infusion electrospray-high-resolution mass spectrometry (FIE-HRMS) to determine potential differences between infected and non-infected badgers based on thoracic blood samples obtained from badgers found dead in Wales. Thoracic blood samples were autoclaved for handling in a containment level 2 (CL2) hazard laboratory. RESULTS: Here we show the major differences associated with with M. bovis infection were changes to folate, pyrimidine, histidine, glycerophospholipid and phosphonate metabolism. CONCLUSIONS: Our studies have indicated differences in the metabolomic signature of badgers found dead in relation to their infection status, suggesting metabolomics could hold potential for developing novel diagnostics for bTB in badgers. As well as highlighting a potential way to handle samples containing a highly pathogenic agent at CL2 for metabolomics studies.