On the undecidability among kinetic models: from model selection to model averaging

Social Stigma and COVID-19 Vaccine Refusal in France

CONTEXT

In 2021, French health authorities strongly promoted vaccination against COVID-19. The authors assumed that refusing this vaccine became a stigma, and they investigated potential public stigma toward unvaccinated people among the French population.

METHODS

A representative sample of the French adult population (N = 2,015) completed an online questionnaire in September 2021. The authors focused on participants who were already vaccinated against COVID-19 or intended to get vaccinated (N = 1,742). A cluster analysis was used to obtain contrasted attitudinal profiles, and the authors investigated associated factors with logistic regressions.

FINDINGS

Regarding attitudes toward unvaccinated people, a majority of respondents supported several pejorative statements, and a significant minority also endorsed social rejection attitudes. The authors found four contrasting attitudinal profiles: moral condemnation only (32% of respondents), full stigma (26%), no stigma (26%), and stigma rejection (16%). Early vaccination, civic motives for it, faith in science, rejection of political extremes, and being aged 65 or older were the main factors associated with stigmatizing attitudes toward unvaccinated people.

CONCLUSIONS

The authors found some evidence of stigmatization toward unvaccinated people, but further research is needed, especially to investigate perceived stigmatization among them. The authors discuss their results with reference to the concept of "folk devils" and from a public health perspective.

Determining the optimal pharmacokinetic modelling and simplified quantification method of [18F]AlF-P16-093 for patients with primary prostate cancer (PPCa)

PURPOSE

This paper discusses the optimization of pharmacokinetic modelling and alternate simplified quantification method for [18F]AlF-P16-093, a novel tracer for in vivo imaging of prostate cancer.

METHODS

Dynamic PET/CT scans were conducted on eight primary prostate cancer patients, followed by a whole-body scan at 60 min post-injection. Time-activity curves (TACs) were obtained by drawing volumes of interest for primary prostatic and metastatic lesions. Optimal kinetic modelling involved evaluating three compartmental models (1T2K, 2T3K, and 2T4K) accounting for fractional blood volume (Vb). The simplified quantification method was then determined based on the correlation between the static uptake measure and total distribution volume (Vt) obtained from the optimal pharmacokinetic analysis.

RESULTS

In total, 17 intraprostatic lesions, 10 lymph nodes, and 36 osseous metastases were evaluated. Visually, the contrast of the tumor increased and showed the steepest incline within the first few minutes, whereas background activity decreased over time. Full pharmacokinetic analysis revealed that a reversible two-compartmental (2T4K) model is the preferred kinetic model for the given tracer. The kinetic parameters K1, k3, Vb, and Vt were all significantly higher in lesions when compared with normal tissue (P < 0.01). Several simplified protocols were tested for approximating comprehensive dynamic quantification in tumors, with image-based SURmean (the ratio of tumor SUVmean to blood SUVmean) within the 28-34 min window found to be sufficient for approximating the total distribution Vt values (R2 = 0.949, P < 0.01). Both Vt and SURmean correlated significantly with the total serum prostate-specific antigen (tPSA) levels (P < 0.01).

CONCLUSIONS

This study introduced an optimized pharmacokinetic modelling approach and a simplified acquisition method for [18F]AlF-P16-093, a novel PSMA-targeted radioligand, highlighting the feasibility of utilizing one static PET imaging (between 30 and 60 min) for the diagnosis of prostate cancer. Note that the image-derived input function in this study may not reflect the true corrected plasma input function, therefore the interpretation of the associated kinetic parameter estimates should be done with caution.

Nor climate nor human impact factors: Chytrid infection shapes the skin bacterial communities of an endemic amphibian in a biodiversity hotspot

The bacterial communities of the amphibian skin (i.e., the bacteriome) are critical to the host's innate immune system. However, it is unclear how different drivers can alter this function by modulating the bacteriome's structure. Our aim was to assess the extent to which different host attributes and extrinsic factors influence the structure of the bacterial communities of the skin. Skin bacterial diversity was examined in 148 individuals of the four-eyed frog (Pleurodema thaul) from 16 localities spanning almost 1800 km in latitude. The richness and beta diversity of bacterial families and the richness and abundance of Bd-inhibitory bacterial genera were used to describe their structure. Predictors associated with the host (developmental stage, genetic lineage, individual Batrachochytrium dendrobatidis [Bd] infection status) and the landscape (current climate, degree of anthropogenic disturbance) were used in the statistical modeling in an information theoretical approach. Bd infection and host developmental stage were the only predictors affecting bacteriome richness, with Bd+ individuals and postmetamorphic stages (adults and juveniles) having higher richness than Bd- ones and tadpoles. High diversity in Bd+ individuals is not driven by bacterial genera with known anti-Bd properties. Beta diversity was not affected by Bd infection and was mainly a consequence of bacterial family turnover rather than nestedness. Finally, for those bacterial genera known to have inhibitory effects on chytrid, Bd+ individuals had a slightly higher diversity than Bd- ones. Our study confirms an association between Bd infection and the host developmental stage with the bacterial communities of the skin of P. thaul. Unexpectedly, macroclimate and human impact factors do not seem to play a role in shaping the amphibian skin microbiome. Our study exemplifies that focusing on a single host-parasite system over a large geographic scale can provide essential insights into the factors driving host-parasite-bacteriome interactions.

Quantification of P-glycoprotein function at the human blood-brain barrier using [18F]MC225 and PET

INTRODUCTION

P-glycoprotein (P-gp) is one of the most studied efflux transporters at the blood-brain barrier. It plays an important role in brain homeostasis by protecting the brain from a variety of endogenous and exogeneous substances. Changes in P-gp function are associated both with the onset of neuropsychiatric diseases, including Alzheimer's disease and Parkinson's disease, and with drug-resistance, for example in treatment-resistant depression. The most widely used approach to measure P-gp function in vivo is (R)-[11C]verapamil PET. (R)-[11C]verapamil is, however, an avid P-gp substrate, which complicates the use of this tracer to measure an increase in P-gp function as its baseline uptake is already very low. [18F]MC225 was developed to measure both increases and decreases in P-gp function.

AIM

The aim of this study was (1) to identify the pharmacokinetic model that best describes [18F]MC225 kinetics in the human brain and (2) to determine test-retest variability.

METHODS

Five (2 male, 3 female) of fourteen healthy subjects (8 male, 6 female, age 67 ± 5 years) were scanned twice (injected dose 201 ± 47 MBq) with a minimum interval of 2 weeks between scans. Each scanning session consisted of a 60-min dynamic [18F]MC225 scan with continuous arterial sampling. Whole brain grey matter data were fitted to a single tissue compartment model, and to reversible and irreversible two tissue-compartment models to obtain various outcome parameters (in particular the volume of distribution (VT), Ki, and the rate constants K1 and k2). In addition, a reversible two-tissue compartment model with fixed k3/k4 was included. The preferred model was selected based on the weighted Akaike Information Criterion (AIC) score. Test-retest variability (TRTV) was determined to assess reproducibility.

RESULTS

Sixty minutes post-injection, the parent fraction was 63.8 ± 4.0%. The reversible two tissue compartment model corrected for plasma metabolites with an estimated blood volume (VB) showed the highest AIC weight score of 34.3 ± 17.6%. The TRVT of the VT for [18F]MC225 PET scans was 28.3 ± 20.4% for the whole brain grey matter region using this preferred model.

CONCLUSION

[18F]MC225 VT, derived using a reversible two-tissue compartment model, is the preferred parameter to describe P-gp function in the human BBB. This outcome parameter has an average test-retest variability of 28%.

TRIAL REGISTRATION

EudraCT 2020-001564-28 . Registered 25 May 2020.

Vaccine hesitancy about the HPV vaccine among French young women and their parents: a telephone survey

BACKGROUND

The human papillomavirus (HPV) vaccine reduces the burden of cervical and other cancers. In numerous countries, a slow uptakeof this vaccine persists, calling for a better understanding of the structural factors leading to vaccine acceptation. We aimed to assess the attitudes toward HPV vaccination among its intended public to explore its specific characteristics.

METHODS

A random cross-sectional telephone survey of the French general population provided data from a sample of 2426 respondents of the target public: the parents of young women and the young women aged 15-25 themselves. We applied cluster analysis to identify contrasting attitudinal profiles, and logistic regressions with a model averaging method to investigate and rank the factors associated with these profiles.

