Nonlinear set-membership estimation: a support vector machine approach

In this paper a support vector machine (SVM) approach for characterizing the feasible parameter set (FPS) in non-linear set-membership estimation problems is presented. It iteratively solves a regression problem from which an approximation of the boundary of the FPS can be determined. To guarantee convergence to the boundary the procedure includes a no-derivative line search and for an appropriate coverage of points on the FPS boundary it is suggested to start with a sequential box pavement procedure. The SVM approach is illustrated on a simple sine and exponential model with two parameters and an agro-forestry simulation model.

Identification of components of the Sigma B Regulon in Listeria monocytogenes that contribute to acid and salt tolerance

Sigma B (σB) is an alternative sigma factor that controls the transcriptional response to stress in Listeria monocytogenes and is also known to play a role in the virulence of this human pathogen. In the present study we investigated the impact of a sigB deletion on the proteome of L. monocytogenes grown in a chemically defined medium both in the presence and in the absence of osmotic stress (0.5 M NaCl). Two new phenotypes associated with the sigB deletion were identified using this medium. (i) Unexpectedly, the strain with the ΔsigB deletion was found to grow faster than the parent strain in the growth medium, but only when 0.5 M NaCl was present. This phenomenon was independent of the carbon source provided in the medium. (ii) The ΔsigB mutant was found to have unusual Gram staining properties compared to the parent, suggesting that σB contributes to the maintenance of an intact cell wall. A proteomic analysis was performed by two-dimensional gel electrophoresis, using cells growing in the exponential and stationary phases. Overall, 11 proteins were found to be differentially expressed in the wild type and the ΔsigB mutant; 10 of these proteins were expressed at lower levels in the mutant, and 1 was overexpressed in the mutant. All 11 proteins were identified by tandem mass spectrometry, and putative functions were assigned based on homology to proteins from other bacteria. Five proteins had putative functions related to carbon utilization (Lmo0539, Lmo0783, Lmo0913, Lmo1830, and Lmo2696), while three proteins were similar to proteins whose functions are unknown but that are known to be stress inducible (Lmo0796, Lmo2391, and Lmo2748). To gain further insight into the role of σB in L. monocytogenes, we deleted the genes encoding four of the proteins, lmo0796, lmo0913, lmo2391, and lmo2748. Phenotypic characterization of the mutants revealed that Lmo2748 plays a role in osmotolerance, while Lmo0796, Lmo0913, and Lmo2391 were all implicated in acid stress tolerance to various degrees. Invasion assays performed with Caco-2 cells indicated that none of the four genes was required for mammalian cell invasion. Microscopic analysis suggested that loss of Lmo2748 might contribute to the cell wall defect observed in the ΔsigB mutant. Overall, this study highlighted two new phenotypes associated with the loss of σB. It also demonstrated clear roles for σB in both osmotic and low-pH stress tolerance and identified specific components of the σB regulon that contribute to the responses observed.

Estimating the volume of tephra deposits: a new simple strategy

Volume determination of tephra deposits is necessary for the assessment of the dynamics and hazards of explosive volcanoes. Several methods have been proposed during the past 40 years that include the analysis of crystal concentration of large pumices, integrations of various thinning relationships, and the inversion of field observations using analytical and computational models. Regardless of their strong dependence on tephra-deposit exposure and distribution of isomass/isopach contours, empirical integrations of deposit thinning trends still represent the most widely adopted strategy due to their practical and fast application. The most recent methods involve the best fitting of thinning data using various exponential seg- ments or a power-law curve on semilog plots of thickness (or mass/area) versus square root of isopach area. The exponential method is mainly sensitive to the number and the choice of straight segments, whereas the power-law method can better reproduce the natural thinning of tephra deposits but is strongly sensitive to the proximal or distal extreme of integration. We analyze a large data set of tephra deposits and propose a new empirical method for the deter- mination of tephra-deposit volumes that is based on the integration of the Weibull function. The new method shows a better agreement with observed data, reconciling the debate on the use of the exponential versus power-law method. In fact, the Weibull best fitting only depends on three free parameters, can well reproduce the gradual thinning of tephra deposits, and does not depend on the choice of arbitrary segments or of arbitrary extremes of integration.

Likelihood-free estimation of model evidence

