A stochastic large deformation model for computational anatomy

In the study of shapes of human organs using computational anatomy, variations are found to arise from inter-subject anatomical differences, disease-specific effects, and measurement noise. This paper introduces a stochastic model for incorporating random variations into the Large Deformation Diffeomorphic Metric Mapping (LDDMM) framework. By accounting for randomness in a particular setup which is crafted to fit the geometrical properties of LDDMM, we formulate the template estimation problem for landmarks with noise and give two methods for efficiently estimating the parameters of the noise fields from a prescribed data set. One method directly approximates the time evolution of the variance of each landmark by a finite set of differential equations, and the other is based on an Expectation-Maximisation algorithm. In the second method, the evaluation of the data likelihood is achieved without registering the landmarks, by applying bridge sampling using a stochastically perturbed version of the large deformation gradient flow algorithm. The method and the estimation algorithms are experimentally validated on synthetic examples and shape data of human corpora callosa

MetAMOS: A modular and open source metagenomic assembly and analysis pipeline

© 2013 Treangen et al. We describe MetAMOS, an open source and modular metagenomic assembly and analysis pipeline. MetAMOS represents an important step towards fully automated metagenomic analysis, starting with next-generation sequencing reads and producing genomic scaffolds, open-reading frames and taxonomic or functional annotations. MetAMOS can aid in reducing assembly errors, commonly encountered when assembling metagenomic samples, and improves taxonomic assignment accuracy while also reducing computational cost. MetAMOS can be downloaded from: https://github.com/treangen/MetAMOS

Gene-Boosted Assembly of a Novel Bacterial Genome from Very Short Reads

Recent improvements in technology have made DNA sequencing dramatically faster and more efficient than ever before. The new technologies produce highly accurate sequences, but one drawback is that the most efficient technology produces the shortest read lengths. Short-read sequencing has been applied successfully to resequence the human genome and those of other species but not to whole-genome sequencing of novel organisms. Here we describe the sequencing and assembly of a novel clinical isolate of Pseudomonas aeruginosa, strain PAb1, using very short read technology. From 8,627,900 reads, each 33 nucleotides in length, we assembled the genome into one scaffold of 76 ordered contiguous sequences containing 6,290,005 nucleotides, including one contig spanning 512,638 nucleotides, plus an additional 436 unordered contigs containing 416,897 nucleotides. Our method includes a novel gene-boosting algorithm that uses amino acid sequences from predicted proteins to build a better assembly. This study demonstrates the feasibility of very short read sequencing for the sequencing of bacterial genomes, particularly those for which a related species has been sequenced previously, and expands the potential application of this new technology to most known prokaryotic species

Feature-by-Feature – Evaluating De Novo Sequence Assembly

The whole-genome sequence assembly (WGSA) problem is among one of the most studied problems in computational biology. Despite the availability of a plethora of tools (i.e., assemblers), all claiming to have solved the WGSA problem, little has been done to systematically compare their accuracy and power. Traditional methods rely on standard metrics and read simulation: while on the one hand, metrics like N50 and number of contigs focus only on size without proportionately emphasizing the information about the correctness of the assembly, comparisons performed on simulated dataset, on the other hand, can be highly biased by the non-realistic assumptions in the underlying read generator. Recently the Feature Response Curve (FRC) method was proposed to assess the overall assembly quality and correctness: FRC transparently captures the trade-offs between contigs' quality against their sizes. Nevertheless, the relationship among the different features and their relative importance remains unknown. In particular, FRC cannot account for the correlation among the different features. We analyzed the correlation among different features in order to better describe their relationships and their importance in gauging assembly quality and correctness. In particular, using multivariate techniques like principal and independent component analysis we were able to estimate the “excess-dimensionality” of the feature space. Moreover, principal component analysis allowed us to show how poorly the acclaimed N50 metric describes the assembly quality. Applying independent component analysis we identified a subset of features that better describe the assemblers performances. We demonstrated that by focusing on a reduced set of highly informative features we can use the FRC curve to better describe and compare the performances of different assemblers. Moreover, as a by-product of our analysis, we discovered how often evaluation based on simulated data, obtained with state of the art simulators, lead to not-so-realistic results

Estrogen regulation of TRPM8 expression in breast cancer cells

