3D time series analysis of cell shape using Laplacian approaches

Background: Fundamental cellular processes such as cell movement, division or food uptake critically depend on cells being able to change shape. Fast acquisition of three-dimensional image time series has now become possible, but we lack efficient tools for analysing shape deformations in order to understand the real three-dimensional nature of shape changes. Results: We present a framework for 3D+time cell shape analysis. The main contribution is three-fold: First, we develop a fast, automatic random walker method for cell segmentation. Second, a novel topology fixing method is proposed to fix segmented binary volumes without spherical topology. Third, we show that algorithms used for each individual step of the analysis pipeline (cell segmentation, topology fixing, spherical parameterization, and shape representation) are closely related to the Laplacian operator. The framework is applied to the shape analysis of neutrophil cells. Conclusions: The method we propose for cell segmentation is faster than the traditional random walker method or the level set method, and performs better on 3D time-series of neutrophil cells, which are comparatively noisy as stacks have to be acquired fast enough to account for cell motion. Our method for topology fixing outperforms the tools provided by SPHARM-MAT and SPHARM-PDM in terms of their successful fixing rates. The different tasks in the presented pipeline for 3D+time shape analysis of cells can be solved using Laplacian approaches, opening the possibility of eventually combining individual steps in order to speed up computations

Decadal predictions with the HiGEM high resolution global coupled climate model: description and basic evaluation

This paper describes the development and basic evaluation of decadal predictions produced using the HiGEM coupled climate model. HiGEM is a higher resolution version of the HadGEM1 Met Office Unified Model. The horizontal resolution in HiGEM has been increased to 1.25◦ × 0.83◦ in longitude and latitude for the atmosphere, and 1/3◦ × 1/3◦ globally for the ocean. The HiGEM decadal predictions are initialised using an anomaly assimilation scheme that relaxes anomalies of ocean temperature and salinity to observed anomalies. 10 year hindcasts are produced for 10 start dates (1960, 1965,..., 2000, 2005). To determine the relative contributions to prediction skill from initial conditions and external forcing, the HiGEM decadal predictions are compared to uninitialised HiGEM transient experiments. The HiGEM decadal predictions have substantial skill for predictions of annual mean surface air temperature and 100 m upper ocean temperature. For lead times up to 10 years, anomaly correlations (ACC) over large areas of the North Atlantic Ocean, the Western Pacific Ocean and the Indian Ocean exceed values of 0.6. Initialisation of the HiGEM decadal predictions significantly increases skill over regions of the Atlantic Ocean,the Maritime Continent and regions of the subtropical North and South Pacific Ocean. In particular, HiGEM produces skillful predictions of the North Atlantic subpolar gyre for up to 4 years lead time (with ACC > 0.7), which are significantly larger than the uninitialised HiGEM transient experiments

The use of phylogeny to interpret cross-cultural patterns in plant use and guide medicinal plant discovery: an example from Pterocarpus (Leguminosae)

The study of traditional knowledge of medicinal plants has led to discoveries that have helped combat diseases and improve healthcare. However, the development of quantitative measures that can assist our quest for new medicinal plants has not greatly advanced in recent years. Phylogenetic tools have entered many scientific fields in the last two decades to provide explanatory power, but have been overlooked in ethnomedicinal studies. Several studies show that medicinal properties are not randomly distributed in plant phylogenies, suggesting that phylogeny shapes ethnobotanical use. Nevertheless, empirical studies that explicitly combine ethnobotanical and phylogenetic information are scarce.In this study, we borrowed tools from community ecology phylogenetics to quantify significance of phylogenetic signal in medicinal properties in plants and identify nodes on phylogenies with high bioscreening potential. To do this, we produced an ethnomedicinal review from extensive literature research and a multi-locus phylogenetic hypothesis for the pantropical genus Pterocarpus (Leguminosae: Papilionoideae). We demonstrate that species used to treat a certain conditions, such as malaria, are significantly phylogenetically clumped and we highlight nodes in the phylogeny that are significantly overabundant in species used to treat certain conditions. These cross-cultural patterns in ethnomedicinal usage in Pterocarpus are interpreted in the light of phylogenetic relationships.This study provides techniques that enable the application of phylogenies in bioscreening, but also sheds light on the processes that shape cross-cultural ethnomedicinal patterns. This community phylogenetic approach demonstrates that similar ethnobotanical uses can arise in parallel in different areas where related plants are available. With a vast amount of ethnomedicinal and phylogenetic information available, we predict that this field, after further refinement of the techniques, will expand into similar research areas, such as pest management or the search for bioactive plant-based compounds

