Adaptive diversity in heterogeneous environments for populations regulated by a mixture of soft and hard selection

The stable co-existence of two haploid genotypes or two species is studied in a spatially heterogeneous environment submitted to a mixture of soft selection (within-patch regulation) and hard selection (outside-patch regulation) and where two kinds of resource are available. This is analysed both at an ecological time-scale (short term) and at an evolutionary time-scale (long term). At an ecological scale, we show that co-existence is very unlikely if the two competitors are symmetrical specialists exploiting different resources. In this case, the most favourable conditions are met when the two resources are equally available, a situation that should favour generalists at an evolutionary scale. Alternatively, low within-patch density dependence (soft selection) enhances the co-existence between two slightly different specialists of the most available resource. This results from the opposing forces that are acting in hard and soft regulation modes. In the case of unbalanced accessibility to the two resources, hard selection favours the most specialized genotype, whereas soft selection strongly favours the less specialized one. Our results suggest that competition for different resources may be difficult to demonstrate in the wild even when it is a key factor in the maintenance of adaptive diversity. At an evolutionary scale, a monomorphic invasive evolutionarily stable strategy (ESS) always exists. When a linear trade-off exists between survival in one habitat versus that in another, this ESS lies between an absolute adjustment of survival to niche size (for mainly soft-regulated populations) and absolute survival (specialization) in a single niche (for mainly hard-regulated populations). This suggests that environments in agreement with the assumptions of such models should lead to an absence of adaptive variation in the long term.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Aggregate monotonic stable single-valued solutions for cooperative games

[cat] En aquest treball caracteritzem les solucions puntuals de jocs cooperatius d'utilitat transferible que compleixen selecció del core i monotonia agregada. També mostrem que aquestes dues propietats són compatibles amb la individualitat racional, la propietat del jugador fals i la propietat de simetria. Finalment, caracteritzem les solucions puntuals que compleixen les cinc propietats a l'hora.

Community-wide plasmid gene mobilization and selection.

Plasmids have long been recognized as an important driver of DNA exchange and genetic innovation in prokaryotes. The success of plasmids has been attributed to their independent replication from the host's chromosome and their frequent self-transfer. It is thought that plasmids accumulate, rearrange and distribute nonessential genes, which may provide an advantage for host proliferation under selective conditions. In order to test this hypothesis independently of biases from culture selection, we study the plasmid metagenome from microbial communities in two activated sludge systems, one of which receives mostly household and the other chemical industry wastewater. We find that plasmids from activated sludge microbial communities carry among the largest proportion of unknown gene pools so far detected in metagenomic DNA, confirming their presumed role of DNA innovators. At a system level both plasmid metagenomes were dominated by functions associated with replication and transposition, and contained a wide variety of antibiotic and heavy metal resistances. Plasmid families were very different in the two metagenomes and grouped in deep-branching new families compared with known plasmid replicons. A number of abundant plasmid replicons could be completely assembled directly from the metagenome, providing insight in plasmid composition without culturing bias. Functionally, the two metagenomes strongly differed in several ways, including a greater abundance of genes for carbohydrate metabolism in the industrial and of general defense factors in the household activated sludge plasmid metagenome. This suggests that plasmids not only contribute to the adaptation of single individual prokaryotic species, but of the prokaryotic community as a whole under local selective conditions.

Sex-linked genetic markers in the shrews of the Sorex araneus group : from free gene flow to speciation

