Comparative performance of supertree algorithms in large data sets using the soapberry family (Sapindaceae) as a case study.

For the last 2 decades, supertree reconstruction has been an active field of research and has seen the development of a large number of major algorithms. Because of the growing popularity of the supertree methods, it has become necessary to evaluate the performance of these algorithms to determine which are the best options (especially with regard to the supermatrix approach that is widely used). In this study, seven of the most commonly used supertree methods are investigated by using a large empirical data set (in terms of number of taxa and molecular markers) from the worldwide flowering plant family Sapindaceae. Supertree methods were evaluated using several criteria: similarity of the supertrees with the input trees, similarity between the supertrees and the total evidence tree, level of resolution of the supertree and computational time required by the algorithm. Additional analyses were also conducted on a reduced data set to test if the performance levels were affected by the heuristic searches rather than the algorithms themselves. Based on our results, two main groups of supertree methods were identified: on one hand, the matrix representation with parsimony (MRP), MinFlip, and MinCut methods performed well according to our criteria, whereas the average consensus, split fit, and most similar supertree methods showed a poorer performance or at least did not behave the same way as the total evidence tree. Results for the super distance matrix, that is, the most recent approach tested here, were promising with at least one derived method performing as well as MRP, MinFlip, and MinCut. The output of each method was only slightly improved when applied to the reduced data set, suggesting a correct behavior of the heuristic searches and a relatively low sensitivity of the algorithms to data set sizes and missing data. Results also showed that the MRP analyses could reach a high level of quality even when using a simple heuristic search strategy, with the exception of MRP with Purvis coding scheme and reversible parsimony. The future of supertrees lies in the implementation of a standardized heuristic search for all methods and the increase in computing power to handle large data sets. The latter would prove to be particularly useful for promising approaches such as the maximum quartet fit method that yet requires substantial computing power.

Going the distance: human population genetics in a clinal world.

Global human genetic variation is greatly influenced by geography, with genetic differentiation between populations increasing with geographic distance and within-population diversity decreasing with distance from Africa. In fact, these 'clines' can explain most of the variation in human populations. Despite this, population genetics inferences often rely on models that do not take geography into account, which could result in misleading conclusions when working at global geographic scales. Geographically explicit approaches have great potential for the study of human population genetics. Here, we discuss the most promising avenues of research in the context of human settlement history and the detection of genomic elements under natural selection. We also review recent technical advances and address the challenges of integrating geography and genetics.

Qualitative and quantitative analysis of thalamo-cortical axons in the somatosensory cortex of normal and barrelless mice

