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.

Comparison of selective genotyping strategies for prediction of breeding values in a population undergoing selection

Genomewide marker information can improve the reliability of breeding value predictions for young selection candidates in genomic selection. However, the cost of genotyping limits its use to elite animals, and how such selective genotyping affects predictive ability of genomic selection models is an open question. We performed a simulation study to evaluate the quality of breeding value predictions for selection candidates based on different selective genotyping strategies in a population undergoing selection. The genome consisted of 10 chromosomes of 100 cM each. After 5,000 generations of random mating with a population size of 100 (50 males and 50 females), generation G(0) (reference population) was produced via a full factorial mating between the 50 males and 50 females from generation 5,000. Different levels of selection intensities (animals with the largest yield deviation value) in G(0) or random sampling (no selection) were used to produce offspring of G(0) generation (G(1)). Five genotyping strategies were used to choose 500 animals in G(0) to be genotyped: 1) Random: randomly selected animals, 2) Top: animals with largest yield deviation values, 3) Bottom: animals with lowest yield deviations values, 4) Extreme: animals with the 250 largest and the 250 lowest yield deviations values, and 5) Less Related: less genetically related animals. The number of individuals in G(0) and G(1) was fixed at 2,500 each, and different levels of heritability were considered (0.10, 0.25, and 0.50). Additionally, all 5 selective genotyping strategies (Random, Top, Bottom, Extreme, and Less Related) were applied to an indicator trait in generation G(0), and the results were evaluated for the target trait in generation G(1), with the genetic correlation between the 2 traits set to 0.50. The 5 genotyping strategies applied to individuals in G(0) (reference population) were compared in terms of their ability to predict the genetic values of the animals in G(1) (selection candidates). Lower correlations between genomic-based estimates of breeding values (GEBV) and true breeding values (TBV) were obtained when using the Bottom strategy. For Random, Extreme, and Less Related strategies, the correlation between GEBV and TBV became slightly larger as selection intensity decreased and was largest when no selection occurred. These 3 strategies were better than the Top approach. In addition, the Extreme, Random, and Less Related strategies had smaller predictive mean squared errors (PMSE) followed by the Top and Bottom methods. Overall, the Extreme genotyping strategy led to the best predictive ability of breeding values, indicating that animals with extreme yield deviations values in a reference population are the most informative when training genomic selection models.

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Strong selection against hybrids at a hybrid zone in the Ensatina ring species complex and its evolutionary implications

The analysis of interactions between lineages at varying levels of genetic divergence can provide insights into the process of speciation through the accumulation of incompatible mutations. Ring species, and especially the Ensatina eschscholtzii system exemplify this approach. The plethodontid salamanders E. eschscholtzii xanthoptica and E. eschscholtzii platensis hybridize in the central Sierran foothills of California. We compared the genetic structure across two transects (southern and northern Calaveras Co.), one of which was resampled over 20 years, and examined diagnostic molecular markers (eight allozyme loci and mitochondrial DNA) and a diagnostic quantitative trait (color pattern). Key results across all studies were: (1) cline centers for all markers were coincident and the zones were narrow, with width estimates of 730 m to 2000 m; (2) cline centers at the northern Calaveras transect were coincident between 1981 and 2001, demonstrating repeatability over five generations; (3) there were very few if any putative F1s, but a relatively high number of backcrossed individuals in the central portion of transects: and (4) we found substantial linkage disequilibrium in all three studies and strong heterozygote deficit both in northern Calaveras, in 2001, and southern Calaveras. Both linkage disequilibrium and heterozygote deficit showed maximum values near the center of the zones. Using estimates of cline width and dispersal, we infer strong selection against hybrids. This is sufficient to promote accumulation of differences at loci that are neutral or under divergent selection, but would still allow for introgression of adaptive alleles. The evidence for strong but incomplete isolation across this centrally located contact is consistent with theory suggesting a gradual increase in postzygotic incompatibility between allopatric populations subject to divergent selection and reinforces the value of Ensatina as a system for the study of divergence and speciation at multiple stages. © 2005 The Society for the Study of Evolution. All rights reserved.

