Geographical and temporal body size variation in a reptile: roles of sex, ecology, phylogeny and ecology structured in phylogeny.

Geographical body size variation has long interested evolutionary biologists, and a range of mechanisms have been proposed to explain the observed patterns. It is considered to be more puzzling in ectotherms than in endotherms, and integrative approaches are necessary for testing non-exclusive alternative mechanisms. Using lacertid lizards as a model, we adopted an integrative approach, testing different hypotheses for both sexes while incorporating temporal, spatial, and phylogenetic autocorrelation at the individual level. We used data on the Spanish Sand Racer species group from a field survey to disentangle different sources of body size variation through environmental and individual genetic data, while accounting for temporal and spatial autocorrelation. A variation partitioning method was applied to separate independent and shared components of ecology and phylogeny, and estimated their significance. Then, we fed-back our models by controlling for relevant independent components. The pattern was consistent with the geographical Bergmann's cline and the experimental temperature-size rule: adults were larger at lower temperatures (and/or higher elevations). This result was confirmed with additional multi-year independent data-set derived from the literature. Variation partitioning showed no sex differences in phylogenetic inertia but showed sex differences in the independent component of ecology; primarily due to growth differences. Interestingly, only after controlling for independent components did primary productivity also emerge as an important predictor explaining size variation in both sexes. This study highlights the importance of integrating individual-based genetic information, relevant ecological parameters, and temporal and spatial autocorrelation in sex-specific models to detect potentially important hidden effects. Our individual-based approach devoted to extract and control for independent components was useful to reveal hidden effects linked with alternative non-exclusive hypothesis, such as those of primary productivity. Also, including measurement date allowed disentangling and controlling for short-term temporal autocorrelation reflecting sex-specific growth plasticity.

The role of phylogeny in the spatial distributions and assembly of mountain plant communities.

A major challenge in community ecology is a thorough understanding of the processes that govern the assembly and composition of communities in time and space. The growing threat of climate change to the vascular plant biodiversity of fragile ecosystems such as mountains has made it equally imperative to develop comprehensive methodologies to provide insights into how communities are assembled. In this perspective, the primary objective of this PhD thesis is to contribute to the theoretical and methodological development of community ecology, by proposing new solutions to better detect the ecological and evolutionary processes that govern community assembly. As phylogenetic trees provide by far, the most advanced tools to integrate the spatial, ecological and evolutionary dynamics of plant communities, they represent the cornerstone on which this work was based. In this thesis, I proposed new solutions to: (i) reveal trends in community assembly on phylogenies, depicted by the transition of signals at the nodes of the different species and lineages responsible for community assembly, (ii) contribute to evidence the importance of evolutionarily labile traits in the distribution of mountain plant species. More precisely, I demonstrated that phylogenetic and functional compositional turnover in plant communities was driven by climate and human land use gradients mostly influenced by evolutionarily labile traits, (iii) predict and spatially project the phylogenetic structure of communities using species distribution models, to identify the potential distribution of phylogenetic diversity, as well as areas of high evolutionary potential along elevation. The altitudinal setting of the Diablerets mountains (Switzerland) provided an appropriate model for this study. The elevation gradient served as a compression of large latitudinal variations similar to a collection of islands within a single area, and allowed investigations on a large number of plant communities. Overall, this thesis highlights that stochastic and deterministic environmental filtering processes mainly influence the phylogenetic structure of plant communities in mountainous areas. Negative density-dependent processes implied through patterns of phylogenetic overdispersion were only detected at the local scale, whereas environmental filtering implied through phylogenetic clustering was observed at both the regional and local scale. Finally, the integration of indices of phylogenetic community ecology with species distribution models revealed the prospects of providing novel and insightful explanations on the potential distribution of phylogenetic biodiversity in high mountain areas. These results generally demonstrate the usefulness of phylogenies in inferring assembly processes, and are worth considering in the theoretical and methodological development of tools to better understand phylogenetic community structure.

