Isolation and characterization of microsatellite loci for the European tree frog (Hyla arborea)

We developed 11 new microsatellite markers for the European tree frog (Hyla arborea), and tested patterns of polymorphism in 54 adults (27 males and 27 females) from two ponds close to Lausanne (Western Switzerland). One marker was sex linked and two pairs displayed linkage disequilibrium. Comparisons of allele numbers with heterozygosity values support a stepwise-mutation model at neutral equilibrium, with mutation rates spanning nearly two orders of magnitude. These markers will prove useful for population genetic studies and fine-scale investigations requiring genetic assignment techniques.

Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

Documenting and preserving the genetic diversity of populations, which conditions their long-term survival, have become a major issue in conservation biology. The loss of diversity often documented in declining populations is usually assumed to result from human disturbances; however, historical biogeographic events, otherwise known to strongly impact diversity, are rarely considered in this context. We apply a multilocus phylogeographic study to investigate the late-Quaternary history of a tree frog (Hyla arborea) with declining populations in the northern and western part of its distribution range. Mitochondrial and nuclear polymorphisms reveal high genetic diversity in the Balkan Peninsula, with a spatial structure moulded by the last glaciations. While two of the main refugial lineages remained limited to the Balkans (Adriatic coast, southern Balkans), a third one expanded to recolonize Northern and Western Europe, loosing much of its diversity in the process. Our findings show that mobile and a priori homogeneous taxa may also display substructure within glacial refugia ('refugia within refugia') and emphasize the importance of the Balkans as a major European biodiversity centre. Moreover, the distribution of diversity roughly coincides with regional conservation situations, consistent with the idea that historically impoverished genetic diversity may interact with anthropogenic disturbances, and increase the vulnerability of populations. Phylogeographic models seem important to fully appreciate the risks of local declines and inform conservation strategies.

Hierarchic species-area relationships and the management of forest habitat islands in intensive farmland

Habitat loss and fragmentation due to land use changes are major threats to biodiversity in forest ecosystems, and they are expected to have important impacts on many taxa and at various spatial scales. Species richness and area relationships (SARs) have been used to assess species diversity patterns and drivers, and thereby in the establishment of conservation and management strategies. Here we propose a hierarchical approach to achieve deeper insights on SARs in small forest islets in intensive farmland and to address the impacts of decreasing naturalness on such relationships. In the intensive dairy landscapes of Northwest Portugal, where small forest stands (dominated by pines, eucalypts or both) represent semi-natural habitat islands, 50 small forest stands were selected and surveyed for vascular plant diversity. A hierarchical analytical framework was devised to determine species richness and inter- and intra-patch SARs for the whole set of forest patches (general patterns) and for each type of forest (specific patterns). Differences in SARs for distinct groups were also tested by considering subsets of species (native, alien, woody, and herbaceous). Overall, values for species richness were confirmed to be different between forest patches exhibiting different levels of naturalness. Whereas higher values of plant diversity were found in pine stands, higher values for alien species were observed in eucalypt stands. Total area of forest (inter-patch SAR) was found not to have a significant impact on species richness for any of the targeted groups of species. However, significant intra-patch SARs were obtained for all groups of species and forest types. A hierarchical approach was successfully applied to scrutinise SARs along a gradient of forest naturalness in intensively managed landscapes. Dominant canopy tree and management intensity were found to reflect differently on distinct species groups as well as to compensate for increasing stand area, buffering SARs among patches, but not within patches. Thus, the maintenance of small semi-natural patches dominated by pines, under extensive practices of forest management, will promote native plant diversity while at the same time contributing to limit the expansion of problematic alien invasive species.

Polymorphism at a sex-linked transcription cofactor in European tree frogs (Hyla arborea): Sex-antagonistic selection or neutral processes?

