Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Multidimensional Deprivation in Contemporary Switzerland Across Social Groups and Time

We have investigated the phenomenon of deprivation in contemporary Switzerland through the adoption of a multidimensional, dynamic approach. By applying Self Organizing Maps (SOM) to a set of 33 non-monetary indicators from the 2009 wave of the Swiss Household Panel (SHP), we identified 13 prototypical forms (or clusters) of well-being, financial vulnerability, psycho-physiological fragility and deprivation within a topological dimensional space. Then new data from the previous waves (2003 to 2008) were classified by the SOM model, making it possible to estimate the weight of the different clusters in time and reconstruct the dynamics of stability and mobility of individuals within the map. Looking at the transition probabilities between year t and year t+1, we observed that the paths of mobility which catalyze the largest number of observations are those connecting clusters that are adjacent on the topological space.

Orthogonal protecting groups in the synthesis of tryptophanyl-hexahydropyrroloindoles

The synthesis of various polycyclic systems containing a C3a-Ni bond between a hexahydropyrrolo[2,3-b]indole and an indole tryptophan is described here. A series of experiments were performed to determine the best combination of five orthogonal protecting groups and the best reaction conditions for formation of said bond, which is a common feature among many recently discovered marine natural products.

Metal dealing at the origin of the Chordata phylum: The metallothionein system and metal overload response in Amphioxus.

Non-vertebrate chordates, specifically amphioxus, are considered of the utmost interest for gaining insight into the evolutionary trends, i.e. differentiation and specialization, of gene/protein systems. In this work, MTs (metallothioneins), the most important metal binding proteins, are characterized for the first time in the cephalochordate subphylum at both gene and protein level, together with the main features defining the amphioxus response to cadmium and copper overload. Two MT genes (BfMT1 and BfMT2) have been identified in a contiguous region of the genome, as well as several ARE (antioxidant response element) and MRE (metal response element) located upstream the transcribed region. Their corresponding cDNAs exhibit identical sequence in the two lancelet species (B. floridae and B. lanceolatum), BfMT2 cDNA resulting from an alternative splicing event. BfMT1 is a polyvalent metal binding peptide that coordinates any of the studied metal ions (Zn, Cd or Cu) rendering complexes stable enough to last in physiological environments, which is fully concordant with the constitutive expression of its gene, and therefore, with a metal homeostasis housekeeping role. On the contrary, BfMT2 exhibits a clear ability to coordinate Cd(II) ions, while it is absolutely unable to fold into stable Cu (I) complexes, even as mixed species. This identifies it as an essential detoxification agent, which is consequently only induced in emergency situations. The cephalochordate MTs are not directly related to vertebrate MTs, neither by gene structure, protein similarity nor metal-binding behavior of the encoded peptides. The closest relative is the echinoderm MT, which confirm proposed phylogenetic relationships between these two groups. The current findings support the existence in most organisms of two types of MTs as for their metal binding preferences, devoted to different biological functions: multivalent MTs for housekeeping roles, and specialized MTs that evolve either as Cd-thioneins or Cu-thioneins, according to the ecophysiological needs of each kind of organisms.

Systematics and phylogeography of the Dysdera erythrina species complex (Araneae, Dysderidae) in Sardinia

Sardinia is the second largest island in the Mediterranean and, together with Corsica and nearby mainland areas, one of the top biodiversity hotspots in the region. The origin of Sardinia traces back to the opening of the western Mediterranean in the late Oligocene. This geological event and the subsequent Messinian Salinity Crisis and Pleistocene glacial cycles have had a major impact on local biodiversity. The Dysdera woodlouse hunter spiders are one of the most diverse ground-dweller groups in the Mediterranean. Here we describe the first two species of this genus endemic to Sardinia: Dysdera jana sp. n. and Dysdera shardana sp. n. The two species show contrasting allopatric distribution: D. jana sp. n. is a narrow endemic while D. shardana sp. n. is distributed throughout most of the island. A multi-gene DNA sequence phylogenetic analys based on mitochondrial and nuclear genes supports the close relationships of the new species to the type species of the genus Dysdera erythrina. Age estimates reject Oligocene origin of the new Dysdera species and identify the Messinian Salinity Crises as the most plausible period for the split between Sardinian endemics and their closest relatives. Phylogeographic analysis reveals deep genetic divergences and population structure in Dysdera shardana sp. n., suggesting that restriction to gene flow probably due to environmental factors could explain local speciation events. Taxonomy, phylogeny, DNA sequencing, Mediterranean biogeography, phylogeography

