Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

Size-dependent predation by Dugesia lugubris (Turbellaria) on Physa acuta (Gastropoda): experiments and model

1. We investigated experimentally predation by the flatworm Dugesia lugubris on the snail Physa acuta in relation to predator body length and to prey morphology [shell length (SL) and aperture width (AW)]. 2. SL and AW correlate strongly in the field, but display significant and independent variance among populations. In the laboratory, predation by Dugesia resulted in large and significant selection differentials on both SL and AW. Analysis of partial effects suggests that selection on AW was indirect, and mediated through its strong correlation with SL. 3. The probability P(ij) for a snail of size category i (SL) to be preyed upon by a flatworm of size category j was fitted with a Poisson-probability distribution, the mean of which increased linearly with predator size (i). Despite the low number of parameters, the fit was excellent (r2 = 0.96). We offer brief biological interpretations of this relationship with reference to optimal foraging theory. 4. The largest size class of Dugesia (>2 cm) did not prey on snails larger than 7 mm shell length. This size threshold might offer Physa a refuge against flatworm predation and thereby allow coexistence in the field. 5. Our results are further discussed with respect to previous field and laboratory observations on P acuta life-history patterns, in particular its phenotypic variance in adult body size.

Vasopressin-dependent coupling between sodium transport and water flow in a mouse cortical collecting duct cell line.

Water balance is achieved through the ability of the kidney to control water reabsorption in the connecting tubule and the collecting duct. In a mouse cortical collecting duct cell line (mCCD(c11)), physiological concentrations of arginine vasopressin increased both electrogenic, amiloride-sensitive, epithelial sodium channel (ENaC)-mediated sodium transport measured by the short-circuit current (Isc) method and water flow (Jv apical to basal) measured by gravimetry with similar activation coefficient K(1/2) (6 and 12 pM, respectively). Jv increased linearly according to the osmotic gradient across the monolayer. A small but highly significant Jv was also measured under isoosmotic conditions. To test the coupling between sodium reabsorption and water flow, mCCD(c11) cells were treated for 24 h under isoosmotic condition with either diluent, amiloride, vasopressin or vasopressin and amiloride. Isc, Jv, and net chemical sodium fluxes were measured across the same monolayers. Around 30% of baseline and 50% of vasopressin-induced water flow is coupled to an amiloride-sensitive, ENaC-mediated, electrogenic sodium transport, whereas the remaining flow is coupled to an amiloride-insensitive, nonelectrogenic sodium transport mediated by an unknown electroneutral transporter. The mCCD(c11) cell line is a first example of a mammalian tight epithelium allowing quantitative study of the coupling between sodium and water transport. Our data are consistent with the 'near isoosmotic' fluid transport model.

Mutations inducing divergent shifts of constitutive activity reveal different modes of binding among catecholamine analogues to the beta(2)-adrenergic receptor.

1. We compared the changes in binding energy generated by two mutations that shift in divergent directions the constitutive activity of the human beta(2) adrenergic receptor (beta(2)AR). 2. A constitutively activating mutant (CAM) and the double alanine replacement (AA mutant) of catechol-binding serines (S204A, S207A) in helix 5 were stably expressed in CHO cell lines, and used to measure the binding affinities of more than 40 adrenergic ligands. Moreover, the efficacy of the same group of compounds was determined as intrinsic activity for maximal adenylyl cyclase stimulation in wild-type beta(2)AR. 3. Although the two mutations had opposite effects on ligand affinity, the extents of change were in both cases largely correlated with the degree of ligand efficacy. This was particularly evident if the extra loss of binding energy due to hydrogen bond deletion in the AA mutant was taken into account. Thus the data demonstrate that there is an overall linkage between the configuration of the binding pocket and the intrinsic equilibrium between active and inactive receptor forms. 4. We also found that AA mutation-induced affinity changes for catecholamine congeners gradually lacking ethanolamine substituents were linearly correlated to the loss of affinity that such modifications of the ligand cause for wild-type receptor. This indicates that the strength of bonds between catechol ring and helix 5 is critically dependent on the rest of interactions of the beta-ethanolamine tail with other residues of the beta(2)-AR binding pocket.

