Phylogeny of subtribe Gentianinae (Gentianaceae): biogeographic inferences despite limitations in temporal calibration points

The subtribe Gentianinae comprises ca. 425 species, most of them within the well-studied genus Gentiana and mainly distributed over the Eurasian continent. Phylogenetic relationships between Gentiana and its closest relatives, the climbing gentians (Crawfurdia, Tripterospermum) and the new genus Metagentiana, remain unclear. All three genera were recently found to be polyphyletic, possibly because of poor sampling of Tripterospermum and Crawfurdia. Highest diversity of Gentianinae occurs in the western Himalaya, but the absence of uncontroversial fossil evidence limits our understanding of its biogeography. In the present study, we generated ITS and atpB-rbcL sequences for 19 species of Tripterospermum, 9 of Crawfurdia and 11 of Metagentiana, together representing about 60 percent of the species diversity of these genera. Our results show that only Metagentiana is polyphyletic and divided into three monophyletic entities. No unambiguous synapomorphies were associated with the three Metagentiana entities. Different combinations of three approximate calibration points were used to generate three divergence time estimation scenarios. Although dating hypotheses were mostly inconsistent, they concurred in associating radiation of Gentiana to an orogenic phase of the Himalaya between 15 and 10 million years ago. Our study illustrates the conceptual difficulties in addressing the time frame of diversification in a group lacking sufficient fossil number and quality.

An empirical hierarchical Bayesian unification of occupational exposure assessment methods.

In occupational exposure assessment of airborne contaminants, exposure levels can either be estimated through repeated measurements of the pollutant concentration in air, expert judgment or through exposure models that use information on the conditions of exposure as input. In this report, we propose an empirical hierarchical Bayesian model to unify these approaches. Prior to any measurement, the hygienist conducts an assessment to generate prior distributions of exposure determinants. Monte-Carlo samples from these distributions feed two level-2 models: a physical, two-compartment model, and a non-parametric, neural network model trained with existing exposure data. The outputs of these two models are weighted according to the expert's assessment of their relevance to yield predictive distributions of the long-term geometric mean and geometric standard deviation of the worker's exposure profile (level-1 model). Bayesian inferences are then drawn iteratively from subsequent measurements of worker exposure. Any traditional decision strategy based on a comparison with occupational exposure limits (e.g. mean exposure, exceedance strategies) can then be applied. Data on 82 workers exposed to 18 contaminants in 14 companies were used to validate the model with cross-validation techniques. A user-friendly program running the model is available upon request.

Therapeutic drug monitoring of imatinib: Bayesian and alternative methods to predict trough levels

Background: The imatinib trough plasma concentration (C(min)) correlates with clinical response in cancer patients. Therapeutic drug monitoring (TDM) of plasma C(min) is therefore suggested. In practice, however, blood sampling for TDM is often not performed at trough. The corresponding measurement is thus only remotely informative about C(min) exposure. Objectives: The objectives of this study were to improve the interpretation of randomly measured concentrations by using a Bayesian approach for the prediction of C(min), incorporating correlation between pharmacokinetic parameters, and to compare the predictive performance of this method with alternative approaches, by comparing predictions with actual measured trough levels, and with predictions obtained by a reference method, respectively. Methods: A Bayesian maximum a posteriori (MAP) estimation method accounting for correlation (MAP-ρ) between pharmacokinetic parameters was developed on the basis of a population pharmacokinetic model, which was validated on external data. Thirty-one paired random and trough levels, observed in gastrointestinal stromal tumour patients, were then used for the evaluation of the Bayesian MAP-ρ method: individual C(min) predictions, derived from single random observations, were compared with actual measured trough levels for assessment of predictive performance (accuracy and precision). The method was also compared with alternative approaches: classical Bayesian MAP estimation assuming uncorrelated pharmacokinetic parameters, linear extrapolation along the typical elimination constant of imatinib, and non-linear mixed-effects modelling (NONMEM) first-order conditional estimation (FOCE) with interaction. Predictions of all methods were finally compared with 'best-possible' predictions obtained by a reference method (NONMEM FOCE, using both random and trough observations for individual C(min) prediction). Results: The developed Bayesian MAP-ρ method accounting for correlation between pharmacokinetic parameters allowed non-biased prediction of imatinib C(min) with a precision of ±30.7%. This predictive performance was similar for the alternative methods that were applied. The range of relative prediction errors was, however, smallest for the Bayesian MAP-ρ method and largest for the linear extrapolation method. When compared with the reference method, predictive performance was comparable for all methods. The time interval between random and trough sampling did not influence the precision of Bayesian MAP-ρ predictions. Conclusion: Clinical interpretation of randomly measured imatinib plasma concentrations can be assisted by Bayesian TDM. Classical Bayesian MAP estimation can be applied even without consideration of the correlation between pharmacokinetic parameters. Individual C(min) predictions are expected to vary less through Bayesian TDM than linear extrapolation. Bayesian TDM could be developed in the future for other targeted anticancer drugs and for the prediction of other pharmacokinetic parameters that have been correlated with clinical outcomes.