Nascent sex chromosomes offer a unique opportunity to investigate the evolutionary fate of genesrecently trapped in non-recombining segments. A housekeeping gene (MED15) was recently shown to lie on the nascent sex-chromosomes of the European tree frog (Hyla arborea), with different alleles fixed on the X and the Y chromosomes. Here we document a polymorphism (glutamine deletion) in the X copy of the gene, and use population surveys and experimental crosses to test whether this polymorphism is neutral or maintained by sex-antagonistic selection. Tadpoles from parents of known genotypes revealed significant discrepancies from Mendelian inheritance, suggesting possible sex-antagonistic effects under laboratory conditions. Quantitatively, however, these effects did not meet the conditions for polymorphism maintenance. Furthermore, field estimates of female genotypic frequencies did not differ from Hardy-Weinberg equilibrium and allelic frequencies on the X chromosome did not differ between sexes. In conclusion, although sex antagonistic effects cannot be excluded given the laboratory conditions, the X-linked polymorphism under study appears neutral in the wild. Alternatively, sex-antagonistic selection might still account for the fixation of a male specific allele on the Y chromosome.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Ever-young sex chromosomes in European tree frogs.

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4-7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.

Tree line shifts in the Swiss Alps: Climate change or land abandonment?

Questions: Did the forest area in the Swiss Alps increase between 1985 and 1997? Does the forest expansion near the tree line represent an invasion into abandoned grasslands (ingrowth) or a true upward shift of the local tree line? What land cover / land use classes did primarily regenerate to forest, and what forest structural types did primarily regenerate? And, what are possible drivers of forest regeneration in the tree line ecotone, climate and/or land use change? Location: Swiss Alps. Methods: Forest expansion was quantified using data from the repeated Swiss land use statistics GEOSTAT. A moving window algorithm was developed to distinguish between forest ingrowth and upward shift. To test a possible climate change influence, the resulting upward shifts were compared to a potential regional tree line. Results: A significant increase of forest cover was found between 1650 to and 2450 m. Above 1650 m, 10% of the new forest areas were identified as true upward shifts whereas 90% represented ingrowth, and we identified both land use and climate change as likely drivers. Most upward shift activities were found to occur within a band of 300 m below the potential regional tree line, indicating land use as the most likely driver. Only 4% of the upward shifts were identified to rise above the potential regional tree line, thus indicating climate change. Conclusions: Land abandonment was the most dominant driver for the establishment of new forest areas, even at the tree line ecotone. However, a small fraction of upwards shift can be attributed to the recent climate warming, a fraction that is likely to increase further if climate continues to warm, and with a longer time-span between warming and measurement of forest cover.

Evolutionary history of almond tree domestication in the Mediterranean basin.

Genetic diversity of contemporary domesticated species is shaped by both natural and human-driven processes. However, until now, little is known about how domestication has imprinted the variation of fruit tree species. In this study, we reconstruct the recent evolutionary history of the domesticated almond tree, Prunus dulcis, around the Mediterranean basin, using a combination of nuclear and chloroplast microsatellites [i.e. simple sequence repeat (SSRs)] to investigate patterns of genetic diversity. Whereas conservative chloroplast SSRs show a widespread haplotype and rare locally distributed variants, nuclear SSRs show a pattern of isolation by distance with clines of diversity from the East to the West of the Mediterranean basin, while Bayesian genetic clustering reveals a substantial longitudinal genetic structure. Both kinds of markers thus support a single domestication event, in the eastern side of the Mediterranean basin. In addition, model-based estimation of the timing of genetic divergence among those clusters is estimated sometime during the Holocene, a result that is compatible with human-mediated dispersal of almond tree out of its centre of origin. Still, the detection of region-specific alleles suggests that gene flow from relictual wild preglacial populations (in North Africa) or from wild counterparts (in the Near East) could account for a fraction of the diversity observed.

Cryptic diversity among Western Palearctic tree frogs: postglacial range expansion, range limits, and secondary contacts of three European tree frog lineages (Hyla arborea group).

