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.

Irrigation performance and gross water productivity in furrow-irrigated ornamental tree production

In the ornamental plant production region of Girona (Spain), which is one of the largest of its kind in southern Europe, most of the surface is irrigated using wide blocked-end furrows. The objectives of this paper were: (1) to evaluate the irrigation scheduling methods used by ornamental plant producers; (2) to analyse different scenarios inorder to assess how they affect irrigation performance; (3) to evaluate the risk of deep percolation; and (4) to calculategross water productivity. A two-year study in a representative commercial field, planted with Prunus cerasifera ‘Nigra’, was carried out. The irrigation dose applied by the farmers was slightly smaller than the required water dose estimated by the use of two different methods: the first based on soil water content, and the second based on evapotranspiration. Distribution uniformity and application efficiency were high, with mean values above 87%. Soil water contentmeasurements revealed that even at the end of the furrow, where the infiltrated water depth was greatest, more than 90% of the infiltrated water was retained in the shallowest 40 cm of the soil; accordingly, the risk of water loss due to deep percolation was minimal. Gross water productivity for ornamental tree production was € 11.70 m–3, approximately 20 times higher than that obtained with maize in the same region

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.

Fast docking using the CHARMM force field with EADock DSS.

The prediction of binding modes (BMs) occurring between a small molecule and a target protein of biological interest has become of great importance for drug development. The overwhelming diversity of needs leaves room for docking approaches addressing specific problems. Nowadays, the universe of docking software ranges from fast and user friendly programs to algorithmically flexible and accurate approaches. EADock2 is an example of the latter. Its multiobjective scoring function was designed around the CHARMM22 force field and the FACTS solvation model. However, the major drawback of such a software design lies in its computational cost. EADock dihedral space sampling (DSS) is built on the most efficient features of EADock2, namely its hybrid sampling engine and multiobjective scoring function. Its performance is equivalent to that of EADock2 for drug-like ligands, while the CPU time required has been reduced by several orders of magnitude. This huge improvement was achieved through a combination of several innovative features including an automatic bias of the sampling toward putative binding sites, and a very efficient tree-based DSS algorithm. When the top-scoring prediction is considered, 57% of BMs of a test set of 251 complexes were reproduced within 2 Å RMSD to the crystal structure. Up to 70% were reproduced when considering the five top scoring predictions. The success rate is lower in cross-docking assays but remains comparable with that of the latest version of AutoDock that accounts for the protein flexibility. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011.

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.

Jasmonate controls polypeptide patterning in undamaged tissue in wounded Arabidopsis leaves.

Wounding initiates a strong and largely jasmonate-dependent remodelling of the transcriptome in the leaf blades of Arabidopsis (Arabidopsis thaliana). How much control do jasmonates exert on wound-induced protein repatterning in leaves? Replicated shotgun proteomic analyses of 2.5-mm-wide leaf strips adjacent to wounds revealed 106 differentially regulated proteins. Many of these gene products have not emerged as being wound regulated in transcriptomic studies. From experiments using the jasmonic acid (JA)-deficient allene oxide synthase mutant we estimated that approximately 95% of wound-stimulated changes in protein levels were deregulated in the absence of JA. The levels of two tonoplast proteins already implicated in defense response regulation, TWO-PORE CHANNEL1 and the calcium-V-ATPase ACA4 increased on wounding, but their transcripts were not wound inducible. The data suggest new roles for jasmonate in controlling the levels of calcium-regulated pumps and transporters, proteins involved in targeted proteolysis, a putative bacterial virulence factor target, a light-dependent catalyst, and a key redox-controlled enzyme in glutathione synthesis. Extending the latter observation we found that wounding increased the proportion of oxidized glutathione in leaves, but only in plants able to synthesize JA. The oxidizing conditions generated through JA signaling near wounds help to define the cellular environment in which proteome remodelling occurs.

Fatty acid ketodienes and fatty acid ketotrienes: Michael addition acceptors that accumulate in wounded and diseased Arabidopsis leaves.

Physical damage and disease are known to lead to changes in the oxylipin signature of plants. We searched for oxylipins produced in response to both wounding and pathogenesis in Arabidopsis leaves. Linoleic acid 9- and 13-ketodienes (KODEs) were found to accumulate in wounded leaves as well as in leaves infected with the pathogen Pseudomonas syringae pv. tomato (Pst). Quantification of the compounds showed that they accumulated to higher levels during the hypersensitive response to Pst avrRpm1 than during infection with a Pst strain lacking an avirulence gene. KODEs are Michael addition acceptors, containing a chemically reactive alpha,beta-unsaturated carbonyl group. When infiltrated into leaves, KODEs were found to induce expression of the GST1 gene, but vital staining indicated that these compounds also damaged plant cells. Several molecules typical of lipid oxidation, including malonaldehyde, also contain the alpha,beta-unsaturated carbonyl reactivity feature, and, when delivered in a volatile form, powerfully induced the expression of GST1. The results draw attention to the potential physiological importance of naturally occurring Michael addition acceptors in plants. In particular, these compounds could act directly, or indirectly via cell damage, as powerful gene activators and might also contribute to host cell death.

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.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

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.