Laboratoire de Biométrie et Biologie Evolutive

UNLABELLED: The package adegenet for the R software is dedicated to the multivariate analysis of genetic markers. It extends the ade4 package of multivariate methods by implementing formal classes and functions to manipulate and analyse genetic markers. Data can be imported from common population genetics software and exported to other software and R packages. adegenet also implements standard population genetics tools along with more original approaches for spatial genetics and hybridization. AVAILABILITY: Stable version is available from CRAN: http://cran.r-project.org/mirrors.html. Development version is available from adegenet website: http://adegenet.r-forge.r-project.org/. Both versions can be installed directly from R. adegenet is distributed under the GNU General Public Licence (v.2).

We present SeaView version 4, a multiplatform program designed to facilitate multiple alignment and phylogenetic tree building from molecular sequence data through the use of a graphical user interface. SeaView version 4 combines all the functions of the widely used programs SeaView (in its previous versions) and Phylo_win, and expands them by adding network access to sequence databases, alignment with arbitrary algorithm, maximum-likelihood tree building with PhyML, and display, printing, and copy-to-clipboard of rooted or unrooted, binary or multifurcating phylogenetic trees. In relation to the wide present offer of tools and algorithms for phylogenetic analyses, SeaView is especially useful for teaching and for occasional users of such software. SeaView is freely available at http://pbil.univ-lyon1.fr/software/seaview.

BACKGROUND: The dramatic progress in sequencing technologies offers unprecedented prospects for deciphering the organization of natural populations in space and time. However, the size of the datasets generated also poses some daunting challenges. In particular, Bayesian clustering algorithms based on pre-defined population genetics models such as the STRUCTURE or BAPS software may not be able to cope with this unprecedented amount of data. Thus, there is a need for less computer-intensive approaches. Multivariate analyses seem particularly appealing as they are specifically devoted to extracting information from large datasets. Unfortunately, currently available multivariate methods still lack some essential features needed to study the genetic structure of natural populations. RESULTS: We introduce the Discriminant Analysis of Principal Components (DAPC), a multivariate method designed to identify and describe clusters of genetically related individuals. When group priors are lacking, DAPC uses sequential K-means and model selection to infer genetic clusters. Our approach allows extracting rich information from genetic data, providing assignment of individuals to groups, a visual assessment of between-population differentiation, and contribution of individual alleles to population structuring. We evaluate the performance of our method using simulated data, which were also analyzed using STRUCTURE as a benchmark. Additionally, we illustrate the method by analyzing microsatellite polymorphism in worldwide human populations and hemagglutinin gene sequence variation in seasonal influenza. CONCLUSIONS: Analysis of simulated data revealed that our approach performs generally better than STRUCTURE at characterizing population subdivision. The tools implemented in DAPC for the identification of clusters and graphical representation of between-group structures allow to unravel complex population structures. Our approach is also faster than Bayesian clustering algorithms by several orders of magnitude, and may be applicable to a wider range of datasets.

Traditionally, measurement of the acute toxicity of pesticides to beneficial arthropods has relied largely on the determination of an acute median lethal dose or concentration. However, the estimated lethal dose during acute toxicity tests may only be a partial measure of the deleterious effects. In addition to direct mortality induced by pesticides, their sublethal effects on arthropod physiology and behavior must be considered for a complete analysis of their impact. An increasing number of studies and methods related to the identification and characterization of these effects have been published in the past 15 years. Review of sublethal effects reported in published literature, taking into account recent data, has revealed new insights into the sublethal effects of pesticides including effects on learning performance, behavior, and neurophysiology. We characterize the different types of sublethal effects on beneficial arthropods, focusing mainly on honey bees and natural enemies, and we describe the methods used in these studies. Finally, we discuss the potential for developing experimental approaches that take into account these sublethal effects in integrated pest management and the possibility of integrating their evaluation in pesticide registration procedures.

SEAVIEW and PHYLO_WIN are two graphic tools for X Windows-Unix computers dedicated to sequence alignment and molecular phylogenetics. SEAVIEW is a sequence alignment editor allowing manual or automatic alignment through an interface with CLUSTALW program. Alignment of large sequences with extensive length differences is made easier by a dot-plot-based routine. The PHYLO_WIN program allows phylogenetic tree building according to most usual methods (neighbor joining with numerous distance estimates, maximum parsimony, maximum likelihood), and a bootstrap analysis with any of them. Reconstructed trees can be drawn, edited, printed, stored, evaluated according to numerous criteria. Taxonomic species groups and sets of conserved regions can be defined by mouse and stored into sequence files, thus avoiding multiple data files. Both tools are entirely mouse driven. On-line help makes them easy to use. They are freely available by anonymous ftp at biom3.univ-lyon1.fr/pub/ mol_phylogeny or http:@acnuc.univ-lyon1.fr/, or by e-mail to galtier@biomserv.univ-lyon1.fr.

