Update on the non-prewhitening model observer in computed tomography for the assessment of the adaptive statistical and model-based iterative reconstruction algorithms.

The state of the art to describe image quality in medical imaging is to assess the performance of an observer conducting a task of clinical interest. This can be done by using a model observer leading to a figure of merit such as the signal-to-noise ratio (SNR). Using the non-prewhitening (NPW) model observer, we objectively characterised the evolution of its figure of merit in various acquisition conditions. The NPW model observer usually requires the use of the modulation transfer function (MTF) as well as noise power spectra. However, although the computation of the MTF poses no problem when dealing with the traditional filtered back-projection (FBP) algorithm, this is not the case when using iterative reconstruction (IR) algorithms, such as adaptive statistical iterative reconstruction (ASIR) or model-based iterative reconstruction (MBIR). Given that the target transfer function (TTF) had already shown it could accurately express the system resolution even with non-linear algorithms, we decided to tune the NPW model observer, replacing the standard MTF by the TTF. It was estimated using a custom-made phantom containing cylindrical inserts surrounded by water. The contrast differences between the inserts and water were plotted for each acquisition condition. Then, mathematical transformations were performed leading to the TTF. As expected, the first results showed a dependency of the image contrast and noise levels on the TTF for both ASIR and MBIR. Moreover, FBP also proved to be dependent of the contrast and noise when using the lung kernel. Those results were then introduced in the NPW model observer. We observed an enhancement of SNR every time we switched from FBP to ASIR to MBIR. IR algorithms greatly improve image quality, especially in low-dose conditions. Based on our results, the use of MBIR could lead to further dose reduction in several clinical applications.

Continuous monitoring and quantification of multiple parameters of daily physical activity in ambulatory Duchenne muscular dystrophy patients.

Multiple motor function and strength assessment tools exist for the evaluation of neuromuscular diseases, but most do not directly assess functional ability in the patients' daily physical activity in their home environment. In this study our aim was to assess: 1) the feasibility and accuracy of physical activity monitoring during two days in a home environment of five DMD patients using a non-commercialized monitor containing a 3D accelerometer and a gyroscope, 2) if a difference in the physical activity parameters could be measured before and one month after starting prednisolone. We reliably quantified the time spend sitting, standing, lying, walking, the number of steps taken, the cadence, the number of walking episodes and their duration as well as how these were distributed over the day. Parameters possibly reflecting endurance, such as the duration of the walking episodes or the succession of two or three walking episodes lasting more than 30 s were the most improved after prednisolone treatment. This degree of detailed determination of physical activity in a home environment has not been previously reported in neuromuscular disorders to our knowledge and some of the reported parameters are potential new outcome measures in clinical trials.

Pharmacokinetic parameters of artesunate and dihydroartemisinin in rats infected with Fasciola hepatica

OBJECTIVES: The pharmacokinetic (PK) parameters of artesunate, recently discovered to possess promising trematocidal activity, and its main metabolite dihydroartemisinin (DHA) were determined in rats infected with hepatic and biliary stages of Fasciola hepatica and compared with uninfected rats after single intragastric and intravenous (iv) doses. METHODS: Rats infected with F. hepatica for 25 and 83 days and uninfected rats were cannulated in the right jugular vein and blood samples were withdrawn at selected timepoints following 10 mg/kg of iv and a single 100 mg/kg oral dose of artesunate. Plasma was analysed for artesunate and DHA by liquid chromatography coupled to tandem mass spectrometry. RESULTS: Rats harbouring juvenile and adult F. hepatica infections revealed considerable changes in PK parameters of artesunate and DHA. Following oral administration, maximum plasma concentrations (C(max)) of artesunate and DHA were 1.8-2.3-fold higher in infected rats [artesunate: 1334 +/- 1404 ng/mL (no infection) versus 2454 +/- 1494 ng/mL (acute infection) and 2768 +/- 538 ng/mL (chronic infection); DHA: 3802 +/- 2149 ng/mL (no infection) versus 6507 +/- 3283 ng/mL (acute infection) and 9093 +/- 884 ng/mL (chronic infection)]. The AUCs of artesunate and DHA were 2.1-4.4-fold greater in infected rats. An opposite trend was observed after iv injection. C(max) and AUC of artesunate and DHA following iv dosing were 5784 +/- 3718 and 140 938 +/- 128 783 ng.min/mL and 3849 +/- 3060 and 86 107 +/- 41 863 ng.min/mL, respectively, in uninfected rats versus 2623 +/- 1554 and 21 617 +/- 12 230 ng.min/mL and 2835 +/- 980 and 64 290 +/- 29 057 ng.min/mL, respectively, in rats harbouring a chronic infection. The elimination half-lives (t(1/2)) of artesunate and DHA were considerably altered in infected rats following oral and iv administration of artesunate. CONCLUSIONS: F. hepatica infections strongly influence the disposition kinetics of artesunate and its metabolite in rats. The clinical implications of this finding need to be carefully studied.