Planarian (Platyhelminthes, Tricladida) Diversity and Molecular Markers: A New View of an Old Group

Planarians are a group of free-living platyhelminths (triclads) best-known largely due to long-standing regeneration and pattern formation research. However, the group"s diversity and evolutionary history has been mostly overlooked. A few taxonomists have focused on certain groups, resulting in the description of many species and the establishment of higher-level groups within the Tricladida. However, the scarcity of morphological features precludes inference of phylogenetic relationships among these taxa. The incorporation of molecular markers to study their diversity and phylogenetic relationships has facilitated disentangling many conundrums related to planarians and even allowed their use as phylogeographic model organisms. Here, we present some case examples ranging from delimiting species in an integrative style, and barcoding them, to analysing their evolutionary history on a lower scale to infer processes affecting biodiversity origin, or on a higher scale to understand the genus level or even higher relationships. In many cases, these studies have allowed proposing better classifications and resulted in taxonomical changes. We also explain shortcomings resulting in a lack of resolution or power to apply the most up-to-date data analyses. Next-generation sequencing methodologies may help improve this situation and accelerate their use as model organisms.

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).

Marked increase in incidence for bloodstream infections due to Escherichia coli, a side effect of previous antibiotic therapy in the elderly.

We conducted a survey including 3334 bloodstream infections (BSIs) due to E. coli diagnosed in 2005-2014 at a stable cohort of hospitals. Marked increases in incidence were observed for community-acquired (CA) BSIs in patients aged >75 years, CA-BSIs of digestive origin in patients aged 60-74 years, healthcare-associated BSIs, and BSIs associated with ESBL (extended-spectrum B-lactamase)-producing E. coli (ESBLEc). Using MLST, we studied the genetic diversity of 412 BSI isolates recovered during the 2014 survey: 7 major sequence type complexes (STCs) were revealed in phylogenetic group B2, 3 in group A/B1 and 2 in group D. Among the 31 ESBLEc isolates, 1/3 belonged to STC 131. We searched for possible associations between clonal groups, clinical determinants and characteristics of BSIs: isolates from groups B2 (except STC 131) and D were susceptible to antibiotics and associated with BSIs of urinary origin in patients <60 years. STC 131 and group A/B1 isolates were multi-drug resistant and associated with CA-BSIs of digestive origin in patients aged 60-74 with a recent history of antibiotic treatment. STC 131 isolates were associated with HCA-BSIs in patients with recent/present hospitalization in a long-stay unit. We provide a unique population-based picture of the epidemiology of E. coli BSI. The aging nature of the population led to an increase in the number of cases caused by the B2 and D isolates generally implicated in BSIs. In addition, the association of a trend toward increasing rates of gut colonization with multi drug-resistant isolates revealed by the rise in the incidence of BSIs of digestive origin caused by STC 131 and A/B1 (STCs 10, 23, and 155) isolates, and a significant increase in the frequency of BSIs in elderly patients with recent antibiotic treatment suggested that antibiotic use may have contributed to the growing incidence of BSI.

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.

Detecting patterns of species diversification in the presence of both rate shifts and mass extinctions.

