Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

The role of basalt replenishment in the generation of basaltic andesites of the ongoing activity at Arenal volcano, Costa Rica: evidence from clinopyroxene and spinel

The bulk composition of magma erupted from Volcan Arenal has remained nearly constant (SiO2 = 53.6-54.9 wt%; MgO = 5.0-4.5 wt%) during almost 30 years of continuous activity (1969-1996). None the less, clinopyroxene (cpx) phenocrysts and their spinel inclusions record a much more complex open-system evolution in which steady-state production of the erupted basaltic andesitic magma is linked to episodic injections of basalt into Arenal's magma conduit/reservoir system. High-resolution major element zoning profiles (electron microprobe) on a large number of phenocrysts (>14,000 analyses), tied to back-scattered electron (BSE) images, have been used to assess the compositional characteristics of the magmatic end members as well as the timing and dynamics of magma replenishment events. No two cpx phenocrysts have exactly the same zoning profile. The vast majority of our analyses record the crystallization of cpx (Cr2O3 < 0.12 wt%; Mg# = 65-79; Al/Ti = 2-7) from a liquid comparable to or more evolved than erupted magma compositions. However, half of all cpx grains are cored by high-Cr cpx (Cr2O3 = 0.2-0.72 wt%) or contain similar basaltic compositions as abrupt growth bands in phenocrysts with and without high-Cr cores; phenocrysts with high-Cr cpx occur throughout the ongoing activity. In a few cases, high-Cr cpx occurs very near the outer margin of the grain without an apparent growth hiatus, particularly in 1968/69 and 1992/93. The main conclusions are: (1) all basaltic andesitic lavas erupted at Arenal during the ongoing activity that began in July, 1968, are the products of magma mixing, (2) clinopyroxenes record multiple replenishment events of basaltic magma in contrast to the near constancy of erupted bulk compositions, (3) some phenocrysts preserve records of multiple interactions with basaltic magmas requiring magmatic processes to operate on time-scales shorter than residence times of some phenocrysts, (4) multiple occurrences of clinopyroxene with high-Cr rims suggest that basalt replenishment events have occurred with sub-decadal frequency and may predate eruption by months or less. From this we infer that Arenal volcano is underlain by a continuously active, small-volume magmatic reservoir maintained in quasi-steady state by basalt recharge over several decades. The monotony of erupting Arenal magmas implies that fractionation, recharge, ascent, and eruption are well balanced in order for magmas to be essentially uniform while containing phenocrysts with vastly different growth histories at the time of eruption.

Stochastic time-concentration activity models for cytotoxicity in 3D brain cell cultures.

BACKGROUND: In vitro aggregating brain cell cultures containing all types of brain cells have been shown to be useful for neurotoxicological investigations. The cultures are used for the detection of nervous system-specific effects of compounds by measuring multiple endpoints, including changes in enzyme activities. Concentration-dependent neurotoxicity is determined at several time points. METHODS: A Markov model was set up to describe the dynamics of brain cell populations exposed to potentially neurotoxic compounds. Brain cells were assumed to be either in a healthy or stressed state, with only stressed cells being susceptible to cell death. Cells may have switched between these states or died with concentration-dependent transition rates. Since cell numbers were not directly measurable, intracellular lactate dehydrogenase (LDH) activity was used as a surrogate. Assuming that changes in cell numbers are proportional to changes in intracellular LDH activity, stochastic enzyme activity models were derived. Maximum likelihood and least squares regression techniques were applied for estimation of the transition rates. Likelihood ratio tests were performed to test hypotheses about the transition rates. Simulation studies were used to investigate the performance of the transition rate estimators and to analyze the error rates of the likelihood ratio tests. The stochastic time-concentration activity model was applied to intracellular LDH activity measurements after 7 and 14 days of continuous exposure to propofol. The model describes transitions from healthy to stressed cells and from stressed cells to death. RESULTS: The model predicted that propofol would affect stressed cells more than healthy cells. Increasing propofol concentration from 10 to 100 μM reduced the mean waiting time for transition to the stressed state by 50%, from 14 to 7 days, whereas the mean duration to cellular death reduced more dramatically from 2.7 days to 6.5 hours. CONCLUSION: The proposed stochastic modeling approach can be used to discriminate between different biological hypotheses regarding the effect of a compound on the transition rates. The effects of different compounds on the transition rate estimates can be quantitatively compared. Data can be extrapolated at late measurement time points to investigate whether costs and time-consuming long-term experiments could possibly be eliminated.

Semi-supervised segmentation of ultrasound images based on patch representation and continuous min cut.

