Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Substitution rates at neutral genes depend on population size under fluctuating demography and overlapping generations.

It is widely accepted that the rate of evolution (substitution rate) at neutral genes is unaffected by population size fluctuations. This result has implications for the analysis of genetic data in population genetics and phylogenetics, and provides, in particular, a justification for the concept of the molecular clock. Here, we show that the substitution rate at neutral genes does depend on population size fluctuations in the presence of overlapping generations. As both population size fluctuations and overlapping generations are expected to be the norm rather than the exception in natural populations, this observation may be relevant for understanding variation in substitution rates within and between lineages.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Evolution of life history traits in a natural population of the shrew Crocidura russula

Abstract Life history traits encompass all the decisions concerning fitness an individual is faced with during his life. The study of these traits is crucial to understand the factors shaping the biology of living organisms. Up until now, most of the information on the evolution of life history traits comes from laboratory studies. While these studies are interesting to test the effect of specific parameters, their conclusions are difficult to extrapolate to natural populations. Investigating the evolution of life history traits in natural populations is of great interest. This may be tricky because it requires information on reproduction, survival and morphology of individuals. Mark-recapture methods allow most of this information to be obtained. However, when direct observations of a species are not possible due to its ecology, indirect methods must be used to infer lifetime reproductive success. In this case, molecular markers are particularly helpful in assessing the genetic relationships between individuals and allow the construction of a pedigree. This thesis focuses on a natural population of a small insectivorous mammal, the greater white-toothed shrew, Crocidura russula. Because of its hidden lifestyle, the two complementary techniques mentioned above were combined to gather information on this population. The data were used to explore diverse aspects of evolutionary biology. We demonstrated that the high genetic variance displayed by the species was not maintained by its mating system because this shrew was less monogamous than previously thought. The large genetic diversity was most likely promoted by gene flow from the neighborhood. Dispersal was thus a central topic in this thesis. We showed that dispersal was not driven by inbreeding avoidance. In addition, we did not find any inbreeding depression in the population. Dispersal was promoted by a high number of vacant territories in the population for both sexes, meaning that territory acquisition played an important role in driving dispersal. Moreover, dispersal propensity was shown to have a genetic basis and, once achieved, to have no effect on individual fitness. Body mass was found to be a life history trait strongly influenced by sexual and viability selection in both sexes. Larger individuals had higher access to reproduction through territory acquisition and defense than lighter ones. By contrast, intermediate size individuals were favored by viability selection presumably because of ecological constraints and metabolic costs. Finally, we demonstrated that the majority of the life history traits in our shrew population has the potential to evolve because they maintained substantial amounts of additive genetic variance. Nonetheless, life history traits had no significant heritability due to their high level of nonadditive or environmental variance. Résumé Les traits d'histoire de vie comprennent toutes les décisions auxquelles un individu est confronté au cours de sa vie et qui concernent sa valeur adaptative. L'étude de ces traits est cruciale pour comprendre les facteurs qui façonnent la biologie des êtres vivants. Jusqu'à ce jour, la majorité des informations sur l'évolution des traits d'histoire de vie provient d'études réalisées en laboratoire. Alors que ces études sont intéressantes pour tester l'effet de paramètres spécifiques, leurs conclusions sont difficilement extrapolables aux populations naturelles. Il est particulièrement intéressant d'étudier l'évolution des traits d'histoire de vie dans des populations naturelles. Toutefois, ces études peuvent se révéler difficiles parce qu'elles requièrent des informations sur la reproduction, la survie et la morphologie des individus. Des méthodes de marquage-recapture permettent d'obtenir ces informations. Cependant, lorsque l'écologie de l'espèce rend les obervations directes impossibles, des méthodes indirectes doivent être utilisées pour obtenir le succès reproducteur des individus. Dans ce cas, les marqueurs moléculaires sont particulièrement utiles pour évaluer les relations génétiques entre individus et permettre la construction d'un pedigree. Cette thèse porte sur une population naturelle d'un petit mammifère insectivore, la musaraigne musette, Crocidura russula. Parce que cette espèce présente un mode de vie souterrain, les deux techniques complémentaires mentionnées ci-dessus ont été combinées pour acquérir les informations nécessaires. Les données ont été utilisées pour explorer divers aspects de biologie evolutive. Nous avons montré que la grande quantité de variance génétique trouvée chez cette espèce n'est pas maintenue par son système d'appariement. Celle-ci s'est en effet avérée être moins monogame que ce qui était admis jusqu'ici. Sa grande diversité génétique est plutôt entretenue par le flux de gènes provenant du voisinage. La dispersion a donc été un sujet phare dans cette thèse. Nous avons montré qu'elle n'est pas provoquée par un évitement de la consanguinité et nous n'avons pas trouvé de dépression de consanguité dans notre population. L'acquisition d'un territoire joue par contre un rôle important dans la dispersion. En outre, la dispersion possède une base génétique chez cette espèce. De plus, une fois qu'ils ont dispersé, les individus n'ont pas une valeur adaptative differente d'individus philopatriques. Le poids s'est avéré être un trait d'histoire de vie fortement influencé par la sélection sexuelle et de viabilité chez les deux sexes. Les gros individus ont accès à la reproduction parce qu'ils acquièrent et défendent un territoire plus facilement que les plus légers. Au contraire, les individus de taille intermédiaire sont favorisés par la sélection de viabilité, certainement à cause de contraintes écologiques et de coûts métaboliques. Finalement, nous avons montré que la majorité des traits d'histoire de vie dans notre population a le potentiel d'évoluer parce qu'elle maintient des quantités considérables de variance génétique additive. Néanmoins, l'héritabilité de ces traits d'histoire de vie n'est pas significative à cause de la grande quantité de variance non-additive ou environmentale associée à ces traits.

