Comparative genomics of metabolic networks of free-living and parasitic eukaryotes

Background Obligate endoparasites often lack particular metabolic pathways as compared to free-living organisms. This phenomenon comprises anabolic as well as catabolic reactions. Presumably, the corresponding enzymes were lost in adaptation to parasitism. Here we compare the predicted core metabolic graphs of obligate endoparasites and non-parasites (free living organisms and facultative parasites) in order to analyze how the parasites' metabolic networks shrunk in the course of evolution. Results Core metabolic graphs comprising biochemical reactions present in the presumed ancestor of parasites and non-parasites were reconstructed from the Kyoto Encyclopedia of Genes and Genomes. While the parasites' networks had fewer nodes (metabolites) and edges (reactions), other parameters such as average connectivity, network diameter and number of isolated edges were similar in parasites and non-parasites. The parasites' networks contained a higher percentage of ATP-consuming reactions and a lower percentage of NAD-requiring reactions. Control networks, shrunk to the size of the parasites' by random deletion of edges, were scale-free but exhibited smaller diameters and more isolated edges. Conclusions The parasites' networks were smaller than those of the non-parasites regarding number of nodes or edges, but not regarding network diameters. Network integrity but not scale-freeness has acted as a selective principle during the evolutionary reduction of parasite metabolism. ATP-requiring reactions in particular have been retained in the parasites' core metabolism while NADH- or NADPH-requiring reactions were lost preferentially.

A Parametric Test of Perfect Ranking in Balanced Ranked Set Sampling

Many techniques based on data which are drawn by Ranked Set Sampling (RSS) scheme assume that the ranking of observations is perfect. Therefore it is essential to develop some methods for testing this assumption. In this article, we propose a parametric location-scale free test for assessing the assumption of perfect ranking. The results of a simulation study in two special cases of normal and exponential distributions indicate that the proposed test performs well in comparison with its leading competitors.

Parameter-free extraction of the functional network topology from intracranial EEG recordings: a time-resolved study of graph properties in focal onset seizures

Rationale: Focal onset epileptic seizures are due to abnormal interactions between distributed brain areas. By estimating the cross-correlation matrix of multi-site intra-cerebral EEG recordings (iEEG), one can quantify these interactions. To assess the topology of the underlying functional network, the binary connectivity matrix has to be derived from the cross-correlation matrix by use of a threshold. Classically, a unique threshold is used that constrains the topology [1]. Our method aims to set the threshold in a data-driven way by separating genuine from random cross-correlation. We compare our approach to the fixed threshold method and study the dynamics of the functional topology. Methods: We investigate the iEEG of patients suffering from focal onset seizures who underwent evaluation for the possibility of surgery. The equal-time cross-correlation matrices are evaluated using a sliding time window. We then compare 3 approaches assessing the corresponding binary networks. For each time window: * Our parameter-free method derives from the cross-correlation strength matrix (CCS)[2]. It aims at disentangling genuine from random correlations (due to finite length and varying frequency content of the signals). In practice, a threshold is evaluated for each pair of channels independently, in a data-driven way. * The fixed mean degree (FMD) uses a unique threshold on the whole connectivity matrix so as to ensure a user defined mean degree. * The varying mean degree (VMD) uses the mean degree of the CCS network to set a unique threshold for the entire connectivity matrix. * Finally, the connectivity (c), connectedness (given by k, the number of disconnected sub-networks), mean global and local efficiencies (Eg, El, resp.) are computed from FMD, CCS, VMD, and their corresponding random and lattice networks. Results: Compared to FMD and VMD, CCS networks present: *topologies that are different in terms of c, k, Eg and El. *from the pre-ictal to the ictal and then post-ictal period, topological features time courses that are more stable within a period, and more contrasted from one period to the next. For CCS, pre-ictal connectivity is low, increases to a high level during the seizure, then decreases at offset. k shows a ‘‘U-curve’’ underlining the synchronization of all electrodes during the seizure. Eg and El time courses fluctuate between the corresponding random and lattice networks values in a reproducible manner. Conclusions: The definition of a data-driven threshold provides new insights into the topology of the epileptic functional networks.

High-density large-scale field emitter arrays for x-ray free electron laser cathodes

High brightness electron sources are of great importance for the operation of the hard X-ray free electron lasers. Field emission cathodes based on the double-gate metallic field emitter arrays (FEAs) can potentially offer higher brightness than the currently used ones. We report on the successful application of electron beam lithography for fabrication of the large-scale single-gate as well as double-gate FEAs. We demonstrate operational high-density single-gate FEAs with sub-micron pitch and total number of tips up to 106 as well as large-scale double-gate FEAs with large collimation gate apertures. The details of design, fabrication procedure and successful measurements of the emission current from the single- and double-gate cathodes are presented.

Influence of topology on the scale setting