SUMMARY: The shrews of the Sorex araneus group are morphologically .very similar, but have undergone a spectacular chromosomal evolution. Altogether, the shrews of this group present a complete array of every possible level of chromosomal and genetic differentiation. In South-Western Europe, four species are recognised: S. antiriorii, S. araneus, S. coronatus and S. granarius, which differ essentially by the amount and the composition of Robertsonian metacentric chromosomés. Additionally, several chromosome races of S. araneus are also present in the same region (i.e. Bretolet, Carlit, Cordon, Jura and Vaud). The objective of this thesis was to examine the genetic relationships between populations, races and /or species of the Sorex araneus group with a special emphasis onsex-specific markers (mtDNA and Y chromosome). We first investigate the evolutionary history of the shrews of the Sorex araneus group distributed in the South-Western Europe. The results of. these analyses confirmed the difficulty to draw a single dichotomic tree within this group. Incongruent mtDNA and Y chromosome phylogenies suggest further that genetic and chromosomal evolution are in this group partially independent processes and that the evolutionary history of the south-western European populations of the S. araneus group can only be understood if we consider secondary contacts between taxa, after their divergence (with genetic exchanges by means of hybridization and / or introgression). Using one male-inherited, one female inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii). Overall there results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps after the last glacial period and came into contact in the Rhône Valley. After the Valais shrew (Sorex antinorii) reached the Swiss Alps, it came into contact with the common shrew (Sorex araneus). When two species come into contact and hybridize, endogenous counter-selection of hybrids is usually first expressed as a reduced fertility or viability in hybrids of the heterogametic sex, a mechanism know as Haldane's rule (Haldane 1922). We first evaluated the extent of introgression for Y chromosome, mtDNA and autosomal markers in a hybrid zone between S. antinoriii and S. araneus. The overall level of genetic and karyotypic differentiation between the two species must be strong .enough to allow the detection asymmetric introgression. Secondly, we compared the levels of gene flow between chromosome common to both species and chromosome differently rearranged in each of them. We detected a significantly stronger genetic structure in rearranged chromosomes. Over a 10-year period, we even observed a decrease of genetic structure for common chromosomes. These results strongly support the role of chromosomal rearrangements in the reproductive barrier between S. araneus and S. anfinorii. Overall, this thesis underlines the need to use different inherited (paternally, maternally and / or biparentally) and chromosomally located (on common vs. on rearranged chromosomes) markers to obtain more accurate pictures of genetic relationships between populations or species. RÉSUMÉ: Les musaraignes du groupe Sorex araneus sont morphologiquement très proches, mais ont connu une spectaculaire évolution chromosomique. Prises dans leur ensemble, les musaraignes de ce groupe présentent tous les nivaux possibles de différenciation génétique et chromosomique. Dans le sud-ouest de l'Europe, quatre espèces appartenant à ce groupe sont présentes : S. antinorii, S. araneus, S. coronatus et S. granarius. Celles-ci diffèrent essentiellement par leur caryotype dont la variabilité est principalement due à des fusions Robertsoniennes. De plus, plusieurs races chromosomiques appartenant à S. araneus sont aussi présentes dans la même région (i.e. les races Bretolet, Carlit, Cordon, Jura et Vaud). L'objectif de cette thèse était d'examiner les relations génétiques entre populations, races et/ou espèces du groupe S. araneus, en utilisant particulièrement des marqueurs liés aux sexes (ADN mitochondrial et Chromosome Y). Nous avons dans un premier temps retracé l'histoire évolutive des musaraignes de ce groupe dans le sud-ouest de l'Europe. Les résultats dé ces analyses confirment qu'il est difficile de tracer un simple arbre dichotomique au sein de ce groupe. Les arbres phylogénétiques obtenus sur l'ADN mitochondrial et le chromosome Y sont incongruents et suggèrent de plus que l'évolution génétique et chromosomique sont des processus indépendants. L'histoire évolutive -des populations de ce groupe ne peut. être comprise qu'en considérant des contacts secondaires entre taxa postérieure à leur divergence et induisant des échanges génétiques par hybridation et/ou introgression. Par la suite, nous avons examiné la structure génétique des populations de la musaraigne du Valais, S. antinorii, en utilisant un marqueur transmis par les mâles, un marqueur transmis par les femelles et huit marqueurs transmis par les 2 sexes. Nos résultats suggèrent que deux lignées génétiquement bien différenciées aient colonisé les Alpes Suisses, après les dernières glaciations et entrent en contact dans là Vallée du Rhône. Après avoir franchi les Alpes Suisses, la musaraigne du Valais est entrée en contact avec là musaraigne commune (S. araneus). Lorsque deux espèces entrent en contact et s'hybrident, la sélection contre les hybrides implique habituellement une baisse de fertilité ou de viabilité des hybrides du sexe hétérogamétique (i.e. les mâles XY chez les mammifères). Ce mécanisme est connu sous le nom de règle de Haldane (Haldane 1922) et implique une plus forte structuration génétique de marqueurs males - spécifiques que des marqueurs femelles spécifiques. Nous avons donc évalué le degré d'introgression des marqueurs situés sur le chromosome Y, sur l'ADN mitochondrial et sur des autosomes dans une zone hybride entre S. araneus et S. antinorii. Le niveau de différenciation chromosomique et génétique entre les 2 espèces doit être suffisamment fort pour ne pas permettre la détection d'une introgression asymétrique entre les sexes. Dans un second temps, nous avons comparé les niveaux de flux de gênes mesurés à l'échelle du chromosome, pour des chromosomes communs aux deux espèces et pour des chromosomes différemment arrangées dans chacune des deux espèces. Nous avons détecté une structure génétique significativement plus forte sur les chromosomes réarrangés et comme la zone hybride a été étudiée à dix années d'intervalle, nous observons même une diminution de la structure génétique pour les chromosomes communs au cours du temps.. Ces résultats soutiennent fortement l'hypothèse d'un rôle des réarrangements chromosomiques dans l'établissement d'une barrière reproductive entre S. araneus et S. antinorii. Ainsi cette thèse souligne l'utilité d'utiliser des marqueurs génétiques avec différents modes de transmission. (par les mâles, par les femelles et/ou par les 2 sexes) ou localisés au niveau du chromosome (chromosomes communs vs chromosomes réarrangés) afin d'obtenir une image plus juste ou du moins plus complète des relations génétiques entre populations ou espèces.