Résumé Les rongeurs utilisent leurs moustaches (vibrisses) pour explorer le milieu environnant. Chaque moustache est mue par un système des muscles. Les récepteurs situés à sa base transmettent les informations au système nerveux central. La transmission vers l'écorce se fait via trois neurones de relais qui se trouvent au niveau du ganglion trigéminé, du tronc cérébral et du thalamus. La représentation corticale d'une vibrisse est une concentration des axones thalamo-corticaux (ATC) autour desquelles s'organisent leurs cibles, les cellules de la couche IV. La structure peut être identifiée histologiquement en coupes tangentielles et porte le nom de « barrel » (« tonneau »). Cette correspondance vibrisse - barrel fait de ce système un model idéal pour étudier l'influence de l'activité périphérique sur l'établissement et le maintien des cartes somatotopiques. Notre laboratoire dispose d'une souche de souris qui a subi une mutation spontanée pour le gène codant l'adenylyl cyclase I (ACI). Cette enzyme membranaire catalyse la formation de l'AMPc et joue un rôle important dans le guidage axonal, la libération des neurotransmetteurs et l'intégration des signaux postsynaptiques. Nous avons démontré dans un premier temps que cette souris adulte ne développe pas de barrels. Cela est dû à un manque d'organisation des ATC et aussi des cellules de la couche IV. De plus, les résultats électrophysiologiques montrent que les informations venant des vibrisses adjacentes ne sont pas intégrées d'une manière normale. Dans ce travail de thèse, j'ai analysé la morphologie des ATC révélés individuellement avec de la biocytine. L'analyse quantitative des ATC a mis en évidence les points suivants: 1. Les axones de la souris normale (NOR) quittent le thalamus, traversent la capsule interne et la substance blanche sous-corticale et pénètrent dans le cortex somato-sensoriel primaire. A l'intérieur de l'écorce ils traversent au maximum 3 colonnes corticales adjacentes dont une contient le barrel cible. En passant à travers les couches VI et V, ces axones arborisent et convergent progressivement vers le barrel dans lequel ils forment une riche arborisation. Un petit nombre des branches « errantes », pleines de boutons synaptiques, pénètrent dans les barrels voisins. Deux axones NOR provenant de corps cellulaires très proches dans le thalamus peuvent avoir un cheminement très divergent lors de la traversée de la capsule interne et de la substance blanche sous-corticale mais, à leur entrée dans le cortex, ils sont distants d'au maximum 2 colonnes corticales de la colonne qui contient le barrel cible et ils convergent progressivement vers ce barrel. 2. Les axones de la souris mutante (BRL) ont le même trajet sous-cortical que les axones NOR, mais leur entrée dans le cortex somato-sensoriel primaire est aléatoire. A l'interface entre la substance blanche sous-corticale et le cortex, l'axone principal se divise rapidement en troncs axonaux qui traversent les couches VI et V d'une manière divergente pour arriver dans la couche IV. Cela contraste beaucoup avec la trajectoire des NOR qui convergent graduellement vers leur barrel cible. Le nombre de branches radiales que les axones BRL utilisent pour entrer dans le cortex et dans la couche IV est double par rapport aux axones NOR. Parmi ces branches, seules quelques-unes donnent des arborisations, les autres ne sont pas développées et leur morphologie est semblable à celle des branches formées par les axones de la souris normale lors du développement. Deux axones BRL issus de corps cellulaires proches dans le thalamus peuvent avoir une trajectoire très divergente jusqu'à leur entrée dans la couche IV, mais à ce niveau ils sont réorientés pour se retrouver et faire un nombre maximal de branches et boutons synaptiques dans la même région corticale. Dans un cas extrême, un des axones observés est entré dans le cortex à la limite entre l'aire somatosensorielle primaire et secondaire et a parcouru une distance de 2 mm pour retrouver son partenaire thalamique et donner avec celui-ci un nombre maximal de branches dans la même région de la couche IV. 3. Les mesures quantitatives ont montré que les arborisations corticales des axones NOR ont une longueur moyenne de 18mm et sont formées par 200 segments qui portent 1200 boutons synaptiques. Par rapport à la souris NOR, les axones BRL ont en moyenne la même longueur, le même nombre de segments et boutons synaptiques, mais donnent deux fois plus de branches radiales. La surface tangentielle occupée par les arborisations BRL dans la couche IV est 2 fois plus grande que celle des NOR. Cela signifie que les 1000 boutons synaptiques qui caractérisent les arborisations NOR et BRL dans la couche IV sont disséminés sur une surface tangentielle double chez les derniers, et donc que la densité des boutons par unité de surface corticale est en moyenne plus faible. En effet, l'augmentation de la surface corticale tangentielle des BRL est due aux surfaces de faible et moyenne densité synaptique (0 - 8 boutons / 400pn2) qui augmentent 2 fois tandis que les surfaces de haute densité synaptiques (8 - 64 boutons / 4001.tm2) sont les mêmes. Nous émettons l'hypothèse selon laquelle, durant le développement, les ATC de la souris BRL divergent et forment un nombre exubérant de branches. Grâce à cette divergence et aux branches supranuméraires, ils trouvent l'endroit de l'écorce où se trouvent leurs voisins thalamiques et arborisent abondamment dans cette région. Cependant, le déficit en AGI ne leurs permet pas par la suite, sous influence de l'activité périphérique, de retirer les branches qui se trouvent dans les endroits inappropriés de l'écorce, avec de possibles conséquences sur la discrimination tactile.

Patterns of pain-free response in 497 cases of classic trigeminal neuralgia treated with Gamma Knife surgery and followed up for least 1 year.