A modified echocardiographic protocol with intrinsic plausibility control to determine intraventricular asynchrony based on TDI and TSI

<p>Abstract</p> <p>Background</p> <p>Established methods to determine asynchrony suffer from high intra- and interobserver variability and failed to improve patient selection for cardiac resynchronization therapy (CRT). Thus, there is a need for easy and robust approaches to reliably assess cardiac asynchrony.</p> <p>Methods and Results</p> <p>We performed echocardiography in 100 healthy subjects and 33 patients with left bundle branch block (LBBB). To detect intraventricular asynchrony, we combined two established methods, i.e., tissue synchronization imaging (TSI) and tissue Doppler imaging (TDI). The time intervals from the onset of aortic valve opening (AVO) to the peak systolic velocity (S') were measured separately in six basal segments in the apical four-, two-, and three-chamber view. Color-coded TSI served as an intrinsic plausibility control and helped to identify the correct S' measuring point in the TDI curves. Next, we identified the segment with the shortest AVO-S' interval. Since this segment most likely represents vital and intact myocardium it served as a reference for other segments. Segments were considered asynchronous when the delay between the segment in question and the reference segment was above the upper limit of normal delays derived from the control population. Intra- and interobserver variability were 7.0% and 7.7%, respectively.</p> <p>Conclusion</p> <p>Our results suggest that combination of TDI and TSI with intrinsic plausibility control improves intra- and interobserver variability and allows easy and reliable assessment of cardiac asynchrony.</p

S100A6 (Calcyclin) is a prostate basal cell marker absent in prostate cancer and its precursors

S100A6 (Calcyclin) is a calcium-binding protein that has been implicated in a variety of biological functions as well as tumorigenesis. The aim of our study was to investigate the involvement of S100A6 during prostate cancer development and progression. Using immunohistochemistry, the expression of S100A6 was examined in benign (n ¼ 66), premalignant (n ¼ 10), malignant (n ¼ 66) and metastatic prostate (n ¼ 5) tissues arranged in a tissue-microarray or whole sections as well as in prostate cancer cell lines. The S100A6 immunostaining pattern in tissues was compared with that of cytokeratin 5 (a basal cell marker) and 18 (a benign luminal cell marker). In all cases of benign epithelium, intense S100A6 expression was seen in the basal cell layer with absent staining in luminal cells. In all cases of prostatic adenocarcinoma (matched), metastatic lesions and 3/10 high-grade prostatic intraepithelial neoplasia lesions, an absence of S100A6 was seen. Western blotting and RT–PCR analysis of cell lines showed S100A6 expression to be absent in LNCaP, LNCaP-LN3 and LNCaP-Pro5 but present in Du145, PC3, PC-3M and PC-3M-LN4. LNCaP cells treated with 5- Azacytidine, caused re-expression of S100A6 mRNA. Sequencing of bisulphite modified DNA showed CpG methylation within the S100A6 promoter region and exon 1 of LNCaP, LNCaP-LN3 and LNCaP-Pro5 cell lines but not in Du145 cells. Our data suggest that loss of S100A6 protein expression is common in prostate cancer development and may occur at an early stage. The mechanism of loss of expression may involve hypermethylation of CpG sites. The finding of intense S100A6 expression in the basal cells of benign glands but loss of expression in cancer could be useful as a novel diagnostic marker for prostate cancer

How to do an evaluation: pitfalls and traps

The recent literature is replete with papers evaluating computational tools (often those operating on 3D structures) for their performance in a certain set of tasks. Most commonly these papers compare a number of docking tools for their performance in cognate re-docking (pose prediction) and/or virtual screening. Related papers have been published on ligand-based tools: pose prediction by conformer generators and virtual screening using a variety of ligand-based approaches. The reliability of these comparisons is critically affected by a number of factors usually ignored by the authors, including bias in the datasets used in virtual screening, the metrics used to assess performance in virtual screening and pose prediction and errors in crystal structures used

How to do an evaluation: pitfalls and traps

The recent literature is replete with papers evaluating computational tools (often those operating on 3D structures) for their performance in a certain set of tasks. Most commonly these papers compare a number of docking tools for their performance in cognate re-docking (pose prediction) and/or virtual screening. Related papers have been published on ligand-based tools: pose prediction by conformer generators and virtual screening using a variety of ligand-based approaches. The reliability of these comparisons is critically affected by a number of factors usually ignored by the authors, including bias in the datasets used in virtual screening, the metrics used to assess performance in virtual screening and pose prediction and errors in crystal structures used

Mechanisms of decadal variability in the Labrador Sea and the wider North Atlantic in a high-resolution climate model

A necessary step before assessing the performance of decadal predictions is the evaluation of the processes that bring memory to the climate system, both in climate models and observations. These mechanisms are particularly relevant in the North Atlantic, where the ocean circulation, related to both the Subpolar Gyre and the Meridional Overturning Circulation (AMOC), is thought to be important for driving significant heat content anomalies. Recently, a rapid decline in observed densities in the deep Labrador Sea has pointed to an ongoing slowdown of the AMOC strength taking place since the mid 90s, a decline also hinted by in-situ observations from the RAPID array. This study explores the use of Labrador Sea densities as a precursor of the ocean circulation changes, by analysing a 300-year long simulation with the state-of-the-art coupled model HadGEM3-GC2. The major drivers of Labrador Sea density variability are investigated, and are characterised by three major contributions. First, the integrated effect of local surface heat fluxes, mainly driven by year-to-year changes in the North Atlantic Oscillation, which accounts for 62% of the total variance. Additionally, two multidecadal-to-centennial contributions from the Greenland-Scotland Ridge outflows are quantified; the first associated with freshwater exports via the East Greenland Current, and the second with density changes in the Denmark Strait Overflow. Finally, evidence is shown that decadal trends in Labrador Sea densities are followed by important atmospheric impacts. In particular, a negative winter NAO response appears to follow the positive Labrador Sea density trends, and provides a phase reversal mechanism

The effects of equine-assisted activities on the social functioning in children and adolescents with autism spectrum disorder

Equine-assisted activities and therapies are increasing in popularity for treatment of ASD symptoms. This research evaluated effects of a 5-week programme of therapeutic riding on social functioning of children/adolescents (N = 15) with ASD. The effectiveness of the programme was evaluated using the autism spectrum quotient, the Vineland Adaptive Behaviour Scale and the empathising and systemising quotient. Results established that the TR intervention increased empathising and reduced maladaptive behaviours. The findings also indicated that specific adaptive behaviours like socialization and communication were not affected by the intervention. Thus, a complex picture of the effects of this intervention emerges: while TR does not change all of the child’s behaviour, it can improve specific aspects of social functioning and also reduce maladaptive ASD traits

Elevational Gradients in Bird Diversity in the Eastern Himalaya: An Evaluation of Distribution Patterns and Their Underlying Mechanisms

BACKGROUND: Understanding diversity patterns and the mechanisms underlying those patterns along elevational gradients is critically important for conservation efforts in montane ecosystems, especially those that are biodiversity hotspots. Despite recent advances, consensus on the underlying causes, or even the relative influence of a suite of factors on elevational diversity patterns has remained elusive. METHODS AND PRINCIPAL FINDINGS: We examined patterns of species richness, density and range size distribution of birds, and the suite of biotic and abiotic factors (primary productivity, habitat variables, climatic factors and geometric constraints) that governs diversity along a 4500-m elevational gradient in the Eastern Himalayan region, a biodiversity hotspot within the world's tallest mountains. We used point count methods for sampling birds and quadrats for estimating vegetation at 22 sites along the elevational gradient. We found that species richness increased to approximately 2000 m, then declined. We found no evidence that geometric constraints influenced this pattern, whereas actual evapotranspiration (a surrogate for primary productivity) and various habitat variables (plant species richness, shrub density and basal area of trees) accounted for most of the variation in bird species richness. We also observed that ranges of most bird species were narrow along the elevation gradient. We find little evidence to support Rapoport's rule for the birds of Sikkim region of the Himalaya. CONCLUSIONS AND SIGNIFICANCE: This study in the Eastern Himalaya indicates that species richness of birds is highest at intermediate elevations along one of the most extensive elevational gradients ever examined. Additionally, primary productivity and factors associated with habitat accounted for most of the variation in avian species richness. The diversity peak at intermediate elevations and the narrow elevational ranges of most species suggest important conservation implications: not only should mid-elevation areas be conserved, but the entire gradient requires equal conservation attention