Constraints on the upper Permian to upper Triassic marine carbon isotope curve : case studies from the Tethys

Secular variations of the seawater carbon isotopic composition provide evidence for paleoceanographic and paleoclimatic changes and may serve for chemiostratigraphic correlations. The present study aimed to improve the current knowledge on the Upper Permian and Triassic segment of the Phanerozoic marine carbon isotope curve, whose Triassic part was poorly constrained by previous studies. Profiles of inorganic carbon isotopes are provided for sections from Himalaya (Salt Range, Kashmir, Spiti and Nepal), Oman and North Dobrogea (Romania) on the basis of whole-rock carbonate analysis. The data acquired, together with a literature compilation confirmed that most of the Upper Permian is characterized by high δ13C values (averaging +40/00) but failed to detect a positive excursion as suggested by recent compilations. In the light of these observations, the large drop in δ13C values associated with the end-Permian mass extinction appears to be driven by a sudden transfer of previously stocked 13C depleted carbon, rather than by the overturn of a Late Permian stratified ocean. The Triassic data-set outlines significant secular variations. The best documented is a carbon isotope positive excursion just across the Lower-Middle Triassic boundary, globally developed since it was detected in various paleogeographic settings. It is interpreted to reflect variations in surface ocean chemistry, possibly related to increased primary productivity, at times when the biotic recovery after the end-Permian mass-extinction began to accelerate significantly and when a sharp rise in seawater δ34S values occurred globally. Strontium isotope data obtained from well preserved biogenic phosphates allow a refinement of the Middle Triassic segment of the seawater strontium isotope curve and show a major inflexion point of the seawater strontium isotope curve also near the Lower Triassic - Middle Triassic boundary. These facts suggest that the transition from the Early to the Middle Triassic was a time of revolutionary global change which represented an important step in the evolution of Mesozoic marine environments. A tentative carbon isotope curve for the Upper Permian to Upper Triassic time interval is proposed. Its major features are: ? high but constant δ13C values during the Late Permian ? a sharp drop in δ13C values in the latest Permian ? subsequent recovery of δ13C values ? a short-lived positive excursion across the Early-Middle Triassic boundary ? a gradual rise in δ13C values starting in the Late Ladinian or in the Early Carnian It is foreseen that these fluctuations of the carbon isotope curve may serve as chronostratigraphic markers and further assist in the correlation of Permian and Triassic carbonate deposits.

Evidence for an episodic model of protein sequence evolution.

The evolution of protein function appears to involve alternating periods of conservative evolution and of relatively rapid change. Evidence for such episodic evolution, consistent with some theoretical expectations, comes from the application of increasingly sophisticated models of evolution to large sequence datasets. We present here some of the recent methods to detect functional shifts, using amino acid or codon models. Both provide evidence for punctual shifts in patterns of amino acid conservation, including the fixation of key changes by positive selection. Although a link to gene duplication, a presumed source of functional changes, has been difficult to establish, this episodic model appears to apply to a wide variety of proteins and organisms.

Evolution of the correlation between expression divergence and protein divergence in mammals.

Divergence of protein sequences and gene expression patterns are two fundamental mechanisms that generate organismal diversity. Here, we have used genome and transcriptome data from eight mammals and one bird to study the positive correlation of these two processes throughout mammalian evolution. We demonstrate that the correlation is stable over time and most pronounced in neural tissues, which indicates that it is the result of strong negative selection. The correlation is not driven by genes with specific functions and may instead best be viewed as an evolutionary default state, which can nevertheless be evaded by certain gene types. In particular, genes with developmental and neural functions are skewed toward changes in gene expression, consistent with selection against pleiotropic effects associated with changes in protein sequences. Surprisingly, we find that the correlation between expression divergence and protein divergence is not explained by between-gene variation in expression level, tissue specificity, protein connectivity, or other investigated gene characteristics, suggesting that it arises independently of these gene traits. The selective constraints on protein sequences and gene expression patterns also fluctuate in a coordinate manner across phylogenetic branches: We find that gene-specific changes in the rate of protein evolution in a specific mammalian lineage tend to be accompanied by similar changes in the rate of expression evolution. Taken together, our findings highlight many new aspects of the correlation between protein divergence and expression divergence, and attest to its role as a fundamental property of mammalian genome evolution.

Patterns of evolution of host proteins involved in retroviral pathogenesis.

BACKGROUND: Evolutionary analysis may serve as a useful approach to identify and characterize host defense and viral proteins involved in genetic conflicts. We analyzed patterns of coding sequence evolution of genes with known (TRIM5alpha and APOBEC3G) or suspected (TRIM19/PML) roles in virus restriction, or in viral pathogenesis (PPIA, encoding Cyclophilin A), in the same set of human and non-human primate species. RESULTS AND CONCLUSION: This analysis revealed previously unidentified clusters of positively selected sites in APOBEC3G and TRIM5alpha that may delineate new virus-interaction domains. In contrast, our evolutionary analyses suggest that PPIA is not under diversifying selection in primates, consistent with the interaction of Cyclophilin A being limited to the HIV-1M/SIVcpz lineage. The strong sequence conservation of the TRIM19/PML sequences among primates suggests that this gene does not play a role in antiretroviral defense.

Biogeochemical evolution of a marine shore tailings deposit during bioremediation