Statistical methods of inference typically require the likelihood function to be computable in a reasonable amount of time. The class of “likelihood-free” methods termed Approximate Bayesian Computation (ABC) is able to eliminate this requirement, replacing the evaluation of the likelihood with simulation from it. Likelihood-free methods have gained in efficiency and popularity in the past few years, following their integration with Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC) in order to better explore the parameter space. They have been applied primarily to estimating the parameters of a given model, but can also be used to compare models. Here we present novel likelihood-free approaches to model comparison, based upon the independent estimation of the evidence of each model under study. Key advantages of these approaches over previous techniques are that they allow the exploitation of MCMC or SMC algorithms for exploring the parameter space, and that they do not require a sampler able to mix between models. We validate the proposed methods using a simple exponential family problem before providing a realistic problem from human population genetics: the comparison of different demographic models based upon genetic data from the Y chromosome.

Bayesian parameter estimation for latent Markov random fields and social networks

Undirected graphical models are widely used in statistics, physics and machine vision. However Bayesian parameter estimation for undirected models is extremely challenging, since evaluation of the posterior typically involves the calculation of an intractable normalising constant. This problem has received much attention, but very little of this has focussed on the important practical case where the data consists of noisy or incomplete observations of the underlying hidden structure. This paper specifically addresses this problem, comparing two alternative methodologies. In the first of these approaches particle Markov chain Monte Carlo (Andrieu et al., 2010) is used to efficiently explore the parameter space, combined with the exchange algorithm (Murray et al., 2006) for avoiding the calculation of the intractable normalising constant (a proof showing that this combination targets the correct distribution in found in a supplementary appendix online). This approach is compared with approximate Bayesian computation (Pritchard et al., 1999). Applications to estimating the parameters of Ising models and exponential random graphs from noisy data are presented. Each algorithm used in the paper targets an approximation to the true posterior due to the use of MCMC to simulate from the latent graphical model, in lieu of being able to do this exactly in general. The supplementary appendix also describes the nature of the resulting approximation.

Characterization of a Listeria monocytogenes Scott A isolate with high tolerance towards high hydrostatic pressure

An isolate of L. monocytogenes Scott A that is tolerant to high hydrostatic pressure (HHP), named AK01, was isolated upon a single pressurization treatment of 400 MPa for 20 min and was further characterized. The survival of exponential- and stationary-phase cells of AK01 in ACES [N-(2-acetamido)-2-aminoethanesulfonic acid] buffer was at least 2 log units higher than that of the wild type over a broad range of pressures (150 to 500 MPa), while both strains showed higher HHP tolerance (piezotolerance) in the stationary than in the exponential phase of growth. In semiskim milk, exponential-phase cells of both strains showed lower reductions upon pressurization than in buffer, but again, AK01 was more piezotolerant than the wild type. The piezotolerance of AK01 was retained for at least 40 generations in rich medium, suggesting a stable phenotype. Interestingly, cells of AK01 lacked flagella, were elongated, and showed slightly lower maximum specific growth rates than the wild type at 8, 22, and 30°C. Moreover, the piezotolerant strain AK01 showed increased resistance to heat, acid, and H2O2 compared with the wild type. The difference in HHP tolerance between the piezotolerant strain and the wild-type strain could not be attributed to differences in membrane fluidity, since strain AK01 and the wild type had identical in situ lipid melting curves as determined by Fourier transform infrared spectroscopy. The demonstrated occurrence of a piezotolerant isolate of L. monocytogenes underscores the need to further investigate the mechanisms underlying HHP resistance of food-borne microorganisms, which in turn will contribute to the appropriate design of safe, accurate, and feasible HHP treatments.

Combined action of S-carvone and mild heat treatment on Listeria monocytogenes Scott A

The combined action of the plant-derived volatile, S-carvone, and mild heat treatment on the food-borne pathogen, Listeria monocytogenes, was evaluated. The viability of exponential phase cultures grown at 8 °C could be reduced by 1·3 log units after exposure to S-carvone (5 mmol l−1) for 30 min at 45 °C, while individual treatment with S-carvone or exposure to 45 °C for 30 min did not result in a loss in viability. Other plant-derived volatiles, namely carvacrol, cinnamaldehyde, thymol and decanal, were also found to reduce the viability of L. monocytogenes in combination with the same mild heat treatment at concentrations of 1·75 mmol l−1, 2·5 mmol l−1, 1·5 mmol l−1 and 2 mmol l−1, respectively. These findings show that essential oil compounds can play an important role in minimally processed foods, and can be used in the concept of Hurdle Technology to reduce the intensity of heat treatment or other individual hurdles.

On the time-dependent analysis of Gamow decay

Gamow's explanation of the exponential decay law uses complex 'eigenvalues' and exponentially growing 'eigenfunctions'. This raises the question, how Gamow's description fits into the quantum mechanical description of nature, which is based on real eigenvalues and square integrable wavefunctions. Observing that the time evolution of any wavefunction is given by its expansion in generalized eigenfunctions, we shall answer this question in the most straightforward manner, which at the same time is accessible to graduate students and specialists. Moreover, the presentation can well be used in physics lectures to students.