Radiosurgery for brain metastases: specific indications for Gamma Knife in lieu of Linac in a single center using both methods

INTRODUCTION: Radiosurgery (RS) is gaining increasing acceptance in the upfront management of brain metastases (BM). It was initially used in so-called radioresistant metastases (melanoma, renal cell, sarcoma) because it allowed delivering higher dose to the tumor. Now, RS is also used for BM of other cancers. The risk of high incidence of new BM questions the need for associated whole-brain radiotherapy (WBRT). Recent evidence suggests that RS alone allows avoiding cognitive impairment related to WBRT, and the latter should be upheld for salvage therapy. Thus the increase use of RS for single and multiple BM raises new technical challenges for treatment delivery and dosimetry. We present our single institution experience focusing on the criteria that led to patients' selection for RS treatment with Gamma Knife (GK) in lieu of Linac. METHODS: Leksell Gamma Knife Perfexion (Elekta, Sweden) was installed in July 2010. Currently, the Swiss federal health care supports the costs of RS for BM with Linac but not with GK. Therefore, in our center, we always consider first the possibility to use Linac for this indication, and only select patients for GK in specific situations. All cases of BM treated with GK were retrospectively reviewed for criteria yielding to GK indication, clinical information, and treatment data. Further work in progress includes a posteriori dosimetry comparison with our Linac planning system (Brainscan V.5.3, Brainlab, Germany). RESULTS: From July 2010 to March 2012, 20 patients had RS for BM with GK (7 patients with single BM, and 13 with multiple BM). During the same period, 31 had Linac-based RS. Primary tumor was melanoma in 9, lung in 7, renal in 2, and gastrointestinal tract in 2 patients. In single BM, the reason for choosing of GK was the anatomical location close to, or in highly functional areas (1 motor cortex, 1 thalamic, 1 ventricular, 1 mesio-temporal, 3 deep cerebellar close to the brainstem), especially since most of these tumors were intended to be treated with high-dose RS (24 Gy at margin) because of their histology (3 melanomas, 1 renal cell). In multiple BM, the reason for choosing GK in relation with the anatomical location of the lesions was either technical (limitations of Linac movements, especially in lower posterior fossa locations) or closeness of multiple lesions to highly functional areas (typically, multiple posterior fossa BM close to the brainstem), precluding optimal dosimetry with Linac. Again, this was made more critical for multiple BM needing high-dose RS (6 melanoma, 2 hypernephroma). CONCLUSION: Radiosurgery for BM may represent some technical challenge in relation with the anatomical location and multiplicity of the lesions. These considerations may be accentuated for so-called radioresistant BM, when higher dose RS in needed. In our experience, Leksell Gamma Knife Perfexion proves to be useful in addressing these challenges for the treatment of BM.

Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