Summary : Mining activities produce enormous amounts of waste material known as tailings which are composed of fine to medium size particles. These tailings often contain sulfides, which oxidation can lead to acid and metal contamination of water; therefore they need to be remediated. In this work a tailings bioremediation approach was investigated by an interdisciplinary study including geochemistry, mineralogy and microbiology. The aim of the work was to study the effect of the implementation of wetland above oxidizing tailings on the hydrogeology and the biogeochemical element cycles, and to assess the system evolution over time. To reach these goals, biogeochemical processes occurring in a marine shore tailings deposit were investigated. The studied tailings deposit is located at the Bahìa de Ite, Pacific Ocean, southern Peru, where between 1940 and 1996 the tailings were discharged from the two porphyry copper mines Cuajone and Toquepala. After the end of deposition, a remediation approach was initiated in 1997 with a wetland implementation above the oxidizing tailings. Around 90% of the tailings deposits (total 16 km2) were thus remediated, except the central delta area and some areas close to the shoreline. The multi-stable isotope study showed that the tailings were saturated with fresh water in spite of the marine setting, due to the high hydraulic gradient resulting from the wetland implementation. Submarine groundwater discharge (SGD) was the major source of SO4 2-, C1-, Na+, Fe2+, and Mn2+ input into the tailings at the original shelf-seawater interface. The geochemical study (aquatic geochemistry and X-Ray diffraction (XRD) and sequential extractions from the solid fraction) showed that iron and sulfur oxidation were the main processes in the non-remediated tailings, which showed a top a low-pH oxidation zone with strong accumulation of efflorescent salts at the surface due to capillary upward transport of heavy metals (Fe, Cu, Zn, Mn, Cd, Co, and Ni) in the arid climate. The study showed also that the implementation of the wetland resulted in very low concentrations of heavy metals in solution (mainly under the detection limit) due to the near neutral pH and more reducing conditions (100-150 mV). The heavy metals, which were taken from solution, precipitated as hydroxides and sulfides or were bound to organic matter. The bacterial community composition analysis by Terminal Restriction Fragment Length Polymorphism (T-RFLP) and cloning and sequencing of 16S rRNA genes combined with a detailed statistical analysis revealed a high correlation between the bacterial distribution and the geochemical variables. Acidophilic autotrophic oxidizing bacteria were dominating the oxidizing tailings, whereas neutrophilic and heterotrophic reducing bacteria were driving the biogeochemical processes in the remediated tailings below the wetland. At the subsurface of the remediated tailings, an iron cycling was highlighted with oxidation and reduction processes due to micro-aerophilic niches provided by the plant rhizosphere in this overall reducing environment. The in situ bioremediation experiment showed that the main parameter to take into account for the effectiveness was the water table and chemistry which controls the system. The constructed remediation cells were more efficient and rapid in metal removal when saturation conditions were available. This study showed that the bioremediation by wetland implementation could be an effective and rapid treatment for some sulfidic mine tailings deposits. However, the water saturation of the tailings has to be managed on a long-term basis in order to guarantee stability. Résumé : L'activité minière produit d'énormes quantités de déchets géologiques connus sous le nom de « tailings » composées de particules de taille fine à moyenne. Ces déchets contiennent souvent des sulfures dont l'oxydation conduit à la formation d'effluents acides contaminés en métaux, d'où la nécessité d'effectuer une remédiation des sites de stockage concernés. Le but de ce travail est dans un premier temps d'étudier l'effet de la bio-remédiation d'un dépôt de tailings oxydés sur l'hydrogéologie du système et les cycles biogéochimiques des éléments et en second lieu, d'évaluer l'évolution du processus de remédiation dans le temps. Le site étudié dans ce travail est situé dans la Bahía de Ite, au sud du Pérou, au bord de l'Océan Pacifique. Les déchets miniers en question sont déposés dans un environnement marin. De 1940 à 1996, les déchets de deux mines de porphyre cuprifère - Cuajone et Toquepala - ont été acheminés sur le site via la rivière Locumba. En 1997, une première remédiation a été initiée avec la construction d'une zone humide sur les tailings. Depuis, environ 90% de la surface du dépôt (16 km2) a été traité, les parties restantes étant la zone centrale du delta du Locumba et certaines zones proches de la plage. Malgré la proximité de l'océan, les études isotopiques menées dans le cadre de ce travail ont montré que les tailings étaient saturés en eau douce. Cette saturation est due à la pression hydraulique résultant de la mise en place des zones humides. Un écoulement d'eau souterrain sous-marin a été à détecté à l'interface entre les résidus et l'ancien fond marin. En raison de la géologie locale, il constitue une source d'entrée de SO4 2-, Cl-, Na+, FeZ+, et Mn2+ dans le système. L'analyse de la géochimie aquatique, la Diffraction aux Rayons X (XRD) et l'extraction séquentielle ont montré que l'oxydation du fer et .des sulfures est le principal processus se produisant dans les déchets non remédiés. Ceci a entraîné le développement d'une zone d'oxydation à pH bas induisant une forte accumulation des sels efflorescents, conséquence de la migration capillaire des métaux lourds (Fe, Cu, Zn, Mn, Cd, Co et Ni) de la solution vers la surface dans ce climat aride. Cette étude a montré également que la construction de la zone humide a eu comme résultats une précipitation des métaux dans des phases minérales en raison du pH neutre et des conditions réductrices (100-150mV). Les métaux lourds ont précipité sous la forme d'hydroxydes et de sulfures ou sont adsorbés à la matière organique. L'analyse de la composition de la communauté bactérienne à l'aide la technique T-RFLP (Terminal Restriction Fragment Length Polymorphism) et par le clonage/séquençage des gènes de l'ARNr 16S a été combinée à une statistique détaillée. Cette dernière a révélé une forte corrélation entre la distribution de bactéries spécifiques et la géochimie : Les bactéries autotrophes acidophiles dominent dans les déchets oxydés non remédiés, tandis que des bactéries hétérotrophes neutrophiles ont mené les processus microbiens dans les déchets remédiés sous la zone humide. Sous la surface de la zone humide, nos analyses ont également mis en évidence un cycle du fer par des processus d'oxydoréduction rendus possibles par la présence de niches micro-aérées par la rhizosphère dans cet environnement réducteur. L'expérience de bio-remédiation in situ a montré que les paramètres clés qui contrôlent l'efficacité du traitement sont le niveau de la nappe aquifère et la chimie de l'eau. Les cellules de remédiation se sont montrées plus efficaces et plus rapides lorsque le système a pu être saturé en eau. Finalement, cette étude a montré que la bio-remédiation de déchets miniers par la construction de zones humides est un moyen de traitement efficace, rapide et peu coûteux. Cependant, la saturation en eau du système doit être gérée sur le long terme afin de garantir la stabilité de l'ensemble du système.

