Molecular dating, evolutionary rates, and the age of the grasses.

Many questions in evolutionary biology require an estimate of divergence times but, for groups with a sparse fossil record, such estimates rely heavily on molecular dating methods. The accuracy of these methods depends on both an adequate underlying model and the appropriate implementation of fossil evidence as calibration points. We explore the effect of these in Poaceae (grasses), a diverse plant lineage with a very limited fossil record, focusing particularly on dating the early divergences in the group. We show that molecular dating based on a data set of plastid markers is strongly dependent on the model assumptions. In particular, an acceleration of evolutionary rates at the base of Poaceae followed by a deceleration in the descendants strongly biases methods that assume an autocorrelation of rates. This problem can be circumvented by using markers that have lower rate variation, and we show that phylogenetic markers extracted from complete nuclear genomes can be a useful complement to the more commonly used plastid markers. However, estimates of divergence times remain strongly affected by different implementations of fossil calibration points. Analyses calibrated with only macrofossils lead to estimates for the age of core Poaceae ∼51-55 Ma, but the inclusion of microfossil evidence pushes this age to 74-82 Ma and leads to lower estimated evolutionary rates in grasses. These results emphasize the importance of considering markers from multiple genomes and alternative fossil placements when addressing evolutionary issues that depend on ages estimated for important groups.

Amphibian albumins as members of the albumin, alpha-fetoprotein, vitamin D-binding protein multigene family.

The Xenopus laevis 68-kd and 74-kd albumin amino acid sequences are examined with respect to their relationship to the other known members of the albumin/alpha-fetoprotein/vitamin D-binding protein gene family. Each of the three members of this family presents a unique pattern of conserved regions indicating a differential selective pressure related to specific functional characteristics. Furthermore, an evolutionary tree of these genes was deduced from the divergence times calculated from direct nucleotide sequence comparisons of individual gene pairs. These calculations indicate that the vitamin D-binding protein/albumin separation occurred 560-600 million years (Myr) ago and the albumin/alpha-fetoprotein divergence 280 Myr ago. This observation leads to the hypothesis according to which the albumin/alpha-fetoprotein gene duplication occurred shortly after the amphibian/reptile separation. Consequently, and unlike mammals, amphibians and fishes should lack an alpha-fetoprotein in their serum at larval stages, which is consistent with a recent analysis of serum proteins in Xenopus laevis larvae. This hypothesis now will have to be tested further in additional lower vertebrates.

Unexpectedly many extinct hominins.

Recent studies indicate that Neanderthal and Denisova hominins may have been separate species, while debate continues on the status of Homo floresiensis. The decade-long debate between "splitters," who recognize over 20 hominin species, and "lumpers," who maintain that all these fossils belong to just a few lineages, illustrates that we do not know how many extinct hominin species to expect. Here, we present probability distributions for the number of speciation events and the number of contemporary species along a branch of a phylogeny. With estimates of hominin speciation and extincton rates, we then show that the expected total number of extinct hominin species is 8, but may be as high as 27. We also show that it is highly unlikely that three very recent species disappeared due to natural, background extinction. This may indicate that human-like remains are too easily considered distinct species. Otherwise, the evidence suggesting that Neanderthal and the Denisova hominin represent distinct species implies a recent wave of extinctions, ostensibly driven by the only survivor, H. sapiens.

Species-specific population structure in rock-specialized sympatric cichlid species in Lake Tanganyika, East Africa.

Species richness and geographical phenotypic variation in East African lacustrine cichlids are often correlated with ecological specializations and limited dispersal. This study compares mitochondrial and microsatellite genetic diversity and structure among three sympatric rock-dwelling cichlids of Lake Tanganyika, Eretmodus cyanostictus, Tropheus moorii, and Ophthalmotilapia ventralis. The species represent three endemic, phylogenetically distinct tribes (Eretmodini, Tropheini, and Ectodini), and display divergent ecomorphological and behavioral specialization. Sample locations span both continuous, rocky shoreline and a potential dispersal barrier in the form of a muddy bay. High genetic diversity and population differentiation were detected in T. moorii and E. cyanostictus, whereas much lower variation and structure were found in O. ventralis. In particular, while a 7-km-wide muddy bay curtails dispersal in all three species to a similar extent, gene flow along mostly continuous habitat appeared to be controlled by distance in E. cyanostictus, further restricted by site philopatry and/or minor habitat discontinuities in T. moorii, and unrestrained in O. ventralis. In contrast to the general pattern of high gene flow along continuous shorelines in rock-dwelling cichlids of Lake Malawi, our study identifies differences in population structure among stenotopic Lake Tanganyika species. The amount of genetic differentiation among populations was not related to the degree of geographical variation of body color, especially since more phenotypic variation is observed in O. ventralis than in the genetically highly structured E. cyanostictus.

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.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins.

