Functional homogenization of bumblebee communities in alpine landscapes under projected climate change

Background: Bumblebees represent an active pollinator group in mountain regions and assure the pollination of many different plant species from low to high elevations. Plant-pollinator interactions are mediated by functional traits. Shift in bumblebee functional structure under climate change may impact plant-pollinator interactions in mountains. Here, we estimated bumblebee upward shift in elevation, community turnover, and change in functional structure under climate change. Method: We sampled bumblebee species at 149 sites along the elevation gradient. We used stacked species distribution models (S-SDMs) forecasted under three climate change scenarios (A2, A1B, RCP3PD) to model the potential distribution of the Bombus species. Furthermore, we used species proboscis length measurements to assess the functional change in bumblebee assemblages along the elevation gradient. Results: We found species-specific response of bumblebee species to climate change. Species differed in their predicted rate of range contraction and expansion. Losers were mainly species currently restricted to high elevation. Under the most severe climate change scenarios (A2), we found a homogenization of proboscis length structure in bumblebee communities along the elevation gradient through the upward colonization of high elevation by species with longer proboscides. Conclusions: Here, we show that in addition to causing the shift in the distribution of bumblebee species, climate change may impact the functional structure of communities. The colonization of high elevation areas by bumblebee species with long proboscides may modify the structure of plant-pollination interaction networks by increasing the diversity of pollination services at high elevation.

Early diagenesis of bone and tooth apatite in fluvial and marine settings: Constraints from combined oxygen isotope, nitrogen and REE analysis

Fossil bones and teeth of Late Pleistocene terrestrial mammals from Rhine River gravels (RS) and the North Sea (NS), that have been exposed to chemically and isotopically distinct diagenetic fluids (fresh water versus seawater), were investigated to study the effects of early diagenesis on biogenic apatite. Changes in phosphate oxygen isotopic composition (delta O-18(PO4)), nitrogen content (wt.% N) and rare earth element (REE) concentrations were measured along profiles within bones that have not been completely fossilized, and in skeletal tissues (bone, dentine, enamel) with different susceptibilities to diagenetic alteration. Early diagenetic changes of elemental and isotopic compositions of apatite in fossil bone are related to the loss of the stabilizing collagen matrix. The REE concentration is negatively correlated with the nitrogen content, and therefore the amount of collagen provides a sensitive proxy for early diagenetic alteration. REE patterns of RS and NS bones indicate initial fossilization in a fresh water fluid with similar REE compositions. Bones from both settings have nearly collagen-free, REE-, U-, F- and Sr-enriched altered outer rims, while the collagen-bearing bone compacta in the central part often display early diagenetic pyrite void-fillings. However, NS bones exposed to Holocene seawater have outer rim delta O-18(PO4) values that are 1.1 to 2.6 parts per thousand higher compared to the central part of the same bones (delta O-18(PO4) = 18.2 +/- 0.9 parts per thousand, n = 19). Surprisingly, even the collagen-rich bone compacta with low REE contents and apatite crystallinity seems altered, as NS tooth enamel (delta O-18(PO4) =15.0 +/- 0.3 parts per thousand, n=4) has about 3%. lower delta O-18(PO4) values, values that are also similar to those of enamel from RS teeth. Therefore, REE concentration, N content and apatite crystallinity are in this case only poor proxies for the alteration of delta O-18(PO4) values. Seawater exposure of a few years up to 8 kyr can change the delta O-18(PO4) values of the bone apatite by > 3 parts per thousand. Therefore, bones fossilized in marine settings must be treated with caution for palaeoclimatic reconstructions. However, enamel seems to preserve pristine delta O-18(PO4) values on this time scale. Using species-specific calibrations for modern mammals, a mean delta O-18(H2O) value can be reconstructed for Late Pleistocene mammalian drinking water of around -9.2 +/- 0.5 parts per thousand, which is similar to that of Late Pleistocene groundwater from central Europe. (c) 2008 Elsevier B.V. All rights reserved.

