Semigroup analysis of structured populations with distributed states-at-birth arising in mathematical aquaculture

Motivated by the modelling of structured parasite populations in aquaculture we consider a class of physiologically structured population models, where individuals may be recruited into the population at different sizes in general. That is, we consider a size-structured population model with distributed states-at-birth. The mathematical model which describes the evolution of such a population is a first order nonlinear partial integro-differential equation of hyperbolic type. First, we use positive perturbation arguments and utilise results from the spectral theory of semigroups to establish conditions for the existence of a positive equilibrium solution of our model. Then we formulate conditions that guarantee that the linearised system is governed by a positive quasicontraction semigroup on the biologically relevant state space. We also show that the governing linear semigroup is eventually compact, hence growth properties of the semigroup are determined by the spectrum of its generator. In case of a separable fertility function we deduce a characteristic equation and investigate the stability of equilibrium solutions in the general case using positive perturbation arguments.

Intestinal content detection in capsule endoscopy using robust features

This work covers two aspects. First, it generally compares and summarizes the similarities and differences of state of the art feature detector and descriptor and second it presents a novel approach of detecting intestinal content (in particular bubbles) in capsule endoscopy images. Feature detectors and descriptors providing invariance to change of perspective, scale, signal-noise-ratio and lighting conditions are important and interesting topics in current research and the number of possible applications seems to be numberless. After analysing a selection of in the literature presented approaches, this work investigates in their suitability for applications information extraction in capsule endoscopy images. Eventually, a very good performing detector of intestinal content in capsule endoscopy images is presented. A accurate detection of intestinal content is crucial for all kinds of machine learning approaches and other analysis on capsule endoscopy studies because they occlude the field of view of the capsule camera and therefore those frames need to be excluded from analysis. As a so called “byproduct” of this investigation a graphical user interface supported Feature Analysis Tool is presented to execute and compare the discussed feature detectors and descriptor on arbitrary images, with configurable parameters and visualized their output. As well the presented bubble classifier is part of this tool and if a ground truth is available (or can also be generated using this tool) a detailed visualization of the validation result will be performed.

Comment on "Streaming potential dependence on water-content in Fontainebleau sand" by V. Allegre, L. Jouniaux, F. Lehmann and P. Sailhac

Allegre et al. recently presented new experimental data regarding the dependence of the streaming potential coupling coefficient with the saturation of the water phase. Such experiments are important to model the self-potential response associated with the flow of water in the vadose zone and the electroseismic/seismoelectric conversions in unsaturated porous media. However, the approach used to interpret the data is questionable and the conclusions reached by Allegre et al. likely incorrect

Mathematical description of bacterial traveling pulses

The Keller-Segel system has been widely proposed as a model for bacterial waves driven by chemotactic processes. Current experiments on E. coli have shown precise structure of traveling pulses. We present here an alternative mathematical description of traveling pulses at a macroscopic scale. This modeling task is complemented with numerical simulations in accordance with the experimental observations. Our model is derived from an accurate kinetic description of the mesoscopic run-and-tumble process performed by bacteria. This model can account for recent experimental observations with E. coli. Qualitative agreements include the asymmetry of the pulse and transition in the collective behaviour (clustered motion versus dispersion). In addition we can capture quantitatively the main characteristics of the pulse such as the speed and the relative size of tails. This work opens several experimental and theoretical perspectives. Coefficients at the macroscopic level are derived from considerations at the cellular scale. For instance the stiffness of the signal integration process turns out to have a strong effect on collective motion. Furthermore the bottom-up scaling allows to perform preliminary mathematical analysis and write efficient numerical schemes. This model is intended as a predictive tool for the investigation of bacterial collective motion.

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

A Comparative Study of the Organic Acid Content of the Hemolymph of Schistosoma mansoni-Resistant and Susceptible Strains of Biomphalaria glabrata

