The biogeography of mitochondrial and nuclear discordance in animals.

Combining nuclear (nuDNA) and mitochondrial DNA (mtDNA) markers has improved the power of molecular data to test phylogenetic and phylogeographic hypotheses and has highlighted the limitations of studies using only mtDNA markers. In fact, in the past decade, many conflicting geographic patterns between mitochondrial and nuclear genetic markers have been identified (i.e. mito-nuclear discordance). Our goals in this synthesis are to: (i) review known cases of mito-nuclear discordance in animal systems, (ii) to summarize the biogeographic patterns in each instance and (iii) to identify common drivers of discordance in various groups. In total, we identified 126 cases in animal systems with strong evidence of discordance between the biogeographic patterns obtained from mitochondrial DNA and those observed in the nuclear genome. In most cases, these patterns are attributed to adaptive introgression of mtDNA, demographic disparities and sex-biased asymmetries, with some studies also implicating hybrid zone movement, human introductions and Wolbachia infection in insects. We also discuss situations where divergent mtDNA clades seem to have arisen in the absence of geographic isolation. For those cases where foreign mtDNA haplotypes are found deep within the range of a second taxon, data suggest that those mtDNA haplotypes are more likely to be at a high frequency and are commonly driven by sex-biased asymmetries and/or adaptive introgression. In addition, we discuss the problems with inferring the processes causing discordance from biogeographic patterns that are common in many studies. In many cases, authors presented more than one explanation for discordant patterns in a given system, which indicates that likely more data are required. Ideally, to resolve this issue, we see important future work shifting focus from documenting the prevalence of mito-nuclear discordance towards testing hypotheses regarding the drivers of discordance. Indeed, there is great potential for certain cases of mitochondrial introgression to become important natural systems within which to test the effect of different mitochondrial genotypes on whole-animal phenotypes.

PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

What animals can teach clinicians about the hippocampus

Abnormalities in hippocampal structure and function have been reported in a number of human neuropathological and neurodevelopmental disorders, including Alzheimer's disease, autism spectrum disorders, Down syndrome, epilepsy, and schizophrenia. Given the complexity of these disorders, animal studies are invaluable and remain to date irreplaceable, providing fundamental knowledge regarding the basic mechanisms underlying normal and pathological human brain structure and function. However, there is a prominent ill-conceived view in current research that scientists should be restricted to using animal models of human diseases that can lead to results applicable to humans within a few years. Although there is no doubt that translational studies of this kind are important and necessary, limiting animal studies to applicable questions is counterproductive and will ultimately lead to a lack of knowledge and an inability to address human health problems. Here, we discuss findings regarding the normal postnatal development of the monkey hippocampal formation, which provide an essential framework to consider the etiologies of different neuropathological disorders affecting human hippocampal structure and function. We focus on studies of gene expression in distinct hippocampal regions that shed light on some basic mechanisms that might contribute to the etiology of schizophrenia. We argue that researchers, as well as clinicians, should not consider the use of animals in research only as 'animal models' of human diseases, as they will continue to need and benefit from a better understanding of the normal structure and functions of the hippocampus in 'model animals'.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Bayesian estimation of speciation and extinction from incomplete fossil occurrence data.

The temporal dynamics of species diversity are shaped by variations in the rates of speciation and extinction, and there is a long history of inferring these rates using first and last appearances of taxa in the fossil record. Understanding diversity dynamics critically depends on unbiased estimates of the unobserved times of speciation and extinction for all lineages, but the inference of these parameters is challenging due to the complex nature of the available data. Here, we present a new probabilistic framework to jointly estimate species-specific times of speciation and extinction and the rates of the underlying birth-death process based on the fossil record. The rates are allowed to vary through time independently of each other, and the probability of preservation and sampling is explicitly incorporated in the model to estimate the true lifespan of each lineage. We implement a Bayesian algorithm to assess the presence of rate shifts by exploring alternative diversification models. Tests on a range of simulated data sets reveal the accuracy and robustness of our approach against violations of the underlying assumptions and various degrees of data incompleteness. Finally, we demonstrate the application of our method with the diversification of the mammal family Rhinocerotidae and reveal a complex history of repeated and independent temporal shifts of both speciation and extinction rates, leading to the expansion and subsequent decline of the group. The estimated parameters of the birth-death process implemented here are directly comparable with those obtained from dated molecular phylogenies. Thus, our model represents a step towards integrating phylogenetic and fossil information to infer macroevolutionary processes.