RESULTS

A third of the respondents had never heard of HPV. However, most of the respondents who had heard of it agreed that it is a severe (93.8%) and frequent (65.1%) infection. Overall, 72.3% of them considered the HPV vaccine to be effective, but 54% had concerns about its side effects. We identified four contrasting profiles based on their perceptions of this vaccine: informed supporters, objectors, uninformed supporters, and those who were uncertain. In multivariate analysis, these attitudinal clusters were the strongest predictors of HPV vaccine uptake, followed by attitudes toward vaccination in general.

CONCLUSIONS

Tailored information campaigns and programs should address the specific and contrasted concerns about HPV vaccination of both young women and of their parents.

Heat tolerance of marine ectotherms in a warming Antarctica

Global warming is affecting the Antarctic continent in complex ways. Because Antarctic organisms are specialized to living in the cold, they are vulnerable to increasing temperatures, although quantitative analyses of this issue are currently lacking. Here we compiled a total of 184 estimates of heat tolerance belonging to 39 marine species and quantified how survival is affected concomitantly by the intensity and duration of thermal stress. Species exhibit thermal limits displaced toward colder temperatures, with contrasting strategies between arthropods and fish that exhibit low tolerance to acute heat challenges, and brachiopods, echinoderms, and molluscs that tend to be more sensitive to chronic exposure. These differences might be associated with mobility. A dynamic mortality model suggests that Antarctic organisms already encounter temperatures that might be physiologically stressful and indicate that these ecological communities are indeed vulnerable to ongoing rising temperatures.

Vaccine hesitancy in health-care providers in Western countries: a narrative review

INTRODUCTION

Vaccine hesitancy (VH) is a leading cause of suboptimal vaccine uptake rates worldwide. The interaction between patients and health-care providers (HCPs) is the keystone in addressing VH. However, significant proportions of HCPs, including those who administer vaccines, are personally and professionally vaccine-hesitant.

AREAS COVERED

This narrative review sought to characterize the nature, extent, correlates, and consequences of VH among HCPs. We included 39 quantitative and qualitative studies conducted in Western countries, published since 2015, that assessed VH among HCPs in general, for several vaccines. Studies were reviewed using the WHO 3Cs model - (lack of) confidence, complacency, and (lack of) convenience.

EXPERT OPINION

Despite the lack of validated tools and substantial heterogeneity in the methods used to measure VH among HCPs, this review confirms its presence in this population, at frequencies that vary by country, profession type, setting, and level of medical education. Lack of knowledge and mistrust in health authorities/pharmaceutical industry/experts were among its principal drivers. Improving the content about vaccination in HCPs' training programs, facilitating access to reliable information for use during consultations, and developing and validating instruments to measure HCPs' VH and its determinants are key to addressing VH among HCPs.

Acceptance of childhood and adolescent vaccination against COVID-19 in France: a national cross-sectional study in May 2021

The French health authorities extended vaccination against COVID-19 to adolescents in June 2021, during the epidemic resurgence linked to the delta variant and because of insufficient vaccination coverage to ensure collective protection. In May 2021, we conducted a national online cross-sectional survey of 2533 adults in France to study their attitudes toward COVID-19 vaccines and their acceptance of child/adolescent vaccination according to targeted age groups (<6 years; 6-11; 12-17) and its determinants. We applied a multi-model averaged logistic regression for each of these age groups to study the determinants of favorability to vaccination. Among the respondents, 62.7% (1597) accepted COVID-19 vaccination for adolescents, 48.3% (1223) for children aged 6-11 years, and only 31% (783) for children under 6 years. Acceptance increased with fear of contracting COVID-19 and trust in institutions and decreased as the COVID-19 vaccine risk perception score increased. People favorable to vaccination in general and those sensitive to social pressure were also more often favorable to vaccinating children/adolescents than those who were not. Drivers of acceptance were ranked differently for the different age groups. Understanding these differences is essential to anticipating obstacles to vaccination of these age groups and designing appropriate information and motivational strategies to support it.

Using Isotopic Data to Evaluate Esox lucius (Linnaeus, 1758) Natal Origins in a Hydrologically Complex River Basin

Otolith microchemistry has emerged as a powerful technique with which to identify the natal origins of fishes, but it relies on differences in underlying geology that may occur over large spatial scales. An examination of how small a spatial scale on which this technique can be implemented, especially in water bodies that share a large proportion of their flow, would be useful for guiding aquatic invasive species control efforts. We examined trace isotopic signatures in northern pike (Esox lucius) otoliths to estimate their provenance between two reservoirs in the Upper Yampa River Basin, Colorado, USA. This is a challenging study area as these reservoirs are only 11-rkm apart on the same river and thus share a high proportion of their inflow. We found that three isotopes (86Sr, 137Ba, and 55Mn) were useful in discriminating between these reservoirs, but their signatures varied annually, and the values overlapped. Strontium isotope ratios (87Sr/86Sr) were different between sites and relatively stable across three years, which made them an ideal marker for determining northern pike provenance. Our study demonstrates the usefulness of otolith microchemistry for natal origin determination within the same river over a relatively small spatial scale when there are geologic differences between sites, especially geologic differences underlying tributaries between sites.

In vivo noninvasive preclinical tumor hypoxia imaging methods: a review

Tumors exhibit areas of decreased oxygenation due to malformed blood vessels. This low oxygen concentration decreases the effectiveness of radiation therapy, and the resulting poor perfusion can prevent drugs from reaching areas of the tumor. Tumor hypoxia is associated with poorer prognosis and disease progression, and is therefore of interest to preclinical researchers. Although there are multiple different ways to measure tumor hypoxia and related factors, there is no standard for quantifying spatial and temporal tumor hypoxia distributions in preclinical research or in the clinic. This review compares imaging methods utilized for the purpose of assessing spatio-temporal patterns of hypoxia in the preclinical setting. Imaging methods provide varying levels of spatial and temporal resolution regarding different aspects of hypoxia, and with varying advantages and disadvantages. The choice of modality requires consideration of the specific experimental model, the nature of the required characterization and the availability of complementary modalities as well as immunohistochemistry.

(+)-[18F]Flubatine as a novel α4β2 nicotinic acetylcholine receptor PET ligand-results of the first-in-human brain imaging application in patients with β-amyloid PET-confirmed Alzheimer's disease and healthy controls

PURPOSES

We present the first in-human brain PET imaging data of the new α4β2 nicotinic acetylcholine receptor (nAChR)-targeting radioligand (+)-[18F]Flubatine. Aims were to develop a kinetic modeling-based approach to quantify (+)-[18F]Flubatine and compare the data of healthy controls (HCs) and patients with Alzheimer's disease (AD); to investigate the partial volume effect (PVE) on regional (+)-[18F]Flubatine binding; and whether (+)-[18F]Flubatine binding and cognitive test data respective β-amyloid radiotracer accumulation were correlated.

METHODS

We examined 11 HCs and 9 mild AD patients. All subjects underwent neuropsychological testing and [11C]PiB PET/MRI examination. (+)-[18F]Flubatine PET data were evaluated using full kinetic modeling and regional as well as voxel-based analyses.

RESULTS

With 270-min p.i., the unchanged parent compound amounted to 97 ± 2%. Adequate fits of the time-activity curves were obtained with the 1 tissue compartment model (1TCM). (+)-[18F]Flubatine distribution volume (binding) was significantly reduced in bilateral mesial temporal cortex in AD patients compared with HCs (right 10.6 ± 1.1 vs 11.6 ± 1.4, p = 0.049; left 11.0 ± 1.1 vs 12.2 ± 1.8, p = 0.046; one-sided t tests each). PVE correction increased not only (+)-[18F]Flubatine binding of approximately 15% but also standard deviation of 0.4-70%. Cognitive test data and (+)-[18F]Flubatine binding were significantly correlated in the left anterior cingulate, right posterior cingulate, and right parietal cortex (r > 0.5, p < 0.05 each). In AD patients, (+)-[18F]Flubatine binding and [11C]PiB standardized uptake value ratios were negatively correlated in several regions; whereas in HCs, a positive correlation between cortical (+)-[18F]Flubatine binding and [11C]PiB accumulation in the white matter was found. No adverse event related to (+)-[18F]Flubatine occurred.

