Safe and complete contig assembly via omnitigs

Contig assembly is the first stage that most assemblers solve when reconstructing a genome from a set of reads. Its output consists of contigs -- a set of strings that are promised to appear in any genome that could have generated the reads. From the introduction of contigs 20 years ago, assemblers have tried to obtain longer and longer contigs, but the following question was never solved: given a genome graph $G$ (e.g. a de Bruijn, or a string graph), what are all the strings that can be safely reported from $G$ as contigs? In this paper we finally answer this question, and also give a polynomial time algorithm to find them. Our experiments show that these strings, which we call omnitigs, are 66% to 82% longer on average than the popular unitigs, and 29% of dbSNP locations have more neighbors in omnitigs than in unitigs.Comment: Full version of the paper in the proceedings of RECOMB 201

A knowledge-driven GIS modeling technique for groundwater potential mapping at the Upper Langat Basin, Malaysia.

The aim of this paper is to use a knowledge-driven expert-based geographical information system (GIS) model coupling with remote-sensing-derived parameters for groundwater potential mapping in an area of the Upper Langat Basin, Malaysia. In this study, nine groundwater storage controlling parameters that affect groundwater occurrences are derived from remotely sensed imagery, available maps, and associated databases. Those parameters are: lithology, slope, lineament, land use, soil, rainfall, drainage density, elevation, and geomorphology. Then the parameter layers were integrated and modeled using a knowledge-driven GIS of weighted linear combination. The weightage and score for each parameter and their classes are based on the Malaysian groundwater expert opinion survey. The predicted groundwater potential map was classified into four distinct zones based on the classification scheme designed by Department of Minerals and Geoscience Malaysia (JMG). The results showed that about 17% of the study area falls under low-potential zone, with 66% on moderate-potential zone, 15% with high-potential zone, and only 0.45% falls under very-high-potential zone. The results obtained in this study were validated with the groundwater borehole wells data compiled by the JMG and showed 76% of prediction accuracy. In addition statistical analysis indicated that hard rock dominant of the study area is controlled by secondary porosity such as distance from lineament and density of lineament. There are high correlations between area percentage of predicted groundwater potential zones and groundwater well yield. Results obtained from this study can be useful for future planning of groundwater exploration, planning and development by related agencies in Malaysia which provide a rapid method and reduce cost as well as less time consuming. The results may be also transferable to other areas of similar hydrological characteristics

A model describing diffusion in prostate cancer

PURPOSE: Quantitative diffusion MRI has frequently been studied as a means of grading prostate cancer. Interpretation of results is complicated by the nature of prostate tissue, which consists of four distinct compartments: vascular, ductal lumen, epithelium, and stroma. Current diffusion measurements are an ill-defined weighted average of these compartments. In this study, prostate diffusion is analyzed in terms of a model that takes explicit account of tissue compartmentalization, exchange effects, and the non-Gaussian behavior of tissue diffusion.  METHOD: The model assumes that exchange between the cellular (ie, stromal plus epithelial) and the vascular and ductal compartments is slow. Ductal and cellular diffusion characteristics are estimated by Monte Carlo simulation and a two-compartment exchange model, respectively. Vascular pseudodiffusion is represented by an additional signal at b = 0. Most model parameters are obtained either from published data or by comparing model predictions with the published results from 41 studies. Model prediction error is estimated using 10-fold cross-validation.  RESULTS: Agreement between model predictions and published results is good. The model satisfactorily explains the variability of ADC estimates found in the literature.  CONCLUSION: A reliable model that predicts the diffusion behavior of benign and cancerous prostate tissue of different Gleason scores has been developed. Magn Reson Med, 2016. © 2016 Wiley Periodicals, Inc

Income in Adult Survivors of Childhood Cancer.

INTRODUCTION: Little is known about the impact of childhood cancer on the personal income of survivors. We compared income between survivors and siblings, and determined factors associated with income. METHODS: As part of the Swiss Childhood Cancer Survivor Study (SCCSS), a questionnaire was sent to survivors, aged ≥18 years, registered in the Swiss Childhood Cancer Registry (SCCR), diagnosed at age 4'500 CHF), even after we adjusted for socio-demographic and educational factors (OR = 0.46, p<0.001). Older age, male sex, personal and parental education, and number of working hours were associated with high income. Survivors of leukemia (OR = 0.40, p<0.001), lymphoma (OR = 0.63, p = 0.040), CNS tumors (OR = 0.22, p<0.001), bone tumors (OR = 0.24, p = 0.003) had a lower income than siblings. Survivors who had cranial irradiation, had a lower income than survivors who had no cranial irradiation (OR = 0.48, p = 0.006). DISCUSSION: Even after adjusting for socio-demographic characteristics, education and working hours, survivors of various diagnostic groups have lower incomes than siblings. Further research needs to identify the underlying causes

FlexOracle: predicting flexible hinges by identification of stable domains

