Correcting for the bias due to expression specificity improves the estimation of constrained evolution of expression between mouse and human.

MOTIVATION: Comparative analyses of gene expression data from different species have become an important component of the study of molecular evolution. Thus methods are needed to estimate evolutionary distances between expression profiles, as well as a neutral reference to estimate selective pressure. Divergence between expression profiles of homologous genes is often calculated with Pearson's or Euclidean distance. Neutral divergence is usually inferred from randomized data. Despite being widely used, neither of these two steps has been well studied. Here, we analyze these methods formally and on real data, highlight their limitations and propose improvements. RESULTS: It has been demonstrated that Pearson's distance, in contrast to Euclidean distance, leads to underestimation of the expression similarity between homologous genes with a conserved uniform pattern of expression. Here, we first extend this study to genes with conserved, but specific pattern of expression. Surprisingly, we find that both Pearson's and Euclidean distances used as a measure of expression similarity between genes depend on the expression specificity of those genes. We also show that the Euclidean distance depends strongly on data normalization. Next, we show that the randomization procedure that is widely used to estimate the rate of neutral evolution is biased when broadly expressed genes are abundant in the data. To overcome this problem, we propose a novel randomization procedure that is unbiased with respect to expression profiles present in the datasets. Applying our method to the mouse and human gene expression data suggests significant gene expression conservation between these species. CONTACT: marc.robinson-rechavi@unil.ch; sven.bergmann@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Tissue-Specific Evolution of Protein Coding Genes in Human and Mouse.

Protein-coding genes evolve at different rates, and the influence of different parameters, from gene size to expression level, has been extensively studied. While in yeast gene expression level is the major causal factor of gene evolutionary rate, the situation is more complex in animals. Here we investigate these relations further, especially taking in account gene expression in different organs as well as indirect correlations between parameters. We used RNA-seq data from two large datasets, covering 22 mouse tissues and 27 human tissues. Over all tissues, evolutionary rate only correlates weakly with levels and breadth of expression. The strongest explanatory factors of purifying selection are GC content, expression in many developmental stages, and expression in brain tissues. While the main component of evolutionary rate is purifying selection, we also find tissue-specific patterns for sites under neutral evolution and for positive selection. We observe fast evolution of genes expressed in testis, but also in other tissues, notably liver, which are explained by weak purifying selection rather than by positive selection.

THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Worldwide Population Structure, Long-Term Demography, and Local Adaptation of Helicobacter pylori.

Helicobacter pylori is an important human pathogen associated with serious gastric diseases. Owing to its medical importance and close relationship with its human host, understanding genomic patterns of global and local adaptation in H. pylori may be of particular significance for both clinical and evolutionary studies. Here we present the first such whole genome analysis of 60 globally distributed strains, from which we inferred worldwide population structure and demographic history and shed light on interesting global and local events of positive selection, with particular emphasis on the evolution of San-associated lineages. Our results indicate a more ancient origin for the association of humans and H. pylori than previously thought. We identify several important perspectives for future clinical research on candidate selected regions that include both previously characterized genes (e.g., transcription elongation factor NusA and tumor necrosis factor alpha-inducing protein Tipα) and hitherto unknown functional genes.

Developmental and Anatomical Constraints on Vertebrate Genome Evolution

