Effects of different mating scenarios on embryo viability in brown trout.

Mating with attractive or dominant males is often predicted to offer indirect genetic benefits to females, but it is still largely unclear how important such non-random mating can be with regard to embryo viability. We sampled a natural population of adult migratory brown trout (Salmo trutta), bred them in vitro in a half-sib breeding design to separate genetic from maternal environmental effects, raised 2098 embryos singly until hatching, and exposed them experimentally to different levels of pathogen stress at a late embryonic stage. We found that the embryos' tolerance to the induced pathogen stress was linked to the major histocompatibility complex (MHC) of their parents, i.e. certain MHC genotypes appeared to provide better protection against infection than others. We also found significant additive genetic variance for stress tolerance. Melanin-based dark skin patterns revealed males with 'good genes', i.e. embryos fathered by dark coloured males had a high tolerance to infection. Mating with large and dominant males would, however, not improve embryo viability when compared to random mating. We used simulations to provide estimates of how mate choice based on MHC or melanin-based skin patterns would influence embryos' tolerance to the experimentally induced pathogen stress.

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.

Consolidació de la Plataforma COMalaWEB com a eina Web d'autoaprenentatge per a l'EEES

La plataforma d’e-learning: COMalaWEB és una eina multimedia de suport a l’estudi, l’experimentació i l’adquisició de tècniques d’autoaprenentatge. COMalaWEB s’ofereix com a punt de trobada entre estudiants, professors i altres professionals relacionats amb el mon de les telecomunicacions i/o de la docència universitària. Mitjançant el present projecte s’ha dut a terme la consolidació de la plataforma COMalaWEB com a eina WWW d'autoaprenentatge per a l'EEES (Estudis de Bachelor i de Màster): Dins de la plataforma s’hi ha integrat un Laboratori Virtual per a comunicacions analògiques i digitals (LaViCAD) que ofereix activitats experimentals amb un gran ventall de possibilitats que van des de les demostracions teòriques fins a l’emulació de sistemes de comunicacions quotidians com per exemple la televisió digital o el sistema Wifi dels sistemes WLAN. L’altre gran component de la plataforma es la base de dades de continguts empaquetada en unitats bàsiques anomenades objectes de coneixement, organitzada en cursos, integren tant continguts teòrics com un conjunt d’exercicis proposats per a aprofundir cadascun dels temes tractats. Amb l’actual projecte s’ha treballat en les següents línies d’actuació: - Integració dels simuladors del laboratori LAVICAD a la plataforma d’autoaprenentatge. - Creació de base de dades de recursos docents basats en paquets SCORM per a oferir materials docents de tipus teòric i col·leccions d’exercicis resolts en el marc de les diferents assignatures participants en el projecte. - Creació de base de dades de tipus qüestionari per a oferir exercicis a treballar en el marc de les diferents assignatures participants en el projecte. - Inserció de metodologies docents basades en els anteriors recursos en diferents assignatures d’estudis d’enginyeria i de màster.

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.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Divorce in the barn owl: securing a compatible or better mate entails the cost of re-pairing with a less ornamented female mate.

Two nonmutually exclusive hypotheses can explain why divorce is an adaptive strategy to improve reproductive success. Under the 'better option hypothesis', only one of the two partners initiates divorce to secure a higher-quality partner and increases reproductive success after divorce. Under the 'incompatibility hypothesis', partners are incompatible and hence they may both increase reproductive success after divorce. In a long-term study of the barn owl (Tyto alba), we address the question of whether one or the two partners derive fitness benefits by divorcing. Our results support the hypothesis that divorce is adaptive: after a poor reproductive season, at least one of the two divorcees increase breeding success up to the level of faithful pairs. By breeding more often together, faithful pairs improve coordination and thereby gain in their efficiency to produce successful fledglings. Males would divorce to obtain a compatible mate rather than a mate of higher quality: a heritable melanin-based signal of female quality did not predict divorce (indicating that female absolute quality may not be the cause of divorce), but the new mate of divorced males was less melanic than their previous mate. This suggests that, at least for males, a cost of divorce may be to secure a lower-quality but compatible mate. The better option hypothesis could not be formally rejected, as only one of the two divorcing partners commonly succeeded in obtaining a higher reproductive success after divorce. In conclusion, incompatible partners divorce to restore reproductive success, and by breeding more often together, faithful partners improve coordination.

Creació d'una eina en Java que a partir d'OWL i JESS generi SCORM

Aquest projecte consisteix en generar un programa escrit en llenguatge Java, que utilitzant un motor d'inferència a través d'una llibreria anomenada JESS, pugui llegir un document en format OWL que és una representació de l'ontologia (representació del coneixement) sobre una assignatura, transformant-lo al format de triples que és capaç d'interpretar JESS.

Monte Carlo simulations of metasomatic enrichment in the lithosphere and implications for the source of alkaline basalts

