Predicting present and future intra-specific genetic structure through niche hindcasting across 24 millennia.

Paleoclimatic reconstructions coupled with species distribution models and identification of extant spatial genetic structure have the potential to provide insights into the demographic events that shape the distribution of intra-specific genetic variation across time. Using the globeflower Trollius europaeus as a case-study, we combined (1) Amplified Fragment Length Polymorphisms, (2) suites of 1000-years stepwise hindcasted species distributions and (3) a model of diffusion through time over the last 24,000 years, to trace the spatial dynamics that most likely fits the species' current genetic structure. We show that the globeflower comprises four gene pools in Europe which, from the dry period preceding the Last Glacial Maximum, dispersed while tracking the conditions fitting its climatic niche. Among these four gene pools, two are predicted to experience drastic range retraction in the near future. Our interdisciplinary approach, applicable to virtually any taxon, is an advance in inferring how climate change impacts species' genetic structures.

Behaviour evaluation for risk-taking adolescents (BERTA): an easy to use and assess instrument to detect adolescent risky behaviours in a clinical setting.

To create an instrument to be used in an outpatient clinic to detect adolescents prone to risk-taking behaviours. Based on previous research, five identified variables (relationship with parents and teachers, liking going to school, average grades, and level of religiosity) were used to create a screening tool to detect at least one of ten risky behaviours (tobacco, alcohol, cannabis and other illegal drugs use; sexual intercourse and sexual risky behaviour; driving while intoxicated, riding with an intoxicated driver, not always using a seat belt, and not always using a helmet). The instrument was tested using the Barcelona Adolescent Health Survey 1993. A Receiver Operating Characteristics curve was used to find the best cut-off point between high and low risk score. Odds ratios and 95% confidence intervals were calculated to detect at least one risky behaviour and for each individual behaviour. In order to assess its predictive value, the analysis was repeated using the Barcelona Adolescent Health Survey 1999. In both cases, analyses were conducted for the whole sample and for younger and older adolescents. Adolescents with a high-risk score were more likely to take at least one risky behaviour both when the whole sample was analysed and by age groups. With very few exceptions, the Behaviour Evaluation for Risk-Taking Adolescents showed significant odds ratios for each individual variable. CONCLUSION: The Behaviour Evaluation for Risk-Taking Adolescents has shown its potential as an easy to use instrument to screen for risk-taking behaviours. Future research must aim towards assessing this instrument's predictive value in the clinical setting and it's application to other populations.

The Paleozoic magmatic history of the Maggia and Sambuco nappes, Lower Penninic, Central Lepontine Alps