Pyrroloquinoline quinone biosynthesis gene pqqC, a novel molecular marker for studying the phylogeny and diversity of phosphate-solubilizing pseudomonads.

Many root-colonizing pseudomonads are able to promote plant growth by increasing phosphate availability in soil through solubilization of poorly soluble rock phosphates. The major mechanism of phosphate solubilization by pseudomonads is the secretion of gluconic acid, which requires the enzyme glucose dehydrogenase and its cofactor pyrroloquinoline quinone (PQQ). The main aim of this study was to evaluate whether a PQQ biosynthetic gene is suitable to study the phylogeny of phosphate-solubilizing pseudomonads. To this end, two new primers, which specifically amplify the pqqC gene of the Pseudomonas genus, were designed. pqqC fragments were amplified and sequenced from a Pseudomonas strain collection and from a natural wheat rhizosphere population using cultivation-dependent and cultivation-independent approaches. Phylogenetic trees based on pqqC sequences were compared to trees obtained with the two concatenated housekeeping genes rpoD and gyrB. For both pqqC and rpoD-gyrB, similar main phylogenetic clusters were found. However, in the pqqC but not in the rpoD-gyrB tree, the group of fluorescent pseudomonads producing the antifungal compounds 2,4-diacetylphloroglucinol and pyoluteorin was located outside the Pseudomonas fluorescens group. pqqC sequences from isolated pseudomonads were differently distributed among the identified phylogenetic groups than pqqC sequences derived from the cultivation-independent approach. Comparing pqqC phylogeny and phosphate solubilization activity, we identified one phylogenetic group with high solubilization activity. In summary, we demonstrate that the gene pqqC is a novel molecular marker that can be used complementary to housekeeping genes for studying the diversity and evolution of plant-beneficial pseudomonads.

A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network.

BACKGROUND: The nuclear receptors are a large family of eukaryotic transcription factors that constitute major pharmacological targets. They exert their combinatorial control through homotypic heterodimerisation. Elucidation of this dimerisation network is vital in order to understand the complex dynamics and potential cross-talk involved. RESULTS: Phylogeny, protein-protein interactions, protein-DNA interactions and gene expression data have been integrated to provide a comprehensive and up-to-date description of the topology and properties of the nuclear receptor interaction network in humans. We discriminate between DNA-binding and non-DNA-binding dimers, and provide a comprehensive interaction map, that identifies potential cross-talk between the various pathways of nuclear receptors. CONCLUSION: We infer that the topology of this network is hub-based, and much more connected than previously thought. The hub-based topology of the network and the wide tissue expression pattern of NRs create a highly competitive environment for the common heterodimerising partners. Furthermore, a significant number of negative feedback loops is present, with the hub protein SHP [NR0B2] playing a major role. We also compare the evolution, topology and properties of the nuclear receptor network with the hub-based dimerisation network of the bHLH transcription factors in order to identify both unique themes and ubiquitous properties in gene regulation. In terms of methodology, we conclude that such a comprehensive picture can only be assembled by semi-automated text-mining, manual curation and integration of data from various sources.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Molecular phylogenetics of soricid shrews (Mammalia) based on mitochondrial cytochrome b gene sequences: with special reference to the Soricinae