The early Miocene rotation of the Corso-Sardinian block. New paleomagnetic constraints for the end of the motion

The paleomagnetic investigations carried out in the 70's on Oligo-Miocene volcanics of Sardinia have demonstrated that the island was turned by 35-30 degrees clockwise from 33 Ma up to 3-1-20.5 Ma and rotated counterclockwise in a few million years [De Jong et al., 1969, 1973; Bobier et Coulon, 1970; Coulon et al., 1974; Manzoni, 1974, 1975; Bellon rr nl.. 1977: Edel et Lortscher, 1977; Edel, 1979, 1980]. Since then, the end of the rotation fixed at 19 Ma by Montigny er al. [1981] was the subject of discussions and several studies associating paleomagnetism and radiometric dating were undertaken [Assorgia er al., 1994: Vigliotti et Langenheim, 1995: Deino et al., 1997; Gattacceca rt Deino, 1999]. This is a contribution to this debate that is hampered by thr important secular variation recorded in the volcanics. The only way to get our of this problem is to sample series of successive flows as completely as possible, and to reduce the effect of secular variation by the calculation of means. Sampling was performed north of Bonorva in 5 pyroclastic flows that belong to the upper ignimbritic series SI2 according to Coulon rr nl. [1974] or LBLS according to Assorgia et al. [1997] (fig. I). Ar-40/Ar-39 dating of biotites from the debris flow (MDF) has yielded an age or 18.35 +/- 0.03 Ma [Dubois, 2000]. Five of the investigated sites are located beneath the debris flow ITV, TVB, TVD, SPM85, SPM86), one site was cured in the matrix of the debris flow (MDF) and one in 4 metric blocks included in the flow (DFC). Another site was sampled in the upper ash flow (PDM) that marks the end of the pyroclastic activity, just before the marine transgression. According to micropaleontological and radiometric dating this transgression has occurred between 18.35 and 17.6 Ma [Dubois, 2000]. After removal of a soft viscous component, the thermal demagnetization generally shows a univectorial behaviour of the remanent magnetization (fig. 2a). The maximum unblocking temperatures of 580-620 degrees (tab. I) and a rapid saturation below 100 mT (fig. 3) indicate that the carrier of the characteristic magnetization is magnetite. The exception comes: from the upper site PDM in which were found two characteristic components, one with a normal polarity and low unblocking temperatures up to 350 degreesC and one with a reversed polarity and maximum unblocking temperatures at 580-600 degreesC of magnetite. After calculation of a mean direction for each flow, the mean << Al >> direction 4 degrees /57 degrees (alpha (95) = 13 degrees) computed with the mean directions for the 5 flows may be considered as weakly affected by secular variation. But the results require a more careful examination. The declinations are N to NNW beneath the debris flow. NNW in the debris flow. and NNE (or SSW) above the debris flow, The elongated distribution of the directions obtained at sites TVB and TVD. scattered from the mean direction of TV to the mean direction of MDF is interpreted as due to partial overprinting during the debris How volcanic episode, The low temperature component PDMa is likely related to the alteration seen on thin sections and is also viewed as an overprint. As NNE/SSW directions occur as well below (mean direction << B >> : 5 degrees /58 degrees) as above the debris flow (PDMb : 200 degrees/-58 degrees). the NNW directions (<< C >> : 337 degrees /64 degrees) associated with the debris flow volcanism may be interpreted as resulting from a magnetic field excursion. According to the polarity scale of Cande and Kent [1992, 1995] and the radiometric age of MDF, the directions with normal polarity (TV, TVB, TVD, SPM85. SPM86a. MDF. DFC) may represent the period 5En. while the directions with reversed polarity PDMb and SPM86b were likely acquired during the period 5Dr. Using the mean << Al >> direction, the mean << B >>, or the PDM direction (tab. I). the deviation in declination with the direction of stable Europe 6.4 degrees /58.7 degrees (alpha (95) = 8 degrees) for a selection of 4 middle Tertiary poles by Besse et Courtillot [1991] or 7 degrees /56 degrees (alpha (95) = 3 degrees) for 19 poles listed by Edel [1980] can be considered as negligible. Using the results from the uppermost ignimbritic layer of Anglona also emplaced around 18.3 Ma [Odin rt al.. 1994]. the mean direction << E >> (3 degrees /51.5 degrees) leads to the same conclusion. On the contrary, when taking into account all dated results available for the period 5En (mean direction << D >> 353 degrees /56 degrees for 45 sites) (tab. II). the deviation 13 degrees is much more significant. As the rotation of Sardinia started around 21-20.5 Ma. the assumption of a constant velocity of rotation and the deviations of the Sardinia directions with respect to the stable Europe direction locate the end of the motion between 18.3 and 17.2 or 16.7 Ma (fig. 4). During the interval 18.35-17.5 Ma, the marine transgression took place. At the same period a NE-SW shortening interpreted as resulting from the collision of Sardinia with Apulia affected different parts of the island [Letouzey et al., 1982]. Consequently, the new paleomagnetic results and the tectono-sedimentary evolution are in favour of an end of the rotation at 17.5-18 Ma.