Object The goal of this study was to establish whether clear patterns of initial pain freedom could be identified when treating patients with classic trigeminal neuralgia (TN) by using Gamma Knife surgery (GKS). The authors compared hypesthesia and pain recurrence rates to see if statistically significant differences could be found. Methods Between July 1992 and November 2010, 737 patients presenting with TN underwent GKS and prospective evaluation at Timone University Hospital in Marseille, France. In this study the authors analyzed the cases of 497 of these patients, who participated in follow-up longer than 1 year, did not have megadolichobasilar artery- or multiple sclerosis-related TN, and underwent GKS only once; in other words, the focus was on cases of classic TN with a single radiosurgical treatment. Radiosurgery was performed with a Leksell Gamma Knife (model B, C, or Perfexion) using both MR and CT imaging targeting. A single 4-mm isocenter was positioned in the cisternal portion of the trigeminal nerve at a median distance of 7.8 mm (range 4.5-14 mm) anterior to the emergence of the nerve. A median maximum dose of 85 Gy (range 70-90 Gy) was delivered. Using empirical methods and assisted by a chart with clear cut-off periods of pain free distribution, the authors were able to divide patients who experienced freedom from pain into 3 separate groups: patients who became pain free within the first 48 hours post-GKS; those who became pain free between 48 hours and 30 days post-GKS; and those who became pain free more than 30 days after GKS. Results The median age in the 497 patients was 68.3 years (range 28.1-93.2 years). The median follow-up period was 43.75 months (range 12-174.41 months). Four hundred fifty-four patients (91.34%) were initially pain free within a median time of 10 days (range 1-459 days) after GKS. One hundred sixty-nine patients (37.2%) became pain free within the first 48 hours (Group PF(≤ 48 hours)), 194 patients (42.8%) between posttreatment Day 3 and Day 30 (Group PF((>48 hours, ≤ 30 days))), and 91 patients (20%) after 30 days post-GKS (Group PF(>30 days)). Differences in postoperative hypesthesia were found: in Group PF(≤ 48 hours) 18 patients (13.7%) developed postoperative hypesthesia, compared with 30 patients (19%) in Group PF((>48 hours, ≤ 30 days)) and 22 patients (30.6%) in Group PF(>30 days) (p = 0.014). One hundred fifty-seven patients (34.4%) who initially became free from pain experienced a recurrence of pain with a median delay of 24 months (range 0.62-150.06 months). There were no statistically significant differences between the patient groups with respect to pain recurrence: 66 patients (39%) in Group PF(≤ 48 hours) experienced pain recurrence, compared with 71 patients (36.6%) in Group PF((>48 hours, ≤ 30 days)) and 27 patients (29.7%) in Group PF(>30 days) (p = 0.515). Conclusions A substantial number of patients (169 cases, 37.2%) became pain free within the first 48 hours. The rate of hypesthesia was higher in patients who became pain free more than 30 days after GKS, with a statistically significant difference between patient groups (p = 0.014).

Étude des sédiments quaternaires, de la molasse et sa tectonique, dans le Grand Lac (Léman) à partir de données sismiques 2D et 3D

