Probabilistic base calling of Solexa sequencing data.

BACKGROUND: Solexa/Illumina short-read ultra-high throughput DNA sequencing technology produces millions of short tags (up to 36 bases) by parallel sequencing-by-synthesis of DNA colonies. The processing and statistical analysis of such high-throughput data poses new challenges; currently a fair proportion of the tags are routinely discarded due to an inability to match them to a reference sequence, thereby reducing the effective throughput of the technology. RESULTS: We propose a novel base calling algorithm using model-based clustering and probability theory to identify ambiguous bases and code them with IUPAC symbols. We also select optimal sub-tags using a score based on information content to remove uncertain bases towards the ends of the reads. CONCLUSION: We show that the method improves genome coverage and number of usable tags as compared with Solexa's data processing pipeline by an average of 15%. An R package is provided which allows fast and accurate base calling of Solexa's fluorescence intensity files and the production of informative diagnostic plots.

Early Variscan I-type pluton in the pre-Alpine basement of the Western Alps: The ca. 360 Ma Cogne diorite (NW-Italy)

Located at the internal border of the Grand-Saint-Bernard Zone, the diorite and its aureole lie on top of intensively studied Alpine eclogitic units but this pluton, poorly studied yet, has kept locally almost undeformed. The pluton intruded, at similar to 360 Ma, country-rocks mostly composed of dark shales with Na2O > K2O and minor mafic intercalations of tholeiitic basalt affinity. This association is characteristic of the Vanoise (France) basement series, where available age determinations suggest an Early Paleozoic age. Parts of the pluton, and of its hornfels aureole that is evidenced here for the first time, in the Punta Bioula section of Valsavaranche valley (NW-Italy), have been well-preserved from the Alpine deformation. Syn-emplacement hardening, dehydration-induced, probably prevented strain-enhanced Alpine recrystallization. Magmatic rock-types range continuously from subordinate mafic types at SiO2 similar to 48%, of hornblendite with cumulative or appinite affinities, to the main body of quartz diorite to quartz monzonite (SiO2 up to 62%). P-T estimates for the pluton emplacement, based on the abundance of garnet in the hornfelses, using also zircon and apatite saturation thermometry and Al-in-hornblende barometry, suggest T similar to 800-950 degrees C and minimum P in the 0.2-0.5 GPa range, with records of higher pressure conditions (up to 1-2 GPa?) in hornblendite phlogopite-cored amphibole. The high-K, Na > K, calcalkaline geochemistry is in line with a destructive plate-margin setting. Based on major element data and radiogenic isotope signature (epsilon Nd-360 Ma from -1.2 to + 0.9, Sr-87/Sr-86(360 MA) from 0.7054 to 0.7063), the parental magmas are interpreted in terms of deep-seated metabasaltic partial melts with limited contamination from shallower sources, the low radiogenic Nd-content excluding a major contribution from Vanoise tholeiites. There is no other preserved evidence for Variscan magmatism of similar age and composition in the Western Alps, but probable analogs are known in the western and northern parts of French Massif Central. Regarding the Alpine tectonics, not only the age of the pluton and its host-rocks (instead of the Permo-Carboniferous age previously believed), but also its upper mylonitic contact, suggest revisions of the Alpine nappe model. The Cogne diorite allegedly constituted the axial part of the E-verging ``pli en retour [backfold] du Valsavaranche'', a cornerstone of popular Alpine structural models: in fact, the alleged fold limbs, as attested here by field and geochemical data, do not belong to the same unit, and the backfold hypothesis is unfounded. (C) 2012 Elsevier B.V. All rights reserved.

Probabilistic evidential assessment of gunshot residue particle evidence (Part II) : Bayesian parameter estimation for experimental count data.