Regional-scale integration of multi-scale hydrological and geophysical data using a two-step Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale for the purpose of improving predictions of groundwater flow and solute transport. However, extending corresponding approaches to the regional scale still represents one of the major challenges in the domain of hydrogeophysics. To address this problem, we have developed a regional-scale data integration methodology based on a two-step Bayesian sequential simulation approach. Our objective is to generate high-resolution stochastic realizations of the regional-scale hydraulic conductivity field in the common case where there exist spatially exhaustive but poorly resolved measurements of a related geophysical parameter, as well as highly resolved but spatially sparse collocated measurements of this geophysical parameter and the hydraulic conductivity. To integrate this multi-scale, multi-parameter database, we first link the low- and high-resolution geophysical data via a stochastic downscaling procedure. This is followed by relating the downscaled geophysical data to the high-resolution hydraulic conductivity distribution. After outlining the general methodology of the approach, we demonstrate its application to a realistic synthetic example where we consider as data high-resolution measurements of the hydraulic and electrical conductivities at a small number of borehole locations, as well as spatially exhaustive, low-resolution estimates of the electrical conductivity obtained from surface-based electrical resistivity tomography. The different stochastic realizations of the hydraulic conductivity field obtained using our procedure are validated by comparing their solute transport behaviour with that of the underlying ?true? hydraulic conductivity field. We find that, even in the presence of strong subsurface heterogeneity, our proposed procedure allows for the generation of faithful representations of the regional-scale hydraulic conductivity structure and reliable predictions of solute transport over long, regional-scale distances.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Use of phylogeny to resolve the taxonomy of the widespread and highly polymorphic African giant shrews (Crocidura olivieri group, Crocidurinae, Mammalia).

The aim of this study is to provide a better understanding of the genetic relationships within the widespread and highly polymorphic group of African giant shrews (Crocidura olivieri group). We sequenced 769 base pairs (bp) of the mitochondrial cytochrome b gene and 472 bp of the mitochondrial control region over the entire geographic range from South Africa to Morocco. The analyses reveal four main clades associated with different biomes. The largest clade occurs over a range covering Northwest and Central Africa and includes samples of C. fulvastra, C. olivieri, and C. viaria. The second clade is composed of C. goliath from Gabon, while South African C. flavescens, and C. hirta form two additional clades. On the basis of these results, the validity of some taxa in the C. olivieri group should be re-evaluated.

A Bayesian network approach to the database serch problem in criminal proceedings

Background The 'database search problem', that is, the strengthening of a case - in terms of probative value - against an individual who is found as a result of a database search, has been approached during the last two decades with substantial mathematical analyses, accompanied by lively debate and centrally opposing conclusions. This represents a challenging obstacle in teaching but also hinders a balanced and coherent discussion of the topic within the wider scientific and legal community. This paper revisits and tracks the associated mathematical analyses in terms of Bayesian networks. Their derivation and discussion for capturing probabilistic arguments that explain the database search problem are outlined in detail. The resulting Bayesian networks offer a distinct view on the main debated issues, along with further clarity. Methods As a general framework for representing and analyzing formal arguments in probabilistic reasoning about uncertain target propositions (that is, whether or not a given individual is the source of a crime stain), this paper relies on graphical probability models, in particular, Bayesian networks. This graphical probability modeling approach is used to capture, within a single model, a series of key variables, such as the number of individuals in a database, the size of the population of potential crime stain sources, and the rarity of the corresponding analytical characteristics in a relevant population. Results This paper demonstrates the feasibility of deriving Bayesian network structures for analyzing, representing, and tracking the database search problem. The output of the proposed models can be shown to agree with existing but exclusively formulaic approaches. Conclusions The proposed Bayesian networks allow one to capture and analyze the currently most well-supported but reputedly counter-intuitive and difficult solution to the database search problem in a way that goes beyond the traditional, purely formulaic expressions. The method's graphical environment, along with its computational and probabilistic architectures, represents a rich package that offers analysts and discussants with additional modes of interaction, concise representation, and coherent communication.