Choosing among techniques for quantifying single-case intervention effectiveness

If single case experimental designs are to be used to establish guidelines for evidence-based interventions in clinical and educational settings, numerical values that reflect treatment effect sizes are required. The present study compares four recently developed procedures for quantifying the magnitude of intervention effect using data with known characteristics. Monte Carlo methods were used to generate AB designs data with potential confounding variables (serial dependence, linear and curvilinear trend, and heteroscedasticity between phases) and two types of treatment effect (level and slope change). The results suggest that data features are important for choosing the appropriate procedure and, thus, inspecting the graphed data visually is a necessary initial stage. In the presence of serial dependence or a change in data variability, the Nonoverlap of All Pairs (NAP) and the Slope and Level Change (SLC) were the only techniques of the four examined that performed adequately. Introducing a data correction step in NAP renders it unaffected by linear trend, as is also the case for the Percentage of Nonoverlapping Corrected Data and SLC. The performance of these techniques indicates that professionals" judgments concerning treatment effectiveness can be readily complemented by both visual and statistical analyses. A flowchart to guide selection of techniques according to the data characteristics identified by visual inspection is provided.

Single cell analysis reveals similar functional competence of dominant and non-dominant CD8 T-cell clonotypes

The properties of CD8 T-cells requiredfor protection from infectiousdisease and cancer are only partiallycharacterized, and only limited data isavailable regarding T-cell clonotypes.It has been proposed that dominantT-cell clonotypes may have higherprotective potential than their nondominantcounterparts. Our objectiveswere to assess memory andeffector functions, stage of differentiationand clonotype selection of tumor-reactive T lymphocytes followingpeptide vaccination in melanomapatients.We also characterized dominantversus non-dominant clonotypesto further understand the in vivo functionof these T-cells based on theirprevalence. Using a novel single-cellapproach for simultaneous ex vivomolecular and functional analysis, wereport the preferential selection andexpansion of several tumor-specificco-dominant clonotypes of intermediateto high frequencies, irrespectiveof whether native or analog peptidewas used for vaccination. Theseclonotypes made up 40 - 95% of thedifferentiated "effector-like" T-cells,but only 25% of the less-differentiated"effector-memory" cells. Bothsubsets also contained non-dominantT-cell clonotypes, but these were significantlymore frequent in the lessdifferentiatedcells. Thus, cell differentiationwas clonotype-dependent.Surprisingly however, the acquisitionof memory and effector T-cell propertieswas clonotype independent, as wefound similar functional profiles indominant and low/ non-dominantT-cell clonotypes. In contrast to analogpeptide vaccination, native peptidevaccination induced T-cell functionsthat were more comprehensive,with more pronounced effector functionscombined with memory cellproperties. In summary, this study revealsthat T-cell functions are determinedprimarily by the antigen andthe stage of T-cell differentiation, butare similar in dominant and non-dominantclonotypes participating in aCD8 T-cell response. The identifiedclonotypic basis of T-cell responsescontributes to the rational developmentof vaccines.

Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location Europe. Methods We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Selection via simulated individual BLUP based on family genotypic effects in sugarcane