The package fitdistrplus provides functions for fitting univariate distributions to different types of data (continuous censored or non-censored data and discrete data) and allowing different estimation methods (maximum likelihood, moment matching, quantile matching and maximum goodness-of-fit estimation). Outputs of fitdist and fitdistcens functions are S3 objects, for which specific methods are provided, including summary, plot and quantile. This package also provides various functions to compare the fit of several distributions to the same data set and can handle to bootstrap parameter estimates. Detailed examples are given in food risk assessment, ecotoxicology and insurance contexts.

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.

BACKGROUND: In patients who have chronic heart failure with reduced left ventricular ejection fraction, severe secondary mitral-valve regurgitation is associated with a poor prognosis. Whether percutaneous mitral-valve repair improves clinical outcomes in this patient population is unknown. METHODS: or a regurgitant volume of >30 ml per beat), a left ventricular ejection fraction between 15 and 40%, and symptomatic heart failure, in a 1:1 ratio, to undergo percutaneous mitral-valve repair in addition to receiving medical therapy (intervention group; 152 patients) or to receive medical therapy alone (control group; 152 patients). The primary efficacy outcome was a composite of death from any cause or unplanned hospitalization for heart failure at 12 months. RESULTS: At 12 months, the rate of the primary outcome was 54.6% (83 of 152 patients) in the intervention group and 51.3% (78 of 152 patients) in the control group (odds ratio, 1.16; 95% confidence interval [CI], 0.73 to 1.84; P=0.53). The rate of death from any cause was 24.3% (37 of 152 patients) in the intervention group and 22.4% (34 of 152 patients) in the control group (hazard ratio, 1.11; 95% CI, 0.69 to 1.77). The rate of unplanned hospitalization for heart failure was 48.7% (74 of 152 patients) in the intervention group and 47.4% (72 of 152 patients) in the control group (hazard ratio, 1.13; 95% CI, 0.81 to 1.56). CONCLUSIONS: Among patients with severe secondary mitral regurgitation, the rate of death or unplanned hospitalization for heart failure at 1 year did not differ significantly between patients who underwent percutaneous mitral-valve repair in addition to receiving medical therapy and those who received medical therapy alone. (Funded by the French Ministry of Health and Research National Program and Abbott Vascular; MITRA-FR ClinicalTrials.gov number, NCT01920698 .).

▪ Abstract In large-herbivore populations, environmental variation and density dependence co-occur and have similar effects on various fitness components. Our review aims to quantify the temporal variability of fitness components and examine how that variability affects changes in population growth rates. Regardless of the source of variation, adult female survival shows little year-to-year variation [coefficient of variation (CV <10%)], fecundity of prime-aged females and yearling survival rates show moderate year-to-year variation (CV <20%), and juvenile survival and fecundity of young females show strong variation (CV >30%). Old females show senescence in both survival and reproduction. These patterns of variation are independent of differences in body mass, taxonomic group, and ecological conditions. Differences in levels of maternal care may fine-tune the temporal variation of early survival. The immature stage, despite a low relative impact on population growth rate compared with the adult stage, may be the critical component of population dynamics of large herbivores. Observed differences in temporal variation may be more important than estimated relative sensitivity or elasticity in determining the relative demographic impact of various fitness components.

We propose a multivariate approach to the study of geographic species distribution which does not require absence data. Building on Hutchinson's concept of the ecological niche, this factor analysis compares, in the multidimensional space of ecological variables, the distribution of the localities where the focal species was observed to a reference set describing the whole study area. The first factor extracted maximizes the marginality of the focal species, defined as the ecological distance between the species optimum and the mean habitat within the reference area. The other factors maximize the specialization of this focal species, defined as the ratio of the ecological variance in mean habitat to that observed for the focal species. Eigenvectors and eigenvalues are readily interpreted and can be used to build habitat-suitability maps. This approach is recommended in situations where absence data are not available (many data banks), unreliable (most cryptic or rare species), or meaningless (invaders). We provide an illustration and validation of the method for the alpine ibex, a species reintroduced in Switzerland which presumably has not yet recolonized its entire range.