<p>Abstract</p> <p>Background</p> <p>Protein motions play an essential role in catalysis and protein-ligand interactions, but are difficult to observe directly. A substantial fraction of protein motions involve hinge bending. For these proteins, the accurate identification of flexible hinges connecting rigid domains would provide significant insight into motion. Programs such as GNM and FIRST have made global flexibility predictions available at low computational cost, but are not designed specifically for finding hinge points.</p> <p>Results</p> <p>Here we present the novel FlexOracle hinge prediction approach based on the ideas that energetic interactions are stronger <it>within </it>structural domains than <it>between </it>them, and that fragments generated by cleaving the protein at the hinge site are independently stable. We implement this as a tool within the Database of Macromolecular Motions, MolMovDB.org. For a given structure, we generate pairs of fragments based on scanning all possible cleavage points on the protein chain, compute the energy of the fragments compared with the undivided protein, and predict hinges where this quantity is minimal. We present three specific implementations of this approach. In the first, we consider only pairs of fragments generated by cutting at a <it>single </it>location on the protein chain and then use a standard molecular mechanics force field to calculate the enthalpies of the two fragments. In the second, we generate fragments in the same way but instead compute their free energies using a knowledge based force field. In the third, we generate fragment pairs by cutting at <it>two </it>points on the protein chain and then calculate their free energies.</p> <p>Conclusion</p> <p>Quantitative results demonstrate our method's ability to predict known hinges from the Database of Macromolecular Motions.</p

Multisensory Perceptual Learning of Temporal Order: Audiovisual Learning Transfers to Vision but Not Audition

Background: An outstanding question in sensory neuroscience is whether the perceived timing of events is mediated by a central supra-modal timing mechanism, or multiple modality-specific systems. We use a perceptual learning paradigm to address this question. Methodology/Principal Findings: Three groups were trained daily for 10 sessions on an auditory, a visual or a combined audiovisual temporal order judgment (TOJ). Groups were pre-tested on a range TOJ tasks within and between their group modality prior to learning so that transfer of any learning from the trained task could be measured by post-testing other tasks. Robust TOJ learning (reduced temporal order discrimination thresholds) occurred for all groups, although auditory learning (dichotic 500/2000 Hz tones) was slightly weaker than visual learning (lateralised grating patches). Crossmodal TOJs also displayed robust learning. Post-testing revealed that improvements in temporal resolution acquired during visual learning transferred within modality to other retinotopic locations and orientations, but not to auditory or crossmodal tasks. Auditory learning did not transfer to visual or crossmodal tasks, and neither did it transfer within audition to another frequency pair. In an interesting asymmetry, crossmodal learning transferred to all visual tasks but not to auditory tasks. Finally, in all conditions, learning to make TOJs for stimulus onsets did not transfer at all to discriminating temporal offsets. These data present a complex picture of timing processes

Engineered hexavalent Fc proteins with enhanced Fc-gamma receptor avidity provide insights into immune-complex interactions

Tania Rowley et al. present multivalent Fc molecules with enhanced avidity for Fc gamma receptors in order to improve the treatment of autoantibody-mediated human diseases. They found several key amino acids involved in Fc receptor binding interactions

Expression of thymidylate synthase in human cells is an early G1 event regulated by CDK4 and p16INK4A but not E2F

Thymidylate synthase (TS) is the enzyme that catalyses the last step in de novo thymidylate synthesis. It is of interest clinically because it is an effective target for drugs such as 5-fluorouracil, often used in combination therapy. Despite a number of earlier reports indicating that TS is a cell cycle-dependent enzyme, this remains equivocal. Here, we show that in HCT116 cells synchronised by serum starvation, there is a clear dissociation between the expression of cyclin E (a well-characterised cell-cycle protein) and TS. Although both cyclin E and TS mRNA and protein increased during G1, TS upregulation was delayed. Moreover, TS levels did not decrease following S-phase completion while cyclin E decreased sharply. Similarly, clear differences were seen between cyclin E and TS as asynchronously growing HCT116 cells were growth-inhibited by low-serum treatment. In contrast to previous reports using rodent cells, adenovirus-mediated over-expression of E2F1 and cyclin E in three human cell lines had no effect on TS. Cell-cycle progression was blocked by treatment of cells with pharmacological inhibitors of CDK2 and CDK4 and by ectopic expression of p16INK4A. Whereas CDK2 inhibition had no effect on TS levels, inhibition of CDK4 was associated with decreased TS protein levels. These results provide the first evidence that drugs targeting CDK4 may be useful with anti-TS drugs as combination therapy for cancer

Exoerythrocytic Plasmodium Parasites Secrete a Cysteine Protease Inhibitor Involved in Sporozoite Invasion and Capable of Blocking Cell Death of Host Hepatocytes

Plasmodium parasites must control cysteine protease activity that is critical for hepatocyte invasion by sporozoites, liver stage development, host cell survival and merozoite liberation. Here we show that exoerythrocytic P. berghei parasites express a potent cysteine protease inhibitor (PbICP, P. berghei inhibitor of cysteine proteases). We provide evidence that it has an important function in sporozoite invasion and is capable of blocking hepatocyte cell death. Pre-incubation with specific anti-PbICP antiserum significantly decreased the ability of sporozoites to infect hepatocytes and expression of PbICP in mammalian cells protects them against peroxide- and camptothecin-induced cell death. PbICP is secreted by sporozoites prior to and after hepatocyte invasion, localizes to the parasitophorous vacuole as well as to the parasite cytoplasm in the schizont stage and is released into the host cell cytoplasm at the end of the liver stage. Like its homolog falstatin/PfICP in P. falciparum, PbICP consists of a classical N-terminal signal peptide, a long N-terminal extension region and a chagasin-like C-terminal domain. In exoerythrocytic parasites, PbICP is posttranslationally processed, leading to liberation of the C-terminal chagasin-like domain. Biochemical analysis has revealed that both full-length PbICP and the truncated C-terminal domain are very potent inhibitors of cathepsin L-like host and parasite cysteine proteases. The results presented in this study suggest that the inhibitor plays an important role in sporozoite invasion of host cells and in parasite survival during liver stage development by inhibiting host cell proteases involved in programmed cell death