Adaptive tuning of perceptual timing to whole body motion.

In a previous work we have shown that sinusoidal whole-body rotations producing continuous vestibular stimulation, affected the timing of motor responses as assessed with a paced finger tapping (PFT) task (Binetti et al. (2010). Neuropsychologia, 48(6), 1842-1852). Here, in two new psychophysical experiments, one purely perceptual and one with both sensory and motor components, we explored the relationship between body motion/vestibular stimulation and perceived timing of acoustic events. In experiment 1, participants were required to discriminate sequences of acoustic tones endowed with different degrees of acceleration or deceleration. In this experiment we found that a tone sequence presented during acceleratory whole-body rotations required a progressive increase in rate in order to be considered temporally regular, consistent with the idea of an increase in "clock" frequency and of an overestimation of time. In experiment 2 participants produced self-paced taps, which entailed an acoustic feedback. We found that tapping frequency in this task was affected by periodic motion by means of anticipatory and congruent (in-phase) fluctuations irrespective of the self-generated sensory feedback. On the other hand, synchronizing taps to an external rhythm determined a completely opposite modulation (delayed/counter-phase). Overall this study shows that body displacements "remap" our metric of time, affecting not only motor output but also sensory input.

Uncovering the genetic basis of adaptive change: on the intersection of landscape genomics and theoretical population genetics

A workshop recently held at the Ecole Polytechnique Federale de Lausanne (EPFL, Switzerland) was dedicated to understanding the genetic basis of adaptive change, taking stock of the different approaches developed in theoretical population genetics and landscape genomics and bringing together knowledge accumulated in both research fields. Indeed, an important challenge in theoretical population genetics is to incorporate effects of demographic history and population structure. But important design problems (e.g. focus on populations as units, focus on hard selective sweeps, no hypothesis-based framework in the design of the statistical tests) reduce their capability of detecting adaptive genetic variation. In parallel, landscape genomics offers a solution to several of these problems and provides a number of advantages (e.g. fast computation, landscape heterogeneity integration). But the approach makes several implicit assumptions that should be carefully considered (e.g. selection has had enough time to create a functional relationship between the allele distribution and the environmental variable, or this functional relationship is assumed to be constant). To address the respective strengths and weaknesses mentioned above, the workshop brought together a panel of experts from both disciplines to present their work and discuss the relevance of combining these approaches, possibly resulting in a joint software solution in the future.

Tissue homing and persistence of defined antigen-specific CD8+ tumor-reactive T-cell clones in long-term melanoma survivors.