<p>Abstract</p> <p>Background</p> <p>The calcium-permeable cation channel TRPM8 (melastatin-related transient receptor potential member 8) is over-expressed in several cancers. The present study aimed at investigating the expression, function and potential regulation of TRPM8 channels by ER alpha (estrogen receptor alpha) in breast cancer.</p> <p>Methods</p> <p>RT-PCR, Western blot, immuno-histochemical, and siRNA techniques were used to investigate TRPM8 expression, its regulation by estrogen receptors, and its expression in breast tissue. To investigate the channel activity in MCF-7 cells, we used the whole cell patch clamp and the calcium imaging techniques.</p> <p>Results</p> <p>TRPM8 channels are expressed at both mRNA and protein levels in the breast cancer cell line MCF-7. Bath application of the potent TRPM8 agonist Icilin (20 μM) induced a strong outwardly rectifying current at depolarizing potentials, which is associated with an elevation of cytosolic calcium concentration, consistent with established TRPM8 channel properties. RT-PCR experiments revealed a decrease in TRPM8 mRNA expression following steroid deprivation for 48 and 72 hours. In steroid deprived medium, addition of 17-beta-estradiol (E₂, 10 nM) increased both TRPM8 mRNA expression and the number of cells which respond to Icilin, but failed to affect the Ca²⁺entry amplitude. Moreover, silencing ERα mRNA expression with small interfering RNA reduced the expression of TRPM8. Immuno-histochemical examination of the expression of TRPM8 channels in human breast tissues revealed an over-expression of TRPM8 in breast adenocarcinomas, which is correlated with estrogen receptor positive (ER⁺) status of the tumours.</p> <p>Conclusion</p> <p>Taken together, these results show that TRPM8 channels are expressed and functional in breast cancer and that their expression is regulated by ER alpha.</p

Bone mineral density and body composition in postmenopausal women with psoriasis and psoriatic arthritis

Introduction: the aim of the present study was to compare bone mineral density (BMD) and body composition (BC) measurements as well as identify risk factors for low BMD and osteoporotic fractures in postmenopausal women with psoriasis (Ps) and psoriatic arthritis (PsA).Methods: A cross-sectional study was carried out in 45 PsA women, 52 Ps women and 98 healthy female controls (HC). Clinical risk factors for low bone density and osteoporotic fracture were evaluated by a specific questionnaire. An X-ray absorptiometry (DXA) at the lumbar spine, total femur and total body was performed on all patients. Skin and joint outcomes were measured by specific tools (PASI, HAQ and DAS28). Morphometric vertebral fractures were evaluated by lumbar and thoracic spine X-ray, according to Genant's method.Results: There were no significant differences in age, body mass index (BMI), total lean mass and bone mineral density among the groups. However, the PsA group had a significantly higher body fat percentage (BF%) than the Ps and HC groups. Osteoporotic fractures were more frequently observed in PsA and Ps groups than in the HC group (P = 0.01). Recurrent falls and a longer duration of disease increased the risk of fracture (odds ratio (OR) = 18.3 and 1.08, respectively) in the PsA group (P = 0.02). Disability was the main factor related to osteoporotic fracture in the Ps group (odds ratio (OR) = 11.1) (P = 0.02).Conclusions: Ps and PsA patients did not present lower BMD. However, they had a higher prevalence of osteoporotic fractures and higher risk of metabolic syndrome. Patients with a longer duration of disease, disability and recurrent falls need preventive measures.Rheumatology Division at UNIFESP/EPMUniversidade Federal de São Paulo, UNIFESP Paulista Sch Med, Div Rheumatol, EPM, BR-04023900 São Paulo, BrazilUniversidade Federal de São Paulo, UNIFESP Paulista Sch Med, Div Rheumatol, EPM, BR-04023900 São Paulo, BrazilWeb of Scienc

Developmentally regulated GTP binding protein 1 (DRG1) controls microtubule dynamics

The mitotic spindle, essential for segregating the sister chromatids into the two evolving daughter cells, is composed of highly dynamic cytoskeletal filaments, the microtubules. The dynamics of microtubules are regulated by numerous microtubule associated proteins. We identify here Developmentally regulated GTP binding protein 1 (DRG1) as a microtubule binding protein with diverse microtubule-associated functions. In vitro, DRG1 can diffuse on microtubules, promote their polymerization, drive microtubule formation into bundles, and stabilize microtubules. HeLa cells with reduced DRG1 levels show delayed progression from prophase to anaphase because spindle formation is slowed down. To perform its microtubule-associated functions, DRG1, although being a GTPase, does not require GTP hydrolysis. However, all domains are required as truncated versions show none of the mentioned activities besides microtubule binding