Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Phylogenetic analyses of C4 photosynthesis evolution in grasses: genetics of convergence and convergence of genes

Summary : During the evolutionary diversification of organisms, similar ecological constraints led to the recurrent appearances of the same traits (phenotypes) in distant lineages, a phenomenon called convergence. In most cases, the genetic origins of the convergent traits remain unknown, but recent studies traced the convergent phenotypes to recurrent alterations of the same gene or, in a few cases, to identical genetic changes. However, these cases remain anecdotal and there is a need for a study system that evolved several times independently and whose genetic determinism is well resolved and straightforward, such as C4 photosynthesis. This adaptation to warm environments, possibly driven by past atmospheric CO2 decreases, consists in a CO2-concentrating pump, created by numerous morphological and biochemical novelties. All genes encoding C4 enzymes already existed in C3 ancestors, and are supposed to have been recruited through gene duplication followed by neo-functionalization, to acquire the cell specific expression pattern and altered kinetic properties that characterize Ca-specific enzymes. These predictions have so far been tested only in species-poor and ecologically marginal C4 dicots. The monocots, and especially the grass family (Poaceae), the most important C4 family in terms of species number, ecological dominance and economical importance, have been largely under-considered as suitable study systems. This thesis aimed at understanding the evolution of the C4 trait in grasses at a molecular level and to use the genetics of C4 photosynthesis to infer the evolutionary history of the C4 phenotype and its driving selective pressures. A molecular phylogeny of grasses and affiliated monocots identified 17 to 18 independent acquisitions of the C4 pathway in the grass family. A relaxed molecular clock was used to date these events and the first C4 evolution was estimated in the Chloridoideae subfamily, between 32-25 million years ago, at a period when atmospheric CO2 abruptly declined. Likelihood models showed that after the COZ decline the probability of evolving the C4 pathway strongly increased, confirming low CO2 as a likely driver of C4 photosynthesis evolution. In order to depict the genetic changes linked to the numerous C4 origins, genes encoding phopshoenolpyruvate carboxylase (PEPC), the key-enzyme responsible for the initial fixation of atmospheric CO2 in the C4 pathway, were isolated from a large sample of C3 and C4 grasses. Phylogenetic analyses were used to reconstruct the evolutionary history of the PEPC multigene family and showed that the evolution of C4-specific PEPC had been driven by positive selection on 21 codons simultaneously in up to eight C4 lineages. These selective pressures led to numerous convergent genetic changes in many different C4 clades, highlighting the repeatability of some evolutionary processes, even at the molecular level. PEPC C4-adaptive changes were traced and used to show multiple appearances of the C, pathway in clades where species tree inferences were unable to differentiate multiple C4 appearances and a single appearance followed by C4 to C3 reversion. Further investigations of genes involved in some of the C4 subtypes only (genes encoding decarboxylating enzymes NADP-malic enzyme and phosphoenolpyruvate carboxykinase) showed that these C4-enzymes also evolved through strong positive selection and underwent parallel genetic changes during the different Ca origins. The adaptive changes on these subtype-specific C4 genes were used to retrace the history of the C4-subtypes phenotypes, which revealed that the evolution of C4-PEPC and C4-decarboxylating enzymes was in several cases disconnected, emphasizing the multiplicity of the C4 trait and the gradual acquisition of the features that create the CO2-pump. Finally, phylogenetic analyses of a gene encoding the Rubisco (the enzyme responsible for the fixation of CO2 into organic compounds in all photosynthetic organisms) showed that C4 evolution switched the selective pressures on this gene. Five codons were recurrently mutated to adapt the enzyme kinetics to the high CO2 concentrations of C4 photosynthetic cells. This knowledge could be used to introgress C4-like Rubisco in C3 crops, which could lead to an increased yield under predicted future high CO2 atmosphere. Globally, the phylogenetic framework adopted during this thesis demonstrated the widespread occurrence of genetic convergence on C4-related enzymes. The genetic traces of C4 photosynthesis evolution allowed reconstructing events that happened during the last 30 million years and proved the usefulness of studying genes directly responsible for phenotype variations when inferring evolutionary history of a given trait. Résumé Durant la diversification évolutive des organismes, des pressions écologiques similaires ont amené à l'apparition récurrente de certains traits (phénotypes) dans des lignées distantes, un phénomène appelé évolution convergente. Dans la plupart des cas, l'origine génétique des traits convergents reste inconnue mais des études récentes ont montré qu'ils étaient dus dans certains cas à des changements répétés du même gène ou, dans de rares cas, à des changements génétiques identiques. Malgré tout, ces cas restent anecdotiques et il y a un réel besoin d'un système d'étude qui ait évolué indépendamment de nombreuses fois et dont le déterminisme génétique soit clairement identifié. La photosynthèse dite en Ça répond à ces critères. Cette adaptation aux environnements chauds, dont l'évolution a pu être encouragé par des baisses passées de la concentration atmosphérique en CO2, est constituée de nombreuses nouveautés morphologiques et biochimiques qui créent une pompe à CO2. La totalité des gènes codant les enzymes Ç4 étaient déjà présents dans les ancêtres C3. Leur recrutement pour la photosynthèse Ç4 est supposé s'être fait par le biais de duplications géniques suivies par une néo-fonctionnalisation pour leur conférer l'expression cellule-spécifique et les propriétés cinétiques qui caractérisent les enzymes C4. Ces prédictions n'ont jusqu'à présent été testées que dans des familles C4 contenant peu d'espèces et ayant un rôle écologique marginal. Les graminées (Poaceae), qui sont la famille C4 la plus importante, tant en termes de nombre d'espèces que de dominance écologique et d'importance économique, ont toujours été considérés comme un système d'étude peu adapté et ont fait le sujet de peu d'investigations évolutives. Le but de cette thèse était de comprendre l'évolution de la photosynthèse en C4 chez les graminées au niveau génétique et d'utiliser les gènes pour inférer l'évolution du phénotype C4 ainsi que les pressions de sélection responsables de son évolution. Une phylogénie moléculaire de la famille des graminées et des monocotylédones apparentés a identifié 17 à 18 acquisitions indépendantes de la photosynthèse chez les graminées. Grâce à une méthode d'horloge moléculaire relâchée, ces évènements ont été datés et la première apparition C4 a été estimée dans la sous-famille des Chloridoideae, il y a 32 à 25 millions d'années, à une période où les concentrations atmosphériques de CO2 ont décliné abruptement. Des modèles de maximum de vraisemblance ont montré qu'à la suite du déclin de CO2, la probabilité d'évoluer la photosynthèse C4 a fortement augmenté, confirmant ainsi qu'une faible concentration de CO2 est une cause potentielle de l'évolution de la photosynthèse C4. Afin d'identifier les mécanismes génétiques responsables des évolutions répétées de la photosynthèse C4, un segment des gènes codant pour la phosphoénolpyruvate carboxylase (PEPC), l'enzyme responsable de la fixation initiale du CO2 atmosphérique chez les plantes C4, ont été séquencés dans une centaine de graminées C3 et C4. Des analyses phylogénétiques ont permis de reconstituer l'histoire évolutive de la famille multigénique des PEPC et ont montré que l'évolution de PEPC spécifiques à la photosynthèse Ça a été causée par de la sélection positive agissant sur 21 codons, et ce simultanément dans huit lignées C4 différentes. Cette sélection positive a conduit à un grand nombre de changements génétiques convergents dans de nombreux clades différents, ce qui illustre la répétabilité de certains phénomènes évolutifs, et ce même au niveau génétique. Les changements sur la PEPC liés au C4 ont été utilisés pour confirmer des évolutions indépendantes du phénotype C4 dans des clades où l'arbre des espèces était incapable de différencier des apparitions indépendantes d'une seule apparition suivie par une réversion de C4 en C3. En considérant des gènes codant des protéines impliquées uniquement dans certains sous-types C4 (deux décarboxylases, l'enzyme malique à NADP et la phosphoénolpyruvate carboxykinase), des études ultérieures ont montré que ces enzymes C4 avaient elles-aussi évolué sous forte sélection positive et subi des changements génétiques parallèles lors des différentes origines de la photosynthèse C4. Les changements adaptatifs sur ces gènes liés seulement à certains sous-types C4 ont été utilisés pour retracer l'histoire des phénotypes de sous-types C4, ce qui a révélé que les caractères formant le trait C4 ont, dans certains cas, évolué de manière déconnectée. Ceci souligne la multiplicité du trait C4 et l'acquisition graduelle de composants participant à la pompe à CO2 qu'est la photosynthèse C4. Finalement, des analyses phylogénétiques des gènes codant pour la Rubisco (l'enzyme responsable de la fixation du CO2 en carbones organiques dans tous les organismes photosynthétiques) ont montré que l'évolution de la photosynthèse Ça a changé les pressions de sélection sur ce gène. Cinq codons ont été mutés de façon répétée afin d'adapter les propriétés cinétiques de la Rubisco aux fortes concentrations de CO2 présentes dans les cellules photosynthétiques des plantes C4. Globalement, l'approche phylogénétique adoptée durant cette thèse de doctorat a permis de démontré des phénomène fréquents de convergence génétique sur les enzymes liées à la photosynthèse C4. Les traces génétiques de l'évolution de la photosynthèse C4 ont permis de reconstituer des évènements qui se sont produits durant les derniers 30 millions d'années et ont prouvé l'utilité d'étudier des gènes directement responsables des variations phénotypiques pour inférer l'histoire évolutive d'un trait donné.

Scale-dependent adaptive evolution and morphological convergence to climatic niche in Californian eriogonoids (Polygonaceae)

Aim Macroevolutionary patterns and processes change substantially depending on levels of taxonomic and ecological organization, and the resolution of environmental and spatial variability. In comparative methods, the resolution of environmental and spatial variability often defines the number of selective regimes used to test whether phenotypic characteristics are adaptively correlated with the environment. Here, we examine how investigator choice of the number of selective regimes, determined by varying the resolution of among-species variability in the species climatic niche (hereafter called ecological scale'), influences trait morphological diversification among Eriogonoideae species. We assess whether adaptive or neutral processes drive the evolution of several morphological traits in these species. Location South-western North America. Methods We applied a phylogenetic framework of three evolutionary models to four morphological traits and the climatic niches of Eriogonoideae (in the buckwheat family, Polygonaceae). We tested whether morphological traits evolve in relation to climate by adaptive or neutral process, and whether the resulting patterns of morphological variability are conserved or convergent across the clade. We inspected adaptive models of evolution under different levels of resolution of among-species variability of the climatic niche. Results We show that morphological traits and climate niches of Eriogonoideae species are not phylogenetically conserved. Further, adaptive evolution of phenotypic traits is specific to climatic niche occupancy across this clade. Finally, the likely evolutionary process and the level of detectable niche conservatism change depending on the resolution of environmental variability of the climatic niche. Main conclusions Our study demonstrates the need to consider both the resolution of environmental variability and alternative evolutionary models to understand the morphological diversification that accompanies divergent adaptive evolution of lineages to climatic conditions.

Formation of aluminosilicate-bearing quartz veins in the Simano nappe (Central Alps): structural, thermobarometric and oxygen isotope constraints

We combined structural analysis, thermobarometry and oxygen isotope geochemistry to constrain the evolution of kyanite and/or andalusite-bearing quartz veins from the amphibolite facies metapelites of the Simano nappe, in the Central Alps of Switzerland. The Simano nappe records a complex polyphase tectonic evolution associated with nappe stacking during Tertiary Alpine collision (D1). The second regional deformation phase (132) is responsible for the main penetrative schistosity and mineral lineation, and formed during top-to-the-north thrusting. During the next stage of deformation (D3) the aluminosilicate-bearing veins formed by crystallization in tension gashes, in tectonic shadows of boudins, as well as along shear bands associated with top-to-the-north shearing. D2 and D3 are coeval with the Early Miocene metamorphic peak, characterised by kyanite + staurolite + garnet + biotite assemblages in metapelites. The peak pressure (P) and temperature (T) conditions recorded are constrained by multiple-equilibrium thermobarometry at 630 +/- 20 degrees C and 8.5 +/- 1 kbar (similar to 27 km depth), which is in agreement with oxygen isotope thermometry indicating isotopic equilibration of quartz-kyanite pairs at 670 +/- 50 degrees C. Quartz-kyanite pairs from the aluminosilicate-bearing quartz veins yield equilibration temperatures of 645 +/- 20 degrees C, confirming that the veins formed under conditions near metamorphic peak. Quartz and kyanite from veins and the surrounding metapelites have comparable isotopic compositions. Local intergranular diffusion in the border of the veins controls the mass-transfer and the growth of the product assemblage, inducing local mobilization of SiO2 and Al2O3. Andalusite is absent from the host rocks, but it is common in quartz veins, where it often pseudomorphs kyanite. For andalusite to be stable at T-max, the pressure in the veins must have been substantially lower than lithostatic. An alternative explanation consistent with structural observations would be inheritance by andalusite of the kyanite isotopic signature during polymorphic transformation after the metamorphic peak.

Stratigraphic, sedimentological and palaeoenvironmental constraints on the rise of the Urgonian platform in the western Swiss Jura

Urgonian-type carbonates are a characteristic feature of many late Early Cretaceous shallow-marine, tropical and subtropical environments. The presence of typical photozoan carbonate-producing communities including corals and rudists indicates the prevalence of warm, transparent and presumably oligotrophic conditions in a period otherwise characterized by the high density of globally occurring anoxic episodes. Of particular interest, therefore, is the exploration of relationships between Urgonian platform growth and palaeoceanographic change. In the French and Swiss Jura Mountains, the onset and evolution of the Urgonian platform have been controversially dated, and a correlation with other, better dated, successions is correspondingly difficult. It is for this reason that the stratigraphy and sedimentology of a series of recently exposed sections (Eclepens, Vaumarcus and Neuchatel) and, in addition, the section of the Gorges de l'Areuse were analysed. Calcareous nannofossil biostratigraphy, the evolution of phosphorus contents of bulk rock, a sequence-stratigraphic interpretation and a correlation of drowning unconformities with better dated sections in the Helvetic Alps were used to constrain the age of the Urgonian platform. The sum of the data and field observations suggests the following evolution: during the Hauterivian, important outward and upward growth of a bioclastic and oolitic carbonate platform is documented in two sequences, separated by a phase of platform drowning during the late Early Hauterivian. Following these two phases of platform growth, a second drowning phase occurred during the latest Hauterivian and Early Barremian, which was accompanied by significant platform erosion and sediment reworking. The Late Barremian witnessed the renewed installation of a carbonate platform, which initiated with a phase of oolite production, and which progressively evolved into a typical Urgonian carbonate platform colonized by corals and rudists. This phase terminated at the latest in the middle Early Aptian, due to a further drowning event. The evolution of this particular platform segment is compatible with that of more distal and well-dated segments of the same northern Tethyan platform preserved in the Helvetic zone of the Alps and in the northern subalpine chains (Chartreuse and Vercors).

Ecology and evolution of social-organization: insights from fire ants and other highly eusocial insects

Social organisms exhibit conspicuous intraspecific variation in all facets of their social organization. A prominent example of such variation in the highly eusocial Hymenoptera is differences in the number of reproductive queens per colony, Differences in queen number in ants are associated with differences in a host of reproductive and social traits, including queen phenotype and breeding strategy, mode of colony reproduction, and pattern of sex allocation. We examine the causes and consequences of changes in colony queen number and associated traits using the fire ant Solenopsis invicta as a principal model. Ecological constraints on mode of colony founding may act as important selective forces causing the evolution of queen number in this and many other ants, with social organization generally perpetuated across generations by means of the social environment molding appropriate queen phenotypes and reproductive strategies. Shifts in colony queen number have profound effects on genetic structure within nests and may also influence genetic structure at higher levels (aggregations of nests or local demes) because of the association of queen number with particular mating and dispersal habits. Divergence of breeding habits between populations with different social organizations has the potential to promote genetic differentiation between these social variants. Thus, evolution of social organization can be important in generating intrinsic selective regimes that channel subsequent social evolution and in initiating the development of significant population genetic structure, including barriers to gene flow important in cladogenesis.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Growth decrease in Alpine whitefish: investigating the relative contribution of fishing induced selection and environmental change and its implication for conservation measures

La pression exercée par les activités humaines menace pratiquement tous les écosystèmes aquatiques du globe. Ainsi, sous l'effet de divers facteurs tels que la pollution, le réchauffement climatique ou encore la pêche industrielle, de nombreuses populations de poissons ont vu leurs effectifs chuter et divers changements morphologiques ont été observés. Dans cette étude, nous nous sommes intéressés à une menace particulière: la sélection induite par la pêche sur la croissance des poissons. En effet, la génétique des populations prédit que la soustraction régulière des individus les plus gros peut entraîner des modifications rapides de certains traits physiques comme la croissance individuelle. Cela a par ailleurs été observé dans de nombreuses populations marines ou lacustres, dont les populations de féras, bondelles et autres corégones des lacs suisses. Toutefois, malgré un nombre croissant d'études décrivant ce phénomène, peu de plans de gestion en tiennent compte, car l'importance des effets génétiques liés à la pêche est le plus souvent négligée par rapport à l'impact des changements environnementaux. Le but premier de cette étude a donc été de quantifier l'importance des facteurs génétiques et environnementaux. Dans le premier chapitre, nous avons étudié la population de palée du lac de Joux (Coregonus palaea). Nous avons déterminé les différentiels de sélection dus à la pêche, c'est-à-dire l'intensité de la sélection sur le taux de croissance, ainsi que les changements nets de croissance au cours du temps. Nous avons observé une baisse marquée de croissance et un différentiel de sélection important indiquant qu'au moins 30% de la diminution de croissance observée était due à la pression de sélection induite par la pêche. Dans le deuxième chapitre, nous avons effectué les mêmes analyses sur deux espèces proches du lac de Brienz (C. albellus et C. fatioi) et avons observé des effets similaires dont l'intensité était spécifique à chaque espèce. Dans le troisième chapitre, nous avons analysé deux autres espèces : C. palaea et C. confusus du lac de Bienne, et avons constaté que le lien entre la pression de sélection et la diminution de croissance était influencé par des facteurs environnementaux. Finalement, dans le dernier chapitre, nous avons étudié les effets potentiels de différentes modifications de la taille des mailles des filets utilisés pour la pêche à l'aide de modèles mathématiques. Nous concluons que la pêche a un effet génétique non négligeable (et donc peu réversible) sur la croissance individuelle dans les populations observée, que cet effet est lié à la compétition pour la nourriture et à la qualité de l'environnement, et que certaines modifications simples de la taille des mailles des filets de pêche pourraient nettement diminuer l'effet de sélection et ainsi ralentir, voir même renverser la diminution de croissance observée.

Ecological-economic optimization of biodiversity conservation under climate change

Substantial investment in climate change research has led to dire predictions of the impacts and risks to biodiversity. The Intergovernmental Panel on Climate Change fourth assessment report(1) cites 28,586 studies demonstrating significant biological changes in terrestrial systems(2). Already high extinction rates, driven primarily by habitat loss, are predicted to increase under climate change(3-6). Yet there is little specific advice or precedent in the literature to guide climate adaptation investment for conserving biodiversity within realistic economic constraints(7). Here we present a systematic ecological and economic analysis of a climate adaptation problem in one of the world's most species-rich and threatened ecosystems: the South African fynbos. We discover a counterintuitive optimal investment strategy that switches twice between options as the available adaptation budget increases. We demonstrate that optimal investment is nonlinearly dependent on available resources, making the choice of how much to invest as important as determining where to invest and what actions to take. Our study emphasizes the importance of a sound analytical framework for prioritizing adaptation investments(4). Integrating ecological predictions in an economic decision framework will help support complex choices between adaptation options under severe uncertainty. Our prioritization method can be applied at any scale to minimize species loss and to evaluate the robustness of decisions to uncertainty about key assumptions.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

Neutral and selection-driven decay of sexual traits in asexual stick insects.

Environmental shifts and lifestyle changes may result in formerly adaptive traits becoming non-functional or maladaptive. The subsequent decay of such traits highlights the importance of natural selection for adaptations, yet its causes have rarely been investigated. To study the fate of formerly adaptive traits after lifestyle changes, we evaluated sexual traits in five independently derived asexual lineages, including traits that are specific to males and therefore not exposed to selection. At least four of the asexual lineages retained the capacity to produce males that display normal courtship behaviours and are able to fertilize eggs of females from related sexual species. The maintenance of male traits may stem from pleiotropy, or from these traits only regressing via drift, which may require millions of years to generate phenotypic effects. By contrast, we found parallel decay of sexual traits in females. Asexual females produced altered airborne and contact signals, had modified sperm storage organs, and lost the ability to fertilize their eggs, impeding reversals to sexual reproduction. Female sexual traits were decayed even in recently derived asexuals, suggesting that trait changes following the evolution of asexuality, when they occur, proceed rapidly and are driven by selective processes rather than drift.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Geochemical evolution of active porphyry copper tailings impoundments : from alkaline deposition towards acidification

J. Smuda: Geochemical evolution of active porphyry copper tailings impoundments Thesis abstract Mine waste is the largest volume of materials handled in the world. The oxidation of sulfidic mine waste may result in the release of acid mine drainage (AMD) rich in heavy metals and arsenic to the environment, one of the major problems the mining industry is facing today. To control and reduce this environmental impact, it is crucial to identify the main geochemical and hydrological processes influencing contaminant liberation, transport, and retention. This thesis presents the results of a geochemical, mineralogical and stable isotope study (δ2H, δ18O, δ34S) from two active porphyry copper tailings impoundments in Mediterranean (Carén tailings impoundment, El Teniente mine, Central Chile) and hyper-arid climate (Talabre tailings impoundment, Chuquicamata, Northern Chile) from the deposition in alkaline environment (pH 10.5) towards acidification after several years of exposure. The major hydrological results were the identification of vertical contaminant and water transport in the uppermost, not water-saturated zone, triggered by capillary rise due to evaporation, and infiltration downwards due to new tailings deposition, and of horizontal transport in the groundwater zone. At the surface of the sedimented tailings, evaporation of pore water led to the precipitation of Na-Ca-Mg sulfates (e.g., gypsum, tenorite), in hyper-arid climate also halite. At the Carén tailings impoundment, renewed deposition in a 4-week interval inhibited a pH decrease below neutral values and the formation of an efflorescent salt crust. At the Talabre tailings impoundment, deposition breaks of several years resulted in the formation of acidic oxidation zones in the timeframe of less than 4 years. This process enabled the transport of liberated Cu, Zn, and Fe via capillary rise to the surface, where these metals precipitated as heavy-metal sulfates (e.g., devilline, krohnkite) and chlorides (eriochalcite, atacamite). Renewed depositing may dissolve efflorescent salts and transport liberated elements towards the groundwater zone. This zone was found to be highly dynamic due to infiltration and mixing with water from different sources, like groundwater, catchment water, and infiltration from superficial waters. There, Cu was found to be partially mobile due to complexation with Cl (in Cl-rich groundwater, Talabre) and dissolved organic matter (in zones with infiltration of catchment water rich in dissolved organic matter, Carén). A laboratory study on the isotopic fractionation of sulfur and oxygen of sulfate in different minerals groups (water-soluble sulfates, low- and high-crystalline Fe(III) oxyhydroxides) contributed to the use of stable isotopes as tracer of geochemical and transport processes for environmental studies. The results highlight that a detailed geochemical, stable isotope and mineralogical study permits the identification of contamination processes and pathways already during the deposition of mine tailings. This knowledge allows the early planning of adequate actions to reduce and control the environmental impact during tailings deposition and after the closing of the impoundment. J. Smuda: Geochemical evolution of active porphyry copper tailings impoundments Résumé de these Les déchets miniers constituent les plus grands volumes de matériel gérés dans le monde. L'oxydation des déchets miniers sulfuriques peut conduire à la libération de drainages miniers acides (DMA) riches en métaux et arsenic dans l'environnement, ce qui est l'un des principaux problèmes de l'industrie minière aujourd'hui. Pour contrôler et réduire ces impacts sur l'environnement, il est crucial d'identifier les principaux processus géochimiques et hydrologiques influençant la libération, le transport et la rétention des contaminants. Cette thèse présente les résultats d'une étude géochimique, minéralogique et des isotopes stables (δ2H, δ18O, δ34S) sur des déchets miniers de 2 sites de dépôt actifs en climat méditerranéen (Dépôt de déchets de Carén, mine de El Teniente, Centre du Chili) et en climat hyper-aride (Dépôt de déchets de Talabre, mine de Chuquicamata, Nord du Chili). L'objectif était d'étudier l'évolution des déchets de la déposition en milieu alcalin (pH = 10.5) vers l'acidification après plusieurs années d'exposition. Le principal résultat hydrologique a été l'identification de 2 types de transport : un transport vertical de l'eau et des contaminants dans la zone non saturée en surface, induit par la montée capillaire due à l'évaporation et par l'infiltration subséquente de la déposition de sédiments frais ; et un transport horizontal dans la zone des eaux souterraines. À la surface des déchets, l'évaporation de l'eau interstitielle conduit à la précipitation de sulfates de Na-Ca-Mg (ex. gypse, ténorite) et halite en climat hyper-aride. Dans le site de Carén, une nouvelle déposition de déchets frais à 4 semaines intervalle a empêché la baise du pH en deçà des valeurs neutres et la formation d'une croûte de sels efflorescentes en surface. Dans le site de Talabre, les fentes de dessiccation des dépôts ont entraîné la formation d'une zone d'oxydation à pH acide en moins de 4 ans. Ce processus a permis la libération et le transport par capillarité de Cu, Zn, Fe vers la surface, où ces éléments précipitent sous forme de sulfates de métaux lourds (ex., dévilline, krohnkite) de chlorures (ex. ériochalcite, atacamite). Une nouvelle déposition de sédiments frais pourrait dissoudre ces sels et les transporter vers la zone des eaux souterraines. Cette dernière zone était très dynamique en raison du mélange d'eaux provenant de différentes sources, comme les eaux souterraines, l'eau de captage et l'infiltration des eaux superficielles. Egalement dans cette zone, le cuivre était partiellement mobile à cause de la formation de complexe avec le chlore (dans les zone riche en Cl, Talabre) et avec la matière organique dissoute (dans les zones où s'infiltre l'eau de captage riche en matière organique, Carén). Une étude en laboratoire sur le fractionnement des isotopes stables de sulfure et d'oxygène des sulfates dans différents groupes de minéraux (sulfates hydrosolubles, sulfures de oxy-hydroxyde de Fe(III) faiblement ou fortement cristallins) a permis d'apporter une contribution à leur utilisation comme traceurs dans l'étude des processus géochimiques et de transport lors d'études environnementales. Les résultats montrent qu'une étude détaillée de la géochimie, des isotopes stables et de la minéralogie permet d'identifier les processus et les voies de contamination déjà pendant la période de dépôt des déchets miniers. Cette connaissance permet de planifier, dès le début de l'exploitation, des mesures adéquates pour réduire et contrôler l'impact sur l'environnement pendant la période de dépôts de déchets miniers et après la fermeture du site.

Late Cretaceous to Eocene geology of the South Central American forearc area (southern Costa Rica and western Panama) : initiation and evolution of an intra-oceanic convergent margin

Résumé : L'arc volcanique du sud de l'Amérique Centrale se situe sur la marge SW de la Plaque Caraïbe, au-dessus des plaques subduites de Cocos et Nazca. Il s'agit de l'un des arcs intra-océaniques les plus étudiés au monde, qui est généralement considéré comme s'étant développé à la fin du Crétacé le long d'un plateau océanique (le Plateau Caraïbe ou CLIP) et se trouvant actuellement dans un régime de subduction érosive. Au cours des dernières décennies, des efforts particuliers ont été faits pour comprendre les processus liés à la subduction sur la base d'études géophysiques et géochimiques. Au sud du Costa Rica et à l'ouest du Panama, des complexes d'accrétions et structures à la base de l'arc volcanique ont été exposés grâce à la subduction de rides asismiques et de failles transformantes. Des affleurements, situés jusqu'à seulement 15 km de la fosse, offrent une possibilité unique de mieux comprendre quelques uns des processus ayant lieu le long de la zone de subduction. Nous présentons de nouvelles contraintes sur l'origine de ces affleurements en alliant une étude de terrain poussée, de nouvelles données géochimiques, sédimentaires et paléontologiques, ainsi que des observations structurales effectuées en télédétection. Une nouvelle stratigraphie tectonique entre le Campanien et l'Éocène est définie pour la région d'avant-arc située entre la Péninsule d'Osa (Costa Rica) et la Péninsule d'Azuero (Panama). Nos résultats montrent que la partie externe de la marge est composée d'un arrangement complexe de roches ignées et de séquences sédimentaires de recouvrement qui comprennent principalement le socle de l'arc, des roches d'arc primitif, des fragments de monts sous-marins accrétés et des mélanges d'accrétion. Des preuves sont données pour le développement de l'arc volcanique du sud de l'Amérique Centrale sur un plateau océanique. Le début de la subduction le long de la marge SW de la Plaque Caraïbe a eu lieu au Campanien et a généré des roches d'arc primitif caractérisées par des affinités géochimiques particulières, globalement intermédiaires entre des affinités de plateau et d'arc insulaire. L'arc était mature au Maastrichtien et formait un isthme essentiellement continu entre l'Amérique du Nord et l'Amérique du Sud. Ceci a permis la migration de faunes terrestres entre les Amériques et pourrait avoir contribué à la crise fin Crétacé -Tertiaire en réduisant les courants océaniques subéquatoriaux entre le Pacifique et l'Atlantique. Plusieurs unités composées de fragments de monts sous-marins accrétés sont définies. La nature et l'arrangement structural de ces unités définissent de nouvelles contraintes sur les modes d'accrétion des monts sous-marins/îles océaniques et sur l'évolution de la marge depuis la formation de la zone de subduction. Entre la fin du Crétacé et l'Éocène moyen, la marge a enregistré plusieurs épisodes ponctuels d'accrétion de monts sous-marins alternant avec de la subduction érosive. A l'Éocène moyen, un événement tectonique régional pourrait avoir causé un fort couplage entre les plaques supérieure et inférieure, menant à des taux plus important d'accrétion de monts sous-marins. Durant cette période, la situation le long de la marge était très semblable à la situation actuelle et caractérisée par la présence de monts sous-marins subductants et l'absence d'accrétion de sédiments. L'enregistrement géologique montre qu'il n'est pas possible d'attribuer une nature érosive ou accrétionnaire à la marge dans le passé ou -par analogie- aujourd'hui, parce que (1) les processus d'accrétion et érosifs varient fortement spatialement et temporellement et (2) il est impossible d'évaluer la quantité exacte de matériel tectoniquement enlevé à la marge depuis le début de la subduction. Au sud du Costa Rica, certains fragments de monts sous-marins accrétés sont représentatifs d'une interaction entre une ride et un point chaud dans le Pacifique au Crétacé terminal/Paléocène. L'existence de ces fragments de monts sous-marins et la morphologie du fond de l'Océan Pacifique indiquent que la formation de la ride de Cocos-Nazca s'est formée au moins ~40 Ma avant l'âge proposé par les modèles tectoniques actuels. Au Panama, nous avons identifié une île océanique d'âge début Éocène qui a été accrétée à l'Éocène moyen. L'accrétion a eu lieu à très faible profondeur par détachement de l'île dans la fosse, et a mené à une exceptionnelle préservation des structures volcaniques. Des affleurement comprenant aussi bien des parties basses et hautes de l'édifice volcanique on été étudiées, depuis la phase sous-marine bouclier jusqu'à la phase subaérienne post-bouclier. La stratigraphie nous a permis de différencier les laves de la phase sous-marine de celles de la phase subaérienne. La composition des laves indique une diminution progressive de l'intensité de la fusion partielle de la source et une diminution de la température des laves produites durant les derniers stades de l'activité volcanique. Nous interprétons ces changements comme étant liés à l'éloignement progressif de l'île océanique de la zone de fusion ou point chaud. Abstract The southern Central American volcanic front lies on the SW edge of the Caribbean Plate, inboard of the subducting Cocos and Nazca Plates. It is one of the most studied intra-oceanic convergent margins around the world, which is generally interpreted to have developed in the late Cretaceous along an oceanic plateau (the Caribbean Large Igneous Province or CLIP) and to be currently undergoing a regime of subduction erosion. In the last decades a particular effort has been made to understand subduction-related processes on the basis of geophysical and geochemical studies. In southern Costa Rica and western Panama accretionary complexes and structures at the base of the volcanic front have been exposed in response to subduction of aseismic ridges and transforms. Onland exposures are located as close as to 15 km from the trench and provide a unique opportunity to better understand some of the processes occurring along the subduction zone. We provide new constraints on the origins of these exposures by integrating a comprehensive field work, new geochemical, sedimentary and paleontological data, as well as structural observations based on remote imaging. A new Campanian to Eocene tectonostratigraphy is defined for the forearc area located between the Osa Peninsula (Costa Rica) and the Azuero Peninsula (Panama). Our results show that the outer margin is composed of a complicated arrangement of igneous complexes and overlapping sedimentary sequences that essentially comprise an arc basement, primitive island-arc rocks, accreted seamount fragments and accretionary mélanges. Evidences are provided for the development of the southern Central American arc on the top an oceanic plateau. The subduction initiation along the SW edge of the Caribbean Plate occurred in the Campanian and led to formation of primitive island-arc rocks characterized by unusual geochemical affinities broadly intermediate between plateau and arc affinities. The arc was mature in the Maastrichtian and was forming a predominantly continuous landbridge between the North and South Americas. This allowed migration of terrestrial fauna between the Americas and may have contributed to the Cretaceous-Tertiary crisis by limiting trans-equatorial oceanic currents between the Pacific and the Atlantic. Several units composed of accreted seamount fragments are defined. The nature of the units and their structural arrangement provide new constraints on the modes of accretion of seamounts/oceanic islands and on the evolution of the margin since subduction initiation. Between the late Cretaceous and the middle Eocene, the margin recorded several local episodes of seamount accretion alternating with tectonic erosion. In the middle Eocene a regional tectonic event may have triggered strong coupling between the overriding and subducting plates, leading to higher rates of seamount accretion. During this period the situation along the margin was very similar to the present and characterized by subducting seamounts and absence of sediment accretion. The geological record shows that it is not possible to ascribe an overall erosive or accretionary nature to the margin in the past and, by analogy, today, because (1) accretionary and erosive processes exhibit significant lateral and temporal variations and (2) it is impossible to estimate the exact amount of material tectonically eroded from the margin since subduction initiation. In southern Costa Rica, accreted seamount fragments point toward a plume-ridge interaction in the Pacific in the late Cretaceous/Paleocene. This occurrence of accreted seamount fragments and morphology of the Pacific Ocean floor is indicative of the formation of the Cocos-Nazca spreading system at least ~40 Ma prior to the age proposed in current tectonic models. In Panama, we identified a remarkably-well preserved early Eocene oceanic island that accreted in the middle Eocene. The accretion probably occurred at very shallow depth by detachment of the island in the trench and led to an exceptional preservation of the volcanic structures. Exposures of both deep and superficial parts of the volcanic edifice have been studied, from the submarine-shield to subaerial-postshield stages. The stratigraphy allowed us to distinguish lavas produced during the submarine and subaerial stages. The lava compositions likely define a progressive diminution of source melting and a decrease in the temperature of erupted melts in the latest stages of volcanic activity. We interpret these changes to primarily reflect the progressive migration of the oceanic island out of the melting region or hotspot.