The freshwater snail Biomphalaria glabrata is an intermediate host of the trematode Schistosoma mansoni. However, some strains of B. glabrata are resistant to successful infection by S. mansoni larvae. The present work examines the profile of organic acids present in S. mansoni-resistant and -susceptible strains of B. glabrata, in order to determine whether the type of organic acid present is related to susceptibility. The organic acids were extracted from the hemolymph of two susceptible B. glabrata strains (PR, Puerto Rico and Ba, Jacobina-Bahia from Brazil), and from the resistant strains 13-16-R1 and 10R2, using solid phase extraction procedures followed by high performance liquid chromatography. The organic acids obtained were analyzed and identified by comparison with known standards. Pyruvate, lactate, succinate, malate, fumarate, acetate, propionate, ß-hydroxybutyrate and acetoacetate were detected in all hemolymph samples. Under standard conditions, the concentration of each of these substances varied among the strains tested and appeared to be specific for each strain. An interesting variation was the low concentration of pyruvate in the hemolymph of PR-snails. Only the concentration of fumarate was consistently different (p£ 0.05) between resistant and susceptible strains

A two-compartment mathematical model of neuroglial metabolism using [1-(11)C] acetate.

The purpose of this study was to develop a two-compartment metabolic model of brain metabolism to assess oxidative metabolism from [1-(11)C] acetate radiotracer experiments, using an approach previously applied in (13)C magnetic resonance spectroscopy (MRS), and compared with an one-tissue compartment model previously used in brain [1-(11)C] acetate studies. Compared with (13)C MRS studies, (11)C radiotracer measurements provide a single uptake curve representing the sum of all labeled metabolites, without chemical differentiation, but with higher temporal resolution. The reliability of the adjusted metabolic fluxes was analyzed with Monte-Carlo simulations using synthetic (11)C uptake curves, based on a typical arterial input function and previously published values of the neuroglial fluxes V(tca)(g), V(x), V(nt), and V(tca)(n) measured in dynamic (13)C MRS experiments. Assuming V(x)(g)=10 × V(tca)(g) and V(x)(n)=V(tca)(n), it was possible to assess the composite glial tricarboxylic acid (TCA) cycle flux V(gt)(g) (V(gt)(g)=V(x)(g) × V(tca)(g)/(V(x)(g)+V(tca)(g))) and the neurotransmission flux V(nt) from (11)C tissue-activity curves obtained within 30 minutes in the rat cortex with a beta-probe after a bolus infusion of [1-(11)C] acetate (n=9), resulting in V(gt)(g)=0.136±0.042 and V(nt)=0.170±0.103 μmol/g per minute (mean±s.d. of the group), in good agreement with (13)C MRS measurements.

Complementarities between mathematical and computational modeling: the case of the repeated Prisoners' Dilemma

We study the properties of the well known Replicator Dynamics when applied to a finitely repeated version of the Prisoners' Dilemma game. We characterize the behavior of such dynamics under strongly simplifying assumptions (i.e. only 3 strategies are available) and show that the basin of attraction of defection shrinks as the number of repetitions increases. After discussing the difficulties involved in trying to relax the 'strongly simplifying assumptions' above, we approach the same model by means of simulations based on genetic algorithms. The resulting simulations describe a behavior of the system very close to the one predicted by the replicator dynamics without imposing any of the assumptions of the mathematical model. Our main conclusion is that mathematical and computational models are good complements for research in social sciences. Indeed, while computational models are extremely useful to extend the scope of the analysis to complex scenarios hard to analyze mathematically, formal models can be useful to verify and to explain the outcomes of computational models.

An approximate mathematical model for solidification of a flowing liquid in a microchannel

A mathematical model is developed to analyse the combined flow and solidification of a liquid in a small pipe or two-dimensional channel. In either case the problem reduces to solving a single equation for the position of the solidification front. Results show that for a large range of flow rates the closure time is approximately constant, and the value depends primarily on the wall temperature and channel width. However, the ice shape at closure will be very different for low and high fluxes. As the flow rate increases the closure time starts to depend on the flow rate until the closure time increases dramatically, subsequently the pipe will never close.