Tumor antigen-specific cytotoxic T cells (CTLs) play a major role in the adaptive immune response to cancers. This CTL response is often insufficient because of functional impairment, tumor escape mechanisms, or inhibitory tumor microenvironment. However, little is known about the fate of given tumor-specific CTL clones in cancer patients. Studies in patients with favorable outcomes may be very informative. In this longitudinal study, we tracked, quantified, and characterized functionally defined antigen-specific T-cell clones ex vivo, in peripheral blood and at tumor sites, in two long-term melanoma survivors. MAGE-A10-specific CD8+ T-cell clones with high avidity to antigenic peptide and tumor lytic capabilities persisted in peripheral blood over more than 10 years, with quantitative variations correlating with the clinical course. These clones were also found in emerging metastases, and, in one patient, circulating clonal T cells displayed a fully differentiated effector phenotype at the time of relapse. Longevity, tumor homing, differentiation phenotype, and quantitative adaptation to the disease phases suggest the contribution of the tracked tumor-reactive clones in the tumor control of these long-term metastatic survivor patients. Focusing research on patients with favorable outcomes may help to identify parameters that are crucial for an efficient antitumor response and to optimize cancer immunotherapy.

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae).

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.

Local adaptation maintains clinal variation in melanin-based coloration of European barn owls (Tyto alba).

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin-based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish-brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.

Cross-frequency coupling in EEG revealed by adaptive oscillation tracking.

Introduction: Neuronal oscillations have been the focus of increasing interest in the neuroscientific community, in part because they have been considered as a possible integrating mechanism through which internal states can influence stimulus processing in a top-down way (Engel et al., 2001). Moreover, increasing evidence indicates that oscillations in different frequency bands interact with one other through coupling mechanisms (Jensen and Colgin, 2007). The existence and the importance of these cross-frequency couplings during various tasks have been verified by recent studies (Canolty et al., 2006; Lakatos et al., 2007). In this study, we measure the strength and directionality of two types of couplings - phase-amplitude couplings and phase-phase couplings - between various bands in EEG data recorded during an illusory contour experiment that were identified using a recently-proposed adaptive frequency tracking algorithm (Van Zaen et al., 2010). Methods: The data used in this study have been taken from a previously published study examining the spatiotemporal mechanisms of illusory contour processing (Murray et al., 2002). The EEG in the present study were from a subset of nine subjects. Each stimulus was composed of 'pac-man' inducers presented in two orientations: IC, when an illusory contour was present, and NC, when no contour could be detected. The signals recorded by the electrodes P2, P4, P6, PO4 and PO6 were averaged, and filtered into the following bands: 4-8Hz, 8-12Hz, 15-25Hz, 35-45Hz, 45-55Hz, 55-65Hz and 65-75Hz. An adaptive frequency tracking algorithm (Van Zaen et al., 2010) was then applied in each band in order to extract the main oscillation and estimate its frequency. This additional step ensures that clean phase information is obtained when taking the Hilbert transform. The frequency estimated by the tracker was averaged over sliding windows and then used to compare the two conditions. Two types of cross-frequency couplings were considered: phase-amplitude couplings and phase-phase couplings. Both types were measured with the phase locking value (PLV, Lachaux et al., 1999) over sliding windows. The phase-amplitude couplings were computed with the phase of the low frequency oscillation and the phase of the amplitude of the high frequency one. Different coupling coefficients were used when measuring phase-phase couplings in order to estimate different m:n synchronizations (4:3, 3:2, 2:1, 3:1, 4:1, 5:1, 6:1, 7:1, 8:1 and 9:1) and to take into account the frequency differences across bands. Moreover, the direction of coupling was estimated with a directionality index (Bahraminasab et al., 2008). Finally, the two conditions IC and NC were compared with ANOVAs with 'subject' as a random effect and 'condition' as a fixed effect. Before computing the statistical tests, the PLV values were transformed into approximately normal variables (Penny et al., 2008). Results: When comparing the mean estimated frequency across conditions, a significant difference was found only in the 4-8Hz band, such that the frequency within this band was significantly higher for IC than NC stimuli starting at ~250ms post-stimulus onset (Fig. 1; solid line shows IC and dashed line NC). Significant differences in phase-amplitude couplings were obtained only when the 4-8 Hz band was taken as the low frequency band. Moreover, in all significant situations, the coupling strength is higher for the NC than IC condition. An example of significant difference between conditions is shown in Fig. 2 for the phase-amplitude coupling between the 4-8Hz and 55-65Hz bands (p-value in top panel and mean PLV values in the bottom panel). A decrease in coupling strength was observed shortly after stimulus onset for both conditions and was greater for the condition IC. This phenomenon was observed with all other frequency bands. The results obtained for the phase-phase couplings were more complex. As for the phase-amplitude couplings, all significant differences were obtained when the 4-8Hz band was considered as the low frequency band. The stimulus condition exhibiting the higher coupling strength depended on the ratio of the coupling coefficients. When this ratio was small, the IC condition exhibited the higher phase-phase coupling strength. When this ratio was large, the NC condition exhibited the higher coupling strength. Fig. 3 shows the phase-phase couplings between the 4-8Hz and 35-45Hz bands for the coupling coefficient 6:1, and the coupling strength was significantly higher for the IC than NC condition. By contrast, for the coupling coefficient 9:1 the NC condition gave the higher coupling strength (Fig. 4). Control analyses verified that it is not a consequence of the frequency difference between the two conditions in the 4-8Hz band. The directionality measures indicated a transfer of information from the low frequency components towards the high frequency ones. Conclusions: Adaptive tracking is a feasible method for EEG analyses, revealing information both about stimulus-related differences and coupling patterns across frequencies. Theta oscillations play a central role in illusory shape processing and more generally in visual processing. The presence vs. absence of illusory shapes was paralleled by faster theta oscillations. Phase-amplitude couplings were decreased more for IC than NC and might be due to a resetting mechanism. The complex patterns in phase-phase coupling between theta and beta/gamma suggest that the contribution of these oscillations to visual binding and stimulus processing are not as straightforward as conventionally held. Causality analyses further suggest that theta oscillations drive beta/gamma oscillations (see also Schroeder and Lakatos, 2009). The present findings highlight the need for applying more sophisticated signal analyses in order to establish a fuller understanding of the functional role of neural oscillations.

