On the fate of sexual traits under asexuality.

Environmental shifts and life-history changes may result in formerly adaptive traits becoming non-functional or maladaptive. In the absence of pleiotropy and other constraints, such traits may decay as a consequence of neutral mutation accumulation or selective processes, highlighting the importance of natural selection for adaptations. A suite of traits are expected to lose their adaptive function in asexual organisms derived from sexual ancestors, and the many independent transitions to asexuality allow for comparative studies of parallel trait maintenance versus decay. In addition, because certain traits, notably male-specific traits, are usually not exposed to selection under asexuality, their decay would have to occur as a consequence of drift. Selective processes could drive the decay of traits associated with costs, which may be the case for the majority of sexual traits expressed in females. We review the fate of male and female sexual traits in 93 animal lineages characterized by asexual reproduction, covering a broad taxon range including molluscs, arachnids, diplopods, crustaceans and eleven different hexapod orders. Many asexual lineages are still able occasionally to produce males. These asexually produced males are often largely or even fully functional, revealing that major developmental pathways can remain quiescent and functional over extended time periods. By contrast, for asexual females, there is a parallel and rapid decay of sexual traits, especially of traits related to mate attraction and location, as expected given the considerable costs often associated with the expression of these traits. The level of decay of female sexual traits, in addition to asexual females being unable to fertilize their eggs, would severely impede reversals to sexual reproduction, even in recently derived asexual lineages. More generally, the parallel maintenance versus decay of different trait types across diverse asexual lineages suggests that neutral traits display little or no decay even after extended periods under relaxed selection, while extensive decay for selected traits occurs extremely quickly. These patterns also highlight that adaptations can fix rapidly in natural populations of asexual organisms, in spite of their mode of reproduction.

A stable and 40Ar/39Ar isotope study of a major thrust in the Helvetic nappes (Swiss Alps) - Evidence for fluid-flow and constraints on nappe kinematics

Stable isotope and Ar-40/Ar-39 measurements,were made on samples associated with a major tectonic discontinuity in the Helvetic Alps, the basal thrust of the Diablerets nappe (external zone of the Alpine Belt) in order to determine both the importance of fluids in this thrust zone and the timing of thrusting. A systematic decrease in the delta(18)O values (up to 6 parts per thousand) of calcite, quartz, and white mica exists within a 10- to 70-m-wide zone over a distance of 37 km along the thrust, and they become more pronounced toward the root of the nappe. A similar decrease in the delta(13)C values of calcite is observed only in the deepest sections (up to 3 parts per thousand). The delta D-SMOW (SMOW = standard mean ocean water) values of white mica are -54 parts per thousand +/- 8 parts per thousand (n = 22) and are independent of the distance from the thrust. These variations are interpreted to reflect syntectonic solution reprecipitation during fluid passage along the thrust. The calculated delta(18)O and delta D values (versus SMOW) for the fluid in equilibrium with the analyzed minerals is 12 parts per thousand to 16 parts per thousand and -30 parts per thousand to +5 parts per thousand, respectively, for assumed temperatures of 250 to 450 degrees C. The isotopic and structural data are consistent with fluids derived from the deep-seated roots of the Helvetic nappes where large volumes of Mesozoic sediments were metamorphosed to the amphibolite facies, It is suggested that connate and metamorphic waters, overpressured by rapid tectonic burial in a subductive system escaped by upward infiltration along moderately dipping pathways until they reached the main shear zone at the base of the moving pile, where they were channeled toward the surface, This model also explains the mechanism by which large amounts of fluid were removed from the Mesozoic sediments during Alpine metamorphism. White mica Ar-49/Ar-39 ages vary from 27 Ma far from the Diablerets thrust to 15 Ma along the thrust. An older component is observed in micas far from the thrust, interpreted as a detrital signature, and indicates that regional metamorphic temperatures were less than about 350 degrees C. The;plateau and near plateau ages nearest the thrust are consistent with either neocrystallization of white mica or argon loss by recrystallization during thrusting, which may have been enhanced in the zones of highest fluid flow. The 15 Ma Ar-40/Ar-39 age plateau measured on white mica sampled exactly on the thrust surface dates the end of both fluid flow and tectonic transport.

Liquid line of descent of a basanitic liquid at 1.5 GPa: Constraints for metasomatic vein formation