Perspective: from molecular endocrinology of mouse to genomic endocrinology of animals.

This short perspective explores some ways in which new genomic methodologies impact the study of endocrine signaling. Emphasis is put on the impact of studying species which are not molecular biology models. This opens the door to using knowledge molecular endocrinology in areas of biology as distant as conservation biology, as well as enriching endocrinology with information from biodiversity and natural variation.

Transmission of MRSA and MSSA between humans and farm animals

Human nares are the main niche of Staphylococcus aureus, but farm animals can be also infected (cows) or colonized (pigs) constituting significant reservoir of this pathogen. Previous studies indicated that human and animal strains are quite distinct but the extent of cross-species specialization and transmission remains largely unknown. However, recent reports from several European countries as well as USA and Canada have indicated that employment in farming is an emerging risk factor for MRSA carriage. Pigs were found to be frequently colonized with MRSA, usually with a strain belonging to CC398. It is not known whether animal-human transmission was specific to this particular MRSA strain. S. aureus isolates from cow mastitis and pig colonization isolates were collected in parallel to nasal swab isolates from the animals' caretakers. The isolates were genotyped by AFLP, spatyping, and when appropriate by MLST. The isolates from cow mastitis were genetically uniform in comparison with human isolates. They were quite distinct from farmers\' carriage isolates, indicating pronounced hostspecialization. However, several cases where an infected cow and a colonized farmer had the same strain were detected, including one farm where two farmers were colonized and two cows were infected with MRSA belonging to CC398. Pig isolates were genetically more diverse than cow isolates. They were different from both human and cow isolates with one notable exception. Large fraction of pigs (20%) and pig caretakers (50%) were colonized with isolates belonging to CC398, majority of which were MSSA (2 cases of MRSA). These results indicate that host specialization in S. aureus is quite pronounced. Transmission between humans and farm animals was consequently quite rare. Both MSSA and MRSA strains belonging to otherwise pig-specific CC398 had increased capacity to colonize humans. Study of the genetic factors responsible for host specialization is underway.

Accessory phases as a tool to investigate magmatic underplating in polymetamorphic metasediments : an example from a fossil lower crust : Ivrea Verbano zne, NW Italy

