Estimating the confidence level of white matter connections obtained with MRI tractography.

BACKGROUND: Since the emergence of diffusion tensor imaging, a lot of work has been done to better understand the properties of diffusion MRI tractography. However, the validation of the reconstructed fiber connections remains problematic in many respects. For example, it is difficult to assess whether a connection is the result of the diffusion coherence contrast itself or the simple result of other uncontrolled parameters like for example: noise, brain geometry and algorithmic characteristics. METHODOLOGY/PRINCIPAL FINDINGS: In this work, we propose a method to estimate the respective contributions of diffusion coherence versus other effects to a tractography result by comparing data sets with and without diffusion coherence contrast. We use this methodology to assign a confidence level to every gray matter to gray matter connection and add this new information directly in the connectivity matrix. CONCLUSIONS/SIGNIFICANCE: Our results demonstrate that whereas we can have a strong confidence in mid- and long-range connections obtained by a tractography experiment, it is difficult to distinguish between short connections traced due to diffusion coherence contrast from those produced by chance due to the other uncontrolled factors of the tractography methodology.

Selectome update: quality control and computational improvements to a database of positive selection.

Selectome (http://selectome.unil.ch/) is a database of positive selection, based on a branch-site likelihood test. This model estimates the number of nonsynonymous substitutions (dN) and synonymous substitutions (dS) to evaluate the variation in selective pressure (dN/dS ratio) over branches and over sites. Since the original release of Selectome, we have benchmarked and implemented a thorough quality control procedure on multiple sequence alignments, aiming to provide minimum false-positive results. We have also improved the computational efficiency of the branch-site test implementation, allowing larger data sets and more frequent updates. Release 6 of Selectome includes all gene trees from Ensembl for Primates and Glires, as well as a large set of vertebrate gene trees. A total of 6810 gene trees have some evidence of positive selection. Finally, the web interface has been improved to be more responsive and to facilitate searches and browsing.

Phenological patterns in Mediterranean pastures and scrubs of Catalonia

On the basis of two parallel data sets refering to phenological events in open scrubs and pastures at two sites (Balaguer and Vic), the authors present a comprehensive report of the phenology of these Mediterranean communities. Four main phenophases (vegetative growth, flowering, fruiting and resting) were recorded monthly in 11 communities over 15 months. The results allow comparisons to be drawn between localities and communities, both at community and species levels, and to consider the effects of contemporary climatic data. This yielded useful information on general trends and on the particular responses of each community type to their corresponding constraints.

Comparing 'visual' effect size indices for single-case designs

Effect size indices are indispensable for carrying out meta-analyses and can also be seen as an alternative for making decisions about the effectiveness of a treatment in an individual applied study. The desirable features of the procedures for quantifying the magnitude of intervention effect include educational/clinical meaningfulness, calculus easiness, insensitivity to autocorrelation, low false alarm and low miss rates. Three effect size indices related to visual analysis are compared according to the aforementioned criteria. The comparison is made by means of data sets with known parameters: degree of serial dependence, presence or absence of general trend, changes in level and/or in slope. The percent of nonoverlapping data showed the highest discrimination between data sets with and without intervention effect. In cases when autocorrelation or trend is present, the percentage of data points exceeding the median may be a better option to quantify the effectiveness of a psychological treatment.

Temporal quantification of MAPK induced expression in single yeast cells.

The quantification of gene expression at the single cell level uncovers novel regulatory mechanisms obscured in measurements performed at the population level. Two methods based on microscopy and flow cytometry are presented to demonstrate how such data can be acquired. The expression of a fluorescent reporter induced upon activation of the high osmolarity glycerol MAPK pathway in yeast is used as an example. The specific advantages of each method are highlighted. Flow cytometry measures a large number of cells (10,000) and provides a direct measure of the dynamics of protein expression independent of the slow maturation kinetics of the fluorescent protein. Imaging of living cells by microscopy is by contrast limited to the measurement of the matured form of the reporter in fewer cells. However, the data sets generated by this technique can be extremely rich thanks to the combinations of multiple reporters and to the spatial and temporal information obtained from individual cells. The combination of these two measurement methods can deliver new insights on the regulation of protein expression by signaling pathways.

Heterogeneity and event dependence in the analysis of sickness absence

Sickness absence (SA) is an important social, economic and public health issue. Identifying and understanding the determinants, whether biological, regulatory or, health services-related, of variability in SA duration is essential for better management of SA. The conditional frailty model (CFM) is useful when repeated SA events occur within the same individual, as it allows simultaneous analysis of event dependence and heterogeneity due to unknown, unmeasured, or unmeasurable factors. However, its use may encounter computational limitations when applied to very large data sets, as may frequently occur in the analysis of SA duration. To overcome the computational issue, we propose a Poisson-based conditional frailty model (CFPM) for repeated SA events that accounts for both event dependence and heterogeneity. To demonstrate the usefulness of the model proposed in the SA duration context, we used data from all non-work-related SA episodes that occurred in Catalonia (Spain) in 2007, initiated by either a diagnosis of neoplasm or mental and behavioral disorders. As expected, the CFPM results were very similar to those of the CFM for both diagnosis groups. The CPU time for the CFPM was substantially shorter than the CFM. The CFPM is an suitable alternative to the CFM in survival analysis with recurrent events,especially with large databases.

Rib fractures at postmortem computed tomography (PMCT) validated against the autopsy

To evaluate the sensitivity of postmortem computed tomography (PMCT) in rib fracture detection validated against autopsy. Fifty-one forensic cases underwent a postmortem CT prior to forensic autopsy. Two image readers (radiologist and forensic pathologist) assessed high resolution CT data sets for rib fractures. Correct recognition rates (CRR), sensitivity and specificity values were calculated over all observations as well as individually for every rib and region. Additionally, for partial rib fractures the sensitivity of autopsy was calculated vice versa. 3876 entries in each study protocol (autopsy, PMCT radiologist and PMCT forensic pathologist) were investigated. A total of 690 fractures (autopsy), 491 (PMCT and radiologist) and 559 (PMCT and forensic pathologist) were detected. The CRR was 0.85. Sensitivity and specificity of PMCT for rib fracture detection were 0.63 (0.58 radiologist, 0.68 forensic pathologist) and 0.97 (both readers 0.97), respectively. Low CRR and sensitivity values were obtained for antero-lateral fractures. Partial rib fractures were better detected by PMCT. PMCT has a rather low sensitivity for rib fracture detection when validated against autopsy and indicates that clinical CT may also demonstrate a reasonable number of false negatives. Partial rib fractures often remain undetected at autopsy.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

A pipeline approach with spatial information for segmenting multiple sclerosis lesions on brain magnetic resonance imaging

Background: Conventional magnetic resonance imaging (MRI) techniques are highly sensitive to detect multiple sclerosis (MS) plaques, enabling a quantitative assessment of inflammatory activity and lesion load. In quantitative analyses of focal lesions, manual or semi-automated segmentations have been widely used to compute the total number of lesions and the total lesion volume. These techniques, however, are both challenging and time-consuming, being also prone to intra-observer and inter-observer variability.Aim: To develop an automated approach to segment brain tissues and MS lesions from brain MRI images. The goal is to reduce the user interaction and to provide an objective tool that eliminates the inter- and intra-observer variability.Methods: Based on the recent methods developed by Souplet et al. and de Boer et al., we propose a novel pipeline which includes the following steps: bias correction, skull stripping, atlas registration, tissue classification, and lesion segmentation. After the initial pre-processing steps, a MRI scan is automatically segmented into 4 classes: white matter (WM), grey matter (GM), cerebrospinal fluid (CSF) and partial volume. An expectation maximisation method which fits a multivariate Gaussian mixture model to T1-w, T2-w and PD-w images is used for this purpose. Based on the obtained tissue masks and using the estimated GM mean and variance, we apply an intensity threshold to the FLAIR image, which provides the lesion segmentation. With the aim of improving this initial result, spatial information coming from the neighbouring tissue labels is used to refine the final lesion segmentation.Results:The experimental evaluation was performed using real data sets of 1.5T and the corresponding ground truth annotations provided by expert radiologists. The following values were obtained: 64% of true positive (TP) fraction, 80% of false positive (FP) fraction, and an average surface distance of 7.89 mm. The results of our approach were quantitatively compared to our implementations of the works of Souplet et al. and de Boer et al., obtaining higher TP and lower FP values.Conclusion: Promising MS lesion segmentation results have been obtained in terms of TP. However, the high number of FP which is still a well-known problem of all the automated MS lesion segmentation approaches has to be improved in order to use them for the standard clinical practice. Our future work will focus on tackling this issue.

Digital forensic osteology--possibilities in cooperation with the Virtopsy project.

The present study was carried out to check whether classic osteometric parameters can be determined from the 3D reconstructions of MSCT (multislice computed tomography) scans acquired in the context of the Virtopsy project. To this end, four isolated and macerated skulls were examined by six examiners. First the skulls were conventionally (manually) measured using 32 internationally accepted linear measurements. Then the skulls were scanned by the use of MSCT with slice thicknesses of 1.25 mm and 0.63 mm, and the 33 measurements were virtually determined on the digital 3D reconstructions of the skulls. The results of the traditional and the digital measurements were compared for each examiner to figure out variations. Furthermore, several parameters were measured on the cranium and postcranium during an autopsy and compared to the values that had been measured on a 3D reconstruction from a previously acquired postmortem MSCT scan. The results indicate that equivalent osteometric values can be obtained from digital 3D reconstructions from MSCT scans using a slice thickness of 1.25 mm, and from conventional manual examinations. The measurements taken from a corpse during an autopsy could also be validated with the methods used for the digital 3D reconstructions in the context of the Virtopsy project. Future aims are the assessment and biostatistical evaluation in respect to sex, age and stature of all data sets stored in the Virtopsy project so far, as well as of future data sets. Furthermore, a definition of new parameters, only measurable with the aid of MSCT data would be conceivable.

FM127: Systemic metal ion levels in patients with modular neck total hip arthroplasties: a prospective cohort study

There is increasing evidence that modular neck stems are prone to corrosion-related complications. Recent studies showed elevated metal ions levels and occasional pseudotumor formation in patients with such implants. The purpose of this study was to compare systemic metal-ion levels in patients after primary THA with modular neck stems to those of patients after non-modular implants. To our knowledge, this is the first cohort study including a control group, THA without CoCr heads and dry-assembled neck-stem connections. Methods: 50 patients after THA at a minimum follow-up of 1 year have been selected for the study. Patients with multiple prosthesis or other implants have been deselected. All received a cementless SPS stem from Symbios (Ti6Al4V). 40 patients have the modular neck (CoCr) version and 10 a monobloc version. All bearings were either ceramic-ceramic or ceramic-polyethylene to minimize other sources of CoCr ion release. In the modular group, the neck was chosen pre-operatively based on a 3D planning, allowing for a dry assembly of the stem and neck on the back table before implantation. A plasma system coupled to mass spectrometry was used for a complete elementary quantification in blood and serum separately. Clinical outcome was measured using the Oxford Hip Score. Results : Complete data sets of 29 patients (24 in the modular neck-group (10male, mean age 63y, 35-84y) and 5 in the monobloc-group (3 male, 69 y, 51-83y) are available to date. Mean Co blood levels were .95 ug/L (.14-12.4) in the modular group vs .27 ug/L (.10-.73) in the monobloc group (p=.2). Respective values for Cr were significantly higher in the modular group (.99 g/L; range .75-1.21) compared to those in the monobloc group (.74 g/L ;.62-.86; p=.001). No significant difference was found when comparing serum levels. 5/24 patients had Co levels above 1 g/L (12/24 for Cr), which is by some considered as a relevant elevation. The maximum Co level was measured in an asymptomatic patient. The Oxford Hip Scores were similar in both groups. Conclusion: Cr levels were significantly elevated in the modular neck group compared to those in the monobloc group. 1/24 patients with a modular prosthesis exhibited Co levels, which are beyond the threshold accepted even for metal-on-metal bearing couples. These results have contributed to our decision to abandon the use of modular neck stems. Routine follow-up including annual measurements of systemic CoCr concentrations should be considered.

Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

A comparison of alpha and gamma spectrometry for environmental natural radioactivity surveys.

An alpha-spectrometry, using automated borate fusion and sequential extraction and exchange chromatography, was used to determine the uranium and thorium based on environmental radioactivity of 20 soil samples. The same set of the samples was analysed using gamma-spectrometry with an HPGe detector. The two data sets were checked for coherence using Z-score and chi2 statistical tests. We show that gamma-spectrometry is a valid alternative to time-consuming alpha-spectrometry for the determination of natural uranium and thorium activity in soil (activity range: 12.5-58.2 Bq/kg). The measured activities were compared with the theoretical activities to ensure secular equilibrium in the 238U and 232Th series. For 226Ra, a special study was made on deconvolution of the 186 keV multiplet with the Levenberg-Marquardt algorithm. Finally, the combined use of Z-score and chi2-tests was found to be a powerful tool for comparing the results obtained with two different methods.

Pulse oximetry and transcutaneous oxygen tension for detection of hypoxemia in critically ill infants and children.

We tested the performance of transcutaneous oxygen monitoring (TcPO2) and pulse oximetry (tcSaO2) in detecting hypoxia in critically ill neonatal and pediatric patients. In 54 patients (178 data sets) with a mean age of 2.4 years (range 1 to 19 years), arterial saturation (SaO2) ranged from 9.5 to 100%, and arterial oxygen tension (PaO2) from 16.4 to 128 mmHg. Linear correlation analysis of pulse oximetry vs measured SaO2 revealed an r value of 0.95 (p less than 0.001) with an equation of y = 21.1 + 0.749x, while PaO2 vs tcPO2 showed a correlation coefficient of r = 0.95 (p less than 0.001) with an equation of y = -1.04 + 0.876x. The mean difference between measured SaO2 and tcSaO2 was -2.74 +/- 7.69% (range +14 to - 29%) and the mean difference between PaO2 and tcPO2 was +7.43 +/- 8.57 mmHg (range -14 to +49 mmHg). Pulse oximetry was reliable at values above 65%, but was inaccurate and overestimated the arterial SaO2 at lower values. TcPO2 tended to underestimate the arterial value with increasing PaO2. Pulse oximetry had the best sensitivity to specificity ratio for hypoxia between 65 and 90% SaO2; for tcPO2 the best results were obtained between 35 and 55 mmHg PaO2.

Software for interpolation of vegetative growth of yerba mate plants in 3D

The objective of this work was to build mock-ups of complete yerba mate plants in several stages of development, using the InterpolMate software, and to compute photosynthesis on the interpolated structure. The mock-ups of yerba-mate were first built in the VPlants software for three growth stages. Male and female plants grown in two contrasting environments (monoculture and forest understory) were considered. To model the dynamic 3D architecture of yerba-mate plants during the biennial growth interval between two subsequent prunings, data sets of branch development collected in 38 dates were used. The estimated values obtained from the mock-ups, including leaf photosynthesis and sexual dimorphism, are very close to those observed in the field. However, this similarity was limited to reconstructions that included growth units from original data sets. The modeling of growth dynamics enables the estimation of photosynthesis for the entire yerba mate plant, which is not easily measurable in the field. The InterpolMate software is efficient for building yerba mate mock-ups.