Part I of this series of articles focused on the construction of graphical probabilistic inference procedures, at various levels of detail, for assessing the evidential value of gunshot residue (GSR) particle evidence. The proposed models - in the form of Bayesian networks - address the issues of background presence of GSR particles, analytical performance (i.e., the efficiency of evidence searching and analysis procedures) and contamination. The use and practical implementation of Bayesian networks for case pre-assessment is also discussed. This paper, Part II, concentrates on Bayesian parameter estimation. This topic complements Part I in that it offers means for producing estimates useable for the numerical specification of the proposed probabilistic graphical models. Bayesian estimation procedures are given a primary focus of attention because they allow the scientist to combine (his/her) prior knowledge about the problem of interest with newly acquired experimental data. The present paper also considers further topics such as the sensitivity of the likelihood ratio due to uncertainty in parameters and the study of likelihood ratio values obtained for members of particular populations (e.g., individuals with or without exposure to GSR).

Auditory motion affects visual biological motion processing.

The processing of biological motion is a critical, everyday task performed with remarkable efficiency by human sensory systems. Interest in this ability has focused to a large extent on biological motion processing in the visual modality (see, for example, Cutting, J. E., Moore, C., & Morrison, R. (1988). Masking the motions of human gait. Perception and Psychophysics, 44(4), 339-347). In naturalistic settings, however, it is often the case that biological motion is defined by input to more than one sensory modality. For this reason, here in a series of experiments we investigate behavioural correlates of multisensory, in particular audiovisual, integration in the processing of biological motion cues. More specifically, using a new psychophysical paradigm we investigate the effect of suprathreshold auditory motion on perceptions of visually defined biological motion. Unlike data from previous studies investigating audiovisual integration in linear motion processing [Meyer, G. F. & Wuerger, S. M. (2001). Cross-modal integration of auditory and visual motion signals. Neuroreport, 12(11), 2557-2560; Wuerger, S. M., Hofbauer, M., & Meyer, G. F. (2003). The integration of auditory and motion signals at threshold. Perception and Psychophysics, 65(8), 1188-1196; Alais, D. & Burr, D. (2004). No direction-specific bimodal facilitation for audiovisual motion detection. Cognitive Brain Research, 19, 185-194], we report the existence of direction-selective effects: relative to control (stationary) auditory conditions, auditory motion in the same direction as the visually defined biological motion target increased its detectability, whereas auditory motion in the opposite direction had the inverse effect. Our data suggest these effects do not arise through general shifts in visuo-spatial attention, but instead are a consequence of motion-sensitive, direction-tuned integration mechanisms that are, if not unique to biological visual motion, at least not common to all types of visual motion. Based on these data and evidence from neurophysiological and neuroimaging studies we discuss the neural mechanisms likely to underlie this effect.

Spontaneous pre-stimulus fluctuations in the activity of right fronto-parietal areas influence inhibitory control performance.

Inhibitory control refers to the ability to suppress planned or ongoing cognitive or motor processes. Electrophysiological indices of inhibitory control failure have been found to manifest even before the presentation of the stimuli triggering the inhibition, suggesting that pre-stimulus brain-states modulate inhibition performance. However, previous electrophysiological investigations on the state-dependency of inhibitory control were based on averaged event-related potentials (ERPs), a method eliminating the variability in the ongoing brain activity not time-locked to the event of interest. These studies thus left unresolved whether spontaneous variations in the brain-state immediately preceding unpredictable inhibition-triggering stimuli also influence inhibitory control performance. To address this question, we applied single-trial EEG topographic analyses on the time interval immediately preceding NoGo stimuli in conditions where the responses to NoGo trials were correctly inhibited [correct rejection (CR)] vs. committed [false alarms (FAs)] during an auditory spatial Go/NoGo task. We found a specific configuration of the EEG voltage field manifesting more frequently before correctly inhibited responses to NoGo stimuli than before FAs. There was no evidence for an EEG topography occurring more frequently before FAs than before CR. The visualization of distributed electrical source estimations of the EEG topography preceding successful response inhibition suggested that it resulted from the activity of a right fronto-parietal brain network. Our results suggest that the fluctuations in the ongoing brain activity immediately preceding stimulus presentation contribute to the behavioral outcomes during an inhibitory control task. Our results further suggest that the state-dependency of sensory-cognitive processing might not only concern perceptual processes, but also high-order, top-down inhibitory control mechanisms.