Ultrasound segmentation is a challenging problem due to the inherent speckle and some artifacts like shadows, attenuation and signal dropout. Existing methods need to include strong priors like shape priors or analytical intensity models to succeed in the segmentation. However, such priors tend to limit these methods to a specific target or imaging settings, and they are not always applicable to pathological cases. This work introduces a semi-supervised segmentation framework for ultrasound imaging that alleviates the limitation of fully automatic segmentation, that is, it is applicable to any kind of target and imaging settings. Our methodology uses a graph of image patches to represent the ultrasound image and user-assisted initialization with labels, which acts as soft priors. The segmentation problem is formulated as a continuous minimum cut problem and solved with an efficient optimization algorithm. We validate our segmentation framework on clinical ultrasound imaging (prostate, fetus, and tumors of the liver and eye). We obtain high similarity agreement with the ground truth provided by medical expert delineations in all applications (94% DICE values in average) and the proposed algorithm performs favorably with the literature.

THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Defining secretion processes in astrocytes

During the last decade, evidence that release of chemical transmitters from astrocytes might modulate neuronal activity (the so-called "gliotransmission") occurs in situ has been extensively provided. Nevertheless, gliotransmission remains a highly debated topic because of the lack of direct morphological and functional evidence. Here we provided new information supporting gliotransmission, by i) deepen knowledge about specific properties of regulated secretion of glutamatergic SLMVs, and ii) investigating the involvement of astrocytes in the transmission of dopamine, a molecule whose interaction with astrocytes is likely to occur, but it's still not proven.¦VGLUT-expressing glutamatergic SLMVs have been previously identified both in situ and in vitro, but description of kinetics of release were still lacking. To elucidate this issue, we took advantage of fluorescent tools (styryl dyes and pHluorin) and adapted experimental paradigms and analysis methods previously developed to study exo-endocytosis and recycling of glutamatergic vesicles at synapses. Parallel use of EPIfluorescence and total internal reflection (TIRF) imaging allowed us to find that exo-endocytosis processes in astrocytes are extremely fast, with kinetics in the order of milliseconds, able to sustain and follow neuronal signalling at synapses. Also, exocytosis of SLMVs is under the control of fast, localized Ca2+ elevations in close proximity of SLMVs and endoplasmatic reticulum (ER) tubules, the intracellular calcium stores. Such complex organization supports the fast stimulus-secretion coupling we described; localized calcium elevations have been recently observed in astrocytes in situ, suggesting that these functional microdomains might be present in the intact tissue. In the second part of the work, we investigated whether astrocytes possess some of the benchmarks of brain dopaminergic cells. It's been known for years that astrocytes are able to metabolize monoamines by the enzymes MAO and COMT, but to date no clear information that glial cells are able to uptake and store monoamines have been provided. Here, we identified a whole apparatus for the storage, degradation and release of monoamines, at the ultrastructural level. Electron microscopy immunohistochemistry allowed us to visualize VMAT2- and dopamine-positive intracellular compartments within astrocytic processes, i.e. dense -core granules and cisterns. These organelles might be responsible for dopamine release and storage, respectively; interestingly, this intracellular distribution is reminiscent of VMAT2 expression in dendrites if neurons, where dopamine release is tonic and plays a role in the regulation of its a basal levels, suggesting that astrocytic VMAT2 is involved in the homeostasis of dopamine in healthy brains of adult mammals.¦Durant cette dernière décennie, de nombreux résultats sur le relâchement des transmetteurs par les astrocytes pouvant modulé l'activité synaptique (gliotransmission) ont été fournis. Néanmoins, la gliotransmission reste un processus encore très débattu, notamment à cause de l'absence de preuves directes, morphologique et fonctionnelle démontrant ce phénomène. Nous présentons dans nos travaux de nombreux résultats confortant l'hypothèse de la gliotransmission, dont i) une étude approfondie sur les propriétés spatiales et temporelles de la sécrétion régulée du glutamate dans les astrocytes, et ii) une étude sur la participation des astrocytes dans la transmission de la dopamine, une neuromodulateur dont l'interaction avec les astrocytes est fortement probable, mais qui n'a encore jamais été prouvée. L'expression des petites vésicules (SLMVs - Synaptic Like Micro Vesicles) glutamatergiques exprimant les transporteurs vésiculaires du glutamate (VGLUTs) dans les astrocytes a déjà été prouvé tant in situ qu'in vitro. Afin de mettre en évidence les propriétés précises de la sécrétion de ces organelles, nous avons adapté à nos études des méthodes expérimentales conçues pour observer les processus de exocytose et endocytose dans les neurones. Les résolutions spatiale et temporelle obtenues, grâce a l'utilisation en parallèle de l'épi fluorescence et de la fluorescence a onde évanescente (TIRF), nous ont permis de montrer que la sécrétion régulée dans les astrocytes est un processus extrêmement rapide (de l'ordre de la milliseconde) et qu'elle est capable de soutenir et de suivre la transmission de signaux entre neurones. Nous avons également découvert que cette sécrétion a lieu dans des compartiments subcellulaires particuliers où nous observons la présence du reticulum endoplasmique (ER) ainsi que des augmentations rapides de calcium. Cette organisation spatiale complexe pourrait être la base morphologique du couplage rapide entre le stimulus et la sécrétion. Par ailleurs, plusieurs études récentes in vivo semblent confirmer l'existence de ces compartiments. Depuis des années nous savons que les astrocytes sont capables de métaboliser les monoamines par les enzymes MAO et COMT. Nous avons donc fourni de nouvelles preuves concernant la présence d'un appareil de stockage dans les astrocytes participant à la dégradation et la libération de monoamines au niveau ultrastructurelle. Grâce à la microscopie électronique, nous avons découvert la présence de compartiments intracellulaires exprimant VMAT2 dans les processus astrocytaires, sous forme de granules et des citernes. Ces organelles pourraient donc être responsables à la fois du relâchement et du stockage de la dopamine. De manière surprenante, cette distribution intracellulaire est similaire aux dendrites des neurones exprimant VMAT2, où la dopamine est libérée de façon tonique permettant d'agir sur la régulation de ses niveaux de base. Ces résultats, suggèrent une certaine participation des VMAT2 présents dans les astrocytes dans le processus d'homéostase de la dopamine dans le cerveau.¦A de nombreuses reprises, dans des émissions scientifiques ou dans des films, il est avancé que les hommes n'utilisent que 10% du potentiel de leur cerveau. Cette légende provient probablement du fait que les premiers chercheurs ayant décrit les cellules du cerveau entre le XIXème et le XXeme siècle, ont montré que les neurones, les cellules les plus connues et étudiées de cet organe, ne représentent seulement que 10% de la totalité des cellules composant du cerveau. Parmi les 90% restantes, les astrocytes sont sans doute les plus nombreuses. Jusqu'au début des années 90, les astrocytes ont été plutôt considérés peu plus que du tissu conjonctif, ayant comme rôles principaux de maintenir certaines propriétés physiques du cerveau et de fournir un support métabolique (énergie, environnement propre) aux neurones. Grace à la découverte que les astrocytes ont la capacité de relâcher des substances neuro-actives, notamment le glutamate, le rôle des astrocytes dans le fonctionnement cérébral a été récemment reconsidérée.¦Le rôle du glutamate provenant des astrocytes et son impact sur la fonctionnalité des neurones n'a pas encore été totalement élucidé, malgré les nombreuses publications démontrant l'importance de ce phénomène en relation avec différentes fonctions cérébrales. Afin de mieux comprendre comment les astrocytes sont impliqués dans la transmission cérébrale, nous avons étudié les propriétés spatio-temporelles de cette libération grâce à l'utilisation des plusieurs marqueurs fluorescents combinée avec différentes techniques d'imagerie cellulaires. Nous avons découvert que la libération du glutamate par les astrocytes (un processus maintenant appelé "gliotransmission") était très rapide et contrôlée par des augmentations locales de calcium. Nous avons relié ces phénomènes à des domaines fonctionnels subcellulaires morphologiquement adaptés pour ce type de transmission. Plus récemment, nous avons concentré nos études sur un autre transmetteur très important dans le fonctionnement du cerveau : la dopamine. Nos résultats morphologiques semblent indiquer que les astrocytes ont la capacité d'interagir avec ce transmetteur, mais d'une manière différente comparée au glutamate, notamment en terme de rapidité de transmission. Ces résultats suggèrent que le astrocytes ont la capacité de modifier leurs caractéristiques et de s'adapter à leur environnement par rapport aux types de transmetteur avec lequel ils doivent interagir.