Pendant ma thèse de doctorat, j'ai utilisé des espèces modèles, comme la souris et le poisson-zèbre, pour étudier les facteurs qui affectent l'évolution des gènes et leur expression. Plus précisément, j'ai montré que l'anatomie et le développement sont des facteurs clés à prendre en compte, car ils influencent la vitesse d'évolution de la séquence des gènes, l'impact sur eux de mutations (i.e. la délétion du gène est-elle létale ?), et leur tendance à se dupliquer. Où et quand il est exprimé impose à un gène certaines contraintes ou au contraire lui donne des opportunités d'évoluer. J'ai pu comparer ces tendances aux modèles classiques d'évolution de la morphologie, que l'on pensait auparavant refléter directement les contraintes s'appliquant sur le génome. Nous avons montré que les contraintes entre ces deux niveaux d'organisation ne peuvent pas être transférées simplement : il n'y a pas de lien direct entre la conservation du génotype et celle de phénotypes comme la morphologie. Ce travail a été possible grâce au développement d'outils bioinformatiques. Notamment, j'ai travaillé sur le développement de la base de données Bgee, qui a pour but de comparer l'expression des gènes entre différentes espèces de manière automatique et à large échelle. Cela implique une formalisation de l'anatomie, du développement et de concepts liés à l'homologie grâce à l'utilisation d'ontologies. Une intégration cohérente de données d'expression hétérogènes (puces à ADN, marqueurs de séquence exprimée, hybridations in situ) a aussi été nécessaire. Cette base de données est mise à jour régulièrement et disponible librement. Elle devrait contribuer à étendre les possibilités de comparaison de l'expression des gènes entre espèces pour des études d'évo-devo (évolution du développement) et de génomique. During my PhD, I used model species of vertebrates, such as mouse and zebrafish, to study factors affecting the evolution of genes and their expression. More precisely I have shown that anatomy and development are key factors to take into account, influencing the rate of gene sequence evolution, the impact of mutations (i.e. is the deletion of a gene lethal?), and the propensity of a gene to duplicate. Where and when genes are expressed imposes constraints, or on the contrary leaves them some opportunity to evolve. We analyzed these patterns in relation to classical models of morphological evolution in vertebrates, which were previously thought to directly reflect constraints on the genomes. We showed that the patterns of evolution at these two levels of organization do not translate smoothly: there is no direct link between the conservation of genotype and phenotypes such as morphology. This work was made possible by the development of bioinformatics tools. Notably, I worked on the development of the database Bgee, which aims at comparing gene expression between different species in an automated and large-scale way. This involves the formalization of anatomy, development, and concepts related to homology, through the use of ontologies. A coherent integration of heterogeneous expression data (microarray, expressed sequence tags, in situ hybridizations) is also required. This database is regularly updated and freely available. It should contribute to extend the possibilities for comparison of gene expression between species in evo-devo and genomics studies.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Developmental constraints on vertebrate genome evolution.

Constraints in embryonic development are thought to bias the direction of evolution by making some changes less likely, and others more likely, depending on their consequences on ontogeny. Here, we characterize the constraints acting on genome evolution in vertebrates. We used gene expression data from two vertebrates: zebrafish, using a microarray experiment spanning 14 stages of development, and mouse, using EST counts for 26 stages of development. We show that, in both species, genes expressed early in development (1) have a more dramatic effect of knock-out or mutation and (2) are more likely to revert to single copy after whole genome duplication, relative to genes expressed late. This supports high constraints on early stages of vertebrate development, making them less open to innovations (gene gain or gene loss). Results are robust to different sources of data -- gene expression from microarrays, ESTs, or in situ hybridizations; and mutants from directed KO, transgenic insertions, point mutations, or morpholinos. We determine the pattern of these constraints, which differs from the model used to describe vertebrate morphological conservation ("hourglass" model). While morphological constraints reach a maximum at mid-development (the "phylotypic" stage), genomic constraints appear to decrease in a monotonous manner over developmental time.

Assessing Recent Selection and Functionality at Long Noncoding RNA Loci in the Mouse Genome.

Long noncoding RNAs (lncRNAs) are one of the most intensively studied groups of noncoding elements. Debate continues over what proportion of lncRNAs are functional or merely represent transcriptional noise. Although characterization of individual lncRNAs has identified approximately 200 functional loci across the Eukarya, general surveys have found only modest or no evidence of long-term evolutionary conservation. Although this lack of conservation suggests that most lncRNAs are nonfunctional, the possibility remains that some represent recent evolutionary innovations. We examine recent selection pressures acting on lncRNAs in mouse populations. We compare patterns of within-species nucleotide variation at approximately 10,000 lncRNA loci in a cohort of the wild house mouse, Mus musculus castaneus, with between-species nucleotide divergence from the rat (Rattus norvegicus). Loci under selective constraint are expected to show reduced nucleotide diversity and divergence. We find limited evidence of sequence conservation compared with putatively neutrally evolving ancestral repeats (ARs). Comparisons of sequence diversity and divergence between ARs, protein-coding (PC) exons and lncRNAs, and the associated flanking regions, show weak, but significantly lower levels of sequence diversity and divergence at lncRNAs compared with ARs. lncRNAs conserved deep in the vertebrate phylogeny show lower within-species sequence diversity than lncRNAs in general. A set of 74 functionally characterized lncRNAs show levels of diversity and divergence comparable to PC exons, suggesting that these lncRNAs are under substantial selective constraints. Our results suggest that, in mouse populations, most lncRNA loci evolve at rates similar to ARs, whereas older lncRNAs tend to show signals of selection similar to PC genes.

Why evolution matters for species conservation: perspectives from three case studies of plant metapopulations