The InterPro database (http://www.ebi.ac.uk/interpro/) is a freely available resource that can be used to classify sequences into protein families and to predict the presence of important domains and sites. Central to the InterPro database are predictive models, known as signatures, from a range of different protein family databases that have different biological focuses and use different methodological approaches to classify protein families and domains. InterPro integrates these signatures, capitalizing on the respective strengths of the individual databases, to produce a powerful protein classification resource. Here, we report on the status of InterPro as it enters its 15th year of operation, and give an overview of new developments with the database and its associated Web interfaces and software. In particular, the new domain architecture search tool is described and the process of mapping of Gene Ontology terms to InterPro is outlined. We also discuss the challenges faced by the resource given the explosive growth in sequence data in recent years. InterPro (version 48.0) contains 36,766 member database signatures integrated into 26,238 InterPro entries, an increase of over 3993 entries (5081 signatures), since 2012.

The nucleic acid sequence bank now contains over 600 protein coding genes of which 107 are from prokaryotic organisms. Codon frequencies in each new prokaryotic gene are given. Analysis of genetic code usage in the 83 sequenced genes of the Escherichia coli genome (chromosome, transposons and plasmids) is presented, taking into account new data on gene expressivity and regulation as well as iso-tRNA specificity and cellular concentration. The codon composition of each gene is summarized using two indexes: one is based on the differential usage of iso-tRNA species during gene translation, the other on choice between Cytosine and Uracil for third base. A strong relationship between codon composition and mRNA expressivity is confirmed, even for genes transcribed in the same operon. The influence of codon use of peptide elongation rate and protein yield is discussed. Finally, the evolutionary aspect of codon selection in mRNA sequences is studied.

BACKGROUND: Under current guidelines, intravenous thrombolysis is used to treat acute stroke only if it can be ascertained that the time since the onset of symptoms was less than 4.5 hours. We sought to determine whether patients with stroke with an unknown time of onset and features suggesting recent cerebral infarction on magnetic resonance imaging (MRI) would benefit from thrombolysis with the use of intravenous alteplase. METHODS: In a multicenter trial, we randomly assigned patients who had an unknown time of onset of stroke to receive either intravenous alteplase or placebo. All the patients had an ischemic lesion that was visible on MRI diffusion-weighted imaging but no parenchymal hyperintensity on fluid-attenuated inversion recovery (FLAIR), which indicated that the stroke had occurred approximately within the previous 4.5 hours. We excluded patients for whom thrombectomy was planned. The primary end point was favorable outcome, as defined by a score of 0 or 1 on the modified Rankin scale of neurologic disability (which ranges from 0 [no symptoms] to 6 [death]) at 90 days. A secondary outcome was the likelihood that alteplase would lead to lower ordinal scores on the modified Rankin scale than would placebo (shift analysis). RESULTS: The trial was stopped early owing to cessation of funding after the enrollment of 503 of an anticipated 800 patients. Of these patients, 254 were randomly assigned to receive alteplase and 249 to receive placebo. A favorable outcome at 90 days was reported in 131 of 246 patients (53.3%) in the alteplase group and in 102 of 244 patients (41.8%) in the placebo group (adjusted odds ratio, 1.61; 95% confidence interval [CI], 1.09 to 2.36; P=0.02). The median score on the modified Rankin scale at 90 days was 1 in the alteplase group and 2 in the placebo group (adjusted common odds ratio, 1.62; 95% CI, 1.17 to 2.23; P=0.003). There were 10 deaths (4.1%) in the alteplase group and 3 (1.2%) in the placebo group (odds ratio, 3.38; 95% CI, 0.92 to 12.52; P=0.07). The rate of symptomatic intracranial hemorrhage was 2.0% in the alteplase group and 0.4% in the placebo group (odds ratio, 4.95; 95% CI, 0.57 to 42.87; P=0.15). CONCLUSIONS: In patients with acute stroke with an unknown time of onset, intravenous alteplase guided by a mismatch between diffusion-weighted imaging and FLAIR in the region of ischemia resulted in a significantly better functional outcome and numerically more intracranial hemorrhages than placebo at 90 days. (Funded by the European Union Seventh Framework Program; WAKE-UP ClinicalTrials.gov number, NCT01525290; and EudraCT number, 2011-005906-32 .).

The expression of most Staphylococcus aureus virulence factors is controlled by the agr locus, which encodes a two-component signaling pathway whose activating ligand is an agr-encoded autoinducing peptide (AIP). A polymorphism in the amino acid sequence of the AIP and of its corresponding receptor divides S. aureus strains into four major groups. Within a given group, each strain produces a peptide that can activate the agr response in the other member strains, whereas the AIPs belonging to different groups are usually mutually inhibitory. We investigated a possible relationship between agr groups and human S. aureus disease by studying 198 S. aureus strains isolated from 14 asymptomatic carriers, 66 patients with suppurative infection, and 114 patients with acute toxemia. The agr group and the distribution of 24 toxin genes were analyzed by PCR, and the genetic background was determined by means of amplified fragment length polymorphism (AFLP) analysis. The isolates were relatively evenly distributed among the four agrgroups, with 61 strains belonging to agr group I, 49 belonging to group II, 43 belonging to group III, and 45 belonging to group IV. Principal coordinate analysis performed on the AFLP distance matrix divided the 198 strains into three main phylogenetic groups, AF1 corresponding to strains of agr group IV, AF2 corresponding to strains of agr groups I and II, and AF3 corresponding to strains of agr group III. This indicated that the agr type was linked to the genetic background. A relationship between genetic background, agr group, and disease type was observed for several toxin-mediated diseases: for instance, agr group IV strains were associated with generalized exfoliative syndromes, and phylogenetic group AF1 strains with bullous impetigo. Among the suppurative infections, endocarditis strains mainly belonged to phylogenetic group AF2 and agr groups I and II. While these results do not show a direct role of the agr type in the type of human disease caused by S. aureus, the agr group may reflect an ancient evolutionary division of S. aureus in terms of this species' fundamental biology.

Summary 1. Phylogenetic signal is the tendency of related species to resemble each other more than species drawn at random from the same tree. This pattern is of considerable interest in a range of ecological and evolutionary research areas, and various indices have been proposed for quantifying it. Unfortunately, these indices often lead to contrasting results, and guidelines for choosing the most appropriate index are lacking. 2. Here, we compare the performance of four commonly used indices using simulated data. Data were generated with numerical simulations of trait evolution along phylogenetic trees under a variety of evolutionary models. We investigated the sensitivity of the approaches to the size of phylogenies, the resolution of tree structure and the availability of branch length information, examining both the response of the selected indices and the power of the associated statistical tests. 3. We found that under a Brownian motion (BM) model of trait evolution, Abouheif’s C mean and Pagel’s λ performed well and substantially better than Moran’s I and Blomberg’s K . Pagel’s λ provided a reliable effect size measure and performed better for discriminating between more complex models of trait evolution, but was computationally more demanding than Abouheif’s C mean . Blomberg’s K was most suitable to capture the effects of changing evolutionary rates in simulation experiments. 4. Interestingly, sample size influenced not only the uncertainty but also the expected values of most indices, while polytomies and missing branch length information had only negligible impacts. 5. We propose guidelines for choosing among indices, depending on (a) their sensitivity to true underlying patterns of phylogenetic signal, (b) whether a test or a quantitative measure is required and (c) their sensitivities to different topologies of phylogenies. 6. These guidelines aim to better assess phylogenetic signal and distinguish it from random trait distributions. They were developed under the assumption of BM, and additional simulations with more complex trait evolution models show that they are to a certain degree generalizable. They are particularly useful in comparative analyses, when requiring a proxy for niche similarity, and in conservation studies that explore phylogenetic loss associated with extinction risks of specific clades.

CONTEXT: Providing accurate estimates of cancer risks is a major challenge in the clinical management of Lynch syndrome. OBJECTIVE: To estimate the age-specific cumulative risks of developing various tumors using a large series of families with mutations of the MLH1, MSH2, and MSH6 genes. DESIGN, SETTING, AND PARTICIPANTS: Families with Lynch syndrome enrolled between January 1, 2006, and December 31, 2009, from 40 French cancer genetics clinics participating in the ERISCAM (Estimation des Risques de Cancer chez les porteurs de mutation des gènes MMR) study; 537 families with segregating mutated genes (248 with MLH1; 256 with MSH2; and 33 with MSH6) were analyzed. MAIN OUTCOME MEASURE: Age-specific cumulative cancer risks estimated using the genotype restricted likelihood (GRL) method accounting for ascertainment bias. RESULTS: Significant differences in estimated cumulative cancer risk were found between the 3 mutated genes (P = .01). The estimated cumulative risks of colorectal cancer by age 70 years were 41% (95% confidence intervals [CI], 25%-70%) for MLH1 mutation carriers, 48% (95% CI, 30%-77%) for MSH2, and 12% (95% CI, 8%-22%) for MSH6. For endometrial cancer, corresponding risks were 54% (95% CI, 20%-80%), 21% (95% CI, 8%-77%), and 16% (95% CI, 8%-32%). For ovarian cancer, they were 20% (95% CI, 1%-65%), 24% (95% CI, 3%-52%), and 1% (95% CI, 0%-3%). The estimated cumulative risks by age 40 years did not exceed 2% (95% CI, 0%-7%) for endometrial cancer nor 1% (95% CI, 0%-3%) for ovarian cancer, irrespective of the gene. The estimated lifetime risks for other tumor types did not exceed 3% with any of the gene mutations. CONCLUSIONS: MSH6 mutations are associated with markedly lower cancer risks than MLH1 or MSH2 mutations. Lifetime ovarian and endometrial cancer risks associated with MLH1 or MSH2 mutations were high but do not increase appreciably until after the age of 40 years.

Microsporidia are obligate intracellular parasites infesting many animal groups. Lacking mitochondria and peroxysomes, these unicellular eukaryotes were first considered a deeply branching protist lineage that diverged before the endosymbiotic event that led to mitochondria. The discovery of a gene for a mitochondrial-type chaperone combined with molecular phylogenetic data later implied that microsporidia are atypical fungi that lost mitochondria during evolution. Here we report the DNA sequences of the 11 chromosomes of the approximately 2.9-megabase (Mb) genome of Encephalitozoon cuniculi (1,997 potential protein-coding genes). Genome compaction is reflected by reduced intergenic spacers and by the shortness of most putative proteins relative to their eukaryote orthologues. The strong host dependence is illustrated by the lack of genes for some biosynthetic pathways and for the tricarboxylic acid cycle. Phylogenetic analysis lends substantial credit to the fungal affiliation of microsporidia. Because the E. cuniculi genome contains genes related to some mitochondrial functions (for example, Fe-S cluster assembly), we hypothesize that microsporidia have retained a mitochondrion-derived organelle.

SUMMARY We present an unconventional procedure (fuzzy coding) to structure biological and environmental information, which uses positive scores to describe the affinity of a species for different modalities (i.e. categories) of a given variable. Fuzzy coding is essential for the synthesis of long‐term ecological data because it enables analysis of diverse kinds of biological information derived from a variety of sources (e.g. samples, literature). A fuzzy coded table can be processed by correspondence analysis. An example using aquatic beetles illustrates the properties of such a fuzzy correspondence analysis. Fuzzy coded tables were used in all articles of this issue to examine relationships between spatial‐temporal habitat variability and species traits, which were obtained from a long‐term study of the Upper Rhône River, France. Fuzzy correspondence analysis can be programmed with the equations given in this paper or can be performed using ADE (Environmental Data Analysis) software that has been adapted to analyse such long‐term ecological data. On Macintosh Apple TM computers, ADE performs simple linear ordination, more recently developed methods (e.g. principal component analysis with respect to instrumental variables, canonical correspondence analysis, co‐inertia analysis, local and spatial analyses), and provides a graphical display of results of these and other types of analysis (e.g. biplot, mapping, modelling curves). ADE consists of a program library that exploits the potential of the HyperCard TM interface. ADE in an open system, which offers the user a variety of facilities to create a specific sequence of programs. The mathematical background of ADE is supported by the algebraic model known as ‘duality diagram’.

The nucleic acid sequence bank now contains 161 mRNAs, 43 new genes are added. One sequence, that of B. mori fibroin, is dropped due to uncertainty on the starting point for translation. Frequencies of all codons are given for each gene added and for each genome type in the total bank. A new series of correspondence analyses on codon use is presented, substantiating the genome hypothesis. Internal regulation of mRNA expression by different third base choices between quartet and duet codons is proposed for bacterial genes.

Programs for Bayesian inference of phylogeny currently implement a unique and fixed suite of models. Consequently, users of these software packages are simultaneously forced to use a number of programs for a given study, while also lacking the freedom to explore models that have not been implemented by the developers of those programs. We developed a new open-source software package, RevBayes, to address these problems. RevBayes is entirely based on probabilistic graphical models, a powerful generic framework for specifying and analyzing statistical models. Phylogenetic-graphical models can be specified interactively in RevBayes, piece by piece, using a new succinct and intuitive language called Rev. Rev is similar to the R language and the BUGS model-specification language, and should be easy to learn for most users. The strength of RevBayes is the simplicity with which one can design, specify, and implement new and complex models. Fortunately, this tremendous flexibility does not come at the cost of slower computation; as we demonstrate, RevBayes outperforms competing software for several standard analyses. Compared with other programs, RevBayes has fewer black-box elements. Users need to explicitly specify each part of the model and analysis. Although this explicitness may initially be unfamiliar, we are convinced that this transparency will improve understanding of phylogenetic models in our field. Moreover, it will motivate the search for improvements to existing methods by brazenly exposing the model choices that we make to critical scrutiny. RevBayes is freely available at http://www.RevBayes.com [Bayesian inference; Graphical models; MCMC; statistical phylogenetics.].