RESUME L'Institut de Géophysique de l'Université de Lausanne a développé au cours de ces dernières années un système d'acquisition de sismique réflexion multitrace à haute résolution 2D et 3D. L'objectif de cette thèse était de poursuivre ce développement tout améliorant les connaissances de la géologie sous le lac Léman, en étudiant en particulier la configuration des grands accidents sous-lacustres dans la Molasse (Tertiaire) qui forme l'essentiel du substratum des formations quaternaires. En configuration 2D, notre système permet d'acquérir des profils sismiques avec une distance inter-CDP de 1,25 m. La couverture varie entre 6 et 18 selon le nombre de traces et la distance inter-tir. Le canon à air (15/15 eu. in.), offre une résolution verticale de 1,25 ni et une pénétration maximale de 300 m sous le fond de l'eau. Nous avons acquis au total plus de 400 km de sections 2D dans le Grand Lac et le Haut Lac entre octobre 2000 et juillet 2004. Une campagne de sismique 3D a fourni des données au large d'Evian sur une surface de 442,5 m sur 1450 m, soit 0,64 km2. La navigation ainsi que le positionnement des hydrophones et de la source ont été réalisés avec des GPS différentiels. Nous avons utilisé un traitement sismique conventionnel, sans appliquer d'AGC et en utilisant une migration post-stack. L'interprétation du substratum antéquaternaire est basée sur l'identification des sismofaciès, sur leurs relations avec les unités géologiques adjacentes au lac, ainsi que sur quelques données de forages. Nous obtenons ainsi une carte des unités géologiques dans le Grand Lac. Nous précisons la position du chevauchement subalpin entre la ville de Lausanne, sur la rive nord, et le bassin de Sciez, sur la rive sud. Dans la Molasse de Plateau, nous avons identifié les décrochements de Pontarlier et de St. Cergue ainsi que plusieurs failles non reconnues jusqu'ici. Nous avons cartographié les accidents qui affectent la Molasse subalpine ainsi que le plan de chevauchement du flysch sur la Molasse près de la rive sud du lac. Une nouvelle carte tectonique de la région lémanique a ainsi pu être dressée. L'analyse du substratum ne montre pas de failles suggérant une origine tectonique de la cuvette lémanique. Par contre, nous suggérons que la forme du creusement glaciaire, donc de la forme du lac Léman, a été influencée par la présence de failles dans le substratum antéquaternaire. L'analyse des sédiments quaternaires nous a permis de tracer des cartes des différentes interfaces ou unités qui les composent. La carte du toit du substratum antéquaternaire montre la présence de chenaux d'origine glaciaire dont la profondeur maximale atteint la cote -200 ni. Leur pente est dirigée vers le nord-est, à l'inverse du sens d'écoulement actuel des eaux. Nous expliquons cette observation par l'existence de circulations sous-glaciaires d'eau artésienne. Les sédiments glaciaires dont l'épaisseur maximale atteint 150 ni au centre du lac ont enregistré les différentes récurrences glaciaires. Dans la zone d'Evian, nous mettons en évidence la présence de lentilles de sédiments glaciolacustres perchées sur le flanc de la cuvette lémanique. Nous avons corrélé ces unités avec des données de forage et concluons qu'il s'agit du complexe inférieur de la pile sédimentaire d'Evian. Celui-ci, âgé de plus de 30 000 ans, serait un dépôt de Kame associé à un lac périglaciaire. La sismique réflexion 3D permet de préciser l'orientation de l'alimentation en matériel détritique de l'unité. La finesse des images obtenues nous permet également d'établir quels types d'érosion ont affecté certaines unités. Les sédiments lacustres, dont l'épaisseur maximale imagée atteint plus de 225 m et sans doute 400 ni sous le delta du Rhône, indiquent plusieurs mécanismes de dépôts. A la base, une mégaturbidite, épaisse d'une trentaine de mètres en moyenne, s'étend entre l'embouchure de la Dranse et le delta du Rhône. Au-dessus, la décantation des particules en suspension d'origine biologique et détritique fournit l'essentiel des sédiments. Dans la partie orientale du lac, les apports détritiques du Rhône forment un delta qui prograde vers l'ouest en s'imbriquant avec les sédiments déposés par décantation. La structure superficielle du delta a brutalement évolué, probablement à la suite de l'évènement catastrophique du Tauredunum (563 A.D.). Sa trace probable se marque par la présence d'une surface érosive que nous avons cartographiée. Le delta a ensuite changé de géométrie, avec notamment un déplacement des chenaux sous-lacustres. Sur l'ensemble de nos sections sismiques, nous n'observons aucune faille dans les sédiments quaternaires qui attesterait d'une tectonique postglaciaire du substratum. ABSTRACT During the last few years the institute of Geophysics of the University of Lausanne cleveloped a 2D and 3D high-resolution multichannel seismic reflection acquisition system. The objective of the present work was to carry on this development white improving our knowledge of the geology under Lake Geneva, in particular by studying the configuration of the large accidents affecting the Tertiary Molasse that makes up the basement of most Quaternary deposits. In its 2D configuration, our system makes it possible to acquire seismic profiles with a CDP interval of 1.25 m. The fold varies from 6 to 18 depending on the number of traces and the shooting interval. Our air gun (15/15 cu. in.) provides a vertical resolution of 1.25 m and a maximum penetration depth of approximately 300 m under water bottom. We acquired more than 400 km of 2D sections in the Grand Lac and the Haut Lac between October 2000 and July 2004. A 3D seismic survey off the city of Evian provided data on a surface of 442.5 m x 1450 m (0.64 km2). Ship's navigation as well as hydrophone- and source positioning were carried out with differential GPS. The seismic data were processed following a conventional sequence without .applying AGC and using post-stack migration. The interpretation of the pre-Quaternary substratum is based on sismofacies, on their relationships with terrestrial geological units and on some borehole data. We thus obtained a map of the geological units in the Grand Lac. We defined the location of the subalpine thrust from Lausanne, on the north shore, to the Sciez Basin, on the south shore. Within the Molasse de Plateau, we identified the already know Pontarlier and St Cergue transforms Fault as well as faults. We mapped faults that affect subalpine Molasse as well as the thrust fault plane between alpine flysch and Molasse near the lake's south shore. A new tectonic map of the Lake Geneva region could thus be drawn up. The substratum does not show faults indicating a tectonic origin for the Lake Geneva Basin. However, we suggest that the orientation of glacial erosion, and thus the shape of Lake Geneva, vas influenced by the presence of faults in the pre-Quaternary basement. The analysis of Quaternary sediments enabled us to draw up maps of various discontinuities or internal units. The top pre-Quaternary basement map shows channels of glacial origin, the deepest of them reaching an altitude of 200 m a.s.l. The channel's slopes are directed to the North-East, in opposite direction of the present water flow. We explain this observation by the presence of artesian subglacial water circulation. Glacial sediments, the maximum thickness of which reaches 150 m in the central part of the lake, record several glacial recurrences. In the Evian area, we found lenses of glacio-lacustrine sediments set high up on the flank of the Lake Geneva Bassin. We correlated these units with on-land borehole data and concluded that they represent the lower complex of the Evian sedimentary pile. The lower complex is aider than 30 000 years, and it could be a Kame deposit associated with a periglacial lake. Our 3D seismic reflexion survey enables us to specify the supply direction of detrital material in this unit. With detailed seismic images we established how some units were affected by different erosion types. The lacustrine sediments we imaged in Lake Geneva are thicker than 225 m and 400 m or more Linder the Rhone Delta. They indicate several depositional mechanisms. Their base is a major turbidite, thirty meters thick on average, that spreads between the Dranse mouth and the Rhone delta. Above this unit, settling of suspended biological and detrital particles provides most of the sediments. In the eastern part of the lake, detrital contribution from the Rhone builds a delta that progrades to the west and imbricates with the settling sediments. The shallow structure of the Rhone delta abruptly evolved, probably after the catastrophic Tauredunum event (563 A.D.). It probably coincides with an erosive surface that we mapped. As a result, the delta geometry changed, in particular associated with a displacement of water bottom channels. In all our seismic sections, we do not observe fault in the Quaternary sediments that would attest postglacial tectonic activity in the basement.