Copper, selenium, zinc, and thiamine balances during continuous venovenous hemodiafiltration in critically ill patients.

BACKGROUND: Acute renal failure is a serious complication in critically ill patients and frequently requires renal replacement therapy, which alters trace element and vitamin metabolism. OBJECTIVE: The objective was to study trace element balances during continuous renal replacement therapy (CRRT) in intensive care patients. DESIGN: In a prospective randomized crossover trial, patients with acute renal failure received CRRT with either sodium bicarbonate (Bic) or sodium lactate (Lac) as a buffering agent over 2 consecutive 24-h periods. Copper, selenium, zinc, and thiamine were measured with highly sensitive analytic methods in plasma, replacement solutions, and effluent during 8-h periods. Balances were calculated as the difference between fluids administered and effluent losses and were compared with the recommended intakes (RI) from parenteral nutrition. RESULTS: Nineteen sessions were conducted in 11 patients aged 65 +/- 10 y. Baseline plasma concentrations of copper were normal, whereas those of selenium and zinc were below reference ranges; glutathione peroxidase was in the lower range of normal. The replacement solutions contained no detectable copper, 0.01 micromol Se/L (Bic and Lac), and 1.42 (Bic) and 0.85 (Lac) micromol Zn/L. Micronutrients were detectable in all effluents, and losses were stable in each patient; no significant differences were found between the Bic and Lac groups. The 24-h balances were negative for selenium (-0.97 micromol, or 2 times the daily RI), copper (-6.54 micromol, or 0.3 times the daily RI), and thiamine (-4.12 mg, or 1.5 times the RI) and modestly positive for zinc (20.7 micromol, or 0.2 times the RI). CONCLUSIONS: CRRT results in significant losses and negative balances of selenium, copper, and thiamine, which contribute to low plasma concentrations. Prolonged CRRT is likely to result in selenium and thiamine depletion despite supplementation at recommended amounts.