Petrology and chemistry of the Chalten Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the Miocene

Résumé Scientific:Pétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe sujet de cette thèse est le Complexe Plutonique de Chaltén (CHPC), situé à la frontière entre le Chili et l'Argentine, en Patagonie (49°15'S). Ce complexe s'est mis en place au début du Miocène, dans un contexte de changements tectoniques importants. La géométrie et la vitesse de migration des plaques en Patagonie a été modifiée suite l'ouverture de la plaque Farallon il y a 25Ma (Pardo-Casas and Molnar 1987) et la subduction de la ride active du Chili sous la plaque sud-américaine il y a 14Ma (Cande and Leslie 1986). Les effets de cette reconfiguration tectonique sur la morphologie et le magmatisme de la plaque supérieure sont encore sujets à discussion. Dans ce contexte, un groupe d'intrusions miocènes - telle que le CHPC - est particulièrement intriguant, car en position transitionnelle entre le batholithe patagonien et l'arc volcanique cénozoïque et récent à l'ouest, et les laves de plateau de Patagonie à l'est (Fig. 1). A cause de leur position tectonique transitoire, ces plutons isolés hors du batholithe représentent un endroit clé pour comprendre les interactions entre la tectonique à large échelle et le magmatisme en Patagonie. Ici, je présente de nouvelles données de terrain, petrologiques, géochimiques et géochronologiques dans le but de caractériser la nature du CHPC, qui était largement inconnu avant cette étude, dans le but de tester l'hypothèse de migration de l'arc et erosion par subduction.Les résultats de l'investigation géochimique (chapitre 2) montrent que le CHPC n'est qu'un exemple parmi les plutons isolés d'arrière arc ave une composition calco-alcaline caractéristique, c-à-d une signature d'arc. La plupart de ces plutons isolés ont une composition alcaline. Le CHPC, contrairement, a une signature calco-alcaline avec Κ intermédiaire, tel que le batholithe patagonien et la plupart des roches volcaniques quaternaires liées à l'arc le long des Andes.De nouvelles données géochronologiques U-Pb de haute précision sur des zircons, acquis par TIMS, sur le CHPC donnent des âges entre 17.0 et 16.4Ma. Les âges absolus sont en accord avec la séquence intrusive déduite des relations de terrain (chapitre 1). Ces données sont les premières contraintes d'âge U-Pb sur le CHPC. Elles montrent clairement que l'histoire magmatique du CHPC n'a pas de lien direct avec la subduction de la ride à cette latitude (Cande and Leslie 1986), car le complexe est au moins 6Ma plus ancien.Une comparaison en profondeur avec les autres intrusions d'âge Miocène en Patagonie révèlent - pour la première fois - une évolution temporelle intéressante. Il y a une tendance E-W distincte au magmatisme calco-alcalin entre 20-16Ma avec une diminution de l'âge vers l'est - le CHPC est l'expression la plus orientale de cette tendance. Je suggère que la relation espace-temps reflète une migration vers l'est (vers le continent) de l'arc magmatique. Je propose que le facteur principal contrôlant cette migration est la subduction rapide suite à la reconfiguration de la vitesse des plaques tectoniques après l'ouverture la plaque Farallon (à ~26Ma) qui résulterait en une déformation importante ainsi qu'à des taux élevés d'érosion dans la fosse de subduction.Les rapports d'isotopes radiogéniques (Pb, Sr, Nd) élevés, une signature 6018 basse et un rapport Th/La élevé sont des paramètres distinctifs pour les roches mafiques du CHPC. Le modèle isotopique présenté (chapitre 2) suggère que cette signature reflète une contamination de la source, dans le coin de manteau, plutôt qu'une contamination crustale. La signature des éléments en trace du CHPC indiquent que le coin de manteau a été contaminé par des composés terrigènes, le plus vraisemblablement par des sédiments paléozoïques.Les travaux de terrain, la pétrographie et la géothermobarométrie ont été utilisés dans le but de comprendre l'histoire interne du CHPC (chapitre 3). Ces données suggèrent deux niveaux distincts de cristallisation : l'un dans la croûte moyenne (6 à 4.5kbar) et l'autre à un niveau peu profond (3.5 à 2kbar). La modélisation isotopique AFC de la contamination crustale indique des taux variables d'assimilation, qui ne sont pas corrélés avec le degré de différenciation. Cela suggère que différents volumes de magma se sont différenciés en profondeur, de façon indépendante. Cela implique que le CHPC se serait formés en plusieurs puises de magmas provenant d'au moins trois sources différentes. Les textures des granodiorites et des granites indiquent des teneurs élevées en cristaux avant la mise en place et, par conséquent, des températures d'emplacement faibles. Les observations de terrain montrent que les roches mafiques sont déformées, alors que ce n'est pas le cas pour les granodiorites et granites (plus jeunes). La déformation des roches mafiques est encore sujet de recherche, afin de savoir si elle est liée à la déformation régionale en régime compressif ou à l'emplacement lui-même. Cependant, la mise en place de grand volume de magma felsique riche en cristaux suggère un régime d'extension.Scientific Abstract:Petrology and chemistry of the Chaltén Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the MioceneThe subject of this thesis is the Chaltén Plutonic Complex (CHPC) located at the frontier between Chile and Argentina in Patagonia (at 49° 15 'Southern latitude). This complex intruded during early Miocene in a context of major tectonics changes. The plate geometry of Patagonia has been modified by changes in the plate motions after the break up of the Farallôn plate at 25Ma (Pardo-Casas and Molnar 1987) and by the subduction of the Chile spreading Ridge beneath South-America at 14 Ma (Cande and Leslie 1986). The effects of this tectonic setting on the morphology and the magmatism of the overriding plate are a matter of on-going discussion. Particularly intriguing in this context is a group of isolated Miocene intrusions - like the CHPC - which are located in a transitional position between the Patagonian Batholith and the Cenozoic and Recent volcanic arc in the West, and the Patagonian plateau lavas in the East (Fig. 1). Due to their transient tectonic position these isolated plutons outside the batholith represent a key to understanding the interaction between global-scale tectonics and magmatism in Patagonia. Here, I present new field, penological, geochemical and geochronological data to characterize the nature of the CHPC, which was largely unknown before this study, in order to test the hypothesis of time- transgressive magmatism.The results of the geochemical investigation (Chapter 2) show that the CHPC is only one among these isolated back-arc plutons with a characteristic calc-alkaline composition, i.e. arc signature. Most of these isolated intrusives have an alkaline character. The CHPC, in contrast, has a medium Κ calc-alkaline signature, like the Patagonian batholith and most of the Quaternary arc-related volcanic rocks along the Andes.New high precision TIMS U-Pb zircon dating of the CHPC yield ages between 17.0 to 16.4 Ma. The absolute ages support the sequence of intrusion relations established in the field (Chapter 1). These data are the first U-Pb age constraints on the CHPC, and clearly show that the magmatic history of CHPC has no direct link to the subduction of the ridge, since this complex is at least 6 Ma older than the time of collision of the Chile ridge at this latitude (Cande and Leslie 1986).An in-depth comparison with other intrusion of Miocene age in Patagonia reveals - for the first time - an interesting temporal pattern. There is a distinct E-W trend of calc-alkaline magmatism between 20-16 Ma with the younging of ages in the East - the CHPC is the easternmost expression of this trend. I suggest that this time-space relation reflects an eastward (landward) migration of the magmatic arc. I propose that main factor controlling this migration is the fast rates of subduction after the major reconfigurations of the plate tectonic motions after the break up of the Farallôn Plate (at -26 ) resulting in strong deformation and high rates of subduction erosion.High radiogenic isotope ratios (Pb, Sr, Nd) ratios, low 5018 signature and high Th/La ratios in mafic rocks are distinctive features of the CHPC. The presented isotopic models (Chapter 2) suggest that this signature reflects source contamination of the mantle wedge rather than crustal contamination. The trace element signature of the CHPC indicates that the mantle wedge was contaminated with a terrigenous component, most likely from Paleozoic sediments.Fieldwork, petrography and geothermobarometry were used to further unravel the internal history of the CHPC (Chapter 3). These data suggest two main levels of crystallization: one a mid crustal levels (6 to 4.5 kbar) and other a shallow level (3.5 to 2 kbar). Isotopic AFC modeling of crustal contamination indicate variable rates of assimilation, which are not correlated with the degree of differentiation. This suggests that different batches of magma differentiate independently at depths. This implies that the CHPC would have formed by several pulses of magmas from at least 3 different sources. Textures of granodiorites and granites indicate a high content of crystals previous to the emplacement and consequently low emplacement temperatures. Field observations show that the mafic rocks are deformed, whereas the (younger) granodiorites and granites are not. It is still subject of investigation whether the deformation of the mafic rocks is related to regional deformation during a compressional regime or to the emplacement it self. However, the emplacement of huge amount of crystal rich felsic magmas suggests an extensional regime.Résumé Grand PublicPétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe Complexe Plutonique de Chaltén (CHPC) est un massif montagneux situé à 49°S à la frontière entre le Chili et l'Argentine, en Patagonie (région la plus au sud de l'Amérique du Sud). Il est composé de montagnes qui peuvent atteindre plus de 3000 mètres d'altitude, telles que le Cerro Fitz Roy (3400m) et le Cerro Torre (3100m). Ces montagnes sont composées de roches plutoniques, c.-à-d. des magmas qui se sont refroidis et ont cristallisés sous la surface terrestre.La composition chimique de ces roches montre que les magmas, qui ont formé ce complexe plutonique, font partie d'un volcanisme d'arc. Celui-ci se forme lorsqu'une plaque océanique plonge sous une plaque continentale. Les géologues appellent ce processus « subduction ». Dans un tel scénario, le manteau terrestre, qui se fait prendre entre ces deux plaques, fond pour former ainsi du magma. Ce magma remonte à travers la plaque continentale vers la surface. Si celui-ci atteint la surface, il forme les roches volcaniques, comme par exemple des laves. S'il n'atteint pas la surface, le magma se refroidit pour former finalement les roches plutoniques.Le long de la marge ouest d'Amérique du Sud, la plaque Nazca - qui se situe au sud-est de la plaque océanique pacifique - passe en dessous de la plaque d'Amérique du Sud. La bordure ouest du sud de la plaque sud-américaine a également été affectée par d'autres processus tectoniques, tels que des changements dramatiques dans les déplacements de plaques (il y a 25Ma) et la collision de la ride du Chili (depuis 15 Ma jusqu'à aujourd'hui). Ces caractéristiques tectoniques et magmatiques font de cette région un haut lieu pour les géologues. La plaque Nazca, s'est formée suite à l'ouverture d'une plaque océanique plus ancienne, il y a 25Ma. Cette ouverture est liée aux vitesses de subduction les plus rapides jamais connues. La ride du Chili est l'endroit où le sol de l'Océan Pacifique s'ouvre, formant deux plaques océaniques : les plaques Nazca et Antarctique. La ride du Chili subducte sous la plaque sud-américaine depuis 15Ma, en association avec la formation de grands volumes de magma ainsi que des changements morphologiques importants. La question de savoir lequel de ces changements tectoniques globaux affecte la géologie et la géographie de Patagonie a été, et est encore, discutée pendant de nombreuses années. De nombreux chercheurs suggèrent que la plupart des caractéristiques morphologiques et magmatiques en Patagonie sont liés à la subduction de la ride du Chili, mais cette suggestion est encore débattue comme le montre notre étude.Le batholithe de Patagonie du sud (SPB) est un énorme massif composé de roches plutoniques et il s'étend tout au long de la côte ouest de Patagonie (au sud de 47°S). Ces roches correspondent certainement aux racines d'un ancien arc volcanique, qui a été soulevé et érodé. Le CHPC, ainsi que d'autres petites intrusions dans la région, se situe dans une position exotique, à 100km à l'est du SPB. Certains chercheurs suggèrent que ces intrusions pourraient être liées à la subduction de la ride du Chili.Afin de débattre de cette problématique, nous avons utilisé différentes méthodes géochronologiques pour déterminer l'âge du CHPC et le comparer (a) à l'âge des roches intrusives similaires du SPB et (b) à l'âge de la collision de la ride du Chili. Dans ce travail, nous prouvons que le CHPC s'est formé au moins 7Ma avant la collision avec la ride du Chili. Sur la base des âges du CHPC et de la composition chimique de ses roches et minéraux, nous proposons que le CHPC fait partie d'un arc volcanique ancien. La migration de l'arc volcanique plus profondément dans le continent résulte de la grande vitesse de subduction entre 25 et lOMa. Des caractéristiques évidentes pour un tel processus - telles qu'une déformation importante et une vitesse d'érosion élevée - peuvent être rencontrées tout au long de la bordure ouest de l'Amérique du sud.