Les Passions filmées: des codes en appropriation, un cinéma en canonisation

Constituant l'un des premiers « genres » de l'histoire du cinéma (dont Burch et Gaudreault ont montré le rôle fondateur dans la standardisation des procédures de montage institutionnalisées), les films mettant en scène la Vie et la Passion du Christ fixent leurs normes en s'appropriant des codes iconographiques préétablis. Dans cet article, Valentine Robert s'attache à déployer le « palimpseste » de ces Passions des premiers temps, à démêler les « séries culturelles » impliquées, à dégager les phénomènes de reprises d'une bande à l'autre, et à replacer certains de ces jeux référentiels dans leur visée de légitimation - ou doit-on dire « canonisation » ? - du médium cinématographique.

Etude des modèles de Whittle markoviens probabilisés

Etude des modèles de Whittle markoviens probabilisés Résumé Le modèle de Whittle markovien probabilisé est un modèle de champ spatial autorégressif simultané d'ordre 1 qui exprime simultanément chaque variable du champ comme une moyenne pondérée aléatoire des variables adjacentes du champ, amortie d'un coefficient multiplicatif ρ, et additionnée d'un terme d'erreur (qui est une variable gaussienne homoscédastique spatialement indépendante, non mesurable directement). Dans notre cas, la moyenne pondérée est une moyenne arithmétique qui est aléatoire du fait de deux conditions : (a) deux variables sont adjacentes (au sens d'un graphe) avec une probabilité 1 − p si la distance qui les sépare est inférieure à un certain seuil, (b) il n'y a pas d'adjacence pour des distances au-dessus de ce seuil. Ces conditions déterminent un modèle d'adjacence (ou modèle de connexité) du champ spatial. Un modèle de Whittle markovien probabilisé aux conditions où p = 0 donne un modèle de Whittle classique qui est plus familier en géographie, économétrie spatiale, écologie, sociologie, etc. et dont ρ est le coefficient d'autorégression. Notre modèle est donc une forme probabilisée au niveau de la connexité du champ de la forme des modèles de Whittle classiques, amenant une description innovante de l'autocorrélation spatiale. Nous commençons par décrire notre modèle spatial en montrant les effets de la complexité introduite par le modèle de connexité sur le pattern de variances et la corrélation spatiale du champ. Nous étudions ensuite la problématique de l'estimation du coefficent d'autorégression ρ pour lequel au préalable nous effectuons une analyse approfondie de son information au sens de Fisher et de Kullback-Leibler. Nous montrons qu'un estimateur non biaisé efficace de ρ possède une efficacité qui varie en fonction du paramètre p, généralement de manière non monotone, et de la structure du réseau d'adjacences. Dans le cas où la connexité du champ est non observée, nous montrons qu'une mauvaise spécification de l'estimateur de maximum de vraisemblance de ρ peut biaiser celui-ci en fonction de p. Nous proposons dans ce contexte d'autres voies pour estimer ρ. Pour finir, nous étudions la puissance des tests de significativité de ρ pour lesquels les statistiques de test sont des variantes classiques du I de Moran (test de Cliff-Ord) et du I de Moran maximal (en s'inspirant de la méthode de Kooijman). Nous observons la variation de puissance en fonction du paramètre p et du coefficient ρ, montrant par cette voie la dualité de l'autocorrélation spatiale entre intensité et connectivité dans le contexte des modèles autorégressifs

