SIZE AS A LINE OF LEAST RESISTANCE II: DIRECT SELECTION ON SIZE OR CORRELATED RESPONSE DUE TO CONSTRAINTS?

Evolutionary change in New World Monkey (NWM) skulls occurred primarily along the line of least resistance defined by size (including allometric) variation (g(max)). Although the direction of evolution was aligned with this axis, it was not clear whether this macroevolutionary pattern results from the conservation of within population genetic covariance patterns (long-term constraint) or long-term selection along a size dimension, or whether both, constraints and selection, were inextricably involved. Furthermore, G-matrix stability can also be a consequence of selection, which implies that both, constraints embodied in g(max) and evolutionary changes observed on the trait averages, would be influenced by selection Here, we describe a combination of approaches that allows one to test whether any particular instance of size evolution is a correlated by-product due to constraints (g(max)) or is due to direct selection on size and apply it to NWM lineages as a case study. The approach is based on comparing the direction and amount of evolutionary change produced by two different simulated sets of net-selection gradients (beta), a size (isometric and allometric size) and a nonsize set. Using this approach it is possible to distinguish between the two hypotheses (indirect size evolution due to constraints or direct selection on size), because although both may produce an evolutionary response aligned with g(max), the amount of change produced by random selection operating through the variance/covariance patterns (constraints hypothesis) will be much smaller than that produced by selection on size (selection hypothesis). Furthermore, the alignment of simulated evolutionary changes with g(max) when selection is not on size is not as tight as when selection is actually on size, allowing a statistical test of whether a particular observed case of evolution along the line of least resistance is the result of selection along it or not. Also, with matrix diagonalization (principal components [PC]) it is possible to calculate directly the net-selection gradient on size alone (first PC [PC1]) by dividing the amount of phenotypic difference between any two populations by the amount of variation in PC1, which allows one to benchmark whether selection was on size or not

Variety matters: adaptive genetic diversity and parasite load in two mouse opossums from the Brazilian Atlantic forest

The adaptive potential of a species to a changing environment and in disease defence is primarily based on genetic variation. Immune genes, such as genes of the major histocompatibility complex (MHC), may thereby be of particular importance. In marsupials, however, there is very little knowledge about natural levels and functional importance of MHC polymorphism, despite their key role in the mammalian evolution. In a previous study, we discovered remarkable differences in the MHC class II diversity between two species of mouse opossums (Gracilinanus microtarsus, Marmosops incanus) from the Brazilian Atlantic forest, which is one of the most endangered hotspots for biodiversity conservation. Since the main forces in generating MHC diversity are assumed to be pathogens, we investigated in this study gastrointestinal parasite burden and functional associations between the individual MHC constitution and parasite load. We tested two contrasting scenarios, which might explain differences in MHC diversity between species. We predicted that a species with low MHC diversity would either be under relaxed selection pressure by low parasite diversity (`Evolutionary equilibrium` scenario), or there was a recent loss in MHC diversity leading to a lack of resistance alleles and increased parasite burden (`Unbalanced situation` scenario). In both species it became apparent that the MHC class II is functionally important in defence against gastrointestinal helminths, which was shown here for the first time in marsupials. On the population level, parasite diversity did not markedly differ between the two host species. However, we did observe considerable differences in the individual parasite load (parasite prevalence and infection intensity): while M. incanus revealed low MHC DAB diversity and high parasite load, G. microtarsus showed a tenfold higher population wide MHC DAB diversity and lower parasite burden. These results support the second scenario of an unbalanced situation.

An Information Theory framework for two-stage binary image operator design

The design of translation invariant and locally defined binary image operators over large windows is made difficult by decreased statistical precision and increased training time. We present a complete framework for the application of stacked design, a recently proposed technique to create two-stage operators that circumvents that difficulty. We propose a novel algorithm, based on Information Theory, to find groups of pixels that should be used together to predict the Output Value. We employ this algorithm to automate the process of creating a set of first-level operators that are later combined in a global operator. We also propose a principled way to guide this combination, by using feature selection and model comparison. Experimental results Show that the proposed framework leads to better results than single stage design. (C) 2009 Elsevier B.V. All rights reserved.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.