CONCLUSION

(+)-[18F]Flubatine is a safe and stable PET ligand. Full kinetic modeling can be realized by 1TCM without metabolite correction. (+)-[18F]Flubatine binding affinity was high enough to detect group differences. Of interest, correlation between white matter β-amyloid PET uptake and (+)-[18F]Flubatine binding indicated an association between white matter integrity and availability of α4β2 nAChRs. Overall, (+)-[18F]Flubatine showed favorable characteristics and has therefore the potential to serve as α4β2 nAChR-targeting PET ligand in further clinical trials.

Attitudes of healthcare workers towards COVID-19 vaccination: a survey in France and French-speaking parts of Belgium and Canada, 2020

In October and November 2020, we conducted a survey of 2,678 healthcare workers (HCWs) involved in general population immunisation in France, French-speaking Belgium and Quebec, Canada to assess acceptance of future COVID-19 vaccines (i.e. willingness to receive or recommend these) and its determinants. Of the HCWs, 48.6% (n = 1,302) showed high acceptance, 23.0% (n = 616) moderate acceptance and 28.4% (n = 760) hesitancy/reluctance. Hesitancy was mostly driven by vaccine safety concerns. These must be addressed before/during upcoming vaccination campaigns.

Spatial heterogeneity of radiolabeled choline positron emission tomography in tumors of patients with non-small cell lung cancer: first-in-patient evaluation of [18F]fluoromethyl-(1,2-2H4)-choline

Purpose: The spatio-molecular distribution of choline and its metabolites in tumors is highly heterogeneous. Due to regulation of choline metabolism by hypoxic transcriptional signaling and other survival factors, we envisage that detection of such heterogeneity in patient tumors could provide the basis for advanced localized therapy. However, non-invasive methods to assess this phenomenon in patients are limited. We investigated such heterogeneity in Non-Small Cell Lung Cancer (NSCLC) with [18F]fluoromethyl-(1,2-2H4) choline ([18F]D4-FCH) and positron emission tomography/computed tomography (PET/CT). Experimental design: [18F]D4-FCH (300.5±72.9MBq [147.60-363.6MBq]) was administered intravenously to 17 newly diagnosed NSCLC patients. PET/CT scans were acquired concurrently with radioactive blood sampling to permit mathematical modelling of blood-tissue transcellular rate constants. Comparisons were made with biopsy-derived choline kinase-α (CHKα) expression and diagnostic [18F]fluorodeoxyglucose ([18F]FDG) scans. Results: Oxidation of [18F]D4-FCH to [18F]D4-fluorobetaine was suppressed (48.58±0.31% parent at 60 min) likely due to the deuterium isotope effect embodied within the design of the radiotracer. Early (5 min) and late (60 min) images showed specific uptake of tracer in all 51 lesions (tumors, lymph nodes and metastases) from 17 patients analyzed. [18F]D4-FCH-derived uptake (SUV60max) in index primary lesions (n=17) ranged between 2.87-10.13; lower than that of [18F]FDG PET [6.89-22.64]. Mathematical modelling demonstrated net irreversible uptake of [18F]D4-FCH at steady-state, and parametric mapping of the entire tumor showed large intratumorally heterogeneity in radiotracer retention, which is likely to have influenced correlations with biopsy-derived CHKα expression. Conclusions: [18F]D4-FCH is detectable in NSCLC with large intratumorally heterogeneity, which could be exploited in the future for targeting localized therapy.

Quantification of 11C-PIB kinetics in cardiac amyloidosis

BACKGROUND

The purpose of this work was to determine the optimal tracer kinetic model of 11C-PIB and to validate the use of the simplified methods retention index (RI) and standardized uptake value (SUV) for quantification of cardiac 11C-PIB uptake in amyloidosis.

METHODS AND RESULTS

Single-tissue, reversible and irreversible two-tissue models were fitted to data from seven cardiac amyloidosis patients who underwent 11C-PIB PET scans and arterial blood sampling for measurement of blood radioactivity and metabolites. The irreversible two-tissue model (2Tirr) best described cardiac 11C-PIB uptake. RI and SUV showed high correlation with the rate of irreversible binding (Ki) from the 2Tirr model (r2 =0.95 and r2 =0.94). Retrospective data from 10 amyloidosis patients and 5 healthy controls were analyzed using RI, SUV, as well as compartment modelling with a population-average metabolite correction. All measures were higher in amyloidosis patients than in healthy controls (p=.001), but with an overlap between groups for Ki.

CONCLUSION

An irreversible two-tissue model best describes the 11C-PIB uptake in cardiac amyloidosis. RI and SUV correlate well with Ki from the 2Tirr model. RI and SUV discriminate better between amyloidosis patients and controls than Ki based on population-average metabolite correction.

Robust nonlinear parameter estimation in tracer kinetic analysis using infinity norm regularization and particle swarm optimization

In positron emission tomography (PET) studies, the voxel-wise calculation of individual rate constants describing the tracer kinetics is quite challenging because of the nonlinear relationship between the rate constants and PET data and the high noise level in voxel data. Based on preliminary simulations using a standard two-tissue compartment model, we can hypothesize that it is possible to reduce errors in the rate constant estimates when constraining the overestimation of the larger of two exponents in the model equation. We thus propose a novel approach based on infinity-norm regularization for limiting this exponent. Owing to the non-smooth cost function of this regularization scheme, which prevents the use of conventional Jacobian-based optimization methods, we examined a proximal gradient algorithm and the particle swarm optimization (PSO) through a simulation study. Because it exploits multiple initial values, the PSO method shows much better convergence than the proximal gradient algorithm, which is susceptible to the initial values. In the implementation of PSO, the use of a Gamma distribution to govern random movements was shown to improve the convergence rate and stability compared to a uniform distribution. Consequently, Gamma-based PSO with regularization was shown to outperform all other methods tested, including the conventional basis function method and Levenberg-Marquardt algorithm, in terms of its statistical properties.

Reserve Flux Capacity in the Pentose Phosphate Pathway by NADPH Binding Is Conserved across Kingdoms

All organisms evolved defense mechanisms to counteract oxidative stress and buildup of reactive oxygen species (ROS). To test whether a potentially conserved mechanism exists for the rapid response, we investigated immediate metabolic dynamics of Escherichia coli, yeast, and human dermal fibroblasts to oxidative stress that we found to be conserved between species. To elucidate the regulatory mechanisms that implement this metabolic response, we developed mechanistic kinetic models for each organism's central metabolism and systematically tested activation and inactivation of each irreversible reaction by each metabolite. This ensemble modeling predicts in vivo relevant metabolite-enzyme interactions based on their ability to quantitatively describe metabolite dynamics. All three species appear to inhibit their oxidative pentose phosphate pathway during normal growth by the redox cofactor NADPH and relieve this inhibition to increase the pathway flux for detoxification of ROS during stress, with the sole exception of yeast when exposed to high levels of stress.

Influx rate of 18F-fluoroaminosuberic acid reflects cystine/glutamate antiporter expression in tumour xenografts

PURPOSE

¹⁸F-fluoroaminosuberic acid (¹⁸F-FASu) is a recently developed amino acid tracer for positron emission tomography (PET) of oxidative stress that may offer improved tumour assessment over the conventional tracer ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG). Our aim was to evaluate and relate dynamic ¹⁸F-FASu and ¹⁸F-FDG uptake with pharmacokinetic modelling to transporter protein expression levels in a panel of diverse tumour xenograft lines.

METHODS

Four different tumour xenograft lines were implanted in female athymic nude mice: MAS98.12 and HBCx3 (breast), TPMX (osteosarcoma) and A549 (lung). Dynamic PET over 60 min was performed on a small animal unit. The time-activity curves (TACs) for ¹⁸F-FASu and ¹⁸F-FDG in individual tumours were used to extract early (SUV_E; 2 min p.i.) and late (SUV_L; 55 min p.i.) standardised uptake values. Pharmacokinetic two-tissue compartment models were applied to the TACs to estimate rate constants K₁-k₄ and blood volume fraction v_B. Relative levels of cystine/glutamate antiporter subunit xCT were assessed by western blotting, and expression of GLUT1 and CD31 by immunohistochemistry.

RESULTS