Evolutionary history and its relevance in understanding and conserving southern African biodiversity

Abstract : Understanding how biodiversity is distributed is central to any conservation effort and has traditionally been based on niche modeling and the causal relationship between spatial distribution of organisms and their environment. More recently, the study of species' evolutionary history and relatedness has permeated the fields of ecology and conservation and, coupled with spatial predictions, provides useful insights to the origin of current biodiversity patterns, community structuring and potential vulnerability to extinction. This thesis explores several key ecological questions by combining the fields of niche modeling and phylogenetics and using important components of southern African biodiversity. The aims of this thesis are to provide comparisons of biodiversity measures, to assess how climate change will affect evolutionary history loss, to ask whether there is a clear link between evolutionary history and morphology and to investigate the potential role of relatedness in macro-climatic niche structuring. The first part of my thesis provides a fine scale comparison and spatial overlap quantification of species richness and phylogenetic diversity predictions for one of the most diverse plant families in the Cape Floristic Region (CFR), the Proteaceae. In several of the measures used, patterns do not match sufficiently to argue that species relatedness information is implicit in species richness patterns. The second part of my thesis predicts how climate change may affect threat and potential extinction of southern African animal and plant taxa. I compare present and future niche models to assess whether predicted species extinction will result in higher or lower V phylogenetic diversity survival than what would be experienced under random extinction processes. l find that predicted extinction will result in lower phylogenetic diversity survival but that this non-random pattern will be detected only after a substantial proportion of the taxa in each group has been lost. The third part of my thesis explores the relationship between phylogenetic and morphological distance in southern African bats to assess whether long evolutionary histories correspond to equally high levels of morphological variation, as predicted by a neutral model of character evolution. I find no such evidence; on the contrary weak negative trends are detected for this group, as well as in simulations of both neutral and convergent character evolution. Finally, I ask whether spatial and climatic niche occupancy in southern African bats is influenced by evolutionary history or not. I relate divergence time between species pairs to climatic niche and range overlap and find no evidence for clear phylogenetic structuring. I argue that this may be due to particularly high levels of micro-niche partitioning. Résumé : Comprendre la distribution de la biodiversité représente un enjeu majeur pour la conservation de la nature. Les analyses se basent le plus souvent sur la modélisation de la niche écologique à travers l'étude des relations causales entre la distribution spatiale des organismes et leur environnement. Depuis peu, l'étude de l'histoire évolutive des organismes est également utilisée dans les domaines de l'écologie et de la conservation. En combinaison avec la modélisation de la distribution spatiale des organismes, cette nouvelle approche fournit des informations pertinentes pour mieux comprendre l'origine des patterns de biodiversité actuels, de la structuration des communautés et des risques potentiels d'extinction. Cette thèse explore plusieurs grandes questions écologiques, en combinant les domaines de la modélisation de la niche et de la phylogénétique. Elle s'applique aux composants importants de la biodiversité de l'Afrique australe. Les objectifs de cette thèse ont été l) de comparer différentes mesures de la biodiversité, 2) d'évaluer l'impact des changements climatiques à venir sur la perte de diversité phylogénétique, 3) d'analyser le lien potentiel entre diversité phylogénétique et diversité morphologique et 4) d'étudier le rôle potentiel de la phylogénie sur la structuration des niches macro-climatiques des espèces. La première partie de cette thèse fournit une comparaison spatiale, et une quantification du chevauchement, entre des prévisions de richesse spécifique et des prédictions de la diversité phylogénétique pour l'une des familles de plantes les plus riches en espèces de la région floristique du Cap (CFR), les Proteaceae. Il résulte des analyses que plusieurs mesures de diversité phylogénétique montraient des distributions spatiales différentes de la richesse spécifique, habituellement utilisée pour édicter des mesures de conservation. La deuxième partie évalue les effets potentiels des changements climatiques attendus sur les taux d'extinction d'animaux et de plantes de l'Afrique australe. Pour cela, des modèles de distribution d'espèces actuels et futurs ont permis de déterminer si l'extinction des espèces se traduira par une plus grande ou une plus petite perte de diversité phylogénétique en comparaison à un processus d'extinction aléatoire. Les résultats ont effectivement montré que l'extinction des espèces liées aux changements climatiques pourrait entraîner une perte plus grande de diversité phylogénétique. Cependant, cette perte ne serait plus grande que celle liée à un processus d'extinction aléatoire qu'à partir d'une forte perte de taxons dans chaque groupe. La troisième partie de cette thèse explore la relation entre distances phylogénétiques et morphologiques d'espèces de chauves-souris de l'Afrique australe. ll s'agit plus précisément de déterminer si une longue histoire évolutive correspond également à des variations morphologiques plus grandes dans ce groupe. Cette relation est en fait prédite par un modèle neutre d'évolution de caractères. Aucune évidence de cette relation n'a émergé des analyses. Au contraire, des tendances négatives ont été détectées, ce qui représenterait la conséquence d'une évolution convergente entre clades et des niveaux élevés de cloisonnement pour chaque clade. Enfin, la dernière partie présente une étude sur la répartition de la niche climatique des chauves-souris de l'Afrique australe. Dans cette étude je rapporte temps de divergence évolutive (ou deux espèces ont divergé depuis un ancêtre commun) au niveau de chevauchement de leurs niches climatiques. Les résultats n'ont pas pu mettre en évidence de lien entre ces deux paramètres. Les résultats soutiennent plutôt l'idée que cela pourrait être I dû à des niveaux particulièrement élevés de répartition de la niche à échelle fine.