Effects of computing parameters and measurement locations on the estimation of 3D NPS in non-stationary MDCT images.

The goal of this study was to investigate the impact of computing parameters and the location of volumes of interest (VOI) on the calculation of 3D noise power spectrum (NPS) in order to determine an optimal set of computing parameters and propose a robust method for evaluating the noise properties of imaging systems. Noise stationarity in noise volumes acquired with a water phantom on a 128-MDCT and a 320-MDCT scanner were analyzed in the spatial domain in order to define locally stationary VOIs. The influence of the computing parameters in the 3D NPS measurement: the sampling distances bx,y,z and the VOI lengths Lx,y,z, the number of VOIs NVOI and the structured noise were investigated to minimize measurement errors. The effect of the VOI locations on the NPS was also investigated. Results showed that the noise (standard deviation) varies more in the r-direction (phantom radius) than z-direction plane. A 25 × 25 × 40 mm(3) VOI associated with DFOV = 200 mm (Lx,y,z = 64, bx,y = 0.391 mm with 512 × 512 matrix) and a first-order detrending method to reduce structured noise led to an accurate NPS estimation. NPS estimated from off centered small VOIs had a directional dependency contrary to NPS obtained from large VOIs located in the center of the volume or from small VOIs located on a concentric circle. This showed that the VOI size and location play a major role in the determination of NPS when images are not stationary. This study emphasizes the need for consistent measurement methods to assess and compare image quality in CT.

MHC Class II Transactivator Is an In Vivo Regulator of Osteoclast Differentiation and Bone Homeostasis Co-opted From Adaptive Immunity.

