Efficient biased estimation of evolutionary distances when substitution rates vary across sites

This paper deals with phylogenetic inference when the variability of substitution rates across sites (VRAS) is modeled by a gamma distribution. We show that underestimating VRAS, which results in underestimates for the evolutionary distances between sequences, usually improves the topological accuracy of phylogenetic tree inference by distance-based methods, especially when the molecular clock holds. We propose a method to estimate the gamma shape parameter value which is most suited for tree topology inference, given the sequences at hand. This method is based on the pairwise evolutionary distances between sequences and allows one to reconstruct the phylogeny of a high number of taxa (>1,000). Simulation results show that the topological accuracy is highly improved when using the gamma shape parameter value given by our method, compared with the true (unknown) value which was used to generate the data. Furthermore, when VRAS is high, the topological accuracy of our distance-based method is better than that of a maximum likelihood approach. Finally, a data set of Maoricicada species sequences is analyzed, which confirms the advantage of our method.     

A phylogenetic analysis of macroevolutionary patterns in fermentative yeasts

When novel sources of ecological opportunity are available, physiological innovations can trigger adaptive radiations. This could be the case of yeasts (Saccharomycotina), in which an evolutionary novelty is represented by the capacity to exploit simple sugars from fruits (fermentation). During adaptive radiations, diversification and morphological evolution are predicted to slow‐down after early bursts of diversification. Here, we performed the first comparative phylogenetic analysis in yeasts, testing the “early burst” prediction on species diversification and also on traits of putative ecological relevance (cell‐size and fermentation versatility). We found that speciation rates are constant during the time‐range we considered (ca., 150 millions of years). Phylogenetic signal of both traits was significant (but lower for cell‐size), suggesting that lineages resemble each other in trait‐values. Disparity analysis suggested accelerated evolution (diversification in trait values above Brownian Motion expectations) in cell‐size. We also found a significant phylogenetic regression between cell‐size and fermentation versatility (R² = 0.10), which suggests correlated evolution between both traits. Overall, our results do not support the early burst prediction both in species and traits, but suggest a number of interesting evolutionary patterns, that warrant further exploration. For instance, we show that the Whole Genomic Duplication that affected a whole clade of yeasts, does not seems to have a statistically detectable phenotypic effect at our level of analysis. In this regard, further studies of fermentation under common‐garden conditions combined with comparative analyses are warranted.     

Tempo and mode of climatic niche evolution in Primates

Climatic niches have increasingly become a nexus in our understanding of a variety of ecological and evolutionary phenomena, from species distributions to latitudinal diversity gradients. Despite the increasing availability of comprehensive datasets on species ranges, phylogenetic histories, and georeferenced environmental conditions, studies on the evolution of climate niches have only begun to understand how niches evolve over evolutionary timescales. Here, using primates as a model system, we integrate recently developed phylogenetic comparative methods, species distribution patterns, and climatic data to explore primate climatic niche evolution, both among clades and over time. In general, we found that simple, constant‐rate models provide a poor representation of how climatic niches evolve. For instance, there have been shifts in the rate of climatic niche evolution in several independent clades, particularly in response to the increasingly cooler climates of the past 10 My. Interestingly, rate accelerations greatly outnumbered rate decelerations. These results highlight the importance of considering more realistic evolutionary models that allow for the detection of heterogeneity in the tempo and mode of climatic niche evolution, as well as to infer possible constraining factors for species distributions in geographical space.     

New molecular marker for phylogenetic reconstruction of black yeast-like fungi (Chaetothyriales) with hypothetical EIF2AK2 kinase gene

In eukaryotes, phosphorylation of the α-subunit of eIF2 is a mechanism to adjust cellular gene expression profiles in response to specific signals. The eIF2α kinases are a group of serine–threonine kinases that perform important functions in response to infection, proteotoxicity, and nutrient scavenging. The conserved nature of eIF2α kinases among fungi makes them potential evolutionary markers, which may contribute to deeper understanding of taxonomy and evolution. To date, only few studies are available of eIF2α kinases in black yeasts, which are members of Chaetothyriales containing potential agents of a gamut of major human diseases, such as chromoblastomycosis, phaeohyphomycosis and mycetoma. To establish the phylogenetic validity of sequences of eIF2α kinases hypothetical genes, we compared these genes between members of different classes of fungi, including black yeasts and allies, aiming at evaluation of the phylogeny of this group using an alternative molecular marker, compared to standard ribosomal genes. Trees generated with eIF2α kinase sequences of fungi were compared with those generated by ribosomal internal transcribed spacers (ITS rDNA) sequences from the same species. Sequences used were obtained from the protein Non-redundant database of NCBI, were aligned using CLUSTALX v1.8 and alignments were analyzed with RAxML v8.2.9 on the CIPRES Science Gateway portal. The trees generated had similar topologies, demonstrating that eIF2α kinases hypothetical gene sequences present a coherent reflection of evolution among fungi, compared to trees reconstructed by the use of ribosomal sequences. Our preliminary findings with a limited dataset strongly suggest that the evolution of kinases among black yeasts follows a similar path as revealed by ribosomal data, which underlines the validity of current taxonomy of black yeasts and relatives.     