Effect of changes in the deuterium content of drinking water on the hydrogen isotope ratio of urinary steroids in the context of sports drug testing.

The hydrogen isotope ratio (HIR) of body water and, therefore, of all endogenously synthesized compounds in humans, is mainly affected by the HIR of ingested drinking water. As a consequence, the entire organism and all of its synthesized substrates will reflect alterations in the isotope ratio of drinking water, which depends on the duration of exposure. To investigate the effect of this change on endogenous urinary steroids relevant to doping-control analysis the hydrogen isotope composition of potable water was suddenly enriched from -50 to 200 0/00 and maintained at this level for two weeks for two individuals. The steroids under investigation were 5β-pregnane-3α,20α-diol, 5α-androst-16-en-3α-ol, 3α-hydroxy-5α-androstan-17-one (ANDRO), 3α-hydroxy-5β-androstan-17-one (ETIO), 5α-androstane-3α,17β-diol, and 5β-androstane-3α,17β-diol (excreted as glucuronides) and ETIO, ANDRO and 3β-hydroxyandrost-5-en-17-one (excreted as sulfates). The HIR of body water was estimated by determination of the HIR of total native urine, to trace the induced changes. The hydrogen in steroids is partly derived from the total amount of body water and cholesterol-enrichment could be calculated by use of these data. Although the sum of changes in the isotopic composition of body water was 150 0/00, shifts of approximately 30 0/00 were observed for urinary steroids. Parallel enrichment in their HIR was observed for most of the steroids, and none of the differences between the HIR of individual steroids was elevated beyond recently established thresholds. This finding is important to sports drug testing because it supports the intended use of this novel and complementary methodology even in cases where athletes have drunk water of different HIR, a plausible and, presumably, inevitable scenario while traveling.

Uncovering hidden duplicated content in public transcriptomics data.

As part of the development of the database Bgee (a dataBase for Gene Expression Evolution), we annotate and analyse expression data from different types and different sources, notably Affymetrix data from GEO and ArrayExpress, and RNA-Seq data from SRA. During our quality control procedure, we have identified duplicated content in GEO and ArrayExpress, affecting ∼14% of our data: fully or partially duplicated experiments from independent data submissions, Affymetrix chips reused in several experiments, or reused within an experiment. We present here the procedure that we have established to filter such duplicates from Affymetrix data, and our procedure to identify future potential duplicates in RNA-Seq data. Database URL: http://bgee.unil.ch/

Stomach content analyses of Simulium perflavum Roubaud 1906 (Diptera: Simuliidae) larvae from streams in Central Amazônia, Brazil

Stomach contents of Simulium perflavum Roubaud larvae were analyzed and compared with plankton and periphyton collected in five streams, in Central Amazonia (Manaus and Presidente Figueiredo counties), in Sep./Oct.1996 (dry season) and Feb./Mar. 1997 (rainy season). A total of 1,400 last-instar larvae were dissected; the stomach contents were analyzed using different methods: fresh and after oxidation. A total of 87 taxa (algae, diatoms and rotifers) were found in the stomach contents. In each stream, qualitative samples of plankton and periphyton were collected; these were mounted between slides and cover slips. A total of 94 taxa of plankton and 54 taxa of periphyton were collected. One species of Rotifera was present in the stomach contents, plankton and periphyton. Cluster analysis based on species composition of the organisms present in the stomach contents grouped the streams into two major groups, each belonging to a different drainage area. Correlations based on presence/absence of species of microalgae in the stomach contents, plankton and periphyton indicated significant associations (p<0.05) between stomach contents and plankton and between plankton and periphyton (z test); the Sorensen coefficient and cluster analysis corroborate the same associations.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.