¹⁸F-FASu showed higher SUV_E, whilst ¹⁸F-FDG exhibited higher SUV_L. Influx rate K₁ for ¹⁸F-FASu was significantly correlated with xCT levels (p = 0.001) and was significantly higher than K₁ for ¹⁸F-FDG (p < 0.001). K₁ for ¹⁸F-FDG was significantly correlated with GLUT1 levels (p = 0.002). v_B estimated from ¹⁸F-FASu and ¹⁸F-FDG TACs was highly consistent and significantly correlated (r = 0.85, p < 0.001). Two qualitatively different ¹⁸F-FASu uptake profiles were identified: type α with low xCT expression and low K₁ (A549 and HBCx3), and type β with high xCT expression and high K₁ (MAS98.12 and TPMX).

CONCLUSION

The influx rate of ¹⁸F-FASu reflects xCT activity in tumour xenografts. Dynamic PET with pharmacokinetic modelling is needed to fully appraise ¹⁸F-FASu distribution routes.

Effects of Land Cover Pattern Along Urban-Rural Gradient on Bird Diversity in Wetlands

Wetlands play an important role in the feeding, breeding, and lives of birds. However, available habitats for bird species are changing due to intensifying human activity, especially in the context of China’s mass urbanization. Urban sprawl has taken over the wetlands along the lakes in the past decades, which places tremendous pressure on wetland ecosystems and, therefore, on bird communities. However, the ways urban land cover pattern along the urban-rural gradient affects bird communities is still unclear. To investigate the influence of land cover pattern on the α and β diversity of birds in the urban-rural gradient we chose 31 sites distributed within the wetlands around the Dianchi Lake in Yunnan, China. We calculated the species richness to indicate α diversity and used the Morisita–Horn index to indicate β diversity. Meanwhile, we assessed the land cover pattern of each site by measuring the proportion of emergent plants, floating plants, submerged plants, ponds, forests, lawns, roads, agricultural lands and built lands in a quadrat of 1 square kilometer. Simple linear regressions, model selection, and an averaging approach based on corrected Akaike information criterion (AICc) were used to test the effects of land cover pattern on bird diversity. Using one-way ANOVA and Tukey’s HSD (honestly significant difference) test, we compared the difference between α and β diversity, respectively, along the urban-rural gradient. Based on our analyses, urban and suburban wetland birds were significantly homogeneous. The community structure in rural wetlands, however, was significantly different from that of the suburban and urban areas. According to our research, the land cover patterns that influenced bird species richness were the built lands acreage, submerged plants acreage, ponds acreage, and the edge density of emergent plants. Meanwhile, of these variables, the built lands acreage, ponds acreage and edge density of emergent plants were significantly different in urban, suburban, and rural wetlands. Therefore, to maintain high biodiversity in wetlands affected by urbanization, we must pay more attention to the land cover patterns.

Seasonal influenza self-vaccination behaviours and attitudes among nurses in Southeastern France

BACKGROUND

Despite seasonal influenza vaccination (SIV) being recommended to healthcare professionals to protect themselves and their patients, uptake is low, especially among nurses. We sought to study self-vaccination behaviours, attitudes and knowledge about SIV among nurses in Southeastern France.

METHODS

A cross-sectional survey with community and hospital-based hospital nurses was conducted with the same standardised questionnaire. Multi-model averaging approaches studied factors associated with the following dependent variables: self-reported SIV uptake; and considering SIV a professional responsibility.

RESULTS

1539 nurses completed the questionnaire (response rate: 85%). SIV was the most frequently cited vaccine (49%) regarding nurses' unfavourable opinions towards specific vaccines. Thirty-four percent of nurses reported being vaccinated at least once during the 2015-2016 or 2016-2017 seasons. A lack of perceived personal vulnerability to influenza, a fear of adverse effects, and a preference for homeopathy constituted the main deterrents of SIV. Nurses held various misconceptions about the SIV, but 69% considered its benefits to be greater than its risks. The multi-model averaging approach showed that considering SIV as a professional responsibility was the main factor associated with SIV uptake among nurses (Nagelkerke's partial R-squared: 15%). This sense of responsibility was strongly associated with trust in various vaccine information sources.

CONCLUSION

Nurses had low SIV uptake rates and held various concerns and a lack of knowledge surrounding the vaccine. This is concerning considering the impact that these factors can have on nurses and patients' health, especially considering the increased role that nurses could have surrounding SIV in the near future.

Studentized bootstrap model-averaged tail area intervals

In many scientific studies, the underlying data-generating process is unknown and multiple statistical models are considered to describe it. For example, in a factorial experiment we might consider models involving just main effects, as well as those that include interactions. Model-averaging is a commonly-used statistical technique to allow for model uncertainty in parameter estimation. In the frequentist setting, the model-averaged estimate of a parameter is a weighted mean of the estimates from the individual models, with the weights typically being based on an information criterion, cross-validation, or bootstrapping. One approach to building a model-averaged confidence interval is to use a Wald interval, based on the model-averaged estimate and its standard error. This has been the default method in many application areas, particularly those in the life sciences. The MA-Wald interval, however, assumes that the studentized model-averaged estimate has a normal distribution, which can be far from true in practice due to the random, data-driven model weights. Recently, the model-averaged tail area Wald interval (MATA-Wald) has been proposed as an alternative to the MA-Wald interval, which only assumes that the studentized estimate from each model has a N(0, 1) or t-distribution, when that model is true. This alternative to the MA-Wald interval has been shown to have better coverage in simulation studies. However, when we have a response variable that is skewed, even these relaxed assumptions may not be valid, and use of these intervals might therefore result in poor coverage. We propose a new interval (MATA-SBoot) which uses a parametric bootstrap approach to estimate the distribution of the studentized estimate for each model, when that model is true. This method only requires that the studentized estimate from each model is approximately pivotal, an assumption that will often be true in practice, even for skewed data. We illustrate use of this new interval in the analysis of a three-factor marine global change experiment in which the response variable is assumed to have a lognormal distribution. We also perform a simulation study, based on the example, to compare the lower and upper error rates of this interval with those for existing methods. The results suggest that the MATA-SBoot interval can provide better error rates than existing intervals when we have skewed data, particularly for the upper error rate when the sample size is small.

Capacity for instantaneous catabolism of preferred and non-preferred carbon sources in Escherichia coli and Bacillus subtilis

Making the right choice for nutrient consumption in an ever-changing environment is a key factor for evolutionary success of bacteria. Here we investigate the regulatory mechanisms that enable dynamic adaptation between non-preferred and preferred carbon sources for the model Gram-negative and -positive species Escherichia coli and Bacillus subtilis, respectively. We focus on the ability for instantaneous catabolism of a gluconeogenic carbon source upon growth on a glycolytic carbon source and vice versa. By following isotopic tracer dynamics on a 1–2 minute scale, we show that flux reversal from the preferred glucose to non-preferred pyruvate as the sole carbon source is primarily transcriptionally regulated. In the opposite direction, however, E. coli can reverse its flux instantaneously by means of allosteric regulation, whereas in B. subtilis this flux reversal is transcriptionally regulated. Upon removal of transcriptional regulation, B. subtilis assumes the ability of instantaneous glucose catabolism. Using an approach that combines quantitative metabolomics and kinetic modelling, we then identify the additionally necessary key metabolite-enzyme interactions that implement the instantaneous flux reversal in the transcriptionally deregulated B. subtilis, and validate the most relevant allosteric interactions.

Reserve Flux Capacity in the Pentose Phosphate Pathway Enables Escherichia coli's Rapid Response to Oxidative Stress

To counteract oxidative stress and reactive oxygen species (ROS), bacteria evolved various mechanisms, primarily reducing ROS through antioxidant systems that utilize cofactor NADPH. Cells must stabilize NADPH levels by increasing flux through replenishing metabolic pathways like pentose phosphate (PP) pathway. Here, we investigate the mechanism enabling the rapid increase in NADPH supply by exposing Escherichia coli to hydrogen peroxide and quantifying the immediate metabolite dynamics. To systematically infer active regulatory interactions governing this response, we evaluated ensembles of kinetic models of glycolysis and PP pathway, each with different regulation mechanisms. Besides the known inactivation of glyceraldehyde 3-phosphate dehydrogenase by ROS, we reveal the important allosteric inhibition of the first PP pathway enzyme by NADPH. This NADPH feedback inhibition maintains a below maximum-capacity PP pathway flux under non-stress conditions. Relieving this inhibition instantly increases PP pathway flux upon oxidative stress. We demonstrate that reducing cells' capacity to rapidly reroute their flux through the PP pathway increases their oxidative stress sensitivity.