Heterogeneous surface expression of EspA translocon filaments by Escherichia coli O157:H7 is controlled at the posttranscriptional level

Type III secretion systems of enteric bacteria enable translocation of effector proteins into host cells. Secreted proteins of verotoxigenic Escherichia coli O157 strains include components of a translocation apparatus, EspA, -B, and -D, as well as "effectors" such as the translocated intimin receptor (Tir) and the mitochondrion-associated protein (Map). This research has investigated the regulation of LEE4 translocon proteins, in particular EspA. EspA filaments could not be detected on the bacterial cell surface when E. coli O157:H7 was cultured in M9 minimal medium but were expressed from only a proportion of the bacterial population when cultured in minimal essential medium modified with 25 mM HEPES. The highest proportions of EspA-filamented bacteria were detected in late exponential phase, after which filaments were lost rapidly from the bacterial cell surface. Our previous research had shown that human and bovine E. coli O157:H7 strains exhibit marked differences in EspD secretion levels. Here it is demonstrated that the proportion of the bacterial population expressing EspA filaments was associated with the level of EspD secretion. The ability of individual bacteria to express EspA filaments was not controlled at the level of LEE1-4 operon transcription, as demonstrated by using both beta-galactosidase and green fluorescent protein (GFP) promoter fusions. All bacteria, whether expressing EspA filaments or not, showed equivalent levels of GFP expression when LEEI-4 translational fusions were used. Despite this, the LEE4-espADB mRNA was more abundant from populations with a high proportion of nonsecreting bacteria (low secretors) than from populations with a high proportion of secreting and therefore filamented bacteria (high secretors). This research demonstrates that while specific environmental conditions are required to induce LEEI-4 expression, a further checkpoint exists before EspA filaments are produced on the bacterial surface and secretion of effector proteins occurs. This checkpoint in E. coli O157:H7 translocon expression is controlled by a posttranscriptional mechanism acting on LEE4-espADB mRNA. The heterogeneity in EspA filamentation could arise from phase-variable expression of regulators that control this posttranscriptional mechanism.

Environmental regulation and colonization attributes of the long polar fimbriae (LPF) of Escherichia coli O157:H7

The ability of Escherichia coli O157:H7 to colonize the intestinal epithelia is dependent on the expression of intimin and other adhesins. The chromosome of E. coli O157:H7 carries two loci encoding long polar fimbriae (LPF). These fimbriae mediate adherence to epithelial cells and are associated with colonization of the intestine. In order to increase our knowledge about the conditions controlling their expression and their role in colonization of an animal model, the environmental cues that promote expression of lpf genes and the role of E. coli O157:H7 LPF in intestinal colonization of lambs were investigated. We found that expression of lpf1 was regulated in response to growth phase, osmolarity, and pH; that lpf2 transcription was stimulated during late exponential growth and iron depletion; and that LPF impacts the ability of E. coli O157:H7 to persist in the intestine of infected 6-week-old lambs.

Experimental and modeled UV erthemal irradiance under overcast conditions: the role of cloud optical depth

This paper evaluates the relationship between the cloud modification factor (CMF) in the ultraviolet erythe- mal range and the cloud optical depth (COD) retrieved from the Aerosol Robotic Network (AERONET) "cloud mode" algorithm under overcast cloudy conditions (confirmed with sky images) at Granada, Spain, mainly for non-precipitating, overcast and relatively homogenous water clouds. Empirical CMF showed a clear exponential dependence on experimental COD values, decreasing approximately from 0.7 for COD=10 to 0.25 for COD=50. In addition, these COD measurements were used as input in the LibRadtran radia tive transfer code allowing the simulation of CMF values for the selected overcast cases. The modeled CMF exhibited a dependence on COD similar to the empirical CMF, but modeled values present a strong underestimation with respect to the empirical factors (mean bias of 22 %). To explain this high bias, an exhaustive comparison between modeled and experimental UV erythemal irradiance (UVER) data was performed. The comparison revealed that the radiative transfer simulations were 8 % higher than the observations for clear-sky conditions. The rest of the bias (~14 %) may be attributed to the substantial underestimation of modeled UVER with respect to experimental UVER under overcast conditions, although the correlation between both dataset was high (R2 ~ 0.93). A sensitive test showed that the main reason responsible for that underestimation is the experimental AERONET COD used as input in the simulations, which has been retrieved from zenith radiances in the visible range. In this sense, effective COD in the erythemal interval were derived from an iteration procedure based on searching the best match between modeled and experimental UVER values for each selected overcast case. These effective COD values were smaller than AERONET COD data in about 80 % of the overcast cases with a mean relative difference of 22 %.