The evolution of parental care and egg size: a comparative analysis in frogs

The evolution of parental care and egg size has attracted considerable attention and theoretical debate. Several different hypotheses have been proposed concerning the trajectories of parental care and egg size evolution and the order of specific evolutionary transitions. Few comparative studies have investigated the predictions of these hypotheses. Here, we investigate the evolutionary association between parental care and egg size in frogs in a phylogenetic context. Data on egg size and presence or absence of parental care in various species of frogs was gathered from the scientific literature. As a basis for our comparative analyses, we developed a phylogenetic supertree, by combining the results of multiple phylogenetic analyses in the literature using matrix representation parsimony. Using phylogenetic pairwise comparisons we demonstrated a significant association between the evolution of parental care and large egg size. We then used recently developed maximum likelihood methods to infer the evolutionary order of specific transitions. This analysis revealed that the evolution of large egg size typically precedes the evolution of parental care, rather than the reverse. We discuss the relevance of our results to previous hypotheses concerning the evolution of parental care and egg size.     

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.     

Relationships within tribe Lycieae (Solanaceae): paraphyly of Lycium and multiple origins of gender dimorphism

We infer phylogenetic relationships among Lycium, Grabowskia, and the monotypic Phrodus microphyllus, using DNA sequence data from the nuclear granule-bound starch synthase gene (GBSSI, waxy) and the chloroplast region trnT-trnF. This is the first comprehensive molecular phylogenetic study of tribe Lycieae (Solanaceae). In addition to providing an understanding of evolutionary relationships, we use the phylogenetic hypotheses to frame our studies of breeding system transitions, floral and fruit evolution, and biogeographical patterns within Lycieae. Whereas Lycium is distributed worldwide, Phrodus and the majority of Grabowskia species are restricted to South America. Tribe Lycieae is strongly supported as monophyletic, but Lycium likely includes both Grabowskia and Phrodus. Results also suggest a single dispersal event from the Americas to the Old World, and frequent dispersal between North and South America. The diversity of fruit types in Lycieae is discussed in light of dispersal patterns and recent work on fruit evolution across Solanaceae. Dimorphic gender expression has been studied previously within Lycium, and results indicate that transitions in sexual expression are convergent, occurring multiple times in North America (a revised estimate from previous studies) and southern Africa.     

Evolution of a membrane protein regulon in Saccharomyces

Expression variation is widespread between species. The ability to distinguish regulatory change driven by natural selection from the consequences of neutral drift remains a major challenge in comparative genomics. In this work, we used observations of mRNA expression and promoter sequence to analyze signatures of selection on groups of functionally related genes in Saccharomycete yeasts. In a survey of gene regulons with expression divergence between Saccharomyces cerevisiae and S. paradoxus, we found that most were subject to variation in trans-regulatory factors that provided no evidence against a neutral model. However, we identified one regulon of membrane protein genes controlled by unlinked cis- and trans-acting determinants with coherent effects on gene expression, consistent with a history of directional, nonneutral evolution. For this membrane protein group, S. paradoxus alleles at regulatory loci were associated with elevated expression and altered stress responsiveness relative to other yeasts. In a phylogenetic comparison of promoter sequences of the membrane protein genes between species, the S. paradoxus lineage was distinguished by a short branch length, indicative of strong selective constraint. Likewise, sequence variants within the S. paradoxus population, but not across strains of other yeasts, were skewed toward low frequencies in promoters of genes in the membrane protein regulon, again reflecting strong purifying selection. Our results support a model in which a distinct expression program for the membrane protein genes in S. paradoxus has been preferentially maintained by negative selection as the result of an increased importance to organismal fitness. These findings illustrate the power of integrating expression- and sequence-based tests of natural selection in the study of evolutionary forces that underlie regulatory change.     

Accommodating natural and sexual selection in butterfly wing pattern evolution