Molecular phylogeny of soricid shrews (Soricidae, Eulipotyphla, Mammalia) based on 1140 bp mitochondrial cytochrome b gene (cytb) sequences was inferred by the maximum likelihood (ML) method. All 13 genera of extant Soricinae and two genera of Crocidurinae were included in the analyses. Anourosorex was phylogenetically distant from the main groupings within Soricinae and Crocidurinae in the ML tree. Thus, it could not be determined to which subfamily Anourosorex should be assigned: Soricinae, Crocidurinae or a new subfamily. Soricinae (excluding Anourosorex) should be divided into four tribes: Neomyini, Notiosoricini, Soricini and Blarinini. However, monophyly of Blarinini was not robust in the present data set. Also, branching orders among tribes of Soricinae and those among genera of Neomyini could not be determined because of insufficient phylogenetic information of the cytb sequences. For water shrews of Neomyini (Chimarrogale, Nectogale and Neomys), monophyly of Neomys and the Chimarrogale-Nectogale group could not be verified, which implies the possibility of multiple origins for the semi-aquatic mode of living among taxa within Neomyini. Episoriculus may contain several separate genera. Blarinella was included in Blarinini not Soricini, based on the cytb sequences, but the confidence level was rather low; hence more phylogenetic information is needed to determine its phylogenetic position. Furthermore, some specific problems of taxonomy of soricid shrews were clarified, for example phylogeny of local populations of Notiosorex crawfordi, Chimarrogale himalayica and Crocidura attenuata.

Phylogenetics of Olea (Oleaceae) based on plastid and nuclear ribosomal DNA sequences: tertiary climatic shifts and lineage differentiation times.

BACKGROUND AND AIMS: The genus Olea (Oleaceae) includes approx. 40 taxa of evergreen shrubs and trees classified in three subgenera, Olea, Paniculatae and Tetrapilus, the first of which has two sections (Olea and Ligustroides). Olive trees (the O. europaea complex) have been the subject of intensive research, whereas little is known about the phylogenetic relationships among the other species. To clarify the biogeographical history of this group, a molecular analysis of Olea and related genera of Oleaceae is thus necessary. METHODS: A phylogeny was built of Olea and related genera based on sequences of the nuclear ribosomal internal transcribed spacer-1 and four plastid regions. Lineage divergence and the evolution of abaxial peltate scales, the latter character linked to drought adaptation, were dated using a Bayesian method. KEY RESULTS: Olea is polyphyletic, with O. ambrensis and subgenus Tetrapilus not sharing a most recent common ancestor with the main Olea clade. Partial incongruence between nuclear and plastid phylogenetic reconstructions suggests a reticulation process in the evolution of subgenus Olea. Estimates of divergence times for major groups of Olea during the Tertiary were obtained. CONCLUSIONS: This study indicates the necessity of revising current taxonomic boundaries in Olea. The results also suggest that main lines of evolution were promoted by major Tertiary climatic shifts: (1) the split between subgenera Olea and Paniculatae appears to have taken place at the Miocene-Oligocene boundary; (2) the separation of sections Ligustroides and Olea may have occurred during the Early Miocene following the Mi-1 glaciation; and (3) the diversification within these sections (and the origin of dense abaxial indumentum in section Olea) was concomitant with the aridification of Africa in the Late Miocene.

BadiRate: Estimating Family Turnover Rates by Likelihood-Based Methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Phylogenetic study of the Grimmia Hedw. (Grimmiales, Bryopsida) based on a combination of morphological and molecular characters