Mapping of Environmental Data Using Kernel-Based Methods

Recently, kernel-based Machine Learning methods have gained great popularity in many data analysis and data mining fields: pattern recognition, biocomputing, speech and vision, engineering, remote sensing etc. The paper describes the use of kernel methods to approach the processing of large datasets from environmental monitoring networks. Several typical problems of the environmental sciences and their solutions provided by kernel-based methods are considered: classification of categorical data (soil type classification), mapping of environmental and pollution continuous information (pollution of soil by radionuclides), mapping with auxiliary information (climatic data from Aral Sea region). The promising developments, such as automatic emergency hot spot detection and monitoring network optimization are discussed as well.

High frequency of functionally active Melan-a-specific T cells in a patient with progressive immunoproteasome-deficient melanoma.

Tumor-reactive T cells play an important role in cancer immunosurveillance. Applying the multimer technology, we report here an unexpected high frequency of Melan-A-specific CTLs in a melanoma patient with progressive lymph node metastases, consisting of 18 and 12.8% of total peripheral blood and tumor-infiltrating CD8+ T cells, respectively. Melan-A-specific CTLs revealed a high cytolytic activity against allogeneic Melan-A-expressing target cells but failed to kill the autologous tumor cells. Loading of the tumor cells with Melan-A peptide reversed the resistance to killing, suggesting impaired function of the MHC class I antigen processing and presentation pathway. Mutations of the coding region of the HLA-A2 binding Melan-A26-35 peptide or down-regulation of the MHC class I heavy chain, the antigenic peptide TAP, and tapasin could be excluded. However, PCR and immunohistochemical analysis revealed a deficiency of the immunoproteasomes low molecular weight protein 2 and low molecular weight protein 7 in the primary tumor cells, which affects the quantity and quality of generated T-cell epitopes and might explain the resistance to killing. This is supported by our data, demonstrating that the resistance to killing can be partially reversed by pre-exposure of the tumor cells to IFN-gamma, which is known to induce the immunoproteasomes. Overall, this is the first report of an extremely high frequency of tumor-specific CTLs that exhibit competent T-cell-effector functions but fail to lyse the autologous tumor cells. Immunotherapeutic approaches should not only focus on the induction of a robust antitumor immune response, but should also have to target tumor immune escape mechanisms.

TCRgamma silencing during alphabeta T cell development depends upon pre-TCR-induced proliferation.

During thymus development, immature T cells become committed to two distinct lineages based upon expression of alphabeta or gammadelta TCR. In the alphabeta lineage, developing thymocytes progressively extinguish transcription of the TCRgamma genes by a poorly understood process known as gamma silencing. We show that alphabeta lineage thymocytes in mice lacking a functional pre-TCR undergo limited proliferation and fail to silence TCRgamma genes during development. Stimulation of pre-TCR-deficient immature thymocytes with anti-CD3 Abs does not directly down-regulate TCRgamma transcription but restores TCRgamma silencing following proliferation. Collectively our data reveal an important role for pre-TCR induced proliferation in activating the TCRgamma silencer in alphabeta lineage thymocytes, a process that may reinforce alphabeta or gammadelta lineage commitment.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Ly49D-mediated ITAM signaling in immature thymocytes impairs development by bypassing the pre-TCR checkpoint.