Visual patterns in animals may serve different functions, such as attracting mates and deceiving predators. If a signal is used for multiple functions, the opportunity arises for conflict among the different functions, preventing optimization for any one visual signal. Here we investigate the hypothesis that spatial separation of different visual signal functions has occurred in Bicyclus butterflies. Using phylogenetic reconstructions of character evolution and comparisons of evolutionary rates, we found dorsal surface characters to evolve at higher rates than ventral characters. Dorsal characters also displayed sex-based differences in evolutionary rates more often than did ventral characters. Thus, dorsal characters corresponded to our predictions of mate signalling while ventral characters appear to play an important role in predator avoidance. Forewing characters also fit a model of mate signalling, and displayed higher rates of evolution than hindwing characters. Our results, as well as the behavioural and developmental data from previous studies of Bicyclus species, support the hypothesis that spatial separation of visual signal functions has occurred in Bicyclus butterflies. This study is the first to demonstrate, in a phylogenetic framework, that spatial separation of signals used for mate signalling and those used for predator avoidance is a viable strategy to accommodate multiple signal functions. This signalling strategy has important ramifications on the developmental evolution of wing pattern elements and diversification of butterfly species.     

Fungal regulatory evolution: cis and trans in the balance

Regulatory divergence is likely a major driving force in evolution. Comparative genomics is being increasingly used to infer the evolution of gene regulation. Ascomycota fungi are uniquely suited among eukaryotes for regulatory evolution studies, due to broad phylogenetic scope, many sequenced genomes, and tractability of genomic analysis. Here we review recent advances in the identification of the contribution of cis- and trans-factors to expression divergence. Whereas current strategies have led to the discovery of surprising signatures and mechanisms, we still understand very little about the adaptive role of regulatory evolution. Empirical studies including experimental evolution, comparative functional genomics and hybrid and engineered strains are showing early promise toward deciphering the contribution of regulatory divergence to adaptation.     

Inferring Phylogenetic Networks with Maximum Pseudolikelihood under Incomplete Lineage Sorting

Phylogenetic networks are necessary to represent the tree of life expanded by edges to represent events such as horizontal gene transfers, hybridizations or gene flow. Not all species follow the paradigm of vertical inheritance of their genetic material. While a great deal of research has flourished into the inference of phylogenetic trees, statistical methods to infer phylogenetic networks are still limited and under development. The main disadvantage of existing methods is a lack of scalability. Here, we present a statistical method to infer phylogenetic networks from multi-locus genetic data in a pseudolikelihood framework. Our model accounts for incomplete lineage sorting through the coalescent model, and for horizontal inheritance of genes through reticulation nodes in the network. Computation of the pseudolikelihood is fast and simple, and it avoids the burdensome calculation of the full likelihood which can be intractable with many species. Moreover, estimation at the quartet-level has the added computational benefit that it is easily parallelizable. Simulation studies comparing our method to a full likelihood approach show that our pseudolikelihood approach is much faster without compromising accuracy. We applied our method to reconstruct the evolutionary relationships among swordtails and platyfishes (Xiphophorus: Poeciliidae), which is characterized by widespread hybridizations.     

Highly variable rates of genome rearrangements between hemiascomycetous yeast lineages

Hemiascomycete yeasts cover an evolutionary span comparable to that of the entire phylum of chordates. Since this group currently contains the largest number of complete genome sequences it presents unique opportunities to understand the evolution of genome organization in eukaryotes. We inferred rates of genome instability on all branches of a phylogenetic tree for 11 species and calculated species-specific rates of genome rearrangements. We characterized all inversion events that occurred within synteny blocks between six representatives of the different lineages. We show that the rates of macro- and microrearrangements of gene order are correlated within individual lineages but are highly variable across different lineages. The most unstable genomes correspond to the pathogenic yeasts Candida albicans and Candida glabrata. Chromosomal maps have been intensively shuffled by numerous interchromosomal rearrangements, even between species that have retained a very high physical fraction of their genomes within small synteny blocks. Despite this intensive reshuffling of gene positions, essential genes, which cluster in low recombination regions in the genome of Saccharomyces cerevisiae, tend to remain syntenic during evolution. This work reveals that the high plasticity of eukaryotic genomes results from rearrangement rates that vary between lineages but also at different evolutionary times of a given lineage.     

Neutral lipid fatty acid composition as trait and constraint in Collembola evolution

Functional traits determine the occurrence of species along environmental gradients and their coexistence with other species. Understanding how traits evolved among coexisting species helps to infer community assembly processes. We propose fatty acid composition in consumer tissue as a functional trait related to both food resources and physiological functions of species. We measured phylogenetic signal in fatty acid profiles of 13 field‐sampled Collembola (springtail) species and then combined the data with published fatty acid profiles of another 24 species. Collembola fatty acid profiles generally showed phylogenetic signal, with related species resembling each other. Long‐chain polyunsaturated fatty acids, related to physiological functions, demonstrated phylogenetic signal. In contrast, most food resource biomarker fatty acids and the ratios between bacterial, fungal, and plant biomarker fatty acids exhibited no phylogenetic signal. Presumably, fatty acids related to physiological functions have been constrained during Collembola evolutionary history: Species with close phylogenetic affinity experienced similar environments during divergence, while niche partitioning in food resources among closely related species favored species coexistence. Measuring phylogenetic signal in ecologically relevant traits of coexisting species provides an evolutionary perspective to contemporary assembly processes of ecological communities. Integrating phylogenetic comparative methods with community phylogenetic and trait‐based approaches may compensate for the limitations of each method when used alone and improve understanding of processes driving and maintaining assembly patterns.     