Phylogenetic surveys on the newt genus Tylototriton sensu lato (Salamandridae, Caudata) reveal cryptic diversity and novel diversification promoted by historical climatic shifts

Global climatic transitions and Tibetan Plateau uplifts are hypothesized to have profoundly impacted biodiversity in southeastern Asia. To further test the hypotheses related to the impacts of these incidents, we investigated the diversification patterns of the newt genus Tylototriton sensu lato, distributed across the mountain ranges of southeastern Asia. Gene-tree and species-tree analyses of two mitochondrial genes and two nuclear genes revealed five major clades in the genus, and suggested several cryptic species. Dating estimates suggested that the genus originated in the early-to-middle Miocene. Under different species delimitating scenarios, diversification analyses with birth-death likelihood tests indicated that the genus held a higher diversification rate in the late Miocene-to-Pliocene era than that in the Pleistocene. Ancestral area reconstructions indicated that the genus originated from the northern Indochina Peninsula. Accordingly, we hypothesized that the Miocene Climatic Transition triggered the diversification of the genus, and the reinforcement of East Asian monsoons associated with the stepwise uplifts of the Tibetan Plateau promoted the radiation of the genus in southeastern Asia during the Miocene-to-Pliocene period. Quaternary glacial cycles likely had limited effects on speciation events in the genus, but mainly had contributions on their intraspecific differentiations.

Parametric mapping using spectral analysis for 11C-PBR28 PET reveals neuroinflammation in mild cognitive impairment subjects

PURPOSE

Neuroinflammation and microglial activation play an important role in amnestic mild cognitive impairment (MCI) and Alzheimer's disease. In this study, we investigated the spatial distribution of neuroinflammation in MCI subjects, using spectral analysis (SA) to generate parametric maps and quantify ¹¹C-PBR28 PET, and compared these with compartmental and other kinetic models of quantification.

METHODS

Thirteen MCI and nine healthy controls were enrolled in this study. Subjects underwent ¹¹C-PBR28 PET scans with arterial cannulation. Spectral analysis with an arterial plasma input function was used to generate ¹¹C-PBR28 parametric maps. These maps were then compared with regional ¹¹C-PBR28 V_T (volume of distribution) using a two-tissue compartment model and Logan graphic analysis. Amyloid load was also assessed with ¹⁸F-Flutemetamol PET.

RESULTS

With SA, three component peaks were identified in addition to blood volume. The ¹¹C-PBR28 impulse response function (IRF) at 90 min produced the lowest coefficient of variation. Single-subject analysis using this IRF demonstrated microglial activation in five out of seven amyloid-positive MCI subjects. IRF parametric maps of ¹¹C-PBR28 uptake revealed a group-wise significant increase in neuroinflammation in amyloid-positive MCI subjects versus HC in multiple cortical association areas, and particularly in the temporal lobe. Interestingly, compartmental analysis detected group-wise increase in ¹¹C-PBR28 binding in the thalamus of amyloid-positive MCI subjects, while Logan parametric maps did not perform well.

CONCLUSIONS

This study demonstrates for the first time that spectral analysis can be used to generate parametric maps of ¹¹C-PBR28 uptake, and is able to detect microglial activation in amyloid-positive MCI subjects. IRF parametric maps of ¹¹C-PBR28 uptake allow voxel-wise single-subject analysis and could be used to evaluate microglial activation in individual subjects.

Evaluation of fitting models for prostate tissue characterization using extended-range b-factor diffusion-weighted imaging

PURPOSE

To compare the fitting and tissue discrimination performance of biexponential, kurtosis, stretched exponential, and gamma distribution models for high b-factor diffusion-weighted images in prostate cancer.

METHODS

Diffusion-weighted images with 15 b-factors ranging from b = 0 to 3500 s/mm² were obtained in 62 prostate cancer patients. Pixel-wise signal decay fits for each model were evaluated with the Akaike Information Criterion (AIC). Parameter values for each model were determined within normal prostate and the index lesion. Their potential to differentiate normal from cancerous tissue was investigated through receiver operating characteristic analysis and comparison with Gleason score.

RESULTS

The biexponential slow diffusion fraction f_slow , the apparent kurtosis diffusion coefficient ADC_K , and the excess kurtosis factor K differ significantly among normal peripheral zone (PZ), normal transition zone (TZ), tumor PZ, and tumor TZ. Biexponential and gamma distribution models result in the lowest AIC, indicating a superior fit. Maximum areas under the curve (AUCs) of all models ranged from 0.93 to 0.96 for the PZ and from 0.95 to 0.97 for the TZ. Similar AUCs also result from the apparent diffusion coefficient (ADC) of a monoexponential fit to a b-factor sub-range up to 1250 s/mm² . For kurtosis and stretched exponential models, single parameters yield the highest AUCs, whereas for the biexponential and gamma distribution models, linear combinations of parameters produce the highest AUCs. Parameters with high AUC show a trend in differentiating low from high Gleason score, whereas parameters with low AUC show no such ability.

CONCLUSION

All models, including a monoexponential fit to a lower-b sub-range, achieve similar AUCs for discrimination of normal and cancer tissue. The biexponential model, which is favored statistically, also appears to provide insight into disease-related microstructural changes. Magn Reson Med 79:2346-2358, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

An integrated network analysis reveals that nitric oxide reductase prevents metabolic cycling of nitric oxide by Pseudomonas aeruginosa

Nitric oxide (NO) is a chemical weapon within the arsenal of immune cells, but is also generated endogenously by different bacteria. Pseudomonas aeruginosa are pathogens that contain an NO-generating nitrite (NO₂^-) reductase (NirS), and NO has been shown to influence their virulence. Interestingly, P. aeruginosa also contain NO dioxygenase (Fhp) and nitrate (NO₃^-) reductases, which together with NirS provide the potential for NO to be metabolically cycled (NO→NO₃^-→NO₂^-→NO). Deeper understanding of NO metabolism in P. aeruginosa will increase knowledge of its pathogenesis, and computational models have proven to be useful tools for the quantitative dissection of NO biochemical networks. Here we developed such a model for P. aeruginosa and confirmed its predictive accuracy with measurements of NO, O₂, NO₂^-, and NO₃^- in mutant cultures devoid of Fhp or NorCB (NO reductase) activity. Using the model, we assessed whether NO was metabolically cycled in aerobic P. aeruginosa cultures. Calculated fluxes indicated a bottleneck at NO₃^-, which was relieved upon O₂ depletion. As cell growth depleted dissolved O₂ levels, NO₃^- was converted to NO₂^- at near-stoichiometric levels, whereas NO₂^- consumption did not coincide with NO or NO₃^- accumulation. Assimilatory NO₂^- reductase (NirBD) or NorCB activity could have prevented NO cycling, and experiments with ΔnirB, ΔnirS, and ΔnorC showed that NorCB was responsible for loss of flux from the cycle. Collectively, this work provides a computational tool to analyze NO metabolism in P. aeruginosa, and establishes that P. aeruginosa use NorCB to prevent metabolic cycling of NO.

Discovery and dissection of metabolic oscillations in the microaerobic nitric oxide response network of Escherichia coli

The virulence of many pathogens depends upon their ability to cope with immune-generated nitric oxide (NO·). In Escherichia coli, the major NO· detoxification systems are Hmp, an NO· dioxygenase (NOD), and NorV, an NO· reductase (NOR). It is well established that Hmp is the dominant system under aerobic conditions, whereas NorV dominates anaerobic conditions; however, the quantitative contributions of these systems under the physiologically relevant microaerobic regime remain ill defined. Here, we investigated NO· detoxification in environments ranging from 0 to 50 μM O2, and discovered a regime in which E. coli NO· defenses were severely compromised, as well as conditions that exhibited oscillations in the concentration of NO·. Using an integrated computational and experimental approach, E. coli NO· detoxification was found to be extremely impaired at low O2 due to a combination of its inhibitory effects on NorV, Hmp, and translational activities, whereas oscillations were found to result from a kinetic competition for O2 between Hmp and respiratory cytochromes. Because at least 777 different bacterial species contain the genetic requirements of this stress response oscillator, we hypothesize that such oscillatory behavior could be a widespread phenomenon. In support of this hypothesis,Pseudomonas aeruginosa, whose respiratory and NO· response networks differ considerably from those of E. coli, was found to exhibit analogous oscillations in low O2 environments. This work provides insight into how bacterial NO· defenses function under the low O2 conditions that are likely to be encountered within host environments.