BACKGROUND: Recent methodological advances allow better examination of speciation and extinction processes and patterns. A major open question is the origin of large discrepancies in species number between groups of the same age. Existing frameworks to model this diversity either focus on changes between lineages, neglecting global effects such as mass extinctions, or focus on changes over time which would affect all lineages. Yet it seems probable that both lineages differences and mass extinctions affect the same groups. RESULTS: Here we used simulations to test the performance of two widely used methods under complex scenarios of diversification. We report good performances, although with a tendency to over-predict events with increasing complexity of the scenario. CONCLUSION: Overall, we find that lineage shifts are better detected than mass extinctions. This work has significance to assess the methods currently used to estimate changes in diversification using phylogenetic trees. Our results also point toward the need to develop new models of diversification to expand our capabilities to analyse realistic and complex evolutionary scenarios.

Phylogenetic relationships of Apodemus Kaup, 1829 (Rodentia: Muridae) species in the Eastern Mediterranean inferred from mitochondrial DNA, with emphasis on iranian species

Wood mice of the genus Apodemus are widely distributed in Eurasia, with the Eastern Mediterranean being considered as a hotspot. Indeed, numerous species have been documented in Iran, including A. witherbyi, A. hyrcanicus, A. uralensis, A. avicennicus, A. hermonensis, and A. arianus. In this study, 129 specimens were collected from different Iranian localities and two specimens from Afghanistan. The animals were identified taxonomically and their phylogenetic relationships were investigated using cytochrome b mitochondrial DNA sequences. Five species of the genus Apodemus were identified in Iran, including A. hyrcanicus, A. witherbyi, A. cf. ponticus, A. uralensis, and A. mystacinus, beside, A. pallipes from Afghanistan. This study found no evidence of A. flavicollis or A. sylvaticus in Iran, despite their occurrence in Turkey, shedding doubt on the status of A. flavicollis in Iran, Asia Minor, and the Levant. Phylogenetic analyses imply that A. witherbyi has priority over A. avicennicus, A. hermonensis, and A. iconicus. Estimation of the divergence time for these taxa suggests a separation at around 7.2 Ma for the subgenera Karstomys (including A. mystacinus and A. epimelas) and Sylvaemus (including A. flavicollis, A. sylvaticus, A. uralensis, A. pallipes, A. hyrcanicus, A. witherbyi, and A. cf. ponticus). Within the subgenus Karstomys, the divergence times for A. mystacinus and A. epimelas were between 3.0 and 6.1 Ma, and divergence times within the subgenus Sylvaemus were between 5.2 and 6.9 Ma for A. witherbyi and other species in this subgenus. It is postulated that vicariance events including the uplifting of the Zagros Mountains and Anatolian Plateau in the middle Miocene and climate oscillations during the Messinian Salinity Crisis besides formation of the Hyrcanian tertiary forests during the Neogene probably played substantial roles in the radiation and distribution of the genus Apodemus in the Eastern Mediterranean.

Coev-web: a web platform designed to simulate and evaluate coevolving positions along a phylogenetic tree.

BACKGROUND: Available methods to simulate nucleotide or amino acid data typically use Markov models to simulate each position independently. These approaches are not appropriate to assess the performance of combinatorial and probabilistic methods that look for coevolving positions in nucleotide or amino acid sequences. RESULTS: We have developed a web-based platform that gives a user-friendly access to two phylogenetic-based methods implementing the Coev model: the evaluation of coevolving scores and the simulation of coevolving positions. We have also extended the capabilities of the Coev model to allow for the generalization of the alphabet used in the Markov model, which can now analyse both nucleotide and amino acid data sets. The simulation of coevolving positions is novel and builds upon the developments of the Coev model. It allows user to simulate pairs of dependent nucleotide or amino acid positions. CONCLUSIONS: The main focus of our paper is the new simulation method we present for coevolving positions. The implementation of this method is embedded within the web platform Coev-web that is freely accessible at http://coev.vital-it.ch/, and was tested in most modern web browsers.