We advocate the advantage of an evolutionary approach to conservation biology that considers evolutionary history at various levels of biological organization. We review work on three separate plant taxa, spanning from one to multiple decades, illustrating extremes in metapopulation functioning. We show how the rare endemics Centaurea corymbosa (Clape Massif, France) and Brassica insularis in Corsica (France) may be caught in an evolutionary trap: disruption of metapopulation functioning due to lack of colonization of new sites may have counterselected traits such as dispersal ability or self-compatibility, making these species particularly vulnerable to any disturbance. The third case study concerns the evolution of life history strategies in the highly diverse genus Leucadendron of the South African fynbos. There, fire disturbance and the recolonization phase after fires are so integral to the functioning of populations that recruitment of new individuals is conditioned by fire. We show how past adaptation to different fire regimes and climatic constraints make species with different life history syndromes more or less vulnerable to global changes. These different case studies suggest that management strategies should promote evolutionary potential and evolutionary processes to better protect extant biodiversity and biodiversification.

Oxygen isotope compositions of iron oxides from high-grade BIF-hosted iron ore deposits of the Central Hamersley Province, Western Australia: Constraints on the evolution of hydrothermal fluids

The Hamersley province of northwest Australia is one of the world's premier iron ore regions with high-grade martite-microplaty hematite iron ore deposits mostly hosted within banded iron formation (BIF) sequences of the Brockman Iron Formations of the Hamersley Group. These high-grade iron ores contain between 60 and 68 wt percent Fe, and formed by the multistage interaction of hydrothermal fluids with the host BIF formation. The oxygen isotope compositions of magnetite and hematite from BIF, hydrothermal alteration assemblages, and high-grade iron Ore were analyzed from the Mount Tom Price, Paraburdoo, and Charmar iron ore deposits. The delta(18)O values of magnetite and hematite from hydrothermal alteration assemblages and high-grade iron ore range from -9.0 to -2.9 per mil, a depletion of 5 to 15 per mil relative to the host BIF. The delta(18)O values are spatially controlled by faults within the deposits, a response to higher fluid flux and larger influence the isotopic compositions by the hydrothermal fluids. The oxygen isotope composition of hydrothermal fluids (delta(18)O(fluid)) indicates that the decrease in the (18)O content of iron oxides was due to the interaction of both basinal brines and meteoric fluids with the original BIF. Late-stage talc-bearing ore at the Mount Tom Price deposit formed in the presence of a pulse of delta(18)O-enriched basinal brine, indicating that hydrothermal fluids may have repeatedly interacted with the BIFs during the Paleoproterozoic.

Origins and functional evolution of Y chromosomes across mammals.

Y chromosomes underlie sex determination in mammals, but their repeat-rich nature has hampered sequencing and associated evolutionary studies. Here we trace Y evolution across 15 representative mammals on the basis of high-throughput genome and transcriptome sequencing. We uncover three independent sex chromosome originations in mammals and birds (the outgroup). The original placental and marsupial (therian) Y, containing the sex-determining gene SRY, emerged in the therian ancestor approximately 180 million years ago, in parallel with the first of five monotreme Y chromosomes, carrying the probable sex-determining gene AMH. The avian W chromosome arose approximately 140 million years ago in the bird ancestor. The small Y/W gene repertoires, enriched in regulatory functions, were rapidly defined following stratification (recombination arrest) and erosion events and have remained considerably stable. Despite expression decreases in therians, Y/W genes show notable conservation of proto-sex chromosome expression patterns, although various Y genes evolved testis-specificities through differential regulatory decay. Thus, although some genes evolved novel functions through spatial/temporal expression shifts, most Y genes probably endured, at least initially, because of dosage constraints.

Constraints on the late thermotectonic evolution of the Western Alps - Evidence for episodic rapid uplift

Forty-two new apatite and zircon fission track ages are presented for samples from the Western Alps in southern Switzerland, northern Italy, and southeastern France. Measured ages plotted against assumed closure temperatures yield cooling patterns for the final cooling, uplift, and exhumation of the Western Alps. Similar fission track zircon ages in the Penninic Gran Paradiso massif, Dent Blanche nappe, Sesia-Lanzo Zone, and Ivrea Zone indicate cooling of all four units to approximately 225-degrees-C by 33 Ma. Differences in apatite ages reveal differential cooling of the four blocks between 33 Ma and the present. In the Sesia-Lanzo Zone, similarity of apatite ages regardless of elevation, together with near-volcanic confined fission track length patterns suggest rapid cooling and uplift at approximately 25 Ma compared with slow cooling of other Western Alps units around 12 Ma. Uplift is thus not continuous but episodic, often over a short time interval beyond the resolution of other methods. Such episodes of uplift, as revealed here in the Sesia-Lanzo Zone, may be the rule rather than the exception.