Numerical simulation of the Filchner overflow

The plume of Ice Shelf Water (ISW) flowing into the Weddell Sea over the Filchner sill contributes to the formation of Antarctic Bottom Water. The Filchner overflow is simulated using a hydrostatic, primitive equation three-dimensional ocean model with a 0.5–2 Sv ISW influx above the Filchner sill. The best fit to mooring temperature observations is found with influxes of 0.5 and 1 Sv, below a previous estimate of 1.6 ± 0.5 Sv based on sparse mooring velocities. The plume first moves north over the continental shelf, and then turns west, along slope of the continental shelf break where it breaks up into subplumes and domes, some of which then move downslope. Other subplumes run into the eastern submarine ridge and propagate along the ridge downslope in a chaotic manner. The next, western ridge is crossed by the plume through several paths. Despite a number of discrepancies with observational data, the model reproduces many attributes of the flow. In particular, we argue that the temporal variability shown by the observations can largely be attributed to the unstable structure of the flow, where the temperature fluctuations are determined by the motion of the domes past the moorings. Our sensitivity studies show that while thermobaricity plays a role, its effect is small for the flows considered. Smoothing the ridges out demonstrate that their presence strongly affects the plume shape around the ridges. An increase in the bottom drag or viscosity leads to slowing down, and hence thickening and widening of the plume

The maximum rate of mammal evolution

How fast can a mammal evolve from the size of a mouse to the size of an elephant? Achieving such a large transformation calls for major biological reorganization. Thus, the speed at which this occurs has important implications for extensive faunal changes, including adaptive radiations and recovery from mass extinctions. To quantify the pace of large-scale evolution we developed a metric, clade maximum rate, which represents the maximum evolutionary rate of a trait within a clade. We applied this metric to body mass evolution in mammals over the last 70 million years, during which multiple large evolutionary transitions occurred in oceans and on continents and islands. Our computations suggest that it took a minimum of 1.6, 5.1, and 10 million generations for terrestrial mammal mass to increase 100-, and 1,000-, and 5,000- fold, respectively. Values for whales were down to half the length (i.e., 1.1, 3, and 5 million generations), perhaps due to the reduced mechanical constraints of living in an aquatic environment. When differences in generation time are considered, we find an exponential increase in maximum mammal body mass during the 35 million years following the Cretaceous–Paleogene (K–Pg) extinction event. Our results also indicate a basic asymmetry in macroevolution: very large decreases (such as extreme insular dwarfism) can happen at more than 10 times the rate of increases. Our findings allow more rigorous comparisons of microevolutionary and macroevolutionary patterns and processes. Keywords: haldanes, biological time, scaling, pedomorphosis

Sparse shape representation using the Laplace-Beltrami eigenfunctions and its application to modeling subcortical structures

We present a new sparse shape modeling framework on the Laplace-Beltrami (LB) eigenfunctions. Traditionally, the LB-eigenfunctions are used as a basis for intrinsically representing surface shapes by forming a Fourier series expansion. To reduce high frequency noise, only the first few terms are used in the expansion and higher frequency terms are simply thrown away. However, some lower frequency terms may not necessarily contribute significantly in reconstructing the surfaces. Motivated by this idea, we propose to filter out only the significant eigenfunctions by imposing l1-penalty. The new sparse framework can further avoid additional surface-based smoothing often used in the field. The proposed approach is applied in investigating the influence of age (38-79 years) and gender on amygdala and hippocampus shapes in the normal population. In addition, we show how the emotional response is related to the anatomy of the subcortical structures.

Femtosecond pulse shaping using differential evolutionary algorithm and wavelet operators

Evolutionary meta-algorithms for pulse shaping of broadband femtosecond duration laser pulses are proposed. The genetic algorithm searching the evolutionary landscape for desired pulse shapes consists of a population of waveforms (genes), each made from two concatenated vectors, specifying phases and magnitudes, respectively, over a range of frequencies. Frequency domain operators such as mutation, two-point crossover average crossover, polynomial phase mutation, creep and three-point smoothing as well as a time-domain crossover are combined to produce fitter offsprings at each iteration step. The algorithm applies roulette wheel selection; elitists and linear fitness scaling to the gene population. A differential evolution (DE) operator that provides a source of directed mutation and new wavelet operators are proposed. Using properly tuned parameters for DE, the meta-algorithm is used to solve a waveform matching problem. Tuning allows either a greedy directed search near the best known solution or a robust search across the entire parameter space.