Positive selection results in frequent reversible amino acid replacements in the G protein gene of human respiratory syncytial virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a flipflop phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dengue Virus Type 3 Adaptive Changes during Epidemics in Sao Jose de Rio Preto, Brazil, 2006-2007

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Worldwide genetic variation at the 3 ' untranslated region of the HLA-G gene: balancing selection influencing genetic diversity

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Kin-bias, breeding site selection and female fitness in a cannibalistic Neotropical frog

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fast fault section estimation in distribution control centers using adaptive genetic algorithm

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evolution of extrafloral nectaries: adaptive process and selective regime changes from forest to savanna

Much effort has been devoted to understanding the function of extrafloral nectaries (EFNs) for antplantherbivore interactions. However, the pattern of evolution of such structures throughout the history of plant lineages remains unexplored. In this study, we used empirical knowledge on plant defences mediated by ants as a theoretical framework to test specific hypotheses about the adaptive role of EFNs during plant evolution. Emphasis was given to different processes (neutral or adaptive) and factors (habitat change and trade-offs with new trichomes) that may have affected the evolution of antplant associations. We measured seven EFN quantitative traits in all 105 species included in a well-supported phylogeny of the tribe Bignonieae (Bignoniaceae) and collected field data on antEFN interactions in 32 species. We identified a positive association between ant visitation (a surrogate of ant guarding) and the abundance of EFNs in vegetative plant parts and rejected the hypothesis of phylogenetic conservatism of EFNs, with most traits presenting K-values < 1. Modelling the evolution of EFN traits using maximum likelihood approaches further suggested adaptive evolution, with static-optimum models showing a better fit than purely drift models. In addition, the abundance of EFNs was associated with habitat shifts (with a decrease in the abundance of EFNs from forest to savannas), and a potential trade-off was detected between the abundance of EFNs and estipitate glandular trichomes (i.e. trichomes with sticky secretion). These evolutionary associations suggest divergent selection between species as well as explains K-values < 1. Experimental studies with multiple lineages of forest and savanna taxa may improve our understanding of the role of nectaries in plants. Overall, our results suggest that the evolution of EFNs was likely associated with the adaptive process which probably played an important role in the diversification of this plant group.

Relaxation of Adaptive Evolution during the HIV-1 Infection Owing to Reduction of CD4+T Cell Counts

Background: The first stages of HIV-1 infection are essential to establish the diversity of virus population within host. It has been suggested that adaptation to host cells and antibody evasion are the leading forces driving HIV evolution at the initial stages of AIDS infection. In order to gain more insights on adaptive HIV-1 evolution, the genetic diversity was evaluated during the infection time in individuals contaminated by the same viral source in an epidemic cluster. Multiple sequences of V3 loop region of the HIV-1 were serially sampled from four individuals: comprising a single blood donor, two blood recipients, and another sexually infected by one of the blood recipients. The diversity of the viral population within each host was analyzed independently in distinct time points during HIV-1 infection. Results: Phylogenetic analysis identified multiple HIV-1 variants transmitted through blood transfusion but the establishing of new infections was initiated by a limited number of viruses. Positive selection (d(N)/d(S)>1) was detected in the viruses within each host in all time points. In the intra-host viruses of the blood donor and of one blood recipient, X4 variants appeared respectively in 1993 and 1989. In both patients X4 variants never reached high frequencies during infection time. The recipient, who X4 variants appeared, developed AIDS but kept narrow and constant immune response against HIV-1 during the infection time. Conclusion: Slowing rates of adaptive evolution and increasing diversity in HIV-1 are consequences of the CD4+ T cells depletion. The dynamic of R5 to X4 shift is not associated with the initial amplitude of humoral immune response or intensity of positive selection.