Little effect of seasonal constraints on population genetic structure in eusocial paper wasps.

Climate has long been suggested to affect population genetic structures of eusocial insect societies. For instance, Hamilton [Journal of Theoretical Biology7 (1964) 17] discusses whether temperate and tropical eusocial insects may show differences in population-level genetic structure and viscosity, and how this might relate to differences in the degree of synchrony in their life cycles or modes of nest founding. Despite the importance of Hamilton's 1964 papers, this specific idea has not been tested in actual populations of wasps, probably due to the paucity of studies on tropical species. Here, we compare colony and population genetic structures in two species of primitively eusocial paper wasps with contrasting ecologies: the tropical species Polistes canadensis and the temperate species P. dominulus. Our results provide important clarifications of Hamilton's discussion. Specifically, we show that the genetic structures of the temperate and tropical species were very similar, indicating that seasonality does not greatly affect population viscosity or inbreeding. For both species, the high genetic differentiation between nests suggests strong selection at the nest level to live with relatives, whereas low population viscosity and low genetic differentiation between nest aggregations might reflect balancing selection to disperse, avoiding competition with relatives. Overall, our study suggests no prevalence of seasonal constraints of the life cycle in affecting the population genetic structure of eusocial paper wasps. These conclusions are likely to apply also to other primitively eusocial insects, such as halictine bees. They also highlight how selection for a kin structure that promotes altruism can override potential effects of ecology in eusocial insects.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.