Polyethylene glycol and polyvinylpirrolidone effect on bacterial rRNA extraction and hybridization from cells exposed to tannins

In order to detect fluctuations in ruminal microbial populations due to forage tannins using 16S ribosomal RNA (rRNA) probes, recovery of intact rRNA is required. The objective of this work was to evaluate the effect of polyethylene glycol (PEG) and polyvinylpirrolidone (PVP) on extraction of bacterial rRNA, in the presence of tannins from tropical legume forages and other sources, that hybridize with oligonucleotide probes. Ruminococcus albus 8 cells were exposed to 8 g/L tannic acid or 1 g/L condensed tannins extracted from Acacia angustissima, banana (Musa sp.) skin, Desmodium ovalifolium, red grape (Vitis vinifera) skin and Inga edulis, or no tannins. Cells were rinsed with Tris buffer pH 7 containing either 8% PEG or 6% PVP prior to cell lysis. Total RNA samples rinsed with either PEG or PVP migrated through denaturing agarose gels. The 16S rRNA bands successfully hybridized with a R. albus species-specific oligonucleotide probe, regardless of tannin source. The effect of rinsing buffers on the density of 16S rRNA bands, as well as on the hybridization signals was compared. There were significant effects (P<0.01) when the controls were compared to either buffer treatments due to tannin type, buffer used and the interaction of tannin type and buffer. The significant interaction indicates the influence of tannin type on the parameters evaluated.

Introduction of snow and geomorphic disturbance variables into predictive models of alpine plant distribution in the Western Swiss Alps

Indirect topographic variables have been used successfully as surrogates for disturbance processes in plant species distribution models (SDM) in mountain environments. However, no SDM studies have directly tested the performance of disturbance variables. In this study, we developed two disturbance variables: a geomorphic index (GEO) and an index of snow redistribution by wind (SNOW). These were developed in order to assess how they improved both the fit and predictive power of presenceabsence SDM based on commonly used topoclimatic (TC) variables for 91 plants in the Western Swiss Alps. The individual contribution of the disturbance variables was compared to TC variables. Maps of models were prepared to spatially test the effect of disturbance variables. On average, disturbance variables significantly improved the fit but not the predictive power of the TC models and their individual contribution was weak (5.6% for GEO and 3.3% for SNOW). However their maximum individual contribution was important (24.7% and 20.7%). Finally, maps including disturbance variables (i) were significantly divergent from TC models in terms of predicted suitable surfaces and connectivity between potential habitats, and (ii) were interpreted as more ecologically relevant. Disturbance variables did not improve the transferability of models at the local scale in a complex mountain system, and the performance and contribution of these variables were highly species-specific. However, improved spatial projections and change in connectivity are important issues when preparing projections under climate change because the future range size of the species will determine the sensitivity to changing conditions.

Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction.

Ionotropic glutamate receptors (iGluRs) are a highly conserved family of ligand-gated ion channels present in animals, plants, and bacteria, which are best characterized for their roles in synaptic communication in vertebrate nervous systems. A variant subfamily of iGluRs, the Ionotropic Receptors (IRs), was recently identified as a new class of olfactory receptors in the fruit fly, Drosophila melanogaster, hinting at a broader function of this ion channel family in detection of environmental, as well as intercellular, chemical signals. Here, we investigate the origin and evolution of IRs by comprehensive evolutionary genomics and in situ expression analysis. In marked contrast to the insect-specific Odorant Receptor family, we show that IRs are expressed in olfactory organs across Protostomia--a major branch of the animal kingdom that encompasses arthropods, nematodes, and molluscs--indicating that they represent an ancestral protostome chemosensory receptor family. Two subfamilies of IRs are distinguished: conserved "antennal IRs," which likely define the first olfactory receptor family of insects, and species-specific "divergent IRs," which are expressed in peripheral and internal gustatory neurons, implicating this family in taste and food assessment. Comparative analysis of drosophilid IRs reveals the selective forces that have shaped the repertoires in flies with distinct chemosensory preferences. Examination of IR gene structure and genomic distribution suggests both non-allelic homologous recombination and retroposition contributed to the expansion of this multigene family. Together, these findings lay a foundation for functional analysis of these receptors in both neurobiological and evolutionary studies. Furthermore, this work identifies novel targets for manipulating chemosensory-driven behaviours of agricultural pests and disease vectors.