Effects of muscle fibre shortening on the characteristics of surface motor unit potentials.

Traditionally, studies dealing with muscle shortening have concentrated on assessing its impact on conduction velocity, and to this end, electrodes have been located between the end-plate and tendon regions. Possible morphologic changes in surface motor unit potentials (MUPs) as a result of muscle shortening have not, as yet, been evaluated or characterized. Using a convolutional MUP model, we investigated the effects of muscle shortening on the shape, amplitude, and duration characteristics of MUPs for different electrode positions relative to the fibre-tendon junction and for different depths of the MU in the muscle (MU-to-electrode distance). It was found that the effects of muscle shortening on MUP morphology depended not only on whether the electrodes were between the end-plate and the tendon junction or beyond the tendon junction, but also on the specific distance to this junction. When the electrodes lie between the end-plate and tendon junction, it was found that (1) the muscle shortening effect is not important for superficial MUs, (2) the sensitivity of MUP amplitude to muscle shortening increases with MU-to-electrode distance, and (3) the amplitude of the MUP negative phase is not affected by muscle shortening. This study provides a basis for the interpretation of the changes in MUP characteristics in experiments where both physiological and geometrical aspects of the muscle are varied.

Long-distance wolf recolonization of France and Switzerland inferred from non-invasive genetic sampling over a period of 10 years

In the early 1900s, the wolf (Canis lupus) was extirpated from France and Switzerland. There is growing evidence that the species is presently recolonizing these countries in the western Alps. By sequencing the mitochondrial DNA (mtDNA) control region of various samples mainly collected in the field (scats, hairs, regurgitates, blood or tissue; n = 292), we could (1) develop a non-invasive method enabling the unambiguous attribution of these samples to wolf, fox (Vulpes vulpes) or dog (Canis familiaris), among others; (2) demonstrate that Italian, French and Swiss wolves share the same mtDNA haplotype, a haplotype that has never been found in any other wolf population world-wide. Combined together, field and genetic data collected over 10 years corroborate the scenario of a natural expansion of wolves from the Italian source population. Furthermore, such a genetic approach is of conservation significance, since it has important consequences for management decisions. This first long-term report using non-invasive sampling demonstrates that long-distance dispersers are common, supporting the hypothesis that individuals may often attempt to colonize far from their native pack, even in the absence of suitable corridors across habitats characterized by intense human activities.