Recently a new method to set the scale in lattice gauge theories, based on the gradient flow generated by the Wilson action, has been proposed, and the systematic errors of the new scales t0 and w0 have been investigated by various groups. The Wilson flow provides also an interesting alternative smoothing procedure particularly useful for the measurement of the topological charge as a pure gluonic observable. We show the viability of this method for N=1 supersymmetric Yang-Mills theory by analysing the configurations produced by the DESY-Muenster Collaboration. The relation between the scale and the topological charge has been investigated showing a strong correlation. We have found that the scale has a linear dependence on the topological charge, the slope of which increases decreasing the volume and the gluino mass. Moreover we have investigated this dependence as a function of the reference parameter used to define the scale: the tuning of this parameter turns out to be fundamental for a more reliable scale setting. Similar conclusions hold for the Sommer parameter r0.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Early molecular and cytogenetic response is predictive for long-term progression-free and overall survival in chronic myeloid leukemia (CML)

In the face of competing first-line treatment options for CML, early prediction of prognosis on imatinib is desirable to assure favorable survival or otherwise consider the use of a second-generation tyrosine kinase inhibitor (TKI). A total of 1303 newly diagnosed imatinib-treated patients (pts) were investigated to correlate molecular and cytogenetic response at 3 and 6 months with progression-free and overall survival (PFS, OS). The persistence of BCR-ABL transcript levels >10% according to the international scale (BCR-ABL(IS)) at 3 months separated a high-risk group (28% of pts; 5-year OS: 87%) from a group with >1-10% BCR-ABL(IS) (41% of pts; 5-year OS: 94%; P=0.012) and from a group with 1% BCR-ABL(IS) (31% of pts; 5-year OS: 97%; P=0.004). Cytogenetics identified high-risk pts by >35% Philadelphia chromosome-positive metaphases (Ph+, 27% of pts; 5-year OS: 87%) compared with 35% Ph+ (73% of pts; 5-year OS: 95%; P=0.036). At 6 months, >1% BCR-ABL(IS) (37% of pts; 5-year OS: 89%) was associated with inferior survival compared with 1% (63% of pts; 5-year OS: 97%; P<0.001) and correspondingly >0% Ph+ (34% of pts; 5-year OS: 91%) compared with 0% Ph+ (66% of pts; 5-year OS: 97%; P=0.015). Treatment optimization is recommended for pts missing these landmarks.

Why is 11beta-hydroxysteroid dehydrogenase type 1 facing the endoplasmic reticulum lumen? Physiological relevance of the membrane topology of 11beta-HSD1

11Beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) is essential for the local activation of glucocorticoid receptors (GR). Unlike unliganded cytoplasmic GR, 11beta-HSD1 is an endoplasmic reticulum (ER)-membrane protein with lumenal orientation. Cortisone might gain direct access to 11beta-HSD1 by free diffusion across membranes, indirectly via intracellular binding proteins or, alternatively, by insertion into membranes. Membranous cortisol, formed by 11beta-HSD1 at the ER-lumenal side, might then activate cytoplasmic GR or bind to ER-lumenal secretory proteins. Compartmentalization of 11beta-HSD1 is important for its regulation by hexose-6-phosphate dehydrogenase (H6PDH), which regenerates cofactor NADPH in the ER lumen and stimulates oxoreductase activity. ER-lumenal orientation of 11beta-HSD1 is also essential for the metabolism of the alternative substrate 7-ketocholesterol (7KC), a major cholesterol oxidation product found in atherosclerotic plaques and taken up from processed cholesterol-rich food. An 11beta-HSD1 mutant adopting cytoplasmic orientation efficiently catalyzed the oxoreduction of cortisone but not 7KC, indicating access to cortisone from both sides of the ER-membrane but to 7KC only from the lumenal side. These aspects may be relevant for understanding the physiological role of 11beta-HSD1 and for developing therapeutic interventions to control glucocorticoid reactivation.

Correspondence establishment in statistical modeling of shapes with arbitrary topology

Correspondence establishment is a key step in statistical shape model building. There are several automated methods for solving this problem in 3D, but they usually can only handle objects with simple topology, like that of a sphere or a disc. We propose an extension to correspondence establishment over a population based on the optimization of the minimal description length function, allowing considering objects with arbitrary topology. Instead of using a fixed structure of kernel placement on a sphere for the systematic manipulation of point landmark positions, we rely on an adaptive, hierarchical organization of surface patches. This hierarchy can be built on surfaces of arbitrary topology and the resulting patches are used as a basis for a consistent, multi-scale modification of the surfaces' parameterization, based on point distribution models. The feasibility of the approach is demonstrated on synthetic models with different topologies.

Changes in species composition in alpine snowbeds with climate change inferred from small-scale spatial patterns