Phylogeny and circumscription of Sapindaceae revisited: molecular sequence data, morphology and biogeography support recognition of a new family, Xanthoceraceae

Background and aims Recent studies have adopted a broad definition of Sapindaceae that includes taxa traditionally placed in Aceraceae and Hippocastanaceae, achieving monophyly but yielding a family difficult to characterize and for which no obvious morphological synapomorphy exists. This expanded circumscription was necessitated by the finding that the monotypic, temperate Asian genus Xanthoceras, historically placed in Sapindaceae tribe Harpullieae, is basal within the group. Here we seek to clarify the relationships of Xanthoceras based on phylogenetic analyses using a dataset encompassing nearly 3/4 of sapindaceous genera, comparing the results with information from morphology and biogeography, in particular with respect to the other taxa placed in Harpullieae. We then re-examine the appropriateness of maintaining the current broad, morphologically heterogeneous definition of Sapindaceae and explore the advantages of an alternative family circumscription. Methods Using 243 samples representing 104 of the 142 currently recognized genera of Sapindaceae s. lat. (including all in Harpullieae), sequence data were analyzed for nuclear (ITS) and plastid (matK, rpoB, trnD-trnT, trnK-matK, trnL-trnF and trnS-trnG) markers, adopting the methodology of a recent family-wide study, performing single-gene and total evidence analyses based on maximum likelihood (ML) and maximum parsimony (MP) criteria, and applying heuristic searches developed for large datasets, viz, a new strategy implemented in RAxML (for ML) and the parsimony ratchet (for MP). Bootstrap analyses were performed for each method to test for congruence between markers. Key results Our findings support earlier suggestions that Harpullieae are polyphyletic: Xanthoceras is confirmed as sister to all other sampled taxa of Sapindaceae s. lat.; the remaining members belong to three other clades within Sapindaceae s. lat., two of which correspond respectively to the groups traditionally treated as Aceraceae and Hippocastanaceae, together forming a clade sister to the largely tropical Sapindaceae s. str., which is monophyletic and morphologically coherent provided Xanthoceras is excluded. Conclusion To overcome the difficulties of a broadly circumscribed Sapindaceae, we resurrect the historically recognized temperate families Aceraceae and Hippocastanaceae, and describe a new family, Xanthoceraceae, thus adopting a monophyletic and easily characterized circumscription of Sapindaceae nearly identical to that used for over a century.

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.

A Bayesian approach for pervasive estimation of breaststroke velocity using a wearable IMU

A ubiquitous assessment of swimming velocity (main metric of the performance) is essential for the coach to provide a tailored feedback to the trainee. We present a probabilistic framework for the data-driven estimation of the swimming velocity at every cycle using a low-cost wearable inertial measurement unit (IMU). The statistical validation of the method on 15 swimmers shows that an average relative error of 0.1 ± 9.6% and high correlation with the tethered reference system (rX,Y=0.91 ) is achievable. Besides, a simple tool to analyze the influence of sacrum kinematics on the performance is provided.

Object-oriented Bayesian networks for evaluating DIP-STR profiling results from unbalanced DNA mixtures.

The genetic characterization of unbalanced mixed stains remains an important area where improvement is imperative. In fact, with current methods for DNA analysis (Polymerase Chain Reaction with the SGM Plus™ multiplex kit), it is generally not possible to obtain a conventional autosomal DNA profile of the minor contributor if the ratio between the two contributors in a mixture is smaller than 1:10. This is a consequence of the fact that the major contributor's profile 'masks' that of the minor contributor. Besides known remedies to this problem, such as Y-STR analysis, a new compound genetic marker that consists of a Deletion/Insertion Polymorphism (DIP), linked to a Short Tandem Repeat (STR) polymorphism, has recently been developed and proposed elsewhere in literature [1]. The present paper reports on the derivation of an approach for the probabilistic evaluation of DIP-STR profiling results obtained from unbalanced DNA mixtures. The procedure is based on object-oriented Bayesian networks (OOBNs) and uses the likelihood ratio as an expression of the probative value. OOBNs are retained in this paper because they allow one to provide a clear description of the genotypic configuration observed for the mixed stain as well as for the various potential contributors (e.g., victim and suspect). These models also allow one to depict the assumed relevance relationships and perform the necessary probabilistic computations.

Implementing statistical learning methods through Bayesian networks (Part 2): bayesian evaluations for results of black toner analyses in forensic document examination