Résumé Les mousses sont la plus ancienne lignée de plantes terrestres et leur longue évolution a été accompagnée par des tendances à la simplification des caractères morphologiques. Ce phénomène a quelque peu compliqué les reconstructions phylogénétiques basées sur la morphologie. Les analyses génétiques ont permis de donner de nouvelles informations dans le cadre des analyses phylogénétiques et une réévaluation de certains caractères morphologiques. La plupart des études combinant les données morphologiques et moléculaires ne concernent que des niveaux systématiques élevés comme l'ordre ou la famille et très peu considèrent le niveau du genre. La présente étude tend à tester les relations phylogénétiques du genre Grimmia à l'aide d'une combinaison de caractères morphologiques et moléculaires. Les 40 espèces de Grimmia utilisées dans la première partie de cette étude représentent la majorité des espèces trouvées en Eurasie, un des centres de diversification du genre. Lors de l'analyse morphologique, 52 caractères morphologiques/anatomiques (33 du gamétophyte et 19 du sporophyte) ont été numérisés. Malgré le peu de support statistique des arbres, la topologie des arbres est stable. Les Grimmia, comme décrit précédemment, sont paraphylétiques. Trois clades, correspondant respectivement aux sous-genres Rhabdogrimmia Limpr, Litoneuron I.Hagen et Gasterogrimmia Schimp. sont présents, tandis que le restant des taxons appartenant aux Grimmia forment un groupe non-résolu et indistinct des autres Grimmiaceae. Les séquences chloroplastiques trnL-trnF et rps4 combinés à la morphologie ont été ensuite utilisés pour reconstruire la phylogénie des Grimmia. Les arbres obtenus soutiennent la monophylie des Grimmiaceae tandis que les Grimmia, sont paraphylétiques. Deux clades principaux correspondant aux "Rhabdogrimmia" et aux "Grimmia" se détachent. Seules les espèces de "Rhabdogrimmia" produisent des gemmules foliaires (reproduction asexuée). Dans une étude considèrant 91 séquences trrIL-trnF les espèces appartenant aux "Rhabdogrimrnia" (reproduction asexuée essentiellement) ont des variabilités intraspécifique très faible et interspécifique relativement élevée tandis que les "Grimmia" possèdent la tendance inverse (plus de reproduction sexuée). Summary The mosses are a very old land plant lineage and their long evolutionary history has been accompanied by a trend of morphological character simplifications. This phenomenon has somewhat complicated morphological based phylogenetic reconstructions. Genetic analyses have provided new insights for phylogenetic studies, and have allowed morphological data to be re¬evaluated. Most of the studies combining morphological and molecular data have concerned the higher systematic levels of order and family and only have few considered the genus. The present study aims to test the phylogenetic relationships of the genus Grimmia using a combination of morphological and molecular characters. The 40 chosen Grimmia species represent the majority of those found in Eurasia, one diversification centers of the genus. For the morphological analysis, 52 morphological/anatomical characters (33 gametophyte and 19 sporophyte characters) were numerized. Although the internal statistical support was relatively low, the tree topologies were stable. Grimmia as currently defined was found to be paraphyletic. Three subclades, corresponding to the subgenera Rhabdogrimmia Limpr., Litoneuron I.Hagen, and Gasterogrimmia Schimp. were observed in the trees, while the reminder of the Grimmia species formed an unresolved group indistinct from other Grimmiaceae. Chloroplast (trnL-trnF and rps4) DNA sequences combined with morphology were used to reconstruct the phylogeny of Grimmia. The resulting trees supported the monophyly of Grimmiaceae and that the genus Grimmia, as currently defined, as paraphyletic. Two main clades were resolved corresponding to "Rhabdogrimmia" and "Grimmia". The species belonging to "Rhabdogrimmia" produce foliar-gemmae (asexual reproduction). In a study using 91 sequences of trnL-trnF,"Rhabdogrimmia" species (mainly asexual reproduction) have very low intraspecific variability and high interspecific variability whereas the "Grimmia" species possess the inverse tendency.

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau.

BACKGROUND: The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section. RESULTS: We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns. CONCLUSION: Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

Study of criteria influencing the success rate of DNA swabs in operational conditions: a contribution to an evidence-based approach to crime scene investigation and triage

DNA is nowadays swabbed routinely to investigate serious and volume crimes, but research remains scarce when it comes to determining the criteria that may impact the success rate of DNA swabs taken on different surfaces and situations. To investigate these criteria in fully operational conditions, DNA analysis results of 4772 swabs taken by the forensic unit of a police department in Western Switzerland over a 2.5-year period (2012-2014) in volume crime cases were considered. A representative and random sample of 1236 swab analyses was extensively examined and codified, describing several criteria such as whether the swabbing was performed at the scene or in the lab, the zone of the scene where it was performed, the kind of object or surface that was swabbed, whether the target specimen was a touch surface or a biological fluid, and whether the swab targeted a single surface or combined different surfaces. The impact of each criterion and of their combination was assessed in regard to the success rate of DNA analysis, measured through the quality of the resulting profile, and whether the profile resulted in a hit in the national database or not. Results show that some situations - such as swabs taken on door and window handles for instance - have a higher success rate than average swabs. Conversely, other situations lead to a marked decrease in the success rate, which should discourage further analyses of such swabs. Results also confirm that targeting a DNA swab on a single surface is preferable to swabbing different surfaces with the intent to aggregate cells deposited by the offender. Such results assist in predicting the chance that the analysis of a swab taken in a given situation will lead to a positive result. The study could therefore inform an evidence-based approach to decision-making at the crime scene (what to swab or not) and at the triage step (what to analyse or not), contributing thus to save resource and increase the efficiency of forensic science efforts.