Arsenic(V) adsorption onto α-Al2O3 between 25 and 70°C

The adsorption of As(V) onto alpha -Al2O3 was investigated at 25, 50 and 70 degreesC using batch adsorption experiments. Results indicate that As is strongly adsorbed at low pH and gets progressively released to the fluid with increasing pH above 7. At any pH, increasing temperature favors aqueous species of As over surface species. Surface complexation constants were determined at the experimental temperatures by fitting the adsorption data. Adsorption reactions were then converted to semi-isocolumbic reactions, i.e, reactions with balanced like-charged aqueous species. Intrinsic adsorption constants of semi-isocolumbic reactions change linearly when plotted against inverse temperature, suggesting that the heat capacity of these reactions remains constant over the temperature range considered. This permitted thermodynamic parameters of intrinsic surface complexation constants to be determined. Changes in surface complexation constants result in a change in the surface speciation with increasing temperature. This change is similar to the one observed for aqueous species, i.e. increasing temperature favors less negatively charged species below a pH of 9 and more negatively charged species above a pH of 10. Comparison with the stability of As surface complexes with Fe suggests that surface complexes with Al are more stable. (C) 2001 Elsevier Science Ltd. All rights reserved.

Demographic effects of extreme winter weather in the barn owl.

Extreme weather events can lead to immediate catastrophic mortality. Due to their rare occurrence, however, the long-term impacts of such events for ecological processes are unclear. We examined the effect of extreme winters on barn owl (Tyto alba) survival and reproduction in Switzerland over a 68-year period (approximately 20 generations). This long-term data set allowed us to compare events that occurred only once in several decades to more frequent events. Winter harshness explained 17 and 49% of the variance in juvenile and adult survival, respectively, and the two harshest winters were associated with major population crashes caused by simultaneous low juvenile and adult survival. These two winters increased the correlation between juvenile and adult survival from 0.63 to 0.69. Overall, survival decreased non-linearly with increasing winter harshness in adults, and linearly in juveniles. In contrast, brood size was not related to the harshness of the preceding winter. Our results thus reveal complex interactions between climate and demography. The relationship between weather and survival observed during regular years is likely to underestimate the importance of climate variation for population dynamics.

Network revenue management with product-specific no-shows

Revenue management practices often include overbooking capacity to account for customerswho make reservations but do not show up. In this paper, we consider the network revenuemanagement problem with no-shows and overbooking, where the show-up probabilities are specificto each product. No-show rates differ significantly by product (for instance, each itinerary andfare combination for an airline) as sale restrictions and the demand characteristics vary byproduct. However, models that consider no-show rates by each individual product are difficultto handle as the state-space in dynamic programming formulations (or the variable space inapproximations) increases significantly. In this paper, we propose a randomized linear program tojointly make the capacity control and overbooking decisions with product-specific no-shows. Weestablish that our formulation gives an upper bound on the optimal expected total profit andour upper bound is tighter than a deterministic linear programming upper bound that appearsin the existing literature. Furthermore, we show that our upper bound is asymptotically tightin a regime where the leg capacities and the expected demand is scaled linearly with the samerate. We also describe how the randomized linear program can be used to obtain a bid price controlpolicy. Computational experiments indicate that our approach is quite fast, able to scale to industrialproblems and can provide significant improvements over standard benchmarks.

Fuzzy coding in constrained ordinations

Canonical correspondence analysis and redundancy analysis are two methods of constrained ordination regularly used in the analysis of ecological data when several response variables (for example, species abundances) are related linearly to several explanatory variables (for example, environmental variables, spatial positions of samples). In this report I demonstrate the advantages of the fuzzy coding of explanatory variables: first, nonlinear relationships can be diagnosed; second, more variance in the responses can be explained; and third, in the presence of categorical explanatory variables (for example, years, regions) the interpretation of the resulting triplot ordination is unified because all explanatory variables are measured at a categorical level.

Phase-based manganese enhanced MRI, a new methodology to enhance brain cytoarchitectural contrast and study manganese uptake.