Arbuscular mycorrhizal fungi mediate below-ground plant-herbivore interactions: a phylogenetic study

Ecological interactions are complex networks, but have typically been studied in a pairwise fashion. Examining how third-party species can modify the outcome of pairwise interactions may allow us to better predict their outcomes in realistic systems. For instance, arbuscular mycorrhizal fungi (AMF) can affect plant interactions with other organisms, including below-ground herbivores, but the mechanisms underlying these effects remain unclear. Here, we use a comparative, phylogenetically controlled approach to test the relative importance of mycorrhizal colonization and plant chemical defences (cardenolides) in predicting plant survival and the abundance of a generalist below-ground herbivore across 14 species of milkweeds (Asclepias spp.). Plants were inoculated with a mixture of four generalist AMF species or left uninoculated. After 1month, larvae of Bradysia sp. (Diptera: Sciaridae), a generalist below-ground herbivore, colonized plant roots. We performed phylogenetically controlled analyses to assess the influence of AMF colonization and toxic cardenolides on plant growth, mortality and infestation by fungus gnats. Overall, plants inoculated with AMF exhibited greater survival than did uninoculated plants. Additionally, surviving inoculated plants had lower numbers of larvae in their roots and fewer non-AM fungi than surviving uninoculated plants. In phylogenetic controlled regressions, gnat density in roots was better predicted by the extent of root colonized by AMF than by root cardenolide concentration. Taken as a whole, AMF modify the effect of below-ground herbivores on plants in a species-specific manner, independent of changes in chemical defence. This study adds to the growing body of literature demonstrating that mycorrhizal fungi may improve plant fitness by conferring protection against antagonists, rather than growth benefits. In addition, we advocate using comparative analyses to disentangle the roles of shared history and ecology in shaping trait expression and to better predict the outcomes of complex multitrophic interactions.

Evolution and functional characterisation of the IR chemosensory receptors in the Drosophila larval gustatory system