Phylogenetic analyses of Cladophora vagabunda (L.) C. Hoek (Cladophorales, Chlorophyta) from Brazil based on SSU rDNA sequences

The complete SSU rDNA was sequenced for 10 individuals of Cladophora vagabunda collected along the coast of Brazil. For C. rupestris (L.) Kütz. a partial SSU rDNA sequence (1634 bp) was obtained. Phylogenetic trees indicate that Cladophora is paraphyletic, but the section Glomeratae sensu lato including C. vagabunda from Brazil, Japan and France, C. albida (Nees) Kütz., C. sericea (Hudson) Kütz., and C. glomerata (L.) Kütz. is monophyletic. Within this group C. vagabunda is paraphyletic. The sequence identity for the SSU rDNA varied from 98.9% to 100% for the Brazilian C. vagabunda, and from 98.3% to 99.7% comparing the Brazilian individuals to the ones from France and Japan. Sequence identity of the Brazilian C. vagabunda to C. albida and C. sericea vary from 98.0% to 98.6%. The SSU rDNA phylogeny support partially the morphological characteristics presented by Brazilian populations of C. vagabunda. On the other hand, C. rupestris from Brazil does not group with C. rupestris from France, both sequences presenting only 96.9% of identity. The inclusion of sequences of individuals from Brazil reinforces the need of taxonomical revision for the genus Cladophora and for the complex C. vagabunda.

Impact of specimen adequacy on the assessment of renal allograft biopsy specimens

The Banff classification was introduced to achieve uniformity in the assessment of renal allograft biopsies. The primary aim of this study was to evaluate the impact of specimen adequacy on the Banff classification. All renal allograft biopsies obtained between July 2010 and June 2012 for suspicion of acute rejection were included. Pre-biopsy clinical data on suspected diagnosis and time from renal transplantation were provided to a nephropathologist who was blinded to the original pathological report. Second pathological readings were compared with the original to assess agreement stratified by specimen adequacy. Cohen's kappa test and Fisher's exact test were used for statistical analyses. Forty-nine specimens were reviewed. Among these specimens, 81.6% were classified as adequate, 6.12% as minimal, and 12.24% as unsatisfactory. The agreement analysis among the first and second readings revealed a kappa value of 0.97. Full agreement between readings was found in 75% of the adequate specimens, 66.7 and 50% for minimal and unsatisfactory specimens, respectively. There was no agreement between readings in 5% of the adequate specimens and 16.7% of the unsatisfactory specimens. For the entire sample full agreement was found in 71.4%, partial agreement in 20.4% and no agreement in 8.2% of the specimens. Statistical analysis using Fisher's exact test yielded a P value above 0.25 showing that - probably due to small sample size - the results were not statistically significant. Specimen adequacy may be a determinant of a diagnostic agreement in renal allograft specimen assessment. While additional studies including larger case numbers are required to further delineate the impact of specimen adequacy on the reliability of histopathological assessments, specimen quality must be considered during clinical decision making while dealing with biopsy reports based on minimal or unsatisfactory specimens.

Suppression of long-branch attraction artefacts in the animal phylogeny using a site-heterogeneous model

Affiliation: Département de Biochimie, Université de Montréal