Paper 3: EUROCAT data quality indicators for population-based registries of congenital anomalies.

The European Surveillance of Congenital Anomalies (EUROCAT) network of population-based congenital anomaly registries is an important source of epidemiologic information on congenital anomalies in Europe covering live births, fetal deaths from 20 weeks gestation, and terminations of pregnancy for fetal anomaly. EUROCAT's policy is to strive for high-quality data, while ensuring consistency and transparency across all member registries. A set of 30 data quality indicators (DQIs) was developed to assess five key elements of data quality: completeness of case ascertainment, accuracy of diagnosis, completeness of information on EUROCAT variables, timeliness of data transmission, and availability of population denominator information. This article describes each of the individual DQIs and presents the output for each registry as well as the EUROCAT (unweighted) average, for 29 full member registries for 2004-2008. This information is also available on the EUROCAT website for previous years. The EUROCAT DQIs allow registries to evaluate their performance in relation to other registries and allows appropriate interpretations to be made of the data collected. The DQIs provide direction for improving data collection and ascertainment, and they allow annual assessment for monitoring continuous improvement. The DQI are constantly reviewed and refined to best document registry procedures and processes regarding data collection, to ensure appropriateness of DQI, and to ensure transparency so that the data collected can make a substantial and useful contribution to epidemiologic research on congenital anomalies.

Continuous monitoring of liver oxygenation with near infrared spectroscopy during naso-gastric tube feeding in neonates.

Near infrared spectroscopy (NIRS) is a non-invasive method of estimating the haemoglobin concentration changes in certain tissues. It is frequently used to monitor oxygenation of the brain in neonates. At present it is not clear whether near infrared spectroscopy of other organs (e.g. the liver as a corresponding site in the splanchnic region, which reacts very sensitively to haemodynamic instability) provides reliable values on their tissue oxygenation. The aim of the study was to test near infrared spectroscopy by measuring known physiologic changes in tissue oxygenation of the liver in newborn infants during and after feeding via a naso-gastric tube. The test-retest variability of such measurements was also determined. On 28 occasions in 25 infants we measured the tissue oxygenation index (TOI) of the liver and the brain continuously before, during and 30 minutes after feeding via a gastric tube. Simultaneously we measured arterial oxygen saturation (SaO2), heart rate (HR) and mean arterial blood pressure (MAP). In 10 other newborn infants we performed a test-retest analysis of the liver tissue oxygenation index to estimate the variability in repeated intra-individual measurements. The tissue oxygenation index of the liver increased significantly from 56.7 +/- 7.5% before to 60.3 +/- 5.6% after feeding (p < 0.005), and remained unchanged for the next 30 minutes. The tissue oxygenation index of the brain (62.1 +/- 9.7%), SaO2 (94.4 +/- 7.1%), heart rate (145 +/- 17.3 min-1) and mean arterial blood pressure (52.8 +/- 10.2 mm Hg) did not change significantly. The test-retest variability for intra-individual measurements was 2.7 +/- 2.1%. After bolus feeding the tissue oxygenation index of the liver increased as expected. This indicates that near infrared spectroscopy is suitable for monitoring changes in tissue oxygenation of the liver in newborn infants.

MAMOT: hidden Markov modeling tool.

Hidden Markov models (HMMs) are probabilistic models that are well adapted to many tasks in bioinformatics, for example, for predicting the occurrence of specific motifs in biological sequences. MAMOT is a command-line program for Unix-like operating systems, including MacOS X, that we developed to allow scientists to apply HMMs more easily in their research. One can define the architecture and initial parameters of the model in a text file and then use MAMOT for parameter optimization on example data, decoding (like predicting motif occurrence in sequences) and the production of stochastic sequences generated according to the probabilistic model. Two examples for which models are provided are coiled-coil domains in protein sequences and protein binding sites in DNA. A wealth of useful features include the use of pseudocounts, state tying and fixing of selected parameters in learning, and the inclusion of prior probabilities in decoding. AVAILABILITY: MAMOT is implemented in C++, and is distributed under the GNU General Public Licence (GPL). The software, documentation, and example model files can be found at http://bcf.isb-sib.ch/mamot