We describe the odorant binding proteins (OBPs) of the red imported fire ant, Solenopsis invicta, obtained from analyses of an EST library and separate 454 sequencing runs of two normalized cDNA libraries. We identified a total of 18 putative functional OBPs in this ant. A third of the fire ant OBPs are orthologs to honey bee OBPs. Another third of the OBPs belong to a lineage-specific expansion, which is a common feature of insect OBP evolution. Like other OBPs, the different fire ant OBPs share little sequence similarity (∼ 20%), rendering evolutionary analyses difficult. We discuss the resulting problems with sequence alignment, phylogenetic analysis, and tests of selection. As previously suggested, our results underscore the importance for careful exploration of the sensitivity to the effects of alignment methods for data comprising widely divergent sequences.

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.

The phylogeography of an alpine leaf beetle: divergence within Oreina elongata spans several ice ages.

The genetic landscape of the European flora and fauna was shaped by the ebb and flow of populations with the shifting ice during Quaternary climate cycles. While this has been well demonstrated for lowland species, less is known about high altitude taxa. Here we analyze the phylogeography of the leaf beetle Oreina elongata from 20 populations across the Alps and Apennines. Three mitochondrial and one nuclear region were sequenced in 64 individuals. Within an mtDNA phylogeny, three of seven subspecies are monophyletic. The species is chemically defended and aposematic, with green and blue forms showing geographic variation and unexpected within-population polymorphism. These warning colors show pronounced east-west geographical structure in distribution, but the phylogeography suggests repeated origin and loss. Basal clades come from the central Alps. Ancestors of other clades probably survived across northern Italy and the northern Adriatic, before separation of eastern, southern and western populations and rapid spread through the western Alps. After reviewing calibrated gene-specific substitution rates in the literature, we use partitioned Bayesian coalescent analysis to date our phylogeography. The major clades diverged long before the last glacial maximum, suggesting that O. elongata persisted many glacial cycles within or at the edges of the Alps and Apennines. When analyzing additional barcoding pairwise distances, we find strong evidence to consider O. elongata as a species complex rather than a single species.

Evidence for adaptive radiation from a phylogenetic study of plant defenses.

One signature of adaptive radiation is a high level of trait change early during the diversification process and a plateau toward the end of the radiation. Although the study of the tempo of evolution has historically been the domain of paleontologists, recently developed phylogenetic tools allow for the rigorous examination of trait evolution in a tremendous diversity of organisms. Enemy-driven adaptive radiation was a key prediction of Ehrlich and Raven's coevolutionary hypothesis [Ehrlich PR, Raven PH (1964) Evolution 18:586-608], yet has remained largely untested. Here we examine patterns of trait evolution in 51 North American milkweed species (Asclepias), using maximum likelihood methods. We study 7 traits of the milkweeds, ranging from seed size and foliar physiological traits to defense traits (cardenolides, latex, and trichomes) previously shown to impact herbivores, including the monarch butterfly. We compare the fit of simple random-walk models of trait evolution to models that incorporate stabilizing selection (Ornstein-Ulenbeck process), as well as time-varying rates of trait evolution. Early bursts of trait evolution were implicated for 2 traits, while stabilizing selection was implicated for several others. We further modeled the relationship between trait change and species diversification while allowing rates of trait evolution to vary during the radiation. Species-rich lineages underwent a proportionately greater decline in latex and cardenolides relative to species-poor lineages, and the rate of trait change was most rapid early in the radiation. An interpretation of this result is that reduced investment in defensive traits accelerated diversification, and disproportionately so, early in the adaptive radiation of milkweeds.

Toxicity of triclosan, penconazole and metalaxyl on Caulobacter crescentus and a freshwater microbial community as assessed by flow cytometry.

Biocides are widely used for domestic hygiene, agricultural and industrial applications. Their widespread use has resulted in their introduction into the environment and raised concerns about potential deleterious effects on aquatic ecosystems. In this study, the toxicity of the biocides triclosan, penconazole and metalaxyl were evaluated with the freshwater bacterium Caulobacter crescentus and with a freshwater microbial community using a combination of single- and double-stain flow cytometric assays. Growth of C.  crescentus and the freshwater community were repressed by triclosan but not by penconazole or metalaxyl at concentrations up to 250 μM. The repressive effect of triclosan was dependent on culture conditions. Caulobacter crescentus was more sensitive to triclosan when grown with high glucose at high cell density than when grown directly in sterilized lake water at low cell density. This suggests that the use of conventional growth conditions may overestimate biocide toxicity. Additional experiments showed that the freshwater community was more sensitive to triclosan than C.  crescentus, with 10 nM of triclosan being sufficient to repress growth and change the phylogenetic composition of the community. These results demonstrate that isolate-based assays may underestimate biocide toxicity and highlight the importance of assessing toxicity directly on natural microbial communities. Because 10 nM of triclosan is within the range of concentrations observed in freshwater systems, these results also raise concerns about the risk of introducing triclosan into the environment.