Chemosensory receptor gene families encode divergent proteins capable of detecting a huge diversity of environmental stimuli that are constantly changing over evolutionary time as organisms adapt to distinct ecological niches. While olfaction is dedicated to the detection of volatile compounds, taste is key to assess food quality for nutritional value and presence of toxic substances. The sense of taste also provides initial signals to mediate endocrine regulation of appetite and food metabolism and plays a role in kin recognition. The fruit fly Drosophila melanogaster is a very good model for studying smell and taste because these senses are very important in insects and because a broad variety of genetic tools are available in Drosophila. Recently, a family of 66 chemosensory receptors, the Ionotropic Receptors (IRs) was described in fruit flies. IRs are distantly related to ionotropic glutamate receptors (iGluRs), but their evolutionary origin from these synaptic receptors is unclear. While 16 IRs are expressed in the olfactory system, nothing is known about the other members of this repertoire. In this thesis, I describe bioinformatic, expression and functional analyses of the IRs aimed at understanding how these receptors have evolved, and at characterising the role of the non-olfactory IRs. I show that these have emerged at the basis of the protostome lineage and probably have acquired their sensory function very early. Moreover, although several IRs are conserved across insects, there are rapid and dramatic changes in the size and divergence of IR repertoires across species. I then performed a comprehensive analysis of IR expression in the larva of Drosophila melanogaster, which is a good model to study taste and feeding mechanisms as it spends most of its time eating or foraging. I found that most of the divergent members of the IR repertoire are expressed in both peripheral and internal gustatory neurons, suggesting that these are involved in taste perception. Finally, through the establishment of a new neurophysiological assay in larvae, I identified for the first time subsets of IR neurons that preferentially detect sugars and amino acids, indicating that IRs might be involved in sensing these compounds. Together, my results indicate that IRs are an evolutionarily dynamic and functionally versatile family of receptors. In contrast to the olfactory IRs that are well-conserved, gustatory IRs are rapidly evolving species-specific receptors that are likely to be involved in detecting a wide variety of tastants. - La plupart des animaux possèdent de grandes familles de récepteurs chimiosensoriels dont la fonction est de détecter l'immense diversité de composés chimiques présents dans l'environnement. Ces récepteurs évoluent en même temps que les organismes s'adaptent à leur écosystème. Il existe deux manières de percevoir ces signaux chimiques : l'olfaction et le goût. Alors que le système olfactif perçoit les composés volatiles, le sens du goût permet d'évaluer, par contact, la qualité de la nourriture, de détecter des substances toxiques et de réguler l'appétit et le métabolisme. L'un des organismes modèles les plus pertinents pour étudier le sens du goût est le stade larvaire de la mouche du vinaigre Drosophila melanogaster. En effet, la principale fonction du stade larvaire est de trouver de la nourriture et de manger. De plus, il est possible d'utiliser tous les outils génétiques développés chez la drosophile. Récemment, une nouvelle famille de 66 récepteurs chimiosensoriels appelés Récepteurs Ionotropiques (IRs) a été découverte chez la drosophile. Bien que leur orogine soit peu claire, ces récepteurs sont similaires aux récepteurs ionotropiques glutamatergiques impliqués dans la transmission synaptique. 16 IRs sont exprimés dans le système olfactif de la mouche adulte, mais pour l'instant on ne connaît rien des autres membres de cette famille. Durant ma thèse, j'ai effectué des recherches sur l'évolution de ces récepteurs ainsi que sur l'expression et la fonction des IRs non olfactifs. Je démontre que les IRs sont apparus chez l'ancêtre commun des protostomiens et ont probablement acquis leur fonction sensorielle très rapidement. De plus, bien qu'un certain nombre d'IRs olfactifs soient conservés chez les insectes, d'importantes variations dans la taille et la divergence des répertoires d'IRs entre les espèces ont été constatées. J'ai également découvert qu'un grand nombre d'IRs non olfactifs sont exprimés dans différents organes gustatifs, ce qui leur confère probablement une fonction dans la perception des goûts. Finalement, pour la première fois, des neurones exprimant des IRs ont été identifiés pour leur fonction dans la perception de sucres et d'acides aminés chez la larve. Mes résultats présentent les IRs comme une famille très dynamique, aux fonctions très variées, qui joue un rôle tant dans l'odorat que dans le goût, et dont la fonction est restée importante tout au long de l'évolution. De plus, l'identification de neurones spécialisés dans la perception de certains composés permettra l'étude des circuits neuronaux impliqués dans le traitement de ces informations.

The genome sequence of taurine cattle: a window to ruminant biology and evolution.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

A matter of time: delayed mate encounter postpones mating window initiation and reduces the strength of female choosiness

Reproductive success is determined by the presence and timing of encounter of mates. The latter depends on species-specific reproductive characteristics (e.g. initiation/duration of the mating window), season, and reproductive strategies (e.g. intensity of choosiness) that may potentially mitigate constraints imposed by mating windows. Despite their potentially crucial role for fitness and population dynamics, limited evidence exists about mating window initiation, duration and reproductive strategies. Here, we experimentally tested the mechanisms of initiation and the duration of the common lizard's Zootoca vivipara mating window, by manipulating the timing of mate encounter and analyzing its effect on (re-)mating probability. We furthermore tested treatment effects on female reproductive strategies, by measuring female choosiness. The timing of mate encounter and season did not significantly affect mating probability. However, a longer delay until mate encounter reduced female choosiness. Re-mating probability decreased with re-mating delay and was independent of mating delay. This indicates that mating window initiation depends on mate encounter, that its duration is fixed, and that plastic reproductive strategies exist. These findings contrast with previous beliefs and shows that mating windows per se may not necessarily constrain reproductive success, which is congruent with rapid range expansion and absence of positive density-effects on reproductive success (Allee effects). In summary, our results show that predicting the effect of mating windows on reproduction is complex and that experimental evidence is essential for evaluating their effect on reproduction and reproductive strategies, both being important determinants of population dynamics and the colonization of new habitats.