Construction and Experimental Validation of a Quantitative Kinetic Model of Nitric Oxide Stress in Enterohemorrhagic Escherichia coli O157:H7

Enterohemorrhagic Escherichia coli (EHEC) are responsible for large outbreaks of hemorrhagic colitis, which can progress to life-threatening hemolytic uremic syndrome (HUS) due to the release of Shiga-like toxins (Stx). The presence of a functional nitric oxide (NO·) reductase (NorV), which protects EHEC from NO· produced by immune cells, was previously found to correlate with high HUS incidence, and it was shown that NorV activity enabled prolonged EHEC survival and increased Stx production within macrophages. To enable quantitative study of EHEC NO· defenses and facilitate the development of NO·-potentiating therapeutics, we translated an existing kinetic model of the E. coli K-12 NO· response to an EHEC O157:H7 strain. To do this, we trained uncertain model parameters on measurements of [NO·] and [O₂] in EHEC cultures, assessed parametric and prediction uncertainty with the use of a Markov chain Monte Carlo approach, and confirmed the predictive accuracy of the model with experimental data from genetic mutants lacking NorV or Hmp (NO· dioxygenase). Collectively, these results establish a methodology for the translation of quantitative models of NO· stress in model organisms to pathogenic sub-species, which is a critical step toward the application of these models for the study of infectious disease.

General practitioners' attitudes and behaviors toward HPV vaccination: A French national survey

OBJECTIVE

General practitioners (GPs) play a crucial role in human papillomavirus (HPV) vaccine acceptance in France. We sought to study: (1) GPs' perceptions of its risks and efficacy and their recommendation behavior; (2) the relative importance of factors associated with the frequency of their recommendations.

METHODS

Cross-sectional observational study in 2014 nested in a national panel of 1712 randomly selected GPs in private practice in France (response rate: 92.4%). We used model averaging to analyze the associations of self-reported frequency of GPs' HPV vaccine recommendations with their perception of its risk-benefit balance and their opinions about the utility of vaccines in general.

RESULTS

Overall, 72% of participants reported frequently recommending HPV vaccination; 60% considered that not enough is known about its risks. The model averaging showed that the factors most associated with infrequent recommendation of this vaccine by GPs were: unfavorable perceptions of its risk-benefit balance (OR=0.13; 95%CI=0.09-0.21; partial R(2)=0.10), a decision not to vaccinate one's own daughter(s) with this vaccine (OR=0.13; 95%CI=0.07-0.24; partial R(2)=0.05), and doubts about vaccine utility in general (OR=0.78; 95%CI=0.71-0.86; partial R(2)=0.03).

CONCLUSION

Although nearly three-quarters of French GPs frequently recommended the HPV vaccine, our findings indicate that a substantial percentage of them are hesitant about it. Doubts about its risks and efficacy strongly influence their recommendation behavior. More research is warranted to help design and evaluate tailored tools and multicomponent intervention strategies to address physician's hesitancy about this vaccine.

Thermal tolerance and climate warming sensitivity in tropical snails

Tropical ectotherms are predicted to be especially vulnerable to climate change because their thermal tolerance limits generally lie close to current maximum air temperatures. This prediction derives primarily from studies on insects and lizards and remains untested for other taxa with contrasting ecologies. We studied the HCT (heat coma temperatures) and ULT (upper lethal temperatures) of 40 species of tropical eulittoral snails (Littorinidae and Neritidae) inhabiting exposed rocky shores and shaded mangrove forests in Oceania, Africa, Asia and North America. We also estimated extremes in animal body temperature at each site using a simple heat budget model and historical (20 years) air temperature and solar radiation data. Phylogenetic analyses suggest that HCT and ULT exhibit limited adaptive variation across habitats (mangroves vs. rocky shores) or geographic locations despite their contrasting thermal regimes. Instead, the elevated heat tolerance of these species (HCT = 44.5 ± 1.8°C and ULT = 52.1 ± 2.2°C) seems to reflect the extreme temperature variability of intertidal systems. Sensitivity to climate warming, which was quantified as the difference between HCT or ULT and maximum body temperature, differed greatly between snails from sunny (rocky shore; Thermal Safety Margin, TSM = -14.8 ± 3.3°C and -6.2 ± 4.4°C for HCT and ULT, respectively) and shaded (mangrove) habitats (TSM = 5.1 ± 3.6°C and 12.5 ± 3.6°C). Negative TSMs in rocky shore animals suggest that mortality is likely ameliorated during extreme climatic events by behavioral thermoregulation. Given the low variability in heat tolerance across species, habitat and geographic location account for most of the variation in TSM and may adequately predict the vulnerability to climate change. These findings caution against generalizations on the impact of global warming across ectothermic taxa and highlight how the consideration of nonmodel animals, ecological transitions, and behavioral responses may alter predictions of studies that ignore these biological details.

A Kinetic Platform to Determine the Fate of Hydrogen Peroxide in Escherichia coli

Hydrogen peroxide (H2O2) is used by phagocytic cells of the innate immune response to kill engulfed bacteria. H2O2 diffuses freely into bacteria, where it can wreak havoc on sensitive biomolecules if it is not rapidly detoxified. Accordingly, bacteria have evolved numerous systems to defend themselves against H2O2, and the importance of these systems to pathogenesis has been substantiated by the many bacteria that require them to establish or sustain infections. The kinetic competition for H2O2 within bacteria is complex, which suggests that quantitative models will improve interpretation and prediction of network behavior. To date, such models have been of limited scope, and this inspired us to construct a quantitative, systems-level model of H2O2 detoxification in Escherichia coli that includes detoxification enzymes, H2O2-dependent transcriptional regulation, enzyme degradation, the Fenton reaction and damage caused by •OH, oxidation of biomolecules by H2O2, and repair processes. After using an iterative computational and experimental procedure to train the model, we leveraged it to predict how H2O2 detoxification would change in response to an environmental perturbation that pathogens encounter within host phagosomes, carbon source deprivation, which leads to translational inhibition and limited availability of NADH. We found that the model accurately predicted that NADH depletion would delay clearance at low H2O2 concentrations and that detoxification at higher concentrations would resemble that of carbon-replete conditions. These results suggest that protein synthesis during bolus H2O2 stress does not affect clearance dynamics and that access to catabolites only matters at low H2O2 concentrations. We anticipate that this model will serve as a computational tool for the quantitative exploration and dissection of oxidative stress in bacteria, and that the model and methods used to develop it will provide important templates for the generation of comparable models for other bacterial species.

Accuracies and Contrasts of Models of the Diffusion-Weighted-Dependent Attenuation of the MRI Signal at Intermediate b-values

The diffusion-weighted-dependent attenuation of the MRI signal E(b) is extremely sensitive to microstructural features. The aim of this study was to determine which mathematical model of the E(b) signal most accurately describes it in the brain. The models compared were the monoexponential model, the stretched exponential model, the truncated cumulant expansion (TCE) model, the biexponential model, and the triexponential model. Acquisition was performed with nine b-values up to 2500 s/mm² in 12 healthy volunteers. The goodness-of-fit was studied with F-tests and with the Akaike information criterion. Tissue contrasts were differentiated with a multiple comparison corrected nonparametric analysis of variance. F-test showed that the TCE model was better than the biexponential model in gray and white matter. Corrected Akaike information criterion showed that the TCE model has the best accuracy and produced the most reliable contrasts in white matter among all models studied. In conclusion, the TCE model was found to be the best model to infer the microstructural properties of brain tissue.

Kinetic Analysis and Quantification of [¹¹C]Martinostat for in Vivo HDAC Imaging of the Brain