We characterize divergence times, intraspecific diversity and distributions for recently recognized lineages within the Hyla arborea species group, based on mitochondrial and nuclear sequences from 160 localities spanning its whole distribution. Lineages of H. arborea, H. orientalis, H. molleri have at least Pliocene age, supporting species level divergence. The genetically uniform Iberian H. molleri, although largely isolated by the Pyrenees, is parapatric to H. arborea, with evidence for successful hybridization in a small Aquitanian corridor (southwestern France), where the distribution also overlaps with H. meridionalis. The genetically uniform H. arborea, spread from Crete to Brittany, exhibits molecular signatures of a postglacial range expansion. It meets different mtDNA clades of H. orientalis in NE-Greece, along the Carpathians, and in Poland along the Vistula River (there including hybridization). The East-European H. orientalis is strongly structured genetically. Five geographic mitochondrial clades are recognized, with a molecular signature of postglacial range expansions for the clade that reached the most northern latitudes. Hybridization with H. savignyi is suggested in southwestern Turkey. Thus, cryptic diversity in these Pliocene Hyla lineages covers three extremes: a genetically poor, quasi-Iberian endemic (H. molleri), a more uniform species distributed from the Balkans to Western Europe (H. arborea), and a well-structured Asia Minor-Eastern European species (H. orientalis).

Tree-based message diffusion for managing replicated data in unreliable and resource-constrained peer-to-peer environment

Abstract This thesis proposes a set of adaptive broadcast solutions and an adaptive data replication solution to support the deployment of P2P applications. P2P applications are an emerging type of distributed applications that are running on top of P2P networks. Typical P2P applications are video streaming, file sharing, etc. While interesting because they are fully distributed, P2P applications suffer from several deployment problems, due to the nature of the environment on which they perform. Indeed, defining an application on top of a P2P network often means defining an application where peers contribute resources in exchange for their ability to use the P2P application. For example, in P2P file sharing application, while the user is downloading some file, the P2P application is in parallel serving that file to other users. Such peers could have limited hardware resources, e.g., CPU, bandwidth and memory or the end-user could decide to limit the resources it dedicates to the P2P application a priori. In addition, a P2P network is typically emerged into an unreliable environment, where communication links and processes are subject to message losses and crashes, respectively. To support P2P applications, this thesis proposes a set of services that address some underlying constraints related to the nature of P2P networks. The proposed services include a set of adaptive broadcast solutions and an adaptive data replication solution that can be used as the basis of several P2P applications. Our data replication solution permits to increase availability and to reduce the communication overhead. The broadcast solutions aim, at providing a communication substrate encapsulating one of the key communication paradigms used by P2P applications: broadcast. Our broadcast solutions typically aim at offering reliability and scalability to some upper layer, be it an end-to-end P2P application or another system-level layer, such as a data replication layer. Our contributions are organized in a protocol stack made of three layers. In each layer, we propose a set of adaptive protocols that address specific constraints imposed by the environment. Each protocol is evaluated through a set of simulations. The adaptiveness aspect of our solutions relies on the fact that they take into account the constraints of the underlying system in a proactive manner. To model these constraints, we define an environment approximation algorithm allowing us to obtain an approximated view about the system or part of it. This approximated view includes the topology and the components reliability expressed in probabilistic terms. To adapt to the underlying system constraints, the proposed broadcast solutions route messages through tree overlays permitting to maximize the broadcast reliability. Here, the broadcast reliability is expressed as a function of the selected paths reliability and of the use of available resources. These resources are modeled in terms of quotas of messages translating the receiving and sending capacities at each node. To allow a deployment in a large-scale system, we take into account the available memory at processes by limiting the view they have to maintain about the system. Using this partial view, we propose three scalable broadcast algorithms, which are based on a propagation overlay that tends to the global tree overlay and adapts to some constraints of the underlying system. At a higher level, this thesis also proposes a data replication solution that is adaptive both in terms of replica placement and in terms of request routing. At the routing level, this solution takes the unreliability of the environment into account, in order to maximize reliable delivery of requests. At the replica placement level, the dynamically changing origin and frequency of read/write requests are analyzed, in order to define a set of replica that minimizes communication cost.

Antigen-secretory IgA immune complexes are retro-transported into mouse Peyer's patches: immunotargeting to dendritic cells