Context-dependent dual role of SKI8 homologs in mRNA synthesis and turnover.

Eukaryotic mRNA transcription and turnover is controlled by an enzymatic machinery that includes RNA polymerase II and the 3' to 5' exosome. The activity of these protein complexes is modulated by additional factors, such as the nuclear RNA polymerase II-associated factor 1 (Paf1c) and the cytoplasmic Superkiller (SKI) complex, respectively. Their components are conserved across uni- as well as multi-cellular organisms, including yeast, Arabidopsis, and humans. Among them, SKI8 displays multiple facets on top of its cytoplasmic role in the SKI complex. For instance, nuclear yeast ScSKI8 has an additional function in meiotic recombination, whereas nuclear human hSKI8 (unlike ScSKI8) associates with Paf1c. The Arabidopsis SKI8 homolog VERNALIZATION INDEPENDENT 3 (VIP3) has been found in Paf1c as well; however, whether it also has a role in the SKI complex remains obscure so far. We found that transgenic VIP3-GFP, which complements a novel vip3 mutant allele, localizes to both nucleus and cytoplasm. Consistently, biochemical analyses suggest that VIP3-GFP associates with the SKI complex. A role of VIP3 in the turnover of nuclear encoded mRNAs is supported by random-primed RNA sequencing of wild-type and vip3 seedlings, which indicates mRNA stabilization in vip3. Another SKI subunit homolog mutant, ski2, displays a dwarf phenotype similar to vip3. However, unlike vip3, it displays neither early flowering nor flower development phenotypes, suggesting that the latter reflect VIP3's role in Paf1c. Surprisingly then, transgenic ScSKI8 rescued all aspects of the vip3 phenotype, suggesting that the dual role of SKI8 depends on species-specific cellular context.

The origins, evolution, and functional potential of alternative splicing in vertebrates.

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

UniPathway: a resource for the exploration and annotation of metabolic pathways.