Activating and inhibitory NK receptors regulate the development and effector functions of NK cells via their ITAM and ITIM motifs, which recruit protein tyrosine kinases and phosphatases, respectively. In the T cell lineage, inhibitory Ly49 receptors are expressed by a subset of activated T cells and by CD1d-restricted NKT cells, but virtually no expression of activating Ly49 receptors is observed. Using mice transgenic for the activating receptor Ly49D and its associated ITAM signaling DAP12 chain, we show in this article that Ly49D-mediated ITAM signaling in immature thymocytes impairs development due to a block in maturation from the double negative (DN) to double positive (DP) stages. A large proportion of Ly49D/DAP12 transgenic thymocytes were able to bypass the pre-TCR checkpoint at the DN3 stage, leading to the appearance of unusual populations of DN4 and DP cells that lacked expression of intracellular (ic) TCRβ-chain. High levels of CD5 were expressed on ic TCRβ(-) DN and DP thymocytes from Ly49D/DAP12 transgenic mice, further suggesting that Ly49D-mediated ITAM signaling mimics physiological ITAM signaling via the pre-TCR. We also observed unusual ic TCRβ(-) single positive thymocytes with an immature CD24(high) phenotype that were not found in the periphery. Importantly, thymocyte development was completely rescued by expression of an Ly49A transgene in Ly49D/DAP12 transgenic mice, indicating that Ly49A-mediated ITIM signaling can fully counteract ITAM signaling via Ly49D/DAP12. Collectively, our data indicate that inappropriate ITAM signaling by activating NK receptors on immature thymocytes can subvert T cell development by bypassing the pre-TCR checkpoint.

EZ-IO(®) intraosseous device implementation in a pre-hospital emergency service: A prospective study and review of the literature.

INTRODUCTION: Intraosseous access is increasingly recognised as an effective alternative vascular access to peripheral venous access. We aimed to prospectively study the patients receiving prehospital intraosseous access with the EZ-IO(®), and to compare our results with those of the available literature. METHODS: Every patient who required an intraosseous access with the EZ-IO from January 1st, 2009 to December 31st, 2011 was included. The main data collected were: age, sex, indication for intraosseous access, localisation of insertion, success rate, drugs and fluids administered, and complications. All published studies concerning the EZ-IO device were systematically searched and reviewed for comparison. RESULTS: Fifty-eight patients representing 60 EZ-IO procedures were included. Mean age was 47 years (range 0.5-91), and the success rate was 90%. The main indications were cardiorespiratory arrest (74%), major trauma (12%), and shock (5%). The anterior tibia was the main route. The main drugs administered were adrenaline (epinephrine), atropine and amiodarone. No complications were reported. We identified 30 heterogeneous studies representing 1603 EZ-IO insertions. The patients' characteristics and success rate were similar to our study. Complications were reported in 13 cases (1.3%). CONCLUSION: The EZ-IO provides an effective way to achieve vascular access in the pre-hospital setting. Our results were similar to the cumulative results of all studies involving the use of the EZ-IO, and that can be used for comparison for further studies.

Visuo-spatial processing in a dynamic and a static working memory paradigm in schizophrenia.

Recent findings suggest that the visuo-spatial sketchpad (VSSP) may be divided into two sub-components processing dynamic or static visual information. This model may be useful to elucidate the confusion of data concerning the functioning of the VSSP in schizophrenia. The present study examined patients with schizophrenia and matched controls in a new working memory paradigm involving dynamic (the Ball Flight Task - BFT) or static (the Static Pattern Task - SPT) visual stimuli. In the BFT, the responses of the patients were apparently based on the retention of the last set of segments of the perceived trajectory, whereas control subjects relied on a more global strategy. We assume that the patients' performances are the result of a reduced capacity in chunking visual information since they relied mainly on the retention of the last set of segments. This assumption is confirmed by the poor performance of the patients in the static task (SPT), which requires a combination of stimulus components into object representations. We assume that the static/dynamic distinction may help us to understand the VSSP deficits in schizophrenia. This distinction also raises questions about the hypothesis that visuo-spatial working memory can simply be dissociated into visual and spatial sub-components.

Projecting the self to past and future in mental time