Epigenetic mechanisms mediated by histone deacetylases (HDACs) have been implicated in a wide-range of CNS disorders and may offer new therapeutic opportunities. In vivo evaluation of HDAC density and drug occupancy has become possible with [(11)C]Martinostat, which exhibits selectivity for a subset of class I/IIb HDAC enzymes. In this study, we characterize the kinetic properties of [(11)C]Martinostat in the nonhuman primate (NHP) brain in preparation for human neuroimaging studies. The goal of this work was to determine whether classic compartmental analysis techniques were appropriate and to further determine if arterial plasma is required for future NHP studies. Using an arterial plasma input function, several analysis approaches were evaluated for robust outcome measurements. [(11)C]Martinostat showed high baseline distribution volume (VT) ranging from 29.9 to 54.4 mL/cm(3) in the brain and large changes in occupancy (up to 99%) with a blocking dose approaching full enzyme saturation. An averaged nondisplaceable tissue uptake (VND) of 8.6 ± 3.7 mL/cm(3) suggests high specific binding of [(11)C]Martinostat. From a two-tissue compartment model, [(11)C]Martinostat exhibits a high K1 (averaged K1 of 0.65 mL/cm(3)/min) and a small k4 (average of 0.0085 min(-1)). Our study supports that [(11)C]Martinostat can be used to detect changes in HDAC density and occupancy in vivo and that simplified analysis not using arterial blood could be appropriate.

Pharmacokinetic Analysis of (64)Cu-ATSM Dynamic PET in Human Xenograft Tumors in Mice

The aim of this study was to evaluate the feasibility to perform voxel-wise kinetic modeling on datasets obtained from tumor-bearing mice that underwent dynamic PET scans with (64)Cu-ATSM and extract useful physiological parameters.

METHODS

Tumor-bearing mice underwent 90-min dynamic PET scans with (64)Cu-ATSM and CT scans with contrast. Irreversible and reversible two-tissue compartment models were fitted to time activity curves (TACs) obtained from whole tumor volumes and compared using the Akaike information criterion (AIC). Based on voxel-wise pharmacokinetic analysis, parametric maps of model rate constants k₁, k₃ and Ki were generated and compared to (64)Cu-ATSM uptake.

RESULTS

Based on the AIC, an irreversible two-tissue compartment model was selected for voxel-wise pharmacokinetic analysis. Of the extracted parameters, k₁ (~perfusion) showed a strong correlation with early tracer uptake (mean spearman R = 0.88) 5 min post injection (pi). Moreover, positive relationships were found between late tracer uptake (90 min pi) and both k₃ and the net influx rate constant, Ki (mean spearman R = 0.56 and R = 0.86; respectively).

CONCLUSION

This study shows the feasibility to extract relevant parameters from voxel-wise pharmacokinetic analysis to be used for preclinical validation of (64)Cu-ATSM as a hypoxia-specific PET tracer.

A spatially explicit metapopulation model and cattle trade analysis suggests key determinants for the recurrent circulation of rift valley Fever virus in a pilot area of madagascar highlands

Rift Valley fever (RVF) is a vector-borne zoonotic disease that causes high morbidity and mortality in ruminants. In 2008-2009, a RVF outbreak affected the whole Madagascar island, including the Anjozorobe district located in Madagascar highlands. An entomological survey showed the absence of Aedes among the potential RVF virus (RVFV) vector species identified in this area, and an overall low abundance of mosquitoes due to unfavorable climatic conditions during winter. No serological nor virological sign of infection was observed in wild terrestrial mammals of the area, suggesting an absence of wild RVF virus (RVFV) reservoir. However, a three years serological and virological follow-up in cattle showed a recurrent RVFV circulation. The objective of this study was to understand the key determinants of this unexpected recurrent transmission. To achieve this goal, a spatial deterministic discrete-time metapopulation model combined with cattle trade network was designed and parameterized to reproduce the local conditions using observational data collected in the area. Three scenarios that could explain the RVFV recurrent circulation in the area were analyzed: (i) RVFV overwintering thanks to a direct transmission between cattle when viraemic cows calve, vectors being absent during the winter, (ii) a low level vector-based circulation during winter thanks to a residual vector population, without direct transmission between cattle, (iii) combination of both above mentioned mechanisms. Multi-model inference methods resulted in a model incorporating both a low level RVFV winter vector-borne transmission and a direct transmission between animals when viraemic cows calve. Predictions satisfactorily reproduced field observations, 84% of cattle infections being attributed to vector-borne transmission, and 16% to direct transmission. These results appeared robust according to the sensitivity analysis. Interweaving between agricultural works in rice fields, seasonality of vector proliferation, and cattle exchange practices could be a key element for understanding RVFV circulation in this area of Madagascar highlands.

HCV seropositivity in inmates and in the general population: an averaging approach to establish priority prevention interventions

OBJECTIVES

Despite the fact that a considerable portion of hepatitis C virus (HCV) positive individuals are viraemic, the risk of transmitting HCV to others is context dependent. Prison is a particularly risky environment as HCV prevention tools are often unavailable. Using data from a cross-sectional study conducted in centres for HCV testing in southeastern France, we aimed to compare the patterns of risk factors in HCV-positive inmates with those in the general population.

SETTING

26 centres for HIV/HCV testing in southeastern France (23 in the general population and 3 in prison).

PRIMARY OUTCOME MEASURE

HCV seropositivity measured with ELISA test.

METHODS

A propensity score method to ensure that the general and inmate populations could be compared and a multimodel averaging to estimate the degree (strong, weak, none) of the association of a number of specific factors with HCV seropositivity in each group.

RESULTS