Application of the phylogenetic informativeness method to chloroplast markers: A test case of closely related species in tribe Hydrangeeae (Hydrangeaceae)

In evolutionary biology appropriate marker selection for the reconstruction of solid phylogenetic hypotheses is fundamental. One of the most challenging tasks addresses the appropriate choice of genomic regions in studies of closely related species. Robust phylogenetic frameworks are central to studies dealing with questions ranging from evolutionary and conservation biology, biogeography to plant breeding. Phylogenetic informativeness profiles provide a quantitative measure of the phylogenetic signal in markers and therefore a method for locus prioritization. The present work profiles phylogenetic informativeness of mostly non-coding chloroplast regions in an angiosperm lineage of closely related species: the popular ornamental tribe Hydrangeeae (Hydrangeaceae, Cornales, Asterids). A recent phylogenetic study denoted a case of resolution contrast between the two strongly supported clades within tribe Hydrangeeae. We evaluate the phylogenetic signal of 13 highly variable plastid markers for estimating relationships within and among the currently recognized monophyletic groups of this tribe. A selection of combined loci based on their phylogenetic informativeness retrieved more robust phylogenetic hypotheses than simply combining individual markers performing best with respect to resolution, nodal support and accuracy or those presenting the highest number of parsimony informative characters. We propose the rpl32–ndhF intergenic spacer (IGS), trnV–ndhC IGS, trnL–rpl32 IGS, psbT–petB region and ndhA intron as the best candidates for future phylogenetic studies in Hydrangeeae and potentially in other Asterids. We also contrasted the phylogenetic informativeness of coded indels against substitutions concluding that, despite their low phylogenetic informativeness, coded indels provide additional phylogenetic signal that is nearly free of noise. Phylogenetic relationships obtained from our total combined analyses showed improved resolution and nodal support with respect to recently published results.     

Repeated evolution and reversibility of self‐fertilization in the volvocine green algae*

Outcrossing and self‐fertilization are fundamental strategies of sexual reproduction, each with different evolutionary costs and benefits. Self‐fertilization is thought to be an evolutionary “dead‐end” strategy, beneficial in the short term but costly in the long term, resulting in self‐fertilizing species that occupy only the tips of phylogenetic trees. Here, we use volvocine green algae to investigate the evolution of self‐fertilization. We use ancestral‐state reconstructions to show that self‐fertilization has repeatedly evolved from outcrossing ancestors and that multiple reversals from selfing to outcrossing have occurred. We use three phylogenetic metrics to show that self‐fertilization is not restricted to the tips of the phylogenetic tree, a finding inconsistent with the view of self‐fertilization as a dead‐end strategy. We also find no evidence for higher extinction rates or lower speciation rates in selfing lineages. We find that self‐fertilizing species have significantly larger colonies than outcrossing species, suggesting the benefits of selfing may counteract the costs of increased size. We speculate that our macroevolutionary results on self‐fertilization (i.e., non‐tippy distribution, no decreased diversification rates) may be explained by the haploid‐dominant life cycle that occurs in volvocine algae, which may alter the costs and benefits of selfing.     

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8–38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group—a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.     

Climatic Niche Evolution in New World Monkeys (Platyrrhini)

Despite considerable interest in recent years on species distribution modeling and phylogenetic niche conservatism, little is known about the way in which climatic niches change over evolutionary time. This knowledge is of major importance to understand the mechanisms underlying limits of species distributions, as well as to infer how different lineages might be affected by anthropogenic climate change. In this study we investigate the tempo and mode climatic niche evolution in New World monkeys (Platyrrhini). Climatic conditions found throughout the distribution of 140 primate species were investigated using a principal component analysis, which indicated that mean temperature (particularly during the winter) is the most important climatic correlate of platyrrhine geographical distributions, accounting for nearly half of the interspecific variation in climatic niches. The effects of precipitation were associated with the second principal component, particularly with respect to the dry season. When models of trait evolution were fit to scores on each of the principal component axes, significant phylogenetic signal was detected for PC1 scores, but not for PC2 scores. Interestingly, although all platyrrhine families occupied comparable regions of climatic space, some aotid species such as Aotus lemurinus, A. jorgehernandezi, and A. miconax show highly distinctive climatic niches associated with drier conditions (high PC2 scores). This shift might have been made possible by their nocturnal habits, which could serve as an exaptation that allow them to be less constrained by humidity during the night. These results underscore the usefulness of investigating explicitly the tempo and mode of climatic niche evolution and its role in determining species distributions.