Institut des Sciences de l'Evolution de Montpellier

M. Raymond, F. Rousset; GENEPOP (Version 1.2): Population Genetics Software for Exact Tests and Ecumenicism, Journal of Heredity, Volume 86, Issue 3, 1 May

UNLABELLED: Analysis of Phylogenetics and Evolution (APE) is a package written in the R language for use in molecular evolution and phylogenetics. APE provides both utility functions for reading and writing data and manipulating phylogenetic trees, as well as several advanced methods for phylogenetic and evolutionary analysis (e.g. comparative and population genetic methods). APE takes advantage of the many R functions for statistics and graphics, and also provides a flexible framework for developing and implementing further statistical methods for the analysis of evolutionary processes. AVAILABILITY: The program is free and available from the official R package archive at http://cran.r-project.org/src/contrib/PACKAGES.html#ape. APE is licensed under the GNU General Public License.

Summary: After more than fifteen years of existence, the R package ape has continuously grown its contents, and has been used by a growing community of users. The release of version 5.0 has marked a leap towards a modern software for evolutionary analyses. Efforts have been put to improve efficiency, flexibility, support for 'big data' (R's long vectors), ease of use and quality check before a new release. These changes will hopefully make ape a useful software for the study of biodiversity and evolution in a context of increasing data quantity. Availability and implementation: ape is distributed through the Comprehensive R Archive Network: http://cran.r-project.org/package=ape. Further information may be found at http://ape-package.ird.fr/.

This note summarizes developments of the genepop software since its first description in 1995, and in particular those new to version 4.0: an extended input format, several estimators of neighbourhood size under isolation by distance, new estimators and confidence intervals for null allele frequency, and less important extensions to previous options. genepop now runs under Linux as well as under Windows, and can be entirely controlled by batch calls.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. For example, a key point that needs to be emphasized is thatthere is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process versus those that measure flux through the autophagy pathway (i.e., the completeprocess including the amount and rate of cargo sequestered and degraded). In particular, a block in macroautophagy that results in autophagosome accumulation must be differentiated from stimuli that increase autophagic activity, defined as increasedautophagy induction coupled with increased delivery to, and degradation within, lysosomes (inmost higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in manycases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. It is worth emphasizing here that lysosomal digestion is a stage of autophagy and evaluating its competence is a crucial part of the evaluation of autophagic flux, or complete autophagy. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as forreviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multipleassays to monitor autophagy. Along these lines, because of the potential for pleiotropic effects due to blocking autophagy through genetic manipulation, it is imperative to target by gene knockout or RNA interference more than one autophagyrelated protein. In addition, some individual Atg proteins, or groups of proteins, are involved in other cellular pathways implying that not all Atg proteins can be used as a specific marker for an autophagic process. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular assays, we hope to encourage technical innovation in the field.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.

Abstract Since the publication of our previous time scale (Berggren and others, 1985c = BKFV85) a large amount of new magneto- and biostratigraphic data and radioisotopic ages have become available. An evaluation of some of the key magnetobiostratigraphic calibration points used in BKFV85, as suggested by high precision 40Ar/ 39Ar dating (e.g., Montanari and others, 1988; Swisher and Prothero, 1990; Prothero and Swisher, 1992; Prothero, 1994), has served as a catalyst for us in developing a revised Cenozoic time scale. For the Neogene Period, astrochronologic data (Shackleton and others, 1990; Hilgen, 1991) required re-evaluation of the calibration of the Pliocene and Pleistocene Epochs. The significantly older ages for the Pliocene-Pleistocene Epochs predicted by astronomical calibrations were soon corroborated by high precision 40Ar/ 39Ar dating (e.g., Baksi and others, 1992; McDougall and others, 1992; Tauxe and others, 1992; Walter and others, 1991; Renne and others, 1993). At the same time, a new and improved definition of the Late Cretaceous and Cenozoic polarity sequence was achieved based on a comprehensive evaluation of global sea-floor magnetic anomaly profiles (Cande and Kent, 1992). This, in turn, led to a revised Cenozoic geomagnetic polarity time scale (GPTS) based on standardization to a model of South Atlantic spreading history (Cande and Kent, 1992/1995 = CK92/95). This paper presents a revised (integrated magnetobiochronologic) Cenozoic time scale (IMBTS) based on an assessment and integration of data from several sources. Biostratigraphic events are correlated to the recently revised global polarity time scale (CK95). The construction of the new GPTS is outlined with emphasis on methodology and newly developed polarity history nomenclature. The radioisotopic calibration points (as well as other relevant data) used to constrain the GPTS are reviewed in their (bio)stratigraphic context. An updated magnetobiostratigraphic (re)assessment of about 150 pre-Pliocene planktonic foraminiferal datum events (including recently available high southern (austral) latitude data) and a new/modified zonal biostratigraphy provides an essentially global biostratigraphic correlation framework. This is complemented by a (re)assessment of nearly 100 calcareous nannofossil datum events. Unrecognized unconformities in the stratigraphic record (and to a lesser extent differences in taxonomic concepts), rather than latitudinal diachrony, is shown to account for discrepancies in magnetobiostratigraphic correlations in many instances, particularly in the Paleogene Period. Claims of diachrony of low amplitude (<2 my) are poorly substantiated, at least in the Paleocene and Eocene Epochs. Finally, we (re)assess the current status of Cenozoic chronostratigraphy and present estimates of the chronology of lower (stage) and higher (system) level units. Although the numerical values of chronostratigraphic units (and their boundaries) have changed in the decade since the previous version of the Cenozoic time scale, the relative duration of these units has remained essentially the same. This is particularly true of the Paleogene Period, where the Paleocene/Eocene and Eocene/Oligocene boundaries have been shifted ~2 my younger and the Cretaceous/Paleogene boundary ~1 my younger. Changes in the Neogene time scale are relatively minor and reflect primarily improved magnetobiostratigraphic calibrations, better understanding of chronostratigraphic and magnetobiostratigraphic relationships, and the introduction of a congruent astronomical/paleomagnetic chronology for the past 6 my (and concomitant adjustments to magnetochron age estimates).

Michel Raymond, François Rousset; AN EXACT TEST FOR POPULATION DIFFERENTIATION, Evolution, Volume 49, Issue 6, 1 December 1995, Pages 1280–1283, https://do

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

To better determine the history of modern birds, we performed a genome-scale phylogenetic analysis of 48 species representing all orders of Neoaves using phylogenomic methods created to handle genome-scale data. We recovered a highly resolved tree that confirms previously controversial sister or close relationships. We identified the first divergence in Neoaves, two groups we named Passerea and Columbea, representing independent lineages of diverse and convergently evolved land and water bird species. Among Passerea, we infer the common ancestor of core landbirds to have been an apex predator and confirm independent gains of vocal learning. Among Columbea, we identify pigeons and flamingoes as belonging to sister clades. Even with whole genomes, some of the earliest branches in Neoaves proved challenging to resolve, which was best explained by massive protein-coding sequence convergence and high levels of incomplete lineage sorting that occurred during a rapid radiation after the Cretaceous-Paleogene mass extinction event about 66 million years ago.

Thanks to the development of high-throughput sequencing technologies, target enrichment sequencing of nuclear ultraconserved DNA elements (UCEs) now allows routine inference of phylogenetic relationships from thousands of genomic markers. Recently, it has been shown that mitochondrial DNA (mtDNA) is frequently sequenced alongside the targeted loci in such capture experiments. Despite its broad evolutionary interest, mtDNA is rarely assembled and used in conjunction with nuclear markers in capture-based studies. Here, we developed MitoFinder, a user-friendly bioinformatic pipeline, to efficiently assemble and annotate mitogenomic data from hundreds of UCE libraries. As a case study, we used ants (Formicidae) for which 501 UCE libraries have been sequenced whereas only 29 mitogenomes are available. We compared the efficiency of four different assemblers (IDBA-UD, MEGAHIT, MetaSPAdes, and Trinity) for assembling both UCE and mtDNA loci. Using MitoFinder, we show that metagenomic assemblers, in particular MetaSPAdes, are well suited to assemble both UCEs and mtDNA. Mitogenomic signal was successfully extracted from all 501 UCE libraries, allowing us to confirm species identification using CO1 barcoding. Moreover, our automated procedure retrieved 296 cases in which the mitochondrial genome was assembled in a single contig, thus increasing the number of available ant mitogenomes by an order of magnitude. By utilizing the power of metagenomic assemblers, MitoFinder provides an efficient tool to extract complementary mitogenomic data from UCE libraries, allowing testing for potential mitonuclear discordance. Our approach is potentially applicable to other sequence capture methods, transcriptomic data and whole genome shotgun sequencing in diverse taxa. The MitoFinder software is available from GitHub (https://github.com/RemiAllio/MitoFinder).

Dendroscope 3 is a new program for working with rooted phylogenetic trees and networks. It provides a number of methods for drawing and comparing rooted phylogenetic networks, and for computing them from rooted trees. The program can be used interactively or in command-line mode. The program is written in Java, use of the software is free, and installers for all 3 major operating systems can be downloaded from www.dendroscope.org. [Phylogenetic trees; phylogenetic networks; software.].

Polymyxins are polycationic antimicrobial peptides that are currently the last-resort antibiotics for the treatment of multidrug-resistant, Gram-negative bacterial infections. The reintroduction of polymyxins for antimicrobial therapy has been followed by an increase in reports of resistance among Gram-negative bacteria. Some bacteria, such as Klebsiella pneumoniae, Pseudomonas aeruginosa, and Acinetobacter baumannii, develop resistance to polymyxins in a process referred to as acquired resistance, whereas other bacteria, such as Proteus spp., Serratia spp., and Burkholderia spp., are naturally resistant to these drugs. Reports of polymyxin resistance in clinical isolates have recently increased, including acquired and intrinsically resistant pathogens. This increase is considered a serious issue, prompting concern due to the low number of currently available effective antibiotics. This review summarizes current knowledge concerning the different strategies bacteria employ to resist the activities of polymyxins. Gram-negative bacteria employ several strategies to protect themselves from polymyxin antibiotics (polymyxin B and colistin), including a variety of lipopolysaccharide (LPS) modifications, such as modifications of lipid A with phosphoethanolamine and 4-amino-4-deoxy-L-arabinose, in addition to the use of efflux pumps, the formation of capsules and overexpression of the outer membrane protein OprH, which are all effectively regulated at the molecular level. The increased understanding of these mechanisms is extremely vital and timely to facilitate studies of antimicrobial peptides and find new potential drugs targeting clinically relevant Gram-negative bacteria.

Wolbachia are a group of intracellular inherited bacteria that infect a wide range of arthropods. They are associated with a number of different reproductive phenotypes in their hosts, such as cytoplasmic incompatibility, parthenogenesis and feminization. While it is known that the bacterial strains responsible for these different host phenotypes form a single clade within the alpha-Proteobacteria, until now it has not been possible to resolve the evolutionary relationships between different Wolbachia strains. To address this issue we have cloned and sequenced a gene encoding a surface protein of Wolbachia (wsp) from a representative sample of 28 Wolbachia strains. The sequences from this gene were highly variable and could be used to resolve the phylogenetic relationships of different Wolbachia strains. Based on the sequence of the wsp gene from different Wolbachia isolates we propose that the Wolbachia pipientis clade be initially divided into 12 groups. As more sequence information becomes available we expect the number of such groups to increase. In addition, we present a method of Wolbachia classification based on the use of group-specific wsp polymerase chain reaction (PGR) primers which will allow Wolbachia isolates to be typed without the need to clone and sequence individual Wolbachia genes. This system should facilitate future studies investigating the distribution and biology of Wolbachia strains from large samples of different host species.

We examine the power of different exact tests of differentiation for diploid populations. Since there is not necessarily random mating within populations, the appropriate hypothesis to construct exact tests is that of independent sampling of genotypes. There are two categories of tests, FST-estimator tests and goodness of fit tests. In this latter category, we distinguish "allelic statistics", which account for the nature of alleles within genotypes, from "genotypic statistics" that do not. We show that the power of FST-estimator tests and of allelic goodness of fit tests are similar when sampling is balanced, and higher than the power of genotypic goodness of fit tests. When sampling is unbalanced, the most powerful tests are shown to belong to the allelic goodness of fit group.

This review proposes ten tentative answers to frequently asked questions about dispersal evolution. I examine methodological issues, model assumptions and predictions, and their relation to empirical data. Study of dispersal evolution points to the many ecological and genetic feedbacks affecting the evolution of this complex trait, which has contributed to our better understanding of life-history evolution in spatially structured populations. Several lines of research are suggested to ameliorate the exchanges between theoretical and empirical studies of dispersal evolution.

Over the last three decades, mitochondrial DNA has been the most popular marker of molecular diversity, for a combination of technical ease-of-use considerations, and supposed biological and evolutionary properties of clonality, near-neutrality and clock-like nature of its substitution rate. Reviewing recent literature on the subject, we argue that mitochondrial DNA is not always clonal, far from neutrally evolving and certainly not clock-like, questioning its relevance as a witness of recent species and population history. We critically evaluate the usage of mitochondrial DNA for species delineation and identification. Finally, we note the great potential of accumulating mtDNA data for evolutionary and functional analysis of the mitochondrial genome.

Pyrethroid-impregnated bednets are playing an increasing role for combating malaria, especially in stable malaria areas. More than 90% of the current annual malaria incidence (c. 500 million clinical cases with up to 2 million deaths) is in Africa where the major vector is Anopheles gambiae s.s. As pyrethroid resistance has been reported in this mosquito, reliable and simple techniques are urgently needed to characterize and monitor this resistance in the field. In insects, an important mechanism of pyrethroid resistance is due to a modification of the voltage-gated sodium channel protein recently shown to be associated with mutations of the para-type sodium channel gene. We demonstrate here that one of these mutations is present in certain strains of pyrethroid resistant A. gambiae s.s. and describe a PCR-based diagnostic test allowing its detection in the genome of single mosquitoes. Using this test, we found this mutation in six out of seven field samples from West Africa, its frequency being closely correlated with survival to pyrethroid exposure. This diagnostic test should bring major improvement for field monitoring of pyrethroid resistance, within the framework of malaria control programmes.

Many dynamical systems, including lakes, organisms, ocean circulation patterns, or financial markets, are now thought to have tipping points where critical transitions to a contrasting state can happen. Because critical transitions can occur unexpectedly and are difficult to manage, there is a need for methods that can be used to identify when a critical transition is approaching. Recent theory shows that we can identify the proximity of a system to a critical transition using a variety of so-called 'early warning signals', and successful empirical examples suggest a potential for practical applicability. However, while the range of proposed methods for predicting critical transitions is rapidly expanding, opinions on their practical use differ widely, and there is no comparative study that tests the limitations of the different methods to identify approaching critical transitions using time-series data. Here, we summarize a range of currently available early warning methods and apply them to two simulated time series that are typical of systems undergoing a critical transition. In addition to a methodological guide, our work offers a practical toolbox that may be used in a wide range of fields to help detect early warning signals of critical transitions in time series data.

The prevention of deforestation and promotion of afforestation have often been cited as strategies to slow global warming. Deforestation releases CO(2) to the atmosphere, which exerts a warming influence on Earth's climate. However, biophysical effects of deforestation, which include changes in land surface albedo, evapotranspiration, and cloud cover also affect climate. Here we present results from several large-scale deforestation experiments performed with a three-dimensional coupled global carbon-cycle and climate model. These simulations were performed by using a fully three-dimensional model representing physical and biogeochemical interactions among land, atmosphere, and ocean. We find that global-scale deforestation has a net cooling influence on Earth's climate, because the warming carbon-cycle effects of deforestation are overwhelmed by the net cooling associated with changes in albedo and evapotranspiration. Latitude-specific deforestation experiments indicate that afforestation projects in the tropics would be clearly beneficial in mitigating global-scale warming, but would be counterproductive if implemented at high latitudes and would offer only marginal benefits in temperate regions. Although these results question the efficacy of mid- and high-latitude afforestation projects for climate mitigation, forests remain environmentally valuable resources for many reasons unrelated to climate.