Among the 52,082 participants, HCV infection prevalence was 1.5% and 5.2% in the general (n=46,125) and inmate (n=5957) populations, respectively. In both populations, 'drug injection without snorting' and 'drug injection with snorting' were very strongly associated with HCV seropositivity. Among inmates, 'drug snorting alone' (OR (95% CI) 2.21 (1.39 to 3.52) was also a strong correlate while tattoos, piercings (OR (95% CI) 1.22 (0.92 to 1.61)) and the sharing of toiletry items (OR (95% CI) 1.44 (0.84 to 2.47)) were weak correlates.

CONCLUSIONS

The pattern of risk factors associated with HCV seropositivity is different between the general and prison populations, injection and snorting practices being more prevalent in the latter. Access to prevention measures in prisons is not only a public health issue but also a human right for inmates who deserve equity of care and prevention.

Cutting the wires: modularization of cellular networks for experimental design

Understanding naturally evolved cellular networks requires the consecutive identification and revision of the interactions between relevant molecular species. In this process, initially often simplified and incomplete networks are extended by integrating new reactions or whole subnetworks to increase consistency between model predictions and new measurement data. However, increased consistency with experimental data alone is not sufficient to show the existence of biomolecular interactions, because the interplay of different potential extensions might lead to overall similar dynamics. Here, we present a graph-based modularization approach to facilitate the design of experiments targeted at independently validating the existence of several potential network extensions. Our method is based on selecting the outputs to measure during an experiment, such that each potential network extension becomes virtually insulated from all others during data analysis. Each output defines a module that only depends on one hypothetical network extension, and all other outputs act as virtual inputs to achieve insulation. Given appropriate experimental time-series measurements of the outputs, our modules can be analyzed, simulated, and compared to the experimental data separately. Our approach exemplifies the close relationship between structural systems identification and modularization, an interplay that promises development of related approaches in the future.

How informative is your kinetic model?: using resampling methods for model invalidation

Background

Kinetic models can present mechanistic descriptions of molecular processes within a cell. They can be used to predict the dynamics of metabolite production, signal transduction or transcription of genes. Although there has been tremendous effort in constructing kinetic models for different biological systems, not much effort has been put into their validation. In this study, we introduce the concept of resampling methods for the analysis of kinetic models and present a statistical model invalidation approach.

Results

We based our invalidation approach on the evaluation of a kinetic model’s predictive power through cross validation and forecast analysis. As a reference point for this evaluation, we used the predictive power of an unsupervised data analysis method which does not make use of any biochemical knowledge, namely Smooth Principal Components Analysis (SPCA) on the same test sets. Through a simulations study, we showed that too simple mechanistic descriptions can be invalidated by using our SPCA-based comparative approach until high amount of noise exists in the experimental data. We also applied our approach on an eicosanoid production model developed for human and concluded that the model could not be invalidated using the available data despite its simplicity in the formulation of the reaction kinetics. Furthermore, we analysed the high osmolarity glycerol (HOG) pathway in yeast to question the validity of an existing model as another realistic demonstration of our method.

Conclusions

With this study, we have successfully presented the potential of two resampling methods, cross validation and forecast analysis in the analysis of kinetic models’ validity. Our approach is easy to grasp and to implement, applicable to any ordinary differential equation (ODE) type biological model and does not suffer from any computational difficulties which seems to be a common problem for approaches that have been proposed for similar purposes. Matlab files needed for invalidation using SPCA cross validation and our toy model in SBML format are provided at http://www.bdagroup.nl/content/Downloads/software/software.php.

Three-dimensional whole-brain perfusion quantification using pseudo-continuous arterial spin labeling MRI at multiple post-labeling delays: accounting for both arterial transit time and impulse response function

Measurement of the cerebral blood flow (CBF) with whole-brain coverage is challenging in terms of both acquisition and quantitative analysis. In order to fit arterial spin labeling-based perfusion kinetic curves, an empirical three-parameter model which characterizes the effective impulse response function (IRF) is introduced, which allows the determination of CBF, the arterial transit time (ATT) and T(1,eff). The accuracy and precision of the proposed model were compared with those of more complicated models with four or five parameters through Monte Carlo simulations. Pseudo-continuous arterial spin labeling images were acquired on a clinical 3-T scanner in 10 normal volunteers using a three-dimensional multi-shot gradient and spin echo scheme at multiple post-labeling delays to sample the kinetic curves. Voxel-wise fitting was performed using the three-parameter model and other models that contain two, four or five unknown parameters. For the two-parameter model, T(1,eff) values close to tissue and blood were assumed separately. Standard statistical analysis was conducted to compare these fitting models in various brain regions. The fitted results indicated that: (i) the estimated CBF values using the two-parameter model show appreciable dependence on the assumed T(1,eff) values; (ii) the proposed three-parameter model achieves the optimal balance between the goodness of fit and model complexity when compared among the models with explicit IRF fitting; (iii) both the two-parameter model using fixed blood T1 values for T(1,eff) and the three-parameter model provide reasonable fitting results. Using the proposed three-parameter model, the estimated CBF (46 ± 14 mL/100 g/min) and ATT (1.4 ± 0.3 s) values averaged from different brain regions are close to the literature reports; the estimated T(1,eff) values (1.9 ± 0.4 s) are higher than the tissue T1 values, possibly reflecting a contribution from the microvascular arterial blood compartment.

Visualization, imaging and new preclinical diagnostics in radiation oncology

Innovative strategies in cancer radiotherapy are stimulated by the growing knowledge on cellular and molecular tumor biology, tumor pathophysiology, and tumor microenvironment. In terms of tumor diagnostics and therapy monitoring, the reliable delineation of tumor boundaries and the assessment of tumor heterogeneity are increasingly complemented by the non-invasive characterization of functional and molecular processes, moving preclinical and clinical imaging from solely assessing tumor morphology towards the visualization of physiological and pathophysiological processes. Functional and molecular imaging techniques allow for the non-invasive characterization of tissues in vivo, using different modalities, including computed tomography (CT), magnetic resonance imaging (MRI), ultrasound, positron emission tomography (PET) and optical imaging (OI). With novel therapeutic concepts combining optimized radiotherapy with molecularly targeted agents focusing on tumor cell proliferation, angiogenesis, and cell death, the non-invasive assessment of tumor microcirculation and tissue water diffusion, together with strategies for imaging the mechanisms of cellular injury and repair is of particular interest. Characterizing the tumor microenvironment prior to and in response to irradiation will help to optimize the outcome of radiotherapy. These novel concepts of personalized multi-modal cancer therapy require careful pre-treatment stratification as well as a timely and efficient therapy monitoring to maximize patient benefit on an individual basis. Functional and molecular imaging techniques are key in this regard to open novel opportunities for exploring and understanding the underlying mechanisms with the perspective to optimize therapeutic concepts and translate them into a personalized form of radiotherapy in the near future.

Phylogenetic analyses: comparing species to infer adaptations and physiological mechanisms

Comparisons among species have been a standard tool in animal physiology to understand how organisms function and adapt to their surrounding environment. During the last two decades, conceptual and methodological advances from different fields, including evolutionary biology and systematics, have revolutionized the way comparative analyses are performed, resulting in the advent of modern phylogenetic statistical methods. This development stems from the realization that conventional analytical methods assume that observations are statistically independent, which is not the case for comparative data because species often resemble each other due to shared ancestry. By taking evolutionary history explicitly into consideration, phylogenetic statistical methods can account for the confounding effects of shared ancestry in interspecific comparisons, improving the reliability of standard approaches such as regressions or correlations in comparative analyses. Importantly, these methods have also enabled researchers to address entirely new evolutionary questions, such as the historical sequence of events that resulted in current patterns of form and function, which can only be studied with a phylogenetic perspective. Here, we provide an overview of phylogenetic approaches and their importance for studying the evolution of physiological processes and mechanisms. We discuss the conceptual framework underlying these methods, and explain when and how phylogenetic information should be employed. We then outline the difficulties and limitations inherent to comparative approaches and discuss potential problems researchers may encounter when designing a comparative study. These issues are illustrated with examples from the literature in which the incorporation of phylogenetic information has been useful, or even crucial, for inferences on how species evolve and adapt to their surrounding environment.

Heterogeneity in lung (18)FDG uptake in pulmonary arterial hypertension: potential of dynamic (18)FDG positron emission tomography with kinetic analysis as a bridging biomarker for pulmonary vascular remodeling targeted treatments

BACKGROUND

Pulmonary arterial hypertension (PAH) is a disease of progressive vascular remodeling, characterized by dysregulated growth of pulmonary vascular cells and inflammation. A prevailing view is that abnormal cellular metabolism, notably aerobic glycolysis that increases glucose demand, underlies the pathogenesis of PAH. Increased lung glucose uptake has been reported in animal models. Few data exist from patients with PAH.

METHODS AND RESULTS

Dynamic positron emission tomography imaging with fluorine-18-labeled 2-fluoro-2-deoxyglucose ((18)FDG) ligand with kinetic analysis demonstrated increased mean lung parenchymal uptake in 20 patients with PAH, 18 with idiopathic PAH (IPAH) (FDG score: 3.27±1.22), and 2 patients with connective tissue disease (5.07 and 7.11) compared with controls (2.02±0.71; P<0.05). Further compartment analysis confirmed increased lung glucose metabolism in IPAH. Lung (18)FDG uptake and metabolism varied within the IPAH population and within the lungs of individual patients, consistent with the recognized heterogeneity of vascular pathology in this disease. The monocrotaline rat PAH model also showed increased lung (18)FDG uptake, which was reduced along with improvements in vascular pathology after treatment with dicholoroacetate and 2 tyrosine kinase inhibitors, imatinib and sunitinib. Hyperproliferative pulmonary vascular fibroblasts isolated from IPAH patients exhibited upregulated glycolytic gene expression, along with increased cellular (18)FDG uptake; both were reduced by dicholoroacetate and imatinib.

CONCLUSIONS

Some patients with IPAH exhibit increased lung (18)FDG uptake. (18)FDG positron emission tomography imaging is a tool to investigate the molecular pathology of PAH and its response to treatment.

Probabilistic inference of short-term synaptic plasticity in neocortical microcircuits

Short-term synaptic plasticity is highly diverse across brain area, cortical layer, cell type, and developmental stage. Since short-term plasticity (STP) strongly shapes neural dynamics, this diversity suggests a specific and essential role in neural information processing. Therefore, a correct characterization of short-term synaptic plasticity is an important step towards understanding and modeling neural systems. Phenomenological models have been developed, but they are usually fitted to experimental data using least-mean-square methods. We demonstrate that for typical synaptic dynamics such fitting may give unreliable results. As a solution, we introduce a Bayesian formulation, which yields the posterior distribution over the model parameters given the data. First, we show that common STP protocols yield broad distributions over some model parameters. Using our result we propose a experimental protocol to more accurately determine synaptic dynamics parameters. Next, we infer the model parameters using experimental data from three different neocortical excitatory connection types. This reveals connection-specific distributions, which we use to classify synaptic dynamics. Our approach to demarcate connection-specific synaptic dynamics is an important improvement on the state of the art and reveals novel features from existing data.