Alpine tectonic evolution of a Jurassic subduction-accretionary complex: Deformation, kinematics and 40Ar/39Ar age constraints on the Mesozoic low-grade schists of the Circum-Rhodope Belt in the eastern Rhodope-Thrace region, Bulgaria-Greece

Deformation of the Circum-Rhodope Belt Mesozoic (Middle Triassic to earliest Lower Cretaceous) low-grade schists underneath an arc-related ophiolitic magmatic suite and associated sedimentary successions in the eastern Rhodope-Thrace region occurred as a two-episode tectonic process: (i) Late Jurassic deformation of arc to margin units resulting from the eastern Rhodope-Evros arc-Rhodope terrane continental margin collision and accretion to that margin, and (ii) Middle Eocene deformation related to the Tertiary crustal extension and final collision resulting in the closure of the Vardar ocean south of the Rhodope terrane. The first deformational event D-1 is expressed by Late Jurassic NW-N vergent fold generations and the main and subsidiary planar-linear structures. Although overprinting, these structural elements depict uniform bulk north-directed thrust kinematics and are geometrically compatible with the increments of progressive deformation that develops in same greenschist-facies metamorphic grade. It followed the Early-Middle Jurassic magmatic evolution of the eastern Rhodope-Evros arc established on the upper plate of the southward subducting Maliac-Meliata oceanic lithosphere that established the Vardar Ocean in a supra-subduction back-arc setting. This first event resulted in the thrust-related tectonic emplacement of the Mesozoic schists in a supra-crustal level onto the Rhodope continental margin. This Late Jurassic-Early Cretaceous tectonic event related to N-vergent Balkan orogeny is well-constrained by geochronological data and traced at a regional-scale within distinct units of the Carpatho-Balkan Belt. Following subduction reversal towards the north whereby the Vardar Ocean was subducted beneath the Rhodope margin by latest Cretaceous times, the low-grade schists aquired a new position in the upper plate, and hence, the Mesozoic schists are lacking the Cretaceous S-directed tectono-metamorphic episode whose effects are widespread in the underlying high-grade basement. The subduction of the remnant Vardar Ocean located behind the colliding arc since the middle Cretaceous was responsible for its ultimate closure, Early Tertiary collision with the Pelagonian block and extension in the region caused the extensional collapse related to the second deformational event D-2. This extensional episode was experienced passively by the Mesozoic schists located in the hanging wall of the extensional detachments in Eocene times. It resulted in NE-SW oriented open folds representing corrugation antiforms of the extensional detachment surfaces, brittle faulting and burial history beneath thick Eocene sediments as indicated by 42.1-39.7 Ma Ar-40/Ar-39 mica plateau ages obtained in the study. The results provide structural constraints for the involvement components of Jurassic paleo-subduction zone in a Late Jurassic arc-continental margin collisional history that contributed to accretion-related crustal growth of the Rhodope terrane. (C) 2011 Elsevier Ltd. All rights reserved.

Constraints on palaeodrainage evolution induced by uplift and exhumation on the southern flank of the Zagros-Iranian Plateau

Foreland sedimentary rocks from the northern Fars region of Iran contain a record of deformation associated with the Cenozoic collision between Arabia and Eurasia that resulted in formation of the Zagros orogen. The timing of the deformation associated with this event is poorly known. To address this we conducted a study of Miocene foreland sedimentary rocks (19.7-14.8 Ma) of the Chahar-Makan syncline using clast composition, clay mineralogy and low-temperature fission-track dating. The results showed that most of the sedimentary rocks were sourced from ophiolitic rocks. Detrital apatite fission-track (AFT) age signatures of Miocene sedimentary rocks record exhumation in the hanging wall of the Main Zagros Thrust and confirm that the change from underthrusting of the stretched Arabian margin to widespread crustal thickening and deformation in the Zagros region is no younger than 19.7 Ma. A transition from Late Oligocene to Mesozoic-Eocene AFT detrital age signatures between 19.7-16.6 Ma and 16.6-13.8 Ma is interpreted to reflect a possible rearrangement of palaeodrainage distribution that resulted from folding and expansion-uplift of the Zagros-Iranian Plateau region.