PURPOSE: As the magnetic susceptibility induced frequency shift increases linearly with magnetic field strength, the present work evaluates manganese as a phase imaging contrast agent and investigates the dose dependence of brain enhancement in comparison to T1 -weighted imaging after intravenous administration of MnCl2 . METHODS: Experiments were carried out on 12 Sprague-Dawley rats. MnCl2 was infused intravenously with the following doses: 25, 75, 125 mg/kg (n=4). Phase, T1 -weighted images and T1 maps were acquired before and 24h post MnCl2 administration at 14.1 Tesla. RESULTS: Manganese enhancement was manifested in phase imaging by an increase in frequency shift differences between regions rich in calcium gated channels and other tissues, together with local increase in signal to noise ratio (from the T1 reduction). Such contrast improvement allowed a better visualization of brain cytoarchitecture. The measured T1 decrease observed across different manganese doses and in different brain regions were consistent with the increase in the contrast to noise ratio (CNR) measured by both T1 -weighted and phase imaging, with the strongest variations being observed in the dentate gyrus and olfactory bulb. CONCLUSION: Overall from its high sensitivity to manganese combined with excellent CNR, phase imaging is a promising alternative imaging protocol to assess manganese enhanced MRI at ultra high field. Magn Reson Med 72:1246-1256, 2014. © 2013 Wiley Periodicals, Inc.

Increase of grain and green matter of corn by liming

Despite the fact that aluminum toxicity to crops is eliminated near soil water pH of 5.5, lime recommendation in many regions aims to increase soil pH up to 6.0 or even higher. For highly buffered soils, high rates of limestone are required to raise the pH from 5.5 to 6.0, resulting in additional, sometimes unnecessary, costs. The objective of this study was to evaluate the effect of soil pH on corn yield in a very acid Hapludox. The experiment was carried out in Lages, Southern Brazil, from 1992 to 1996. The soil had water pH of 4.7, Al3+ of 33 mmol c kg-1, O.M. of 45 g kg-1 and lime requirement to pH 6.0 of 9.0 t ha-1. Dolomitic limestone at rates of 0, 4.5, 9.0, 13.5 and 18.0 t ha-1 (equivalent to pure CaCO3) was incorporated into the soil down to 17 cm depth, in 1992. Liming increased linearly the values of soil pH (from 4.7 to 6.6) and Ca and Mg, eliminated Al3+ with rates of 9.0 t ha-1 or higher, decreased slightly Al-CuCl2, Fe and Cu, and did not affect Zn and Mn. Maximum average corn yield for grain (7.9 t ha-1) and for green matter for silage (GM) (59 t ha-1) was obtained, respectively, at soil pH of 6.0 (12 t ha-1 of limestone) and of 6.1 (14 t ha-1 of limestone); maximum economic efficiency for grain was obtained at pH 5.6 (7.5 t ha-1 of limestone). Maximum yield increments due to liming were 17% for grain and 20% for GM.

An Efficient Nominal Unification Algorithm

Nominal Unification is an extension of first-order unification where terms can contain binders and unification is performed modulo α equivalence. Here we prove that the existence of nominal unifiers can be decided in quadratic time. First, we linearly-reduce nominal unification problems to a sequence of freshness and equalities between atoms, modulo a permutation, using ideas as Paterson and Wegman for first-order unification. Second, we prove that solvability of these reduced problems may be checked in quadràtic time. Finally, we point out how using ideas of Brown and Tarjan for unbalanced merging, we could solve these reduced problems more efficiently

Saturable metabolism of continuous high-dose ifosfamide with mesna and GM-CSF: a pharmacokinetic study in advanced sarcoma patients. Swiss Group for Clinical Cancer Research (SAKK).