This paper presents and discusses the use of Bayesian procedures - introduced through the use of Bayesian networks in Part I of this series of papers - for 'learning' probabilities from data. The discussion will relate to a set of real data on characteristics of black toners commonly used in printing and copying devices. Particular attention is drawn to the incorporation of the proposed procedures as an integral part in probabilistic inference schemes (notably in the form of Bayesian networks) that are intended to address uncertainties related to particular propositions of interest (e.g., whether or not a sample originates from a particular source). The conceptual tenets of the proposed methodologies are presented along with aspects of their practical implementation using currently available Bayesian network software.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Phylogeny of shrews (Soricidae) : taxonomic and biogeographic implications

Résumé Les Soricidae sont l'une des plus grandes familles de mammifères avec plus de 300 espèces décrites. Elle a été récemment divisée en trois sous-familles, les Soricidae, qui sont distribuées dans la région Holarctique, les Crocidurinae en Afrique et en Eurasie, et les Myosoricinae en Afrique. La diversité spécifique de cette famille a conduit à des interprétations taxonomiques multiples, qui sont à l'origine de polémiques entre spécialistes, et même les premiers résultats moléculaires ont été fortement contradictoires. Le but de cette thèse est donc d'appliquer des meilleures techniques sur des échantillons mieux ciblés, afin de résoudre les contradictions taxonomiques et comprendre l'histoire de cette famille. Par le biais de marqueurs génétiques mitochondriaux et nucléaires, j'ai étudié: (i) Les relations taxonomiques à différent niveaux hiérarchiques au sein des Soricidae, c'est-à dire, entre les sous-familles, tribus, et genres, ainsi qu'au sein de deux complexes d'espèces largement distribués, et d'une espèce européenne, le but étant d'établir la congruence entre les données génétiques et les interprétations morphologiques classiques. (ii) Les relations biogéographiques, soit l'origine potentielle des différentes sous-familles, tribus, et genres, le nombre d'échanges intercontinentaux, ainsi que la structure phylogéographique à un niveau (péri)-spécifique, afin d'établir l'histoire de la diversification de cette famille. Les analyses combinées d'ADN mitochondrial et nucléaire ont montré un rapport clair entre les taxa à un niveau taxonomique élevé, mettant en évidence les rapports entre les sous-familles, les tribus, et les genres. Bien que Myosorex constitue un groupe monophylétique distinct, sa définition en tant que sous-famille séparée ne peut pas être reconnue. Ainsi, nous proposons d'attribuer un niveau de tribu pour ce clade (inclus dans les Crocidurinae). Nous avons également montré l'inclusion du genre Anourosorex dans les Soricinae et non en position basale dans les Soricidae. Au sein des Crocidurinae, Suncus s'est révélé être paraphylétique, et le genre Diplomesodon devrait être considéré d'un point de vue génétique comme invalide, puisque il se trouve au sein du clade du genre Crocidura. À un niveau taxonomique plus bas, nous avons montré la monophylie de deux complexes d'espèces largement distribués, le groupe de C. suaveolens et de C. olivieri. Néanmoins à l'intérieur de ceux-ci, des différences majeures avec la classification morphologique se sont révélées. Par exemples, C. sibirica n'est pas une espèce valide, les analyses de phylogénie moléculaire ne montrant pas de variations génétiques entre celle-ci et un échantillon de la localité type de C. suaveolens. D'un point de vue biogéographique, les fluctuations climatiques et les activités tectoniques des 20 derniers millions d'années ont fortement influencé la diversité actuelle des Soricidae. À un niveau taxonomique élevé, l'apparition de connexions de terre temporaires entre le Vieux et le Nouveau Monde au Miocène moyen ont mené à plusieurs colonisations indépendantes de l'Amérique par les Soricinae. Celles-ci ónt conduit à une diversification d'une tribu (Notiosoricini), ainsi que de genres (par ex: Cryptotis, Blarina) et d'un sous-genre (Otisorex) endémique au Néarctique. Dans le Vieux Monde, les barrières entre l'Afrique et Eurasie étaient plus perméables, menant à plusieurs échanges bidirectionnels de Crocidurinae. La diversification des clades principaux s'est produite au Miocène, certains clades étant endémiques d'Afrique ou d'Eurasie, tandis que d'autres se sont diversifiés à travers le Vieux Monde. À un niveau spécifique ou péri-spécifique, la fluctuation climatique du Pliocène et les glaciations du Pléistocène ont fortement divisé les populations dans tout le Paléarctique, menant à des entités génétiques distinctes. En Europe, les populations du groupe de C. suaveolens ont été divisées en une lignée Sud-Ouest et une Sud-Est, alors qu'au Proche-Orient et au Moyen-Orient, la diversité de clades est plus importante. En conclusion, mes études ont révélé que du Miocène à nos jours, la diversification des Soricidae a été provoquée par la colonisation de nouveaux habitats (dispersion), ainsi que par l'isolement des populations par diverses barrières (vicariance). Abstract The Soricidae is one of the largest mammalian families with more than 300 species described. It has been recently divided into three subfamilies, the Soricinae, which are distributed in the Holartic region, the Crocidurinae in Africa and Eurasia, and the Myosoricinae in Africa. The specific diversity of this family have led to multiple systematic interpretations and controversies between authors. Fortunately, today, cytotaxonomic, allozymic and molecular studies have permitted to clarify some uncertainties. Nevertheless, the Soricidae remains still poorly known. In this thesis, we aim at understanding with the use of mitochondrial and nuclear markers: (i) the taxonomic relationships at different hierarchical levels within Soricidae, i.e., between the subfamilies, tribes, and genera, as well as within two largely distributed species complexes, and within a European species, the goal being to establish congruence between the genetic data and traditional morphological interpretations; (ii) the biogeographic relationships, especially the potential origin of the different subfamilies, tribes, and genera, the number of transcontinental exchanges, as well as the phylogeographic structure at a (peri)-specific level, in order to establish the history of the genetic diversification of this family. The combined analyses of mitochondrial and nuclear DNA highlight for the first time a clear relationship between taxa at a high taxonomical level, permitting to distinguish the relationships between subfamilies, tribes, and genera. Although Myosorex formed a distinct monophyletic group, its definition as a distinct sub-family cannot be advocated. Thus, we propose to attribute a tribe level for this Glade (included within the Crocidurinae). Additionally, this combination of genes pleads in favour of the inclusion of the genus Anourosorex within the Soricinae and not in a basal position within the Soricidae. Within the Crocidurinae, Suncus appeared to be paraphyletic, and Diplomesodon should be considered from a genetic point of view as invalid, and is presently considered as Crocidura. At a lower taxonomic level, we showed the monophyly of two widely distributed species complexes, the C. suaveolens group and the C. olivieri group. Nevertheless within those, we showed major differences compared to morphological classification. For examples, C. sibirica revealed to not be a valid species, the molecular phylogenetic analyses failed to evidence genetical variations between it and samples of the type locality of C. suaveolens. In a biogeographic point of view, the climatic fluctuations and the tectonic plate activities of the last 20 Myr have strongly influenced the actual diversity of the family. At a high taxonomic level, the successive land bridge connections between the Old and the New World, which occurred during the Middle Miocene, have led to several independent colonisations of America by Soricinae, and a subsequent diversification of endemic Nearctic's tribe (Notiosoricini), genera (e.g. Cryptotis, Blaring) and sub-genus (Otisorex) within the Soricinae. Within the Old World, the barriers between Africa and Eurasia were more permeable, leading to several bidirectional exchanges within the Crocidurinae. The diversification of major clades occurred through the Miocene, some clades being endemic to Africa or Eurasia, whereas others diversified through the Old World. At a species level or a peri-specific level, the Pliocene climatic fluctuation and the Pleistocene glaciations have strongly divided the populations throughout the Palaearctic, leading to well defined genetic entities. In Europe, populations of the C. suaveolens group were split in a classical south-western and south-eastern lineage. In contrast, the Near East and the Middle East reveal many differentiated clades. In conclusion, our studies revealed that, from the Miocene to present, the diversification and speciation events within the Soricidae were caused by natural colonisation of new habitats (dispersion) and isolation of populations by various barriers (vicariance).

Using a Bayesian method to assign individuals to karyotypic taxa in shrew hybrid zones.

Individuals sampled in hybrid zones are usually analysed according to their sampling locality, morphology, behaviour or karyotype. But the increasing availability of genetic information more and more favours its use for individual sorting purposes and numerous assignment methods based on the genetic composition of individuals have been developed. The shrews of the Sorex araneus group offer good opportunities to test the genetic assignment on individuals identified by their karyotype. Here we explored the potential and efficiency of a Bayesian assignment method combined or not with a reference dataset to study admixture and individual assignment in the difficult context of two hybrid zones between karyotypic species of the Sorex araneus group. As a whole, we assigned more than 80% of the individuals to their respective karyotypic categories (i.e. 'pure' species or hybrids). This assignment level is comparable to what was obtained for the same species away from hybrid zones. Additionally, we showed that the assignment result for several individuals was strongly affected by the inclusion or not of a reference dataset. This highlights the importance of such comparisons when analysing hybrid zones. Finally, differences between the admixture levels detected in both hybrid zones support the hypothesis of an impact of chromosomal rearrangements on gene flow.