The metasomatism observed in the oceanic and continental lithosphere is generally interpreted to represent a continuous differentiation process forming anhydrous and hydrous veins plus a cryptic enrichment in the surrounding peridotite. In order to constrain the mechanisms of vein formation and potentially clarify the nature and origin of the initial metasomatic agent, we performed a series of high-pressure experiments simulating the liquid line of descent of a basanitic magma differentiating within continental or mature oceanic lithosphere. This series of experiments has been conducted in an end-loaded piston cylinder apparatus starting from an initial hydrous ne-normative basanite at 1.5 GPa and temperature varying between 1,250 and 980°C. Near-pure fractional crystallization process was achieved in a stepwise manner in 30°C temperature steps and starting compositions corresponding to the liquid composition of the previous, higher-temperature glass composition. Liquids evolve progressively from basanite to peralkaline, aluminum-rich compositions without significant SiO2 variation. The resulting cumulates are characterized by an anhydrous clinopyroxene + olivine assemblage at high temperature (1,250-1,160°C), while at lower temperature (1,130-980°C), hydrous cumulates with dominantly amphibole + minor clinopyroxene, spinel, ilmenite, titanomagnetite and apatite (1,130-980°C) are formed. This new data set supports the interpretation that anhydrous and hydrous metasomatic veins could be produced during continuous differentiation processes of primary, hydrous alkaline magmas at high pressure. However, the comparison between the cumulates generated by the fractional crystallization from an initial ne-normative liquid or from hy-normative initial compositions (hawaiite or picrobasalt) indicates that for all hydrous liquids, the different phases formed upon differentiation are mostly similar even though the proportions of hydrous versus anhydrous minerals could vary significantly. This suggests that the formation of amphibole-bearing metasomatic veins observed in the lithospheric mantle could be linked to the differentiation of initial liquids ranging from ne-normative to hy-normative in composition. The present study does not resolve the question whether the metasomatism observed in lithospheric mantle is a precursor or a consequence of alkaline magmatism; however, it confirms that the percolation and differentiation of a liquid produced by a low degree of partial melting of a source similar or slightly more enriched than depleted MORB mantle could generate hydrous metasomatic veins interpreted as a potential source for alkaline magmatism by various authors.

The MIGCLIM R package - seamless integration of dispersal constraints into projections of species distribution models

The MIGCLIM R package is a function library for the open source R software that enables the implementation of species-specific dispersal constraints into projections of species distribution models under environmental change and/or landscape fragmentation scenarios. The model is based on a cellular automaton and the basic modeling unit is a cell that is inhabited or not. Model parameters include dispersal distance and kernel, long distance dispersal, barriers to dispersal, propagule production potential and habitat invasibility. The MIGCLIM R package has been designed to be highly flexible in the parameter values it accepts, and to offer good compatibility with existing species distribution modeling software. Possible applications include the projection of future species distributions under environmental change conditions and modeling the spread of invasive species.

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.

Developmental changes in lateralized inhibition of symmetric movements in children with and without Developmental Coordination Disorder.

The present study investigates developmental changes in selective inhibition of symmetric movements with a lateralized switching task from bimanual to unimanual tapping in typically developing (TD) children and with Developmental Coordination Disorder (DCD) from 7 to 10 years old. Twelve right-handed TD children and twelve gender-matched children with DCD and probable DCD produce a motor switching task in which they have (1) to synchronize with the beat of an auditory metronome to produce bimanual symmetrical tapping and (2) to selectively inhibit their left finger's tapping while continuing their right finger's tapping and conversely. We assess (1) the development of the capacity to inhibit the stopping finger (number of supplementary taps after the stopping instruction) and (2) the development of the capacity to maintain the continuing finger (changes in the mean tempo and its variability for the continuing finger's tapping) and (3) the evolution of performance through trials. Results indicate that (1) TD children present an age-related increase in the capacity to inhibit and to maintain the left finger's tapping, (2) DCD exhibits persistent difficulties to inhibit the left finger's tapping, and (3) both groups improve their capacity to inhibit the left finger's movements through trials. In conclusion, the lateralized switching task provides a simple and fine tool to reveal differences in selective inhibition of symmetric movements in TD children and children with DCD. More theoretically, the specific improvement in selective inhibition of the left finger suggests a progressive development of inter-hemispheric communication during typical development that is absent or delayed in children with DCD.

Constraints on the permeability structure of alluvial aquifers from P-wave sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

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.

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.

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.

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.