The molecular networks controlling bone homeostasis are not fully understood. The common evolution of bone and adaptive immunity encourages the investigation of shared regulatory circuits. MHC Class II Transactivator (CIITA) is a master transcriptional co-activator believed to be exclusively dedicated for antigen presentation. CIITA is expressed in osteoclast precursors, and its expression is accentuated in osteoporotic mice. We thus asked whether CIITA plays a role in bone biology. To this aim, we fully characterized the bone phenotype of two mouse models of CIITA overexpression, respectively systemic and restricted to the monocyte-osteoclast lineage. Both CIITA-overexpressing mouse models revealed severe spontaneous osteoporosis, as assessed by micro-computed tomography and histomorphometry, associated with increased osteoclast numbers and enhanced in vivo bone resorption, whereas osteoblast numbers and in vivo bone-forming activity were unaffected. To understand the underlying cellular and molecular bases, we investigated ex vivo the differentiation of mutant bone marrow monocytes into osteoclasts and immune effectors, as well as osteoclastogenic signaling pathways. CIITA-overexpressing monocytes differentiated normally into effector macrophages or dendritic cells but showed enhanced osteoclastogenesis, whereas CIITA ablation suppressed osteoclast differentiation. Increased c-fms and receptor activator of NF-κB (RANK) signaling underlay enhanced osteoclast differentiation from CIITA-overexpressing precursors. Moreover, by extending selected phenotypic and cellular analyses to additional genetic mouse models, namely MHC Class II deficient mice and a transgenic mouse line lacking a specific CIITA promoter and re-expressing CIITA in the thymus, we excluded MHC Class II expression and T cells from contributing to the observed skeletal phenotype. Altogether, our study provides compelling genetic evidence that CIITA, the molecular switch of antigen presentation, plays a novel, unexpected function in skeletal homeostasis, independent of MHC Class II expression and T cells, by exerting a selective and intrinsic control of osteoclast differentiation and bone resorption in vivo. © 2014 American Society for Bone and Mineral Research.

Clypeal patterning in the paper wasp Polistes dominulus: no evidence of adaptive value in the wild

Status signals function in a number of species to communicate competitive ability to conspecific rivals during competition for resources. In the paper wasp Polistes dominulus, variable black clypeal patterns are thought to be important in mediating competition among females. Results of previous behavioral experiments in the lab indicate that P dominulus clypeal patterns provide information about an individual's competitive ability to rivals during agonistic interactions. To date, however, there has been no detailed examination of the adaptive value of clypeal patterns in the wild. To address this, we looked for correlations between clypeal patterning and various fitness measures, including reproductive success, hierarchical rank, and survival, in a large, free-living population of P. dominulus in southern Spain. Reproductive success over the nesting season was not correlated with clypeal patterning. Furthermore, there was no relationship between a female's clypeal patterning and the rank she achieved within the hierarchy or her survival during nest founding. Overall, we found no evidence that P dominulus clypeal patterns are related to competitive ability or other aspects of quality in our population. This result is consistent with geographical variation in the adaptive value of clypeal patterns between P. dominulus populations; however, data on the relationship between patterning and fitness from other populations are required to test this hypothesis.

Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

The kinetic parameters and energy cost of the Hsp70 chaperone as a polypeptide unfoldase.

Hsp70-Hsp40-NEF and possibly Hsp100 are the only known molecular chaperones that can use the energy of ATP to convert stably pre-aggregated polypeptides into natively refolded proteins. However, the kinetic parameters and ATP costs have remained elusive because refolding reactions have only been successful with a molar excess of chaperones over their polypeptide substrates. Here we describe a stable, misfolded luciferase species that can be efficiently renatured by substoichiometric amounts of bacterial Hsp70-Hsp40-NEF. The reactivation rates increased with substrate concentration and followed saturation kinetics, thus allowing the determination of apparent V(max)' and K(m)' values for a chaperone-mediated renaturation reaction for the first time. Under the in vitro conditions used, one Hsp70 molecule consumed five ATPs to effectively unfold a single misfolded protein into an intermediate that, upon chaperone dissociation, spontaneously refolded to the native state, a process with an ATP cost a thousand times lower than expected for protein degradation and resynthesis.