One hypothesis for the origin of alkaline lavas erupted on oceanic islands and in intracontinental settings is that they represent the melts of amphibole-rich veins in the lithosphere (or melts of their dehydrated equivalents if metasomatized lithosphere is recycled into the convecting mantle). Amphibole-rich veins are interpreted as cumulates produced by crystallization of low-degree melts of the underlying asthenosphere as they ascend through the lithosphere. We present the results of trace-element modelling of the formation and melting of veins formed in this way with the goal of testing this hypothesis and for predicting how variability in the formation and subsequent melting of such cumulates (and adjacent cryptically and modally metasomatized lithospheric peridotite) would be manifested in magmas generated by such a process. Because the high-pressure phase equilibria of hydrous near-solidus melts of garnet lherzolite are poorly constrained and given the likely high variability of the hypothesized accumulation and remelting processes, we used Monte Carlo techniques to estimate how uncertainties in the model parameters (e.g. the compositions of the asthenospheric sources, their trace-element contents, and their degree of melting; the modal proportions of crystallizing phases, including accessory phases, as the asthenospheric partial melts ascend and crystallize in the lithosphere; the amount of metasomatism of the peridotitic country rock; the degree of melting of the cumulates and the amount of melt derived from the metasomatized country rock) propagate through the process and manifest themselves as variability in the trace-element contents and radiogenic isotopic ratios of model vein compositions and erupted alkaline magma compositions. We then compare the results of the models with amphibole observed in lithospheric veins and with oceanic and continental alkaline magmas. While the trace-element patterns of the near-solidus peridotite melts, the initial anhydrous cumulate assemblage (clinopyroxene +/- garnet +/- olivine +/- orthopyroxene), and the modelled coexisting liquids do not match the patterns observed in alkaline lavas, our calculations show that with further crystallization and the appearance of amphibole (and accessory minerals such as rutile, ilmenite, apatite, etc.) the calculated cumulate assemblages have trace-element patterns that closely match those observed in the veins and lavas. These calculated hydrous cumulate assemblages are highly enriched in incompatible trace elements and share many similarities with the trace-element patterns of alkaline basalts observed in oceanic or continental setting such as positive Nb/La, negative Ce/Pb, and similiar slopes of the rare earth elements. By varying the proportions of trapped liquid and thus simulating the cryptic and modal metasomatism observed in peridotite that surrounds these veins, we can model the variations in Ba/Nb, Ce/Pb, and Nb/U ratios that are observed in alkaline basalts. If the isotopic compositions of the initial low-degree peridotite melts are similar to the range observed in mid-ocean ridge basalt, our model calculations produce cumulates that would have isotopic compositions similar to those observed in most alkaline ocean island basalt (OIB) and continental magmas after similar to 0 center dot 15 Gyr. However, to produce alkaline basalts with HIMU isotopic compositions requires much longer residence times (i.e. 1-2 Gyr), consistent with subduction and recycling of metasomatized lithosphere through the mantle. such as a heterogeneous asthenosphere. These modelling results support the interpretation proposed by various researchers that amphibole-bearing veins represent cumulates formed during the differentiation of a volatile-bearing low-degree peridotite melt and that these cumulates are significant components of the sources of alkaline OIB and continental magmas. The results of the forward models provide the potential for detailed tests of this class of hypotheses for the origin of alkaline magmas worldwide and for interpreting major and minor aspects of the geochemical variability of these magmas.

El joc del TERCETI : una unitat didàctica multimèdia sobre GNU/LINUX amb suport DVB-MHP i SCORM.

Fent ús del programari lliure s'ha generat la unitat didàctica del joc del terceti per a la televisió digital interactiva i s'ha desada en format SCORM. La facilitat d'ús del comandament a distància fa possible practicar un tipus de t-learning fonamentat en l'edutainment (education+entertaiment) i provoca que el teleespectador passiu passi a ser un "teleusuari" més actiu.

'Good-genes' and 'compatible-genes' effects in an Alpine whitefish and the information content of breeding tubercles over the course of the spawning season.

Some models of sexual selection predict that individuals vary in their genetic quality and reveal some of this variation in their secondary sexual characteristics. Alpine whitefish (Coregonus sp.) develop breeding tubercles shortly before their spawning season. These tubercles are epidermal structures that are distributed regularly along the body sides of both males and females. There is still much unexplained variation in the size of breeding tubercles within both sexes and with much overlap between the sexes. It has been suggested that breeding tubercles function to maintain body contact between the mating partners during spawning, act as weapons for defence of spawning territories, or are sexual signals that reveal aspects of genetic quality. We took two samples of whitefish from their spawning place, one at the beginning and one around the peak of spawning season. We found that females have on average smaller breeding tubercles than males, and that tubercle size partly reveals the stage of gonad maturation. Two independent full-factorial breeding experiments revealed that embryo mortality was significantly influenced by male and female effects. This finding demonstrates that the males differed in their genetic quality (because offspring get nothing but genes from their fathers). Tubercle size was negatively linked to some aspects of embryo mortality in the first breeding experiment but not significantly so in the second. This lack of consistency adds to inconsistent results that were reported before and suggests that (i) some aspects of genetic quality are not revealed in breeding tubercles while others are, or (ii) individuals vary in their signaling strategies and the information content of breeding tubercles is not always reliable. Moreover, the fact that female whitefish have breeding tubercles of significant size while males seem to have few reasons to be choosy suggests that the tubercles might also serve some functions that are not linked to sexual signaling.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.