Alpine snowbeds are characterised by a very short growing season. However, the length of the snow-free period is increasingly prolonged due to climate change, so that snowbeds become susceptible to invasions from neighbouring alpine meadow communities. We hypothesised that spatial distribution of species generated by plant interactions may indicate whether snowbed species will coexist with or will be out-competed by invading alpine species – spatial aggregation or segregation will point to coexistence or competitive exclusion, respectively. We tested this hypothesis in snowbeds of the Swiss Alps using the variance ratio statistics. We focused on the relationships between dominant snowbed species, subordinate snowbed species, and potentially invading alpine grassland species. Subordinate snowbed species were generally spatially aggregated with each other, but were segregated from alpine grassland species. Competition between alpine grassland and subordinate snowbed species may have caused this segregation. Segregation between these species groups increased with earlier snowmelt, suggesting an increasing importance of competition with climate change. Further, a dominant snowbed species (Alchemilla pentaphyllea) was spatially aggregated with subordinate snowbed species, while two other dominants (Gnaphalium supinum and Salix herbacea) showed aggregated patterns with alpine grassland species. These dominant species are known to show distinct microhabitat preferences suggesting the existence of hidden microhabitats with different susceptibility to invaders. These results allow us to suggest that alpine snowbed areas are likely to be reduced as a consequence of climate change and that invading species from nearby alpine grasslands could outcompete subordinate snowbed species. On the other hand, microhabitats dominated by Gnaphalium or Salix seem to be particularly prone to invasions by non-snowbed species.

Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale

The development of susceptibility maps for debris flows is of primary importance due to population pressure in hazardous zones. However, hazard assessment by process-based modelling at a regional scale is difficult due to the complex nature of the phenomenon, the variability of local controlling factors, and the uncertainty in modelling parameters. A regional assessment must consider a simplified approach that is not highly parameter dependant and that can provide zonation with minimum data requirements. A distributed empirical model has thus been developed for regional susceptibility assessments using essentially a digital elevation model (DEM). The model is called Flow-R for Flow path assessment of gravitational hazards at a Regional scale (available free of charge under http://www.flow-r.org) and has been successfully applied to different case studies in various countries with variable data quality. It provides a substantial basis for a preliminary susceptibility assessment at a regional scale. The model was also found relevant to assess other natural hazards such as rockfall, snow avalanches and floods. The model allows for automatic source area delineation, given user criteria, and for the assessment of the propagation extent based on various spreading algorithms and simple frictional laws. We developed a new spreading algorithm, an improved version of Holmgren's direction algorithm, that is less sensitive to small variations of the DEM and that is avoiding over-channelization, and so produces more realistic extents. The choices of the datasets and the algorithms are open to the user, which makes it compliant for various applications and dataset availability. Amongst the possible datasets, the DEM is the only one that is really needed for both the source area delineation and the propagation assessment; its quality is of major importance for the results accuracy. We consider a 10 m DEM resolution as a good compromise between processing time and quality of results. However, valuable results have still been obtained on the basis of lower quality DEMs with 25 m resolution.

Reduced bone breakage and increased bone strength in free range laying hens fed omega-3 polyunsaturated fatty acid supplemented diets

INTRODUCTION The omega-3 and omega-6 polyunsaturated fatty acids (PUFAs) are the immediate precursors to a number of important mediators of immunity, inflammation and bone function, with products of omega-6 generally thought to promote inflammation and favour bone resorption. Western diets generally provide a 10 to 20-fold deficit in omega-3 PUFAs compared with omega-6, and this is thought to have contributed to the marked rise in incidence of disorders of modern human societies, such as heart disease, colitis and perhaps osteoporosis. Many of our food production animals, fed on grains rich in omega-6, are also exposed to a dietary deficit in omega-3, with perhaps similar health consequences. Bone fragility due to osteoporotic changes in laying hens is a major economic and welfare problem, with our recent estimates of breakage rates indicating up to 95% of free range hens suffer breaks during lay. METHODS Free range hens housed in full scale commercial systems were provided diets supplemented with omega-3 alpha linolenic acid, and the skeletal benefits were investigated by comparison to standard diets rich in omega-6. RESULTS There was a significant 40-60% reduction in keel bone breakage rate, and a corresponding reduction in breakage severity in the omega-3 supplemented hens. There was significantly greater bone density and bone mineral content, alongside increases in total bone and trabecular volumes. The mechanical properties of the omega-3 supplemented hens were improved, with strength, energy to break and stiffness demonstrating significant increases. Alkaline phosphatase (an osteoblast marker) and tartrate-resistant acid phosphatase (an osteoclast marker) both showed significant increases with the omega-3 diets, indicating enhanced bone turnover. This was corroborated by the significantly lower levels of the mature collagen crosslinks, hydroxylysyl pyridinoline, lysyl pyridinoline and histidinohydroxy-lysinonorleucine, with a corresponding significant shift in the mature:immature crosslink ratio. CONCLUSIONS The improved skeletal health in laying hens corresponds to as many as 68million fewer hens suffering keel fractures in the EU each year. The biomechanical and biochemical evidence suggests that increased bone turnover has enhanced the bone mechanical properties, and that this may suggest potential benefits for human osteoporosis.