BACKGROUND: The aim of this study was to assess the pharmacology, toxicity and activity of high-dose ifosfamide mesna +/- GM-CSF administered by a five-day continuous infusion at a total ifosfamide dose of 12-18 g/m2 in adult patients with advanced sarcomas. PATIENTS AND METHODS: Between January 1991 and October 1992 32 patients with advanced or metastatic sarcoma were entered the study. Twenty-seven patients were pretreated including twenty-three with prior ifosfamide at less than 8 g/m2 total dose/cycle. In 25 patients (27 cycles) extensive pharmacokinetic analyses were performed. RESULTS: The area under the plasma concentration-time curve (AUC) for ifosfamide increased linearly with dose while the AUC's of the metabolites measured in plasma by thin-layer chromatography did not increase with dose, particularly that of the active metabolite isophosphoramide mustard. Furthermore the AUC of the inactive carboxymetabolite did not increase with dose. Interpatient variability of pharmacokinetic parameters was high. Dose-limiting toxicity was myelosuppression at 18 g/m2 total dose with grade 4 neutropenia in five of six patients and grade 4 thrombocytopenia in four of six patients. Therefore the maximum tolerated dose was considered to be 18 g/m2 total dose. There was one CR and eleven PR in twenty-nine evaluable patients (overall response rate 41%). CONCLUSION: Both the activation and inactivation pathways of ifosfamide are non-linear and saturable at high-doses although the pharmacokinetics of the parent drug itself are dose linear. Ifosfamide doses greater than 14-16 g/m2 per cycle appear to result in a relative decrease of the active metabolite isophosphoramide mustard. These data suggest a dose-dependent saturation or even inhibition of ifosfamide metabolism by increasing high dose ifosfamide and suggest the need for further metabolic studies.

A note on quadrics through an algebraic curve

In this note we describe the intersection of all quadric hypersur- faces containing a given linearly normal smooth projective curve of genus n and degree 2n + 1

Oscillatory magnetic tweezers based on ferromagnetic beads and simple coaxial coils

We report the design and validation of simple magnetic tweezers for oscillating ferromagnetic beads in the piconewton and nanometer scales. The system is based on a single pair of coaxial coils operating in two sequential modes: permanent magnetization of the beads through a large and brief pulse of magnetic field and generation of magnetic gradients to produce uniaxial oscillatory forces. By using this two step method, the magnetic moment of the beads remains constant during measurements. Therefore, the applied force can be computed and varies linearly with the driving signal. No feedback control is required to produce well defined force oscillations over a wide bandwidth. The design of the coils was optimized to obtain high magnetic fields (280 mT) and gradients (2 T/m) with high homogeneity (5% variation) within the sample. The magnetic tweezers were implemented in an inverted optical microscope with a videomicroscopy-based multiparticle tracking system. The apparatus was validated with 4.5 ¿m magnetite beads obtaining forces up to ~2 pN and subnanometer resolution. The applicability of the device includes microrheology of biopolymer and cell cytoplasm, molecular mechanics, and mechanotransduction in living cells.

Effect of expiratory flow limitation on respiratory mechanical impedance: a model study

Large phasic variations of respiratory mechanical impedance (Zrs) have been observed during induced expiratory flow limitation (EFL) (M. Vassiliou, R. Peslin, C. Saunier, and C. Duvivier. Eur. Respir. J. 9: 779-786, 1996). To clarify the meaning of Zrs during EFL, we have measured from 5 to 30 Hz the input impedance (Zin) of mechanical analogues of the respiratory system, including flow-limiting elements (FLE) made of easily collapsible rubber tubing. The pressures upstream (Pus) and downstream (Pds) from the FLE were controlled and systematically varied. Maximal flow (Vmax) increased linearly with Pus, was close to the value predicted from wave-speed theory, and was obtained for Pus-Pds of 4-6 hPa. The real part of Zin started increasing abruptly with flow (V) >85%Vmax and either further increased or suddenly decreased in the vicinity of V¿max. The imaginary part of Zin decreased markedly and suddenly above 95%Vmax. Similar variations of Zin during EFL were seen with an analogue that mimicked the changes of airway transmural pressure during breathing. After pressure andV measurements upstream and downstream from the FLE were combined, the latter was analyzed in terms of a serial (Zs) and a shunt (Zp) compartment. Zs was consistent with a large resistance and inertance, and Zp with a mainly elastic element having an elastance close to that of the tube walls. We conclude that Zrs data during EFL mainly reflect the properties of the FLE.