The objective of this work was to propose a new selection strategy for the initial stages of sugarcane improvement, based on the methodology 'simulated individual BLUP (BLUPIS)', which promotes a dynamic allocation of individuals selected in each full-sib family, using BLUP as a base for both the genotypic effects of the referred families and plot effects. The method proposed applies to single full-sib families or those obtained from unbalanced or balanced diallel crosses, half-sib families and self-pollinated families. BLUPIS indicates the number of individuals to be selected within each family, the total number of clones to be advanced, and the number of families to contribute with selected individuals. Correlation between BLUPIS and true BLUP was 0.96, by method validation. Additionally, BLUPIS allows the identification of which replication contains the best individuals of each family.

Divergence and genetic variability among superior rubber tree genotypes

The objective of this work was to estimate the genetic variability and divergence among 22 superior rubber tree (Hevea sp.) genotypes of the IAC 400 series. Univariate and multivariate analyses were performed using eight quantitative traits (descriptors), including yield. In the univariate analyses, the estimated parameters were: genetic and environmental variances; genetic and environmental coefficients of variation; and the variation index. The Mahalanobis generalized distance, the Tocher agglomerative method and canonical variables were used for the multivariate analyses. In the univariate analyses, variability was verified among the genotypes for all the variables evaluated. The Tocher method grouped the genotypes into 11 clusters of dissimilarity. The first four canonical variables explained 87.93% of the cumulative variation. The highest genetic variability was found in rubber yield-related traits, which contributed the most to the genetic divergence. The most divergent pairs of genotypes are suggested for crossbreeding. The genotypes evaluated are suitable for breeding and may be used to continue the IAC rubber tree breeding program.

Assessment of rubber tree panels under crowns resistant to South American leaf blight

The objective of this work was to assess the performance of panel clones under crowns resistant to South American leaf blight (Microcyclus ulei). The experiment was carried out with 18 panel clones crown-budded with Hevea pauciflora x H. guianensis, in a Xanthic Ferralsol (Oxisol) in Manaus, AM, Brazil. The following parameters were evaluated: dry rubber yield, plant nutritional status, and anatomical and physiological characteristics of the latex vessels. In the first three years of evaluation, the panel clones IAN 2878, IAN 2903, CNS AM 7905, CNS AM 7905 P1, and PB 28/59 showed the highest dry rubber yield potential, while the clones IAN 6158, IAN 6590, and IAN 6515 should not be recommended for crown budding. Higher potassium and copper foliar content in panel clones were associated to an increase in dry rubber yield. The simultaneous evaluation of anatomical and physiological characteristics of latex is fundamental for the selection of panel clones in the Amazon region. Crown budding is an efficient technology for South American leaf blight management in endemic regions.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Genetic parameters and estimated genetic gains in young rubber tree progenies

The objective of this work was to assess the genetic parameters and to estimate genetic gains in young rubber tree progenies. The experiments were carried out during three years, in a randomized block design, with six replicates and ten plants per plot, in three representative Hevea crop regions of the state of São Paulo, Brazil. Twenty-two progenies were evaluated, from three to five years old, for rubber yield and annual girth growth. Genetic gain was estimated with the multi-effect index (MEI). Selection by progenies means provided greater estimated genetic gain than selection based on individuals, since heritability values of progeny means were greater than the ones of individual heritability, for both evaluated variables, in all the assessment years. The selection of the three best progenies for rubber yield provided a selection gain of 1.28 g per plant. The genetic gains estimated with MEI using data from early assessments (from 3 to 5-year-old) were generally high for annual girth growth and rubber yield. The high genetic gains for annual girth growth in the first year of assessment indicate that progenies can be selected at the beginning of the breeding program. Population effective size was consistent with the three progenies selected, showing that they were not related and that the population genetic variability is ensured. Early selection with the genetic gains estimated by MEI can be made on rubber tree progenies.