Résumé : Les corps magmatiques sont des indicateurs essentiels dans toute reconstitution paléogéographique et/ou géodynamique d'un cycle orogénique, en particulier en contexte polycyclique, où la plupart des autres indices ont été oblitérés. Ils sont aisément datables et leurs caractéristiques géochimiques permettent de contraindre leur contexte tectonique de mise en place. Cette approche a été appliquée aux socles pré-mésozoïques des nappes penniques inférieures de Sambuco et de la Maggia, dans les Alpes centrales lepontines. Plusieurs événements magmatiques ont été identifiés dans le socle de Sambuco et datés par la méthode U-Pb sur zircon couplée à la technique LA-ICPMS. La suite calco-alcaline mafique rubanée de Scheggia est datée du Cambrien inférieur à 540-530 Ma ; le métagranite alumineux oeillé de Sasso Nero a un âge de 480-470 Ma, tout comme bien d'autres «older orthogneisses» des socles alpins. Il contient des zircons hérités d'âge panafricain à 630-610 Ma, indicateur d'une affiliation gondwanienne de ces terrains. Le pluton calco-alcalin du Matorello est daté à environ 300-310 Ma, et les filons lamprophyriques qu'il abrite à 300 Ma. La granodiorite de Cocco et le leucogranite de Ruscada, tous deux intrudés dans le socle de la nappe adjacente de la Maggia, ont des âges similaires à celui du Matorello. Ceci ajouté aux similitudes magmatiques observées entre Cocco et Matorello suggère une proximité paléogéographique des deux nappes au Permien-Carbonifère. Or ces dernières sont actuellement considérées appartenir à deux domaines paléogéographiques mésozoïques distincts : helvétique pour Sambuco et briançonnais pour Maggia, séparés par un bassin océanique. Si tel fut le cas, aucun mouvement décrochant ne doit avoir décalé les marges continentales de l'océan, retrouvées en parfaite coïncidence lors de sa fermeture. Le Matorello est un pluton recristallisé en faciès amphibolite et plissé par cinq phases successives de déformation non-coaxiales, qui ont conduit à son renversement complet, attesté par des indicateurs de paléogravité. Il préserve de spectaculaires phénomènes de coexistence liquide de magmas (essaims d'enclaves et Bills composites). Ce pluton était originellement tabulaire, construit par l'accumulation de multiples injections de magma en feuillets d'épaisseur métrique à décamétrique. Suivant le rythme de mise en place, les injections successives ont rapidement cristallisé avec des contours nets et bien définis (Bills composites) ou se sont mélangées avec les précédentes pour former une couche non consolidée de plusieurs dizaines de mètres d'épaisseur (granodiorite principale). Les injections individuelles sont délimitées par de subtils contrastes en granulométrie, proportions modales ou ségrégation de minéraux (schlieren), ou par des phénomènes d'érosion le long des surfaces de contact. Deux couches métriques à contour sinueux consistent en une accumulation compacte d'enclaves mafiques arrondies dans une matrice granodioritique fine. Le granoclassement des enclaves, la présence de figures de charge et de phénomènes érosifs en base de couche, ainsi que des schlieren de biotite entrecroisés évoquent l'injection de coulées de magma chargé d'enclaves et de faible viscosité en régime hydrodynamique turbulent dans un encaissant granodioritique encore largement liquide. La nature hybride des roches implique une chambre magmatique sous-jacente, en cours de différenciation et périodiquement réalimentée. Les magmas sont des liquides mafiques dérivés du manteau et des liquides anatectiques d'origine crustale, comme l'indique la gamme mesurée des rapports isotopiques initiaux du Sr (0.704 à 0.709) et des valeurs epsilon Nd (-2.1 à -4.7). Ces données montrent également que la contribution crustale est dominante, en accord avec les isotopes du plomb. Les phénomènes d'hybridation ont vraisemblablement eu lieu en base de croûte et dans la chambre magmatique sous-jacente au laccolite du Matorello. Les indicateurs de paléogravité du Matorello contribuent accessoirement à la compréhension de l'architecture actuelle de la nappe de Sambuco. Des plis isoclinaux à surface axiale verticale peuvent être mis en évidence par le contact entre les faciès dioritique et granodioritique. L'antiforme dont le Matorello forme le coeur est un synclinal, ce qui le positionne dans le Flanc inverse du grand pli couché que forme la nappe de Sambuco. Par ailleurs, des blocs de gneiss retrouvés dans le wildflysch sommital de la couverture de la nappe d'Antigorio ont été affiliés dans cette étude au pluton du Matorello. Ceci implique que le front de la nappe de Sambuco chevauchait déjà la partie est du bassin d'Antigorio au moment de sa fermeture. Par conséquent, ce n'est qu'en position externe que la nappe du Lebendun chevauche directement la nappe d'Antigorio. Abstract Magmatic bodies are important markers in paleo-geographic or geodynamic reconstructions of orogenic cycles, even more so in the case of polycyclic events where many of the other markers have been overwritten or destroyed. Plutons are relatively easy to date and their geochemical properties help constrain the tectonic context in which they were emplaced. This study focuses on the pre-mesozoic basement in the Sambuco and Maggia lower Penninic nappes located in the central Lepontine domain of the Alps. A number of magmatic events have been identified in the Sambuco basement. These events were dated using LA-ICPMS U/Pb on zircon grains. The mafic calc-alkaline banded Scheggia suite is dated as lower Cambrian, 540-530 Ma. The Al-rich Sasso-Nero lenticular gneiss is 480-470 Ma old (similarly to many older orfhogneisses of the Alpine basement) and contains 630-610 Ma old pan-African inherited zircons that illustrate the Gondwanian origin of these terranes.The calc-alkaline Matorello pluton is dated as 310-300 Ma whereas the lamprophyric bodies it contains are of 300 Ma. The Cocco granodiorite and the Ruscada leucogranite both intrude the basement of the adjacent Maggia nappe and are of similar ages to the Matorello. The ages as well as the geochemical similarities between the Cocco, Rucada and Matorello plutons suggest their paleo-geographic proximity at the Permian-Carboniferous boundary. However, these nappes are currently considered as belonging to two different Mesozoic paleo-geographic domains. Indeed, the Sambuco is considered as Helvetic whereas the Maggia is said to be Briançonnais, both separated by an oceanic basin. If this is the case, then it is essential that nostrike-slip movement has misaligned both continental margins since these coincide perfectly now that the oceanic domain closed. The Matorello pluton was originally a tabular intrusion, built up by the accumulation of multiple, several meter-thick, subhorizontal sheet-like injections of magma. Depending on their emplacement rate, the successive magma injections either solidified rapidly with sharp and rather well-defined boundaries (like the composite sills) or mingled with previous injections generating a thick molten layer up to several tens to hundred meters thick, like in the main granodioritic facies. These coalesced injections are hardly distinguishable, however subtle contrasts in granulometry, mineral modal proportions or mineral sorting (cross-bedded biotite-rich schlieren), as well as erosional features and/or crystal entrapment along contact surfaces allow to distinguish between the different injections. Two exceptional meter-thick layers display sinuous boundaries with the host granodiorite and consist of a densely packed accumulation of mafic enclaves in a granodioritic matrix. Gravitational sorting of the enclaves with load cast features at the base of the layers and sinuous biotite schlieren point to injection of low viscosity turbulent composite magma flows in the still largely molten granodiorite host. The hybrid nature of these rocks implies the existence of á periodically replenished and differentiated underlying magma chamber. Magmas are mafic liquids derived from the mantle and anatectic liquids of crustal origin, as shown by the (87Sr/86Sr), and epsilon Nd values (0.704-0.709 and -2.1 to -4.7 respectively. These data show that the crustal contribution is important, as confirmed by the Pb isotopes. The hybridisation processes seem to have occurred in the lower crust in magma chambers underlying the Matorello laccolith. The paleo-gravity markers in the Matorello help understand the architecture of the Sambuco nappe. Isoclinal folds with a vertical axial plane can be seen at the contact between dioritic and granodioritic facies. The antiform structure of which the Matorello is the heart is in fact a syncline. This places it in the inverse flanc of the large recumbent fold that constitutes the Sambuco nappe. The gneiss blocs found in the summital wildflysh cover of the Antigorio nappe have been linked to the Matorello pluton. This means that the front of the Sambuco nappe already overlapped the Antigorio basin when it closed. This implies that the Lebendun nappe can only overlap the Antigorio nappe in it's external position. Résumé grand public La chaîne alpine est la conséquence de la collision tertiaire entre deux masses continentales, l'Europe au nord et la péninsule apulienne africaine au sud, originellement séparées par l'océan mésozoïque téthysien. Cette collision a fermé un espace large de plusieurs centaines de km avec pour résultat l'écaillage de la croûte terrestre en unités tectoniques de dimensions variables, qui se sont empilées, imbriquées, éventuellement replissées en nappes de géométrie complexe. Cet amoncellement de 40 km d'épaisseur a vu sa température et sa pression lithostatique internes augmenter jusqu'à des valeurs de l'ordre de 680 °C et 6000 bars, induisant une recristallisation métamorphique des roches. L'un des objectifs de la géologie alpine est de reconstituer la géographie de la région aux temps mésozoïques de l'océan téthysien, en d'autres termes, de replacer chacune des unités tectoniques identifiées au sein de l'empilement alpin dans sa position originelle. Le défi est de taille et peut être comparé à celui de la reconstitution d'un vaste puzzle, dont certaines pièces seraient endommagées au niveau de leur contour ou leurs couleurs (métamorphisme), dissimulées par d'autres (enfouissement), voire tombées de la table de jeu (subduction, échappement latéral). Plusieurs approches ont été mises en oeuvre au cours du siècle écoulé. On citera en particulier la stratigraphie, la tectonique et le paléomagnétisme. Dans ce travail, nous avons essentiellement utilisé des techniques de datation isotopique absolue des roches (U/Pb sur zircon) qui, sur la base des connaissances acquises par l'ensemble des autres disciplines géologiques, nous ont permis de mieux contraindre ta paléogéographie mésozoïque du domaine «pennique inférieur » des Alpes centrales lépontines. Et au-delà? Nous savons tous que la disposition des continents à la surface de la Terre évolue constamment. Il est donc tentant d'essayer de remonter plus loin encore dans le temps et de reconstituer la physionomie de la marge sud européenne, tout au moins certains éléments de son histoire, au cours de l'ère paléozoïque. Les traces de ces événements très anciens sont naturellement ténues et dans ce contexte, les techniques de datation mentionnées ci-dessus deviennent les outils les plus performants. Ainsi, des datations u/Pb sur zircon nous ont permis de recenser plusieurs intrusions magmatiques, attribuées à quatre événements orogéniques anté-alpins. Des âges néoprotérozoïques (630-610 millions d'années ou Ma), cambrien inférieur (540-530 Ma), ordovicien inférieur (480-470 Ma) et carbonifère supérieur-permien inférieur (310-285 Ma) ont été obtenus dans le socle de la nappe de Sambuco. Des âges similaires à 300 Ma ont été obtenus dans la nappe voisine de la Maggia, qui permettent de relier ces deux unités. Aujourd'hui côte à côte, ces deux nappes devaient également se trouver proches l'une de l'autre il y a 300 Ma, lors de l'extension post-varisque. Les structures magmatiques spectaculaires préservées dans le pluton du Matorello (300 Ma) contraignent la géométrie actuelle de la nappe de Sambuco dans laquelle l'intrusion s'est mise en place. La forme originelle du pluton, aujourd'hui retourné et replissé plusieurs fois, s'avère être tabulaire, faite d'intrusions de faible épaisseur (1-300 m) s'étalant en forme de disque (30m à 2 km de diamètre). Les injections successives de magma se sont accumulées sous un toit dioritique précoce; elles sont issues, par le refais de fractures, d'une chambre magmatique plus profonde, périodiquement réalimentée par des magmas calco-alcalins d'origine mantellique contaminés parla croûte continentale profonde (εNd = -2.1 à -4.7). Des accumulations d'enclaves magmatiques arrondies et granoclassées dans des paléo-chenaux à fond érosif témoignent de conditions de mise en place hydrodynamiques à haute énergie. Ces enclaves sont emmenées de la chambre magmatique sous-jacente à la faveur d'épisodes de fracturation hydraulique liés à l'injection de magmas matelliques chauds dans des liquides différenciés riches en eau. Cette hypothèse est étayée par l'existence de filons composites. Une paléohorizontale a pu être déduite au sein du pluton, indiquant que cette partie de la nappe de Sambuco est verticalisée et isoclinalement replissée par la déformation alpine. Finalement, des blocs érodés du socle Sambuco ont été retrouvés dans le wildflysch sommital de la couverture sédimentaire mésozoïque de la nappe d'Antigorio sous-jacente. Ceci suggère que les blocs ont été fournis parle front de la nappe de Sambuco en train de chevaucher sur la nappe d'Antigorio au moment de la fermeture du bassin sédimentaire de cette dernière.

A study on the expressive power of some fragments of the modal µ-Calculus

Résumé Le μ-calcul est une extension de la logique modale par des opérateurs de point fixe. Dans ce travail nous étudions la complexité de certains fragments de cette logique selon deux points de vue, différents mais étroitement liés: l'un syntaxique (ou combinatoire) et l'autre topologique. Du point de vue syn¬taxique, les propriétés définissables dans ce formalisme sont classifiées selon la complexité combinatoire des formules de cette logique, c'est-à-dire selon le nombre d'alternances des opérateurs de point fixe. Comparer deux ensembles de modèles revient ainsi à comparer la complexité syntaxique des formules as¬sociées. Du point de vue topologique, les propriétés définissables dans cette logique sont comparées à l'aide de réductions continues ou selon leurs positions dans la hiérarchie de Borel ou dans celle projective. Dans la première partie de ce travail nous adoptons le point de vue syntax¬ique afin d'étudier le comportement du μ-calcul sur des classes restreintes de modèles. En particulier nous montrons que: (1) sur la classe des modèles symétriques et transitifs le μ-calcul est aussi expressif que la logique modale; (2) sur la classe des modèles transitifs, toute propriété définissable par une formule du μ-calcul est définissable par une formule sans alternance de points fixes, (3) sur la classe des modèles réflexifs, il y a pour tout η une propriété qui ne peut être définie que par une formule du μ-calcul ayant au moins η alternances de points fixes, (4) sur la classe des modèles bien fondés et transitifs le μ-calcul est aussi expressif que la logique modale. Le fait que le μ-calcul soit aussi expressif que la logique modale sur la classe des modèles bien fondés et transitifs est bien connu. Ce résultat est en ef¬fet la conséquence d'un théorème de point fixe prouvé indépendamment par De Jongh et Sambin au milieu des années 70. La preuve que nous donnons de l'effondrement de l'expressivité du μ-calcul sur cette classe de modèles est néanmoins indépendante de ce résultat. Par la suite, nous étendons le langage du μ-calcul en permettant aux opérateurs de point fixe de lier des occurrences négatives de variables libres. En montrant alors que ce formalisme est aussi ex¬pressif que le fragment modal, nous sommes en mesure de fournir une nouvelle preuve du théorème d'unicité des point fixes de Bernardi, De Jongh et Sambin et une preuve constructive du théorème d'existence de De Jongh et Sambin. RÉSUMÉ Pour ce qui concerne les modèles transitifs, du point de vue topologique cette fois, nous prouvons que la logique modale correspond au fragment borélien du μ-calcul sur cette classe des systèmes de transition. Autrement dit, nous vérifions que toute propriété définissable des modèles transitifs qui, du point de vue topologique, est une propriété borélienne, est nécessairement une propriété modale, et inversement. Cette caractérisation du fragment modal découle du fait que nous sommes en mesure de montrer que, modulo EF-bisimulation, un ensemble d'arbres est définissable dans la logique temporelle Ε F si et seulement il est borélien. Puisqu'il est possible de montrer que ces deux propriétés coïncident avec une caractérisation effective de la définissabilité dans la logique Ε F dans le cas des arbres à branchement fini donnée par Bojanczyk et Idziaszek [24], nous obtenons comme corollaire leur décidabilité. Dans une deuxième partie, nous étudions la complexité topologique d'un sous-fragment du fragment sans alternance de points fixes du μ-calcul. Nous montrons qu'un ensemble d'arbres est définissable par une formule de ce frag¬ment ayant au moins η alternances si et seulement si cette propriété se trouve au moins au n-ième niveau de la hiérarchie de Borel. Autrement dit, nous vérifions que pour ce fragment du μ-calcul, les points de vue topologique et combina- toire coïncident. De plus, nous décrivons une procédure effective capable de calculer pour toute propriété définissable dans ce langage sa position dans la hiérarchie de Borel, et donc le nombre d'alternances de points fixes nécessaires à la définir. Nous nous intéressons ensuite à la classification des ensembles d'arbres par réduction continue, et donnons une description effective de l'ordre de Wadge de la classe des ensembles d'arbres définissables dans le formalisme considéré. En particulier, la hiérarchie que nous obtenons a une hauteur (ωω)ω. Nous complétons ces résultats en décrivant un algorithme permettant de calculer la position dans cette hiérarchie de toute propriété définissable.

Event-free survival at 24 months is a robust end point for disease-related outcome in diffuse large B-cell lymphoma treated with immunochemotherapy.

PURPOSE: Studies of diffuse large B-cell lymphoma (DLBCL) are typically evaluated by using a time-to-event approach with relapse, re-treatment, and death commonly used as the events. We evaluated the timing and type of events in newly diagnosed DLBCL and compared patient outcome with reference population data. PATIENTS AND METHODS: Patients with newly diagnosed DLBCL treated with immunochemotherapy were prospectively enrolled onto the University of Iowa/Mayo Clinic Specialized Program of Research Excellence Molecular Epidemiology Resource (MER) and the North Central Cancer Treatment Group NCCTG-N0489 clinical trial from 2002 to 2009. Patient outcomes were evaluated at diagnosis and in the subsets of patients achieving event-free status at 12 months (EFS12) and 24 months (EFS24) from diagnosis. Overall survival was compared with age- and sex-matched population data. Results were replicated in an external validation cohort from the Groupe d'Etude des Lymphomes de l'Adulte (GELA) Lymphome Non Hodgkinien 2003 (LNH2003) program and a registry based in Lyon, France. RESULTS: In all, 767 patients with newly diagnosed DLBCL who had a median age of 63 years were enrolled onto the MER and NCCTG studies. At a median follow-up of 60 months (range, 8 to 116 months), 299 patients had an event and 210 patients had died. Patients achieving EFS24 had an overall survival equivalent to that of the age- and sex-matched general population (standardized mortality ratio [SMR], 1.18; P = .25). This result was confirmed in 820 patients from the GELA study and registry in Lyon (SMR, 1.09; P = .71). Simulation studies showed that EFS24 has comparable power to continuous EFS when evaluating clinical trials in DLBCL. CONCLUSION: Patients with DLBCL who achieve EFS24 have a subsequent overall survival equivalent to that of the age- and sex-matched general population. EFS24 will be useful in patient counseling and should be considered as an end point for future studies of newly diagnosed DLBCL.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.

Phylogenetic approach reveals that virus genotype largely determines HIV set-point viral load.

HIV virulence, i.e. the time of progression to AIDS, varies greatly among patients. As for other rapidly evolving pathogens of humans, it is difficult to know if this variance is controlled by the genotype of the host or that of the virus because the transmission chain is usually unknown. We apply the phylogenetic comparative approach (PCA) to estimate the heritability of a trait from one infection to the next, which indicates the control of the virus genotype over this trait. The idea is to use viral RNA sequences obtained from patients infected by HIV-1 subtype B to build a phylogeny, which approximately reflects the transmission chain. Heritability is measured statistically as the propensity for patients close in the phylogeny to exhibit similar infection trait values. The approach reveals that up to half of the variance in set-point viral load, a trait associated with virulence, can be heritable. Our estimate is significant and robust to noise in the phylogeny. We also check for the consistency of our approach by showing that a trait related to drug resistance is almost entirely heritable. Finally, we show the importance of taking into account the transmission chain when estimating correlations between infection traits. The fact that HIV virulence is, at least partially, heritable from one infection to the next has clinical and epidemiological implications. The difference between earlier studies and ours comes from the quality of our dataset and from the power of the PCA, which can be applied to large datasets and accounts for within-host evolution. The PCA opens new perspectives for approaches linking clinical data and evolutionary biology because it can be extended to study other traits or other infectious diseases.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

What Was I Thinking? Eye-Tracking Experiments Underscore the Bias that Architecture Exerts on Nuclear Grading in Prostate Cancer.

We previously reported that nuclear grade assignment of prostate carcinomas is subject to a cognitive bias induced by the tumor architecture. Here, we asked whether this bias is mediated by the non-conscious selection of nuclei that "match the expectation" induced by the inadvertent glance at the tumor architecture. 20 pathologists were asked to grade nuclei in high power fields of 20 prostate carcinomas displayed on a computer screen. Unknown to the pathologists, each carcinoma was shown twice, once before a background of a low grade, tubule-rich carcinoma and once before the background of a high grade, solid carcinoma. Eye tracking allowed to identify which nuclei the pathologists fixated during the 8 second projection period. For all 20 pathologists, nuclear grade assignment was significantly biased by tumor architecture. Pathologists tended to fixate on bigger, darker, and more irregular nuclei when those were projected before kigh grade, solid carcinomas than before low grade, tubule-rich carcinomas (and vice versa). However, the morphometric differences of the selected nuclei accounted for only 11% of the architecture-induced bias, suggesting that it can only to a small part be explained by the unconscious fixation on nuclei that "match the expectation". In conclusion, selection of « matching nuclei » represents an unconscious effort to vindicate the gravitation of nuclear grades towards the tumor architecture.

Three-dimensional magnetic resonance myocardial motion tracking from a single image plane.

Three-dimensional imaging for the quantification of myocardial motion is a key step in the evaluation of cardiac disease. A tagged magnetic resonance imaging method that automatically tracks myocardial displacement in three dimensions is presented. Unlike other techniques, this method tracks both in-plane and through-plane motion from a single image plane without affecting the duration of image acquisition. A small z-encoding gradient is subsequently added to the refocusing lobe of the slice-selection gradient pulse in a slice following CSPAMM acquisition. An opposite polarity z-encoding gradient is added to the orthogonal tag direction. The additional z-gradients encode the instantaneous through plane position of the slice. The vertical and horizontal tags are used to resolve in-plane motion, while the added z-gradients is used to resolve through-plane motion. Postprocessing automatically decodes the acquired data and tracks the three-dimensional displacement of every material point within the image plane for each cine frame. Experiments include both a phantom and in vivo human validation. These studies demonstrate that the simultaneous extraction of both in-plane and through-plane displacements and pathlines from tagged images is achievable. This capability should open up new avenues for the automatic quantification of cardiac motion and strain for scientific and clinical purposes.

ZHARP: three-dimensional motion tracking from a single image plane.

Three-dimensional imaging and quantification of myocardial function are essential steps in the evaluation of cardiac disease. We propose a tagged magnetic resonance imaging methodology called zHARP that encodes and automatically tracks myocardial displacement in three dimensions. Unlike other motion encoding techniques, zHARP encodes both in-plane and through-plane motion in a single image plane without affecting the acquisition speed. Postprocessing unravels this encoding in order to directly track the 3-D displacement of every point within the image plane throughout an entire image sequence. Experimental results include a phantom validation experiment, which compares zHARP to phase contrast imaging, and an in vivo study of a normal human volunteer. Results demonstrate that the simultaneous extraction of in-plane and through-plane displacements from tagged images is feasible.

Monte Carlo simulation of a clearance box monitor used for nuclear power plant decommissioning.

When decommissioning a nuclear facility it is important to be able to estimate activity levels of potentially radioactive samples and compare with clearance values defined by regulatory authorities. This paper presents a method of calibrating a clearance box monitor based on practical experimental measurements and Monte Carlo simulations. Adjusting the simulation for experimental data obtained using a simple point source permits the computation of absolute calibration factors for more complex geometries with an accuracy of a bit more than 20%. The uncertainty of the calibration factor can be improved to about 10% when the simulation is used relatively, in direct comparison with a measurement performed in the same geometry but with another nuclide. The simulation can also be used to validate the experimental calibration procedure when the sample is supposed to be homogeneous but the calibration factor is derived from a plate phantom. For more realistic geometries, like a small gravel dumpster, Monte Carlo simulation shows that the calibration factor obtained with a larger homogeneous phantom is correct within about 20%, if sample density is taken as the influencing parameter. Finally, simulation can be used to estimate the effect of a contamination hotspot. The research supporting this paper shows that activity could be largely underestimated in the event of a centrally-located hotspot and overestimated for a peripherally-located hotspot if the sample is assumed to be homogeneously contaminated. This demonstrates the usefulness of being able to complement experimental methods with Monte Carlo simulations in order to estimate calibration factors that cannot be directly measured because of a lack of available material or specific geometries.

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence-absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size (n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

Tracking the ghost of coevolution in specific plant-insect interactions : the case of the pollination system between the globeflower Trollius europaeus and seven european Chiastocheta flies

This study aims at understanding the evolutionary processes at work in specialized species interactions. Prom the macroevolutionary perspective, coevolution among specialized taxa was proposed to be one of the major processes generating biodiversity. We challenge this idea from the theoretical and practical perspective and through a literature review and show that the major hypotheses linking coevolutionary process with macroevolutionary patterns do not necessarily predict lineage co diversification and parallel speciation, limit¬ing the utility of the comparative phylogenenetic approach for investigating coevolution¬ary processes. We also point to the rarity of observed long-term coevolutionary dynamics among lineages and propose that coevolution rather occurs in shorter timescales, followed by ecological fitting. Prom the empirical point, we focus on the nursery pollination interaction between the European globeflower Trollius europaeus (Ranunculaceae) and its associated Chiastocheta flies (Anthomyiidae; Diptera) as a model system of evolution and maintenance of special¬ized interactions. The flies are obligate parasites of the seeds, but also pollinate the plant - it was thus proposed that both species are mutually dependent. Contrasting with the paradigm used for two decades of research on this system, we show that the female fitness component of the plant is similar in the populations with and without Chiastocheta. The plant is thus not exclusively dependent on the flies for reproduction. We discuss this result in the context of the factors responsible for the evolution of mutualistic systems. Understanding the evolution of a biological system requires understanding of its phylo- genetic context. Previous studies showed large mismatch between mtDNA phylogeny and morphological taxonomy in Chiastocheta. By using a large set of RAD-sequencing loci, we delineate the species limits that are congruent with morphology, and show that the discordance is best explained by the scenario of mitochondrial capture among fly species. Finally, we examine this system from a phylogeographic perspective, and identify the lack of congruence in spatial genetic structures of the plant and associated insects across their whole geographic range. The flies show lower numbers of spatial genetic groups than the plant, indicating that not all of the plant réfugia were shared by all the fly species or that the migration dynamics homogenized some of the groups. The incongruence in spatial genetic patterns indicates that fly migrations were largely independent from the genetic background of the plant, following rather a scenario of resource tracking, without the signature of coevolutionary process at this scale. Indeed, while the flies require the plant to survive climatic oscillations, the opposite is not true. Eventually, we show that there is no phylogenetic signal of spatial genetic structures, meaning that neither histories nor life- history traits are shared among closely related species and that species are characterized by unique trajectories of their genes. -- Cette étude vise à comprendre les processus évolutifs à l'oeuvre au sein d'interactions en¬tre espèces spécialisées. Du point de vue macroévolutif, la coévolution entre les taxons spécialisée a été considérée comme l'un des principaux processus générateur de biodiversité. Nous contestons cette idée du point de vue théorique et pratique à travers une revue de la littérature. Nous montrons que les hypothèses majeures reliant les processus coévolutifs avec les patterns de diversité au niveau macroévolutif ne prédisent pas nécessairement la co- diversification des lignées et leur spéciation parallèle, ce qui limite l'utilité de l'approche de phylogénie comparative pour étudier les processus coévolutifs . Nous rappelons également le peu d'exemples de dynamique coévolutive à long terme et proposons que la coévolution se produit plutôt dans des intervalles courts, suivis d'ajustements écologiques. Du point empirique, nous nous concentrons sur l'interaction de pollinisation entre le Trolle d'Europe Trollius europaeus (Ranunculaceae) et ses pollinisateurs associés, du genre Chiastocheta (Anthomyiidae; Diptera) en tant que système-modèle pour étudier l'évolution et le maintien des interactions spécialisées. Les mouches sont des parasites obligatoires des semences, mais pollinisent également la plante. Il a donc été proposé que les deux espèces soient mutuellement dépendantes. Contrastant avec le paradigme utilisé pendant deux décennies de recherche sur ce système, nous montrons, que la composante de fitness femelle de la plante est similaire dans les populations avec et sans Chiastocheta. La plante ne dépend donc pas exclusivement de son interaction avec les mouches pour la reproduction. Nous discutons de ce résultat dans le contexte des facteurs responsables de l'évolution des systèmes mutualistes. Comprendre l'évolution d'un système biologique nécessite la compréhension de son con- texte phylogénétique. Des études antérieures ont montré, chez Chiastocheta, de grandes disparités entre les phylogénies obtenues à partir d'ADN mitochondrial et la taxonomie basée sur les critères morphologiques. En utilisant un grand nombre de loci obtenus par RAD-sequencing, nous traçons les limites des espèces, qui concordent avec les car¬actéristiques morphologies, et montrons que la discordance s'explique en fait par un scénario de capture mitochondriale entre espèces de mouches. Enfin, nous examinons le système d'un point de vue phylogéographique, et identi¬fions les incohérences entre structurations génétiques spatiales de la plante et des insectes associés dans toute leur aire de distribution géographique. Les mouches présentent un nombre de groupes génétiques inférieur à la plante, indiquant que tous les refuges de la plante n'étaient pas partagés par toutes les espèces de mouches ou que les dynamiques migratoires ont homogénéisés certains des groupes chez les mouches. Les différences ob¬servées dans les patrons de structuration génétique spatiale indique que les migrations et dispersions des mouches ont été indépendantes du contexte génétique de la plante, et ces dernières ont été uniquement tributaires de la disponibilité des ressources, sans qu'il n'y ait de signature du processus de coévolution à cette échelle. En effet, tandis que les mouches ont besoin de la plante pour survivre aux oscillations climatiques, le contraire n'est pas exact. Finalement, nous montrons qu'il n'y a pas de signal phylogénétique des structurations génétiques spatiales chez les mouches, ce qui signifie que ni l'histoire, ni les traits d'histoire de vie ne sont partagés entre les espèces phylogénétiquement proches et que les espèces sont caractérisées par des trajectoires uniques de leurs gènes.