RÉSUMÉ Les plaques de Peyer (PP) représentent le site d'entrée majeur des pathogènes au niveau des muqueuses intestinales. Après avoir traversé la cellule M, l'antigène est pris en charge par les cellules dendritiques (DC) de la région sub-épithéliale du dôme des PP. Ces dernières activent une réponse immunitaire qui conduit à la production de l'IgA de sécrétion (SIgA), l'anticorps majeur au niveau muqueux. Des études précédentes dans notre laboratoire ont démontré qu'après administration de SIgA dans des anses intestinales de souris, les SIgA se lient spécifiquement aux cellules M, entrent dans les PP, et sont éventuellement internalisées par les DC. Le but de ce travail est de comprendre la relevance biologique de l'entrée des SIgA dans les PP et leur relevance physiologique dans l'homéostasie mucosale. Dans un premier temps, nous avons montré en utilisant une méthode de purification optimisée basée sur une isolation magnétique, que, en plus des DC myéloïdes (CD11c+/CD11b+) et des DC lymphoïdes (CD11c+/CD8+), les PP de souris contiennent un nouveau sous-type de DC exprimant les marqueurs CD11c et CD19. L'utilisation de la microscopie confocale nous a permis de démontrer que les DC myéloïdes internalisent des SIgA, contrairement aux DC lymphoïdes qui n'interagissent pas avec les SIgA, alors que le nouveau sous-type de DC exprimant CD19 lie les SIgA. En plus, nous avons démontré qu'aucune des DC de rate, de ganglion bronchique ou de ganglion inguinal interagit avec les SIgA. Dans le but d'explorer si les SIgA peuvent délivrer des antigènes aux DC des PP in vivo, nous avons administré des complexes immunitaires formés de Shigella flexneri complexées à des SIgA, dans des anses intestinales de souris. Nous avons observé une entrée dans les PP, suivie d'une migration vers les ganglions mésentériques drainants, contrairement aux Shigella flexneri seules, qui n'infectent pas la souris par la voie intestinale. Shigella flexneri délivrée par SIgA n'induit pas de destruction tissulaire au niveau de l'intestin. En plus de l'exclusion immunitaire, ces résultats suggèrent un nouveau rôle des SIgA, qui consiste à transporter des antigènes à l'intérieur des PP dans un contexte non-inflammatoire. RÉSUMÉ DESTINÉ À UN LARGE PUBLIC L'intestin a pour rôle principal d'absorber les nutriments digérés tout au long du tube digestif, et de les faire passer dans le compartiment intérieur sanguin. Du fait de son exposition chronique avec un monde extérieur constitué d'aliments et de bactéries, l'intestin est un endroit susceptible aux infections et a donc besoin d'empêcher l'entrée de microbes. Pour cela, l'intestin est tapissé de "casernes" appelées les plaques de Peyer, qui appartiennent à un système de défense appelé système immunitaire muqueux. Les plaques de Peyer sont composées de différents types de cellules, ayant pour rôle de contrôler l'entrée de microbes et de développer une réaction immunitaire lors d'infection. Cette réaction immunitaire contre les microbes (antigènes) débute par la prise en charge de l'antigène par des sentinelles, les cellules dendritiques. L'antigène est préparé de façon à être reconnu par d'autres cellules appelées lymphocytes T capables d'activer d'autres cellules, les lymphocytes B. La réaction immunitaire résulte dans la production par les lymphocytes B d'un anticorps spécifique appelé IgA de sécrétion (SIgA) au niveau de la lumière intestinale. De manière classique, le rôle de SIgA au niveau de la lumière intestinale consiste à enrober les microbes et donc exclure leur entrée dans le compartiment intérieur. Dans ce travail, nous avons découvert une nouvelle fonction des SIgA qui consiste à introduire des antigènes dans les plaques de Peyer, et de les diriger vers les cellules dendritiques. Sachant que les SIgA sont des anticorps qui ne déclenchent pas de réactions de défense violentes dites inflammatoires, l'entrée des antigènes via SIgA serait en faveur d'une défense intestinale maîtrisée sans qu'il y ait d'inflammation délétère. Ces résultats nous laissent supposer que l'entrée d'antigènes via SIgA pourrait conduire le système immunitaire muqueux à reconnaître ces antigènes de manière appropriée. Ce mécanisme pourrait expliquer les désordres immunitaires de types allergiques et maladies auto-immunitaires que l'on rencontre chez certaines personnes déficientes en IgA, chez qui cette lecture d'antigènes de manière correcte serait inadéquate. ABSTRACT Peyer's patches (PP) represent the primary site for uptake and presentation of ingested antigens in the intestine. Antigens are sampled by M cells, which pass them to underlying antigen-presenting cells including dendritic cells (DC). This leads to the induction of mucosal T cell response that is important for the production of secretory IgA (SIgA), the chief antibody at mucosal surfaces. Previous studies in the laboratory have shown that exogenous SIgA administrated into mouse intestinal loop binds specifically to M cells, enter into PP, and is eventually internalized by DC. The aim of this work is to understand the biological significance of the SIgA uptake by PP DC and its physiological relevance for mucosal homeostasis. As a first step, we have shown by using an optimized MACS method that, in addition to the CD11c+/CD11b+ (myeloid DC) and CD11c+/CD8+ (lymphoid DC) subtypes, mouse PP contain a novel DC subtype exhibiting both CD11c and CD19 markers. By using a combination of MACS isolation and confocal microscopy, we have demonstrated that in contrast to the lymphoid DC which do not interact with SIgA, the myeloid DC internalize SIgA, while the CD19+ subtype binds SIgA on its surface. Neither spleen DC, nor bronchial-lymph node DC, nor inguinal lymph node DC exhibit such a binding specificity. To test whether SIgA could deliver antigens to PP DC in vivo, we administered SIgA-Shigella flexneri immune complexes into mouse intestinal loop containing a PP. We found that (i) SIgA-Shigella flexneri immune complexes enter the PP and are internalized by sub-epithelial dome PP DC, in contrast to Shigella flexneri alone that does not penetrate the intestinal epithelia in mice, (ii) immune complexes migrate to the draining mesenteric lymph node, (iii) Shigella flexneri carried via SIgA do not induce intestinal tissue destruction. Our results suggest that in addition to immune exclusion, SIgA transports antigens back to the PP under non-inflammatory conditions.