C4 Photosynthesis evolved in grasses via parallel adaptive genetic changes.

Phenotypic convergence is a widespread and well-recognized evolutionary phenomenon. However, the responsible molecular mechanisms remain often unknown mainly because the genes involved are not identified. A well-known example of physiological convergence is the C4 photosynthetic pathway, which evolved independently more than 45 times [1]. Here, we address the question of the molecular bases of the C4 convergent phenotypes in grasses (Poaceae) by reconstructing the evolutionary history of genes encoding a C4 key enzyme, the phosphoenolpyruvate carboxylase (PEPC). PEPC genes belong to a multigene family encoding distinct isoforms of which only one is involved in C4 photosynthesis [2]. By using phylogenetic analyses, we showed that grass C4 PEPCs appeared at least eight times independently from the same non-C4 PEPC. Twenty-one amino acids evolved under positive selection and converged to similar or identical amino acids in most of the grass C4 PEPC lineages. This is the first record of such a high level of molecular convergent evolution, illustrating the repeatability of evolution. These amino acids were responsible for a strong phylogenetic bias grouping all C4 PEPCs together. The C4-specific amino acids detected must be essential for C4 PEPC enzymatic characteristics, and their identification opens new avenues for the engineering of the C4 pathway in crops.

Molecular evidence for ancient asexuality in timema stick insects.

Asexuality is rare in animals in spite of its apparent advantage relative to sexual reproduction, indicating that it must be associated with profound costs [1-9]. One expectation is that reproductive advantages gained by new asexual lineages will be quickly eroded over time [3, 5-7]. Ancient asexual taxa that have evolved and adapted without sex would be "scandalous" exceptions to this rule, but it is often difficult to exclude the possibility that putative asexuals deploy some form of "cryptic" sex, or have abandoned sex more recently than estimated from divergence times to sexual relatives [10]. Here we provide evidence, from high intraspecific divergence of mitochondrial sequence and nuclear allele divergence patterns, that several independently derived Timema stick-insect lineages have persisted without recombination for more than a million generations. Nuclear alleles in the asexual lineages displayed significantly higher intraindividual divergences than in related sexual species. In addition, within two asexuals, nuclear allele phylogenies suggested the presence of two clades, with sequences from the same individual appearing in both clades. These data strongly support ancient asexuality in Timema and validate the genus as an exceptional opportunity to attack the question of how asexual reproduction can be maintained over long periods of evolutionary time.

Loss of egg yolk genes in mammals and the origin of lactation and placentation

Embryonic development in nonmammalian vertebrates depends entirely on nutritional reserves that are predominantly derived from vitellogenin proteins and stored in egg yolk. Mammals have evolved new resources, such as lactation and placentation, to nourish their developing and early offspring. However, the evolutionary timing and molecular events associated with this major phenotypic transition are not known. By means of sensitive comparative genomics analyses and evolutionary simulations, we here show that the three ancestral vitellogenin-encoding genes were progressively lost during mammalian evolution (until around 30-70 million years ago, Mya) in all but the egg-laying monotremes, which have retained a functional vitellogenin gene. Our analyses also provide evidence that the major milk resource genes, caseins, which have similar functional properties as vitellogenins, appeared in the common mammalian ancestor approximately 200-310 Mya. Together, our data are compatible with the hypothesis that the emergence of lactation in the common mammalian ancestor and the development of placentation in eutherian and marsupial mammals allowed for the gradual loss of yolk-dependent nourishment during mammalian evolution

Identification and characterization of pathogenic Pestalotiopsis species to pecan tree in Brazil

The objective of this work was to characterize and cluster isolates of Pestalotiopsis species and to identify those that are pathogenic to pecan, based on morphological and molecular characters. Pestalotiopsis spp. isolates were identified by sequencing the internal transcribed spacer (ITS) and β?tubulin regions. Identification methods were compared to indicate the key morphological characters for species characterization. Thirteen isolates were used for the pathogenicity tests. Morphological characterization was performed using the following variables: mycelial growth rate, sporulation, colony pigmentation, and conidial length and width. Ten pathogenic isolates were identified, three as -tubulin regions. Identification methods were compared to indicate the key morphological characters for species characterization. Thirteen isolates were used for the pathogenicity tests. Morphological characterization was performed using the following variables: mycelial growth rate, sporulation, colony pigmentation, and conidial length and width. Ten pathogenic isolates were identified, three as Pestalotiopsis clavispora and three as P. cocculi. The other isolates remained as an undefined species. The morphological characters were efficient for an initial separation of the isolates, which were grouped according to differences at species level, mainly colony diameter, which was identified as an important morphological describer. Beta-tubulin gene sequencing was less informative than the ITS region sequencing for species identification.