Introduction: Human experience takes place in the line of mental-time (MT) created through imagination of oneself in different time-points in past or future (self-projection in time). Here we manipulated self-projection in MT not only with respect to one's life-events but also with respect to one's faces from different past and future time-points. Methods: We here compared MTT with respect to one's facial images from different time points in past and future (study 1: MT-faces) as well as with respect to different past and future life events (study 2: MT-events). Participants were asked to make judgments about past and future face images and past and future events from three different time-points: the present (Now), eight years earlier (Past) or eight years later (Future). In addition, as a control task participants were asked to make recognition judgments with respect to faces and memory-related judgments with respect to events without changing their habitual self-location in time. Behavioral measures and functional magnetic resonance imaging (fMRI) activity after subtraction of recognition and memory related activities show both absolute MT and relative MT effects for faces and events, signifying a fundamental brain mechanism of MT, disentangled from episodic memory functions. Results: Behavioural and event-related fMRI activity showed three independent effects characterized by (1) similarity between past recollection and future imagination, (2) facilitation of judgments related to the future as compared to the past, and (3) facilitation of judgments related to time-points distant from the present. These effects were found with respect to faces and events suggesting that the brain mechanisms of MT are independent of whether actual life episodes have to be re-/pre-experienced and recruited a common cerebral network including the medial-temporal, precuneus, inferior-frontal, temporo-parietal, and insular cortices. Conclusions: These behavioural and neural data suggest that self-projection in time is a crucial aspect of MT, relying on neural structures encoding memory, mental imagery, and self. Furthermore our results emphasize the idea that mental temporal processing is more strongly directed to future prediction than to past recollection.

Automatic Mapping and Classification of Spatial Environmental Data

The paper deals with the development and application of the generic methodology for automatic processing (mapping and classification) of environmental data. General Regression Neural Network (GRNN) is considered in detail and is proposed as an efficient tool to solve the problem of spatial data mapping (regression). The Probabilistic Neural Network (PNN) is considered as an automatic tool for spatial classifications. The automatic tuning of isotropic and anisotropic GRNN/PNN models using cross-validation procedure is presented. Results are compared with the k-Nearest-Neighbours (k-NN) interpolation algorithm using independent validation data set. Real case studies are based on decision-oriented mapping and classification of radioactively contaminated territories.

A fast 3D approach for coronary MRA.

Two-dimensional (2D)-breath-hold coronary magnetic resonance angiography (MRA) has been shown to be a fast and reliable method to depict the proximal coronary arteries. Recent developments, however, allow for free-breathing navigator gated and navigator corrected three-dimensional (3D) coronary MRA. These 3D approaches have potential for improved signal-to-noise ratio (SNR) and allow for the acquisition of adjacent thin slices without the misregistration problems known from 2D approaches. Still, a major impediment of a 3D acquisition is the increased scan time. The purpose of this study was the implementation of a free-breathing navigator gated and corrected ultra-fast 3D coronary MRA technique, which allows for scan times of less than 5 minutes. Twelve healthy adult subjects were examined in the supine position using a navigator gated and corrected ECG triggered ultra-fast 3D interleaved gradient echo planar imaging sequence (TFE-EPI). A 3D slab, consisting of 20 slices with a reconstructed slice thickness of 1.5 mm, was acquired with free-breathing. The diastolic TFE-EPI acquisition block was preceded by a T2prep pre-pulse, a diaphragmatic navigator pulse, and a fat suppression pre-pulse. With a TR of 19 ms and an effective TE of 5.4 ms, the duration of the data acquisition window duration was 38 ms. The in-plane spatial resolution was 1.0-1.3 mm*1.5-1.9 mm. In all cases, the entire left main (LM) and extensive portions of the left anterior descending (LAD) and right coronary artery (RCA) could be visualized with an average scan time for the entire 3D-volume data set of 2:57 +/- 0:51 minutes. Average contiguous vessel length visualized was 53 +/- 11 mm (range: 42 to 75 mm) for the LAD and 84 +/- 14 mm (range: 62 to 112 mm) for the RCA. Contrast-to-noise between coronary blood and myocardium was 5.0 +/- 2.3 for the LM/LAD and 8.0 +/- 2.9 for the RCA, resulting in an excellent suppression of myocardium. We present a new approach for free-breathing 3D coronary MRA, which allows for scan times superior to corresponding 2D coronary MRA approaches, and which takes advantage of the enhanced SNR of 3D acquisitions and the post-processing benefits of thin adjacent slices. The robust image quality and the short average scanning time suggest that this approach may be useful for screening the major coronary arteries or identification of anomalous coronary arteries. J. Magn. Reson. Imaging 1999;10:821-825.