Climatic extremes improve predictions of spatial patterns of tree species.

Understanding niche evolution, dynamics, and the response of species to climate change requires knowledge of the determinants of the environmental niche and species range limits. Mean values of climatic variables are often used in such analyses. In contrast, the increasing frequency of climate extremes suggests the importance of understanding their additional influence on range limits. Here, we assess how measures representing climate extremes (i.e., interannual variability in climate parameters) explain and predict spatial patterns of 11 tree species in Switzerland. We find clear, although comparably small, improvement (+20% in adjusted D(2), +8% and +3% in cross-validated True Skill Statistic and area under the receiver operating characteristics curve values) in models that use measures of extremes in addition to means. The primary effect of including information on climate extremes is a correction of local overprediction and underprediction. Our results demonstrate that measures of climate extremes are important for understanding the climatic limits of tree species and assessing species niche characteristics. The inclusion of climate variability likely will improve models of species range limits under future conditions, where changes in mean climate and increased variability are expected.

Heterozygosity-fitness correlations among wild populations of European tree frogs (Hyla arborea) detect fixation load.

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.

The contribution of patch topology and demographic parameters to PVA predictions: the case of the European tree frog

Population viability analyses (PVA) are increasingly used in metapopulation conservation plans. Two major types of models are commonly used to assess vulnerability and to rank management options: population-based stochastic simulation models (PSM such as RAMAS or VORTEX) and stochastic patch occupancy models (SPOM). While the first set of models relies on explicit intrapatch dynamics and interpatch dispersal to predict population levels in space and time, the latter is based on spatially explicit metapopulation theory where the probability of patch occupation is predicted given the patch area and isolation (patch topology). We applied both approaches to a European tree frog (Hyla arborea) metapopulation in western Switzerland in order to evaluate the concordances of both models and their applications to conservation. Although some quantitative discrepancies appeared in terms of network occupancy and equilibrium population size, the two approaches were largely concordant regarding the ranking of patch values and sensitivities to parameters, which is encouraging given the differences in the underlying paradigms and input data.