The Ivrea and the Strona-Ceneri zones, NW italy and S Switzerland, offer the possibility to study the continental crust of the Southern Alps. Because of its high metamorphic degree and the abundant Permo- Carboniferous mafic intrusions, the Ivrea Zone is classically interpreted an exposed section trough the Permian lower crust. The present work is focused here on metasedimentary slices (septa) intercalated within Permian gabbro (mafic complex). In particular I studied the evolution of accessory phases such as rutile and zircon and the chemistry of the metasediments. The septa build an irregular and discontinuous band that cut obliquely the mafic complex from its deepest part (N) to its roof (S). The chemistry of the metasediments evolves along the band and the chemical evolution can be compared with that observed in the country-rock surrounding the mafic intrusion to the NE and overprinted by a main regional metamorphic event. This suggests that the degree of chemical depletion of the septa was mainly established during the same regional metamorphic event. Moreover it suggests that incorporation of the septa within the gabbro did not modify their original stratigraphie distribution within the crust. It implies that the mafic complex has been emplaced following a dynamic substantially different from the classic model of « gabbro glacier » (Quick et al., 1992; Quick et al., 1994). It is more likely that it has been emplaced by repeated injections of sills at different depths during a protracted period of time. Zircon trace elements and U-Pb ages suggest that regional metamorphism occurred 330-320Ma, the first sills in the deepest part of the Mafic Complex are injected at ~300Ma, the mafic magmas reached higher levels in the crust at 285Ma and the magmatic activity continued locally until 275Ma. The ages of detrital cores in zircons fix the maximal sedimentation age at ~370Ma, this age corresponds therefore with the maximal age of the incorporation of the Ivrea zone within the lower crust. I propose that the Ivrea zone has been accreted to the lower crust during the Hercynian orogeny sensu lato. The analysis of detrital ages suggests that the source terrains for the Ivrea zone and those for the Strona-Ceneri zone have a completely different Palaeozoic history. The systematic analysis of rutile in partially molten metasediments of the Ivrea zone reveals the occurrence of two generations. The two generations are characterized by a different chemistry and textural distribution. A first generation is formed during pro-grade metamorphism in the restitic counterpart. The second generation is formed in the melts during cooling at the same time that part of the first generation re-equilibrate. Re-equilibration of the first generation seems to be spatially controlled by the presence of fluids. Locally the second generation forms overgrowths on the first generation. Considered the different diffusivity of U and Pb in rutile, U heterogeneities have important implication for U-Pb dating of rutile. ID-TIMS and LA-ICPMS dating coupled with a careful textural investigation (SEM) suggest that rutile grains are characterized by multiple path along which Pb diffusion can occur: volume diffusion is an important process, but intragrain and subgrain boundaries provide additional high diffusivity pathways for Pb escape and reduce drastically the effective diffusion length. -- La zone d'Ivrea et la zone de Strona-Ceneri, en Italie nord-occidentale et Suisse méridionale, offrent la possibilité d'étudier la croûte continentale des Alpes du Sud. En raison du haut degré métamorphique et l'abondance d'intrusions mafiques d'âge Permo-Carbonifère [complexe mafique), la zone d'Ivrea est interprétée classiquement comme de la croûte inférieure permienne. Ce travail ce concentre sur des bandes metasédimentaires (septa) incorporées dans les magmas mafiques lors de l'intrusion. Les septa forment une bande irrégulière qui coupe obliquement le complexe mafique du bas (N) vers le haut (S). La chimie des septa évolue du bas vers le haut et l'évolution chimique se rapproche de l'évolution observé dans la roche encaissante l'intrusion affecté par un événement métamorphique régionale. Cette relation suggère que le degré d'appauvrissement chimique des septa a été établit principalement lors de l'événement métamorphique régional. De plus l'incorporation dans les gabbros n'a pas perturbée la distribution stratigraphique originelle des septa. Ces deux observations impliquent que le métamorphisme dans la roche encaissante précède la mise en place du gabbro et que cette dernière ne se fait pas selon le modèle classique (« gabbro glacier » de Quick et al., 1992, 1994), mais se fait plutôt par injections répétées de sills a différentes profondeurs. Les âges U-Pb et les éléments traces des zircons suggèrent que le métamorphisme régionale a eu lieu 330-320Ma, alors que les premiers sills dans la partie profonde du Mafic Complex s'injectent à ~300Ma, le magmatisme mafique atteigne des niveaux supérieurs à 285Ma et continue localement jusqu'à 270Ma. Les âges des coeurs détritiques des zircons permettent de fixer l'âge maximale de sédimentation à ~370Ma ce qui correspond donc à l'âge maximale de l'incorporation de la zone d'Ivrea dans la croûte inférieur. L'analyse systématique des rutiles, nous a permit de montrer l'existence de plusieurs générations qui ont une répartition texturale et une chimie différente. Une génération se forme lors de l'événement UHT dans les restites, une autre génération se forme dans les liquides lors du refroidissement, au même temps qu'une partie de la première génération se rééquilibre au niveau du Zr. Localement la deuxième génération peut former des surcroissances autour de la première génération. Dans ces cas, des fortes différences en uranium entre les deux générations ont des importantes implications pour la datation U-Pb sur rutile. Classiquement les ratios Pb/U dans le rutile sont interprétés comme indiquant l'âges du refroidissement du minéral sous une température à la quelle la diffusion du Pb dans le minéral n'est plus détectable et la diffusion à plus hautes températures est assumée se faire par «volume diffusion» dans le grain (Mezger et al., 1989). Par des datations ID-TIMS (sur grain entier) et LA-ICPMS (in-situ) et une analyse texturale (MEB) approfondie nous montrons que cette supposition est trop simpliste et que le rutile est repartie en sous-domaines. Chacun de ces domaines a ça propre longueur ou chemin de diffusion spécifique. Nous proposons donc une nouvelle approche plus cohérente pour l'interprétation des âges U-Pb sur rutile.