Bioclimatic constraints to Andean cat distribution: a modelling application for rare species

AimTo identify the bioclimatic niche of the endangered Andean cat (Leopardus jacobita), one of the rarest and least known felids in the world, by developing a species distribution model.LocationSouth America, High Andes and Patagonian steppe. Peru, Bolivia, Chile, Argentina.MethodsWe used 108 Andean cat records to build the models, and 27 to test them, applying the Maxent algorithm to sets of uncorrelated bioclimatic variables from global databases, including elevation. We based our biogeographical interpretations on the examination of the predicted geographic range, the modelled response curves and latitudinal variations in climatic variables associated with the locality data.ResultsSimple bioclimatic models for Andean cats were highly predictive with only 3-4 explanatory variables. The climatic niche of the species was defined by extreme diurnal variations in temperature, cold minimum and moderate maximum temperatures, and aridity, characteristic not only of the Andean highlands but also of the Patagonian steppe. Argentina had the highest representation of suitable climates, and Chile the lowest. The most favourable conditions were centrally located and spanned across international boundaries. Discontinuities in suitable climatic conditions coincided with three biogeographical barriers associated with climatic or topographic transitions.Main conclusionsSimple bioclimatic models can produce useful predictions of suitable climatic conditions for rare species, including major biogeographical constraints. In our study case, these constraints are also known to affect the distribution of other Andean species and the genetic structure of Andean cat populations. We recommend surveys of areas with suitable climates and no Andean cat records, including the corridor connecting two core populations. The inclusion of landscape variables at finer scales, crucially the distribution of Andean cat prey, would contribute to refine our predictions for conservation applications.