UniPathway (http://www.unipathway.org) is a fully manually curated resource for the representation and annotation of metabolic pathways. UniPathway provides explicit representations of enzyme-catalyzed and spontaneous chemical reactions, as well as a hierarchical representation of metabolic pathways. This hierarchy uses linear subpathways as the basic building block for the assembly of larger and more complex pathways, including species-specific pathway variants. All of the pathway data in UniPathway has been extensively cross-linked to existing pathway resources such as KEGG and MetaCyc, as well as sequence resources such as the UniProt KnowledgeBase (UniProtKB), for which UniPathway provides a controlled vocabulary for pathway annotation. We introduce here the basic concepts underlying the UniPathway resource, with the aim of allowing users to fully exploit the information provided by UniPathway.

Development and application of computational methodologies in qualitative modeling

La présente thèse s'intitule "Développent et Application des Méthodologies Computationnelles pour la Modélisation Qualitative". Elle comprend tous les différents projets que j'ai entrepris en tant que doctorante. Plutôt qu'une mise en oeuvre systématique d'un cadre défini a priori, cette thèse devrait être considérée comme une exploration des méthodes qui peuvent nous aider à déduire le plan de processus regulatoires et de signalisation. Cette exploration a été mue par des questions biologiques concrètes, plutôt que par des investigations théoriques. Bien que tous les projets aient inclus des systèmes divergents (réseaux régulateurs de gènes du cycle cellulaire, réseaux de signalisation de cellules pulmonaires) ainsi que des organismes (levure à fission, levure bourgeonnante, rat, humain), nos objectifs étaient complémentaires et cohérents. Le projet principal de la thèse est la modélisation du réseau de l'initiation de septation (SIN) du S.pombe. La cytokinèse dans la levure à fission est contrôlée par le SIN, un réseau signalant de protéines kinases qui utilise le corps à pôle-fuseau comme échafaudage. Afin de décrire le comportement qualitatif du système et prédire des comportements mutants inconnus, nous avons décidé d'adopter l'approche de la modélisation booléenne. Dans cette thèse, nous présentons la construction d'un modèle booléen étendu du SIN, comprenant la plupart des composantes et des régulateurs du SIN en tant que noeuds individuels et testable expérimentalement. Ce modèle utilise des niveaux d'activité du CDK comme noeuds de contrôle pour la simulation d'évènements du SIN à différents stades du cycle cellulaire. Ce modèle a été optimisé en utilisant des expériences d'un seul "knock-out" avec des effets phénotypiques connus comme set d'entraînement. Il a permis de prédire correctement un set d'évaluation de "knock-out" doubles. De plus, le modèle a fait des prédictions in silico qui ont été validées in vivo, permettant d'obtenir de nouvelles idées de la régulation et l'organisation hiérarchique du SIN. Un autre projet concernant le cycle cellulaire qui fait partie de cette thèse a été la construction d'un modèle qualitatif et minimal de la réciprocité des cyclines dans la S.cerevisiae. Les protéines Clb dans la levure bourgeonnante présentent une activation et une dégradation caractéristique et séquentielle durant le cycle cellulaire, qu'on appelle communément les vagues des Clbs. Cet évènement est coordonné avec la courbe d'activation inverse du Sic1, qui a un rôle inhibitoire dans le système. Pour l'identification des modèles qualitatifs minimaux qui peuvent expliquer ce phénomène, nous avons sélectionné des expériences bien définies et construit tous les modèles minimaux possibles qui, une fois simulés, reproduisent les résultats attendus. Les modèles ont été filtrés en utilisant des simulations ODE qualitatives et standardisées; seules celles qui reproduisaient le phénotype des vagues ont été gardées. L'ensemble des modèles minimaux peut être utilisé pour suggérer des relations regulatoires entre les molécules participant qui peuvent ensuite être testées expérimentalement. Enfin, durant mon doctorat, j'ai participé au SBV Improver Challenge. Le but était de déduire des réseaux spécifiques à des espèces (humain et rat) en utilisant des données de phosphoprotéines, d'expressions des gènes et des cytokines, ainsi qu'un réseau de référence, qui était mis à disposition comme donnée préalable. Notre solution pour ce concours a pris la troisième place. L'approche utilisée est expliquée en détail dans le dernier chapitre de la thèse. -- The present dissertation is entitled "Development and Application of Computational Methodologies in Qualitative Modeling". It encompasses the diverse projects that were undertaken during my time as a PhD student. Instead of a systematic implementation of a framework defined a priori, this thesis should be considered as an exploration of the methods that can help us infer the blueprint of regulatory and signaling processes. This exploration was driven by concrete biological questions, rather than theoretical investigation. Even though the projects involved divergent systems (gene regulatory networks of cell cycle, signaling networks in lung cells), as well as organisms (fission yeast, budding yeast, rat, human), our goals were complementary and coherent. The main project of the thesis is the modeling of the Septation Initiation Network (SIN) in S.pombe. Cytokinesis in fission yeast is controlled by the SIN, a protein kinase signaling network that uses the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this thesis, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN. Another cell cycle related project that is part of this thesis was to create a qualitative, minimal model of cyclin interplay in S.cerevisiae. CLB proteins in budding yeast present a characteristic, sequential activation and decay during the cell cycle, commonly referred to as Clb waves. This event is coordinated with the inverse activation curve of Sic1, which has an inhibitory role in the system. To generate minimal qualitative models that can explain this phenomenon, we selected well-defined experiments and constructed all possible minimal models that, when simulated, reproduce the expected results. The models were filtered using standardized qualitative ODE simulations; only the ones reproducing the wave-like phenotype were kept. The set of minimal models can be used to suggest regulatory relations among the participating molecules, which will subsequently be tested experimentally. Finally, during my PhD I participated in the SBV Improver Challenge. The goal was to infer species-specific (human and rat) networks, using phosphoprotein, gene expression and cytokine data and a reference network provided as prior knowledge. Our solution to the challenge was selected as in the final chapter of the thesis.

The evolution of alternative splicing patterns and regulatory mechanisms

Alternative splicing produces multiple isoforms from the same gene, thus increasing the number of transcripts of the species. Alternative splicing is a virtually ubiquitous mechanism in eukaryotes, for example more than 90% of protein-coding genes in human are alternatively spliced. Recent evolutionary studies showed that alternative splicing is a fast evolving and highly species- specific mechanism. The rapid evolution of alternative splicing was considered as a contribution to the phenotypic diversity between species. However, the function of many isoforms produced by alternative splicing remains unclear and they might be the result of noisy splicing. Thus, the functional relevance of alternative splicing and the evolutionary mechanisms of its rapid divergence among species are still poorly understood. During my thesis, I performed a large-scale analysis of the regulatory mechanisms that drive the rapid evolution of alternative splicing. To study the evolution of alternative splicing regulatory mechanisms, I used an extensive RNA-sequencing dataset comprising 12 tetrapod species (human, chimpanzee and bonobo, gorilla, orangutan, macaque, marmoset, mouse, opossum, platypus, chicken and frog) and 8 tissues (cerebellum, brain, heart, kidney, liver, testis, placenta and ovary). To identify the catalogue of alternative splicing eis-acting regulatory elements in the different tetrapod species, I used a previously defined computational approach. This approach is a statistical analysis of exons/introns and splice sites composition and relies on a principle of compensation between splice sites strength and the presence of additional regulators. With an evolutionary comparative analysis of the exonic eis-acting regulators, I showed that these regulatory elements are generally shared among primates and more conserved than non-regulatory elements. In addition, I showed that the usage of these regulatory elements is also more conserved than expected by chance. In addition to the identification of species- specific eis-acting regulators, these results may explain the rapid evolution of alternative splicing. I also developed a new approach based on evolutionary sequence changes and corresponding alternative splicing changes to identify potential splicing eis-acting regulators in primates. The identification of lineage-specific substitutions and corresponding lineage-specific alternative splicing changes, allowed me to annotate the genomic sequences that might have played a role in the alternative splicing pattern differences among primates. Finally, I showed that the identified splicing eis-acting regulator datasets are enriched in human disease-causing mutations, thus confirming their biological relevance.

Past climate-driven range shifts and population genetic diversity in arctic plants

AimHigh intra-specific genetic diversity is necessary for species adaptation to novel environments under climate change, but species tracking suitable conditions are losing alleles through successive founder events during range shift. Here, we investigated the relationship between range shift since the Last Glacial Maximum (LGM) and extant population genetic diversity across multiple plant species to understand variability in species responses. LocationThe circumpolar Arctic and northern temperate alpine ranges. MethodsWe estimated the climatic niches of 30 cold-adapted plant species using range maps coupled with species distribution models and hindcasted species suitable areas to reconstructions of the mid-Holocene and LGM climates. We computed the species-specific migration distances from the species glacial refugia to their current distribution and correlated distances to extant genetic diversity in 1295 populations. Differential responses among species were related to life-history traits. ResultsWe found a negative association between inferred migration distances from refugia and genetic diversities in 25 species, but only 11 had statistically significant negative slopes. The relationships between inferred distance and population genetic diversity were steeper for insect-pollinated species than wind-pollinated species, but the difference among pollination system was marginally independent from phylogenetic autocorrelation. Main conclusionThe relationships between inferred migration distances and genetic diversities in 11 species, independent from current isolation, indicate that past range shifts were associated with a genetic bottleneck effect with an average of 21% loss of genetic diversity per 1000km(-1). In contrast, the absence of relationship in many species also indicates that the response is species specific and may be modulated by plant pollination strategies or result from more complex historical contingencies than those modelled here.