Diversity Dynamics in Nymphalidae Butterflies: Effect of Phylogenetic Uncertainty on Diversification Rate Shift Estimates

The species rich butterfly family Nymphalidae has been used to study evolutionary interactions between plants and insects. Theories of insect-hostplant dynamics predict accelerated diversification due to key innovations. In evolutionary biology, analysis of maximum credibility trees in the software MEDUSA (modelling evolutionary diversity using stepwise AIC) is a popular method for estimation of shifts in diversification rates. We investigated whether phylogenetic uncertainty can produce different results by extending the method across a random sample of trees from the posterior distribution of a Bayesian run. Using the MultiMEDUSA approach, we found that phylogenetic uncertainty greatly affects diversification rate estimates. Different trees produced diversification rates ranging from high values to almost zero for the same clade, and both significant rate increase and decrease in some clades. Only four out of 18 significant shifts found on the maximum clade credibility tree were consistent across most of the sampled trees. Among these, we found accelerated diversification for Ithomiini butterflies. We used the binary speciation and extinction model (BiSSE) and found that a hostplant shift to Solanaceae is correlated with increased net diversification rates in Ithomiini, congruent with the diffuse cospeciation hypothesis. Our results show that taking phylogenetic uncertainty into account when estimating net diversification rate shifts is of great importance, as very different results can be obtained when using the maximum clade credibility tree and other trees from the posterior distribution.     

Nodes in phylogenetic trees: the relation between imbalance and number of descendent species

The imbalance of a node in a phylogenetic tree can be defined in terms of the relative numbers of species (or higher taxa) on the branches that originate at the node. Empirically, imbalance also turns out to depend on the absolute total number of species on the branches: in a sample of large trees, nodes with more descendent species tend to be more unbalanced. Subsidiary analyses suggest that this pattern is not a result of errors in tree estimation. Instead, the increase in imbalance with species is consistent with a cumulative effect of differences in diversification rates between branches. [Equal-rates Markov model; imbalance; phylogeny shape; proportional-to-distinguishable-arrangements model.].     

Masting in whitebark pine (Pinus albicaulis) depletes stored nutrients

? In masting trees, synchronized, heavy reproductive events are thought to deplete stored resources and to impose a replenishment period before subsequent masting. However, direct evidence of resource depletion in wild, masting trees is very rare. Here, we examined the timing and magnitude (local vs individual-level) of stored nutrient depletion after a heavy mast event in Pinus albicaulis. ? In 2005, the mast year, we compared seasonal changes in leaf and sapwood nitrogen (N) and phosphorus (P) concentrations and leaf photosynthetic rates in cone-bearing branches, branches that never produced cones, and branches with experimentally removed cones. We also compared nutrient concentrations in cone branches and branches that had never had cones between 2005 and 2006, and measured tree ring width and new shoot growth during 2005. ? During the mast year, N or P depletion occurred only in tissue fractions of reproductive branches, where photosynthetic rates were reduced. However, by the end of the following year, nutrients were depleted in all branches, indicating individual-level resource depletion. New shoot and radial growth were not affected by masting. ? We provide direct evidence that mast events in wild trees deplete stored nutrients. Our results highlight the importance of evaluating reproductive costs over time and at the individual level.     

Exploring the tempo of species diversification in legumes

Whatever criteria are used to measure evolutionary success – species numbers, geographic range, ecological abundance, ecological and life history diversity, background diversification rates, or the presence of rapidly evolving clades – the legume family is one of the most successful lineages of flowering plants. Despite this, we still know rather little about the dynamics of lineage and species diversification across the family through the Cenozoic, or about the underlying drivers of diversification. There have been few attempts to estimate net species diversification rates or underlying speciation and extinction rates for legume clades, to test whether among-lineage variation in diversification rates deviates from null expectations, or to locate species diversification rate shifts on specific branches of the legume phylogenetic tree. In this study, time-calibrated phylogenetic trees for a set of species-rich legume clades – Calliandra, Indigofereae, Lupinus, Mimosa and Robinieae – and for the legume family as a whole, are used to explore how we might approach these questions. These clades are analysed using recently developed maximum likelihood and Bayesian methods to detect species diversification rate shifts and test for among-lineage variation in speciation, extinction and net diversification rates. Possible explanations for rate shifts in terms of extrinsic factors and/or intrinsic trait evolution are discussed. In addition, several methodological issues and limitations associated with these analyses are highlighted emphasizing the potential to improve our understanding of the evolutionary dynamics of legume diversification by using much more densely sampled phylogenetic trees that integrate information across broad taxonomic, geographical and temporal levels.     

Comparative analysis reveals that polyploidy does not decelerate diversification in fish

While the proliferation of the species‐rich teleost fish has been ascribed to an ancient genome duplication event at the base of this group, the broader impact of polyploidy on fish evolution and diversification remains poorly understood. Here, we investigate the association between polyploidy and diversification in several fish lineages: the sturgeons (Acipenseridae: Acipenseriformes), the botiid loaches (Botiidae: Cypriniformes), Cyprininae fishes (Cyprinidae: Cypriniformes) and the salmonids (Salmonidae: Salmoniformes). Using likelihood‐based evolutionary methodologies, we co‐estimate speciation and extinction rates associated with polyploid vs. diploid fish lineages. Family‐level analysis of Acipenseridae and Botiidae revealed no significant difference in diversification rates between polyploid and diploid relatives, while analysis of the subfamily Cyprininae revealed higher polyploid diversification. Additionally, order‐level analysis of the polyploid Salmoniformes and its diploid sister clade, the Esociformes, did not support a significantly different net diversification rate between the two groups. Taken together, our results suggest that polyploidy is generally not associated with decreased diversification in fish – a pattern that stands in contrast to that previously observed in plants. While there are notable differences in the time frame examined in the two studies, our results suggest that polyploidy is associated with different diversification patterns in these two major branches of the eukaryote tree of life.     

Multilocus analyses of an Antarctic fish species flock (Teleostei, Notothenioidei, Trematominae): phylogenetic approach and test of the early-radiation event

Clades that have undergone episodes of rapid cladogenesis are challenging from a phylogenetic point of view. They are generally characterised by short or missing internal branches in phylogenetic trees and by conflicting topologies among individual gene trees. This may be the case of the subfamily Trematominae, a group of marine teleosts of coastal Antarctic waters, which is considered to have passed through a period of rapid diversification. Despite much phylogenetic attention, the relationships among Trematominae species remain unclear. In contrast to previous studies that were mostly based on concatenated datasets of mitochondrial and/or single nuclear loci, we applied various single-locus and multilocus phylogenetic approaches to sequences from 11 loci (eight nuclear) and we also used several methods to assess the hypothesis of a radiation event in Trematominae evolution. Diversification rate analyses support the hypothesis of a period of rapid diversification during Trematominae history and only a few nodes in the hypothetical species tree were consistently resolved with various phylogenetic methods. We detected significant discrepancies among trees from individual genes of these species, most probably resulting from incomplete lineage sorting, suggesting that concatenation of loci is not the most appropriate way to investigate Trematominae species interrelationships. These data also provide information about the possible effects of historic climate changes on the diversification rate of this group of fish.     

The impact of gene-tree/species-tree discordance on diversification-rate estimation

Molecular phylogenies are often used to test hypotheses about the tempo and mode of speciation and extinction. One commonly used statistic is Pybus and Harvey's γ, which measures the density of ordered internode distances on an ultrametric tree to infer earlier (negative γ) or later (positive γ) bursts of diversification. However, coalescent theory predicts that γ might be biased toward negative values (inferring early bursts of diversification) when using gene trees rather than species trees. Gene divergences predate species divergences, increasingly so at higher effective population sizes (N(e)), and proportionally more so toward the tips of the tree. Thus, gene trees will have a higher density of older nodes in many cases (particularly at higher N(e)), due to the disproportionate lengthening of terminal branches. This will yield an artifactual signature of early bursts of diversification when estimating γ from gene trees. We simulate gene trees within species trees under both Yule (pure-birth) and birth-death processes, and demonstrate support for these predictions. However, for most realistic estimates of θ in natural populations, gene trees provide relatively good estimates of γ, despite the disproportionate overestimation of younger node ages. This is corroborated with an empirical dataset of North American fence lizards (Sceloporus).     

Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal-vicariance analysis of the thrushes (Aves: Turdus)

The phylogeny of the thrushes (Aves: Turdus) has been difficult to reconstruct due to short internal branches and lack of node support for certain parts of the tree. Reconstructing the biogeographic history of this group is further complicated by the fact that current implementations of biogeographic methods, such as dispersal-vicariance analysis (DIVA; Ronquist, 1997), require a fully resolved tree. Here, we apply a Bayesian approach to dispersal-vicariance analysis that accounts for phylogenetic uncertainty and allows a more accurate analysis of the biogeographic history of lineages. Specifically, ancestral area reconstructions can be presented as marginal distributions, thus displaying the underlying topological uncertainty. Moreover, if there are multiple optimal solutions for a single node on a certain tree, integrating over the posterior distribution of trees often reveals a preference for a narrower set of solutions. We find that despite the uncertainty in tree topology, ancestral area reconstructions indicate that the Turdus clade originated in the eastern Palearctic during the Late Miocene. This was followed by an early dispersal to Africa from where a worldwide radiation took place. The uncertainty in tree topology and short branch lengths seems to indicate that this radiation took place within a limited time span during the Late Pliocene. The results support the role of Africa as a probable source area for intercontinental dispersals as suggested for other passerine groups, including basal diversification within the songbird tree.     

Simulated herbivory advances autumn phenology in Acer rubrum

To determine the degree to which herbivory contributes to phenotypic variation in autumn phenology for deciduous trees, red maple (Acer rubrum) branches were subjected to low and high levels of simulated herbivory and surveyed at the end of the season to assess abscission and degree of autumn coloration. Overall, branches with simulated herbivory abscised ～7 % more leaves at each autumn survey date than did control branches within trees. While branches subjected to high levels of damage showed advanced phenology, abscission rates did not differ from those of undamaged branches within trees because heavy damage induced earlier leaf loss on adjacent branch nodes in this treatment. Damaged branches had greater proportions of leaf area colored than undamaged branches within trees, having twice the amount of leaf area colored at the onset of autumn and having ～16 % greater leaf area colored in late October when nearly all leaves were colored. When senescence was scored as the percent of all leaves abscised and/or colored, branches in both treatments reached peak senescence earlier than did control branches within trees: dates of 50 % senescence occurred 2.5 days earlier for low herbivory branches and 9.7 days earlier for branches with high levels of simulated damage. These advanced rates are of the same time length as reported delays in autumn senescence and advances in spring onset due to climate warming. Thus, results suggest that should insect damage increase as a consequence of climate change, it may offset a lengthening of leaf life spans in some tree species.     

Understanding angiosperm diversification using small and large phylogenetic trees

How will the emerging possibility of inferring ultra-large phylogenies influence our ability to identify shifts in diversification rate? For several large angiosperm clades (Angiospermae, Monocotyledonae, Orchidaceae, Poaceae, Eudicotyledonae, Fabaceae, and Asteraceae), we explore this issue by contrasting two approaches: (1) using small backbone trees with an inferred number of extant species assigned to each terminal clade and (2) using a mega-phylogeny of 55473 seed plant species represented in GenBank. The mega-phylogeny approach assumes that the sample of species in GenBank is at least roughly proportional to the actual species diversity of different lineages, as appears to be the case for many major angiosperm lineages. Using both approaches, we found that diversification rate shifts are not directly associated with the major named clades examined here, with the sole exception of Fabaceae in the GenBank mega-phylogeny. These agreements are encouraging and may support a generality about angiosperm evolution: major shifts in diversification may not be directly associated with major named clades, but rather with clades that are nested not far within these groups. An alternative explanation is that there have been increased extinction rates in early-diverging lineages within these clades. Based on our mega-phylogeny, the shifts in diversification appear to be distributed quite evenly throughout the angiosperms. Mega-phylogenetic studies of diversification hold great promise for revealing new patterns, but we will need to focus more attention on properly specifying null expectation.     

Crown alterations induced by decline: a study of relationships between growth rate and crown morphology in beech (Fagus sylvatica L.)

One of the first symptoms expressed by declining trees is reduced growth in stem diameter and length increment. The possibility of a relationship between length increment and crown thinning in beech (Fagus sylvatica L.) was investigated by developing a computer model to simulate first order branching patterns of the apical 2 m of monopodially branching beech trees, 70–100 years old, for a range of length increment rates. The model was based on values for branching angle, main axis and branch length increment, number of branches produced per year and branch mortality rates for six healthy and declining trees. Shoot growth rates in the apical 2 m of the sample trees ranged from about 5 cm/year (decline class 3) to 43 cm/ year (healthy). Simulations of branching patterns in the apical 2 m of trees growing at different rates indicated that, when growth rate exceeded about 20 cm/year, total first order branch length and area explored were independent of growth rate. When growth rates fell below this value there was a reduction in total area explored and first order branch length due primarily to the formation of fewer branches. More acute branching angles contributed to a reduction in the area explored. Growth rate-related crown thinning could increase the risk of bark necrosis and secondary pathogen infection during dry and/or hot spells.     

Bats, clocks, and rocks: diversification patterns in Chiroptera

Identifying nonrandom clade diversification is a critical first step toward understanding the evolutionary processes underlying any radiation and how best to preserve future phylogenetic diversity. However, differences in diversification rates have not been quantitatively assessed for the majority of groups because of the lack of necessary analytical tools (e.g., complete species-level phylogenies, estimates of divergence times, and robust statistics which incorporate phylogenetic uncertainty and test appropriate null models of clade growth). Here, for the first time, we investigate diversification rate heterogeneity in one of the largest groups studied thus far, the bats (Mammalia: Chiroptera). We use a recent, robust statistical approach (whole-tree likelihood-based relative rate tests) on complete dated species-level supertree phylogenies. As has been demonstrated previously for most other groups, among-lineage diversification rate within bats has not been constant. However, we show that bat diversification is more heterogeneous than in other mammalian clades thus far studied. The whole-tree likelihood-based relative rates tests suggest that clades within the families Phyllostomidae and Molossidae underwent a number of significant changes in relative diversification rate. There is also some evidence for rate shifts within Pteropodidae, Emballonuridae, Rhinolophidae, Hipposideridae, and Vespertilionidae, but the significance of these shifts depends on polytomy resolution within each family. Diversification rate in bats has also not been constant, with the largest diversification rate shifts occurring 30-50 million years ago, a time overlapping with the greatest number of shifts in flowering plant diversification rates.     

Incorporating phylogenetic uncertainty on phylogeny‐based palaeontological dating and the timing of turtle diversification

Methods improving the performance of molecular dating of divergence time of clades have improved dramatically in recent years. The calibration of molecular dating using the first appearance of a clade in the fossil record is a crucial step towards inferring the minimal diversification time of various groups and the choice of extinct taxa can strongly influence the molecular dates. Here, we evaluate the uncertainty on the phylogenetic position of extinct taxa through non‐parametric bootstrapping. The recognition of phylogenetic uncertainty resulted in the definition of the Bootstrap Uncertainty Range (BUR) for the age of first appearance of a given clade. The BUR is calculated as the interval of geological time in which the diversification of a given clade can be inferred to have occurred, based on the temporal information of the fossil record and the topologies of the bootstrap trees. Divergence times based on BUR analyses were calculated for three clades of turtles: Testudines, Pleurodira and Cryptodira. This resulted in extensive uncertainty ranges of topology‐dependent minimal divergence dates for these clades.     

A phylogenetic analysis of genes of the influenza virus. Phylogenetic trees and fixation rates

The phylogenetic trees of influenza virus genes of hemagglutinins, neuraminidases, and of NS genes were composed. Considering properties of synonimic replacements to be neutral and their rates constant at each tree, the dates of ancestor branch points were calculated, and the rates of fixation of synonimic (Ks) and non-synonimic (Kns) replacements estimated. The epidemic branches were mostly shown to be "deadlocks", non-epidemic ones being internal or "roots." The ratios of the numbers of synonimic to non-synonimic replacements (vs/vns) were correspondingly 1.32+/-0.42 and 4.78+/-1.28 for all trees, the difference being significant. It was shown that the dated branch points for hemagglutinins are non-randomly clustered around the initial points of the main genetic shifts of the A-type virus, corresponding to the influenza pandemics. It seems that these ancestor forms of virus behave similar to the "train" of these shifts, reproducing together with the pandemic forms under conditions of decreased immune resistance of host population. The rates of fixation of non-synonimic replacements in the epidemic branches of this tree are 4 times increased, as compared to non-epidemic ones.     

Trait evolution rates shape continental patterns of species richness in North America's most diverse angiosperm genus (Carex,Cyperaceae)

Ecological opportunity has been associated with increases in diversification rates across the tree of life. Under an ecological diversification model, the emergence of novel environments is hypothesized to promote morpho- and ecospace evolution. Whether this model holds at the clade level within the most species-rich angiosperm genus found in North America (Carex, Cyperaceae) is yet to be tested. Recent works demonstrate a temporal coupling of climate cooling and widespread colonization of Carex in North America, implicating ecological diversification. In addition, research has consistently found asymmetric patterns of lineage-level diversification in the genus. Why does variation in clade sizes exist in the genus? Is ecological diversification involved? In this study, we tested whether rates of morphological and ecological trait evolution are correlated with clade-level species richness in Carex of North America north of Mexico. We constructed a phylogeny of 477 species—an almost complete regional sample. We estimated rates of evolution of morphological traits, habitat, and climatic niche and assessed whether differences in rates of evolution correlate with species richness differences in replicate non-nested sister clades. Our work demonstrates significant positive correlations between climatic niche rates, habitat and reproductive morphological evolution, and species richness. This coupling of trait and niche evolution and species richness in a diverse, continental clade sample strongly suggests that the ability of clades to explore niche and functional space has shaped disparities in richness and functional diversity across the North American flora region. Our findings highlight the importance of the evolutionary history of trait and niche evolution in shaping continental and regional floras.     

Whole-tree methods for detecting differential diversification rates

Prolific cladogenesis, adaptive radiation, species selection, key innovations, and mass extinctions are a few examples of biological phenomena that lead to differential diversification among lineages. Central to the study of differential diversification rates is the ability to distinguish chance variation from that which requires deterministic explanation. To detect diversification rate variation among lineages, we propose a number of methods that incorporate information on the topological distribution of species diversity from all internal nodes of a phylogenetic tree. These whole-tree methods (M(Pi), M(Sigma), and M(R)) are explicitly connected to a null model of random diversification--the equal-rates Markov (ERM) random branching model--and an alternative model of differential diversification: M(Pi) is based on the product of individual nodal ERM probabilities; M(Sigma) is based on the sum of individual nodal ERM probabilities, and M(R) is based on a transformation of ERM probabilities that corresponds to a formalized system that orders trees by their relative symmetry. These methods have been implemented in a freely available computer program, SYMMETREE, to detect clades with variable diversification rates, thereby allowing the study of biological processes correlated with and possibly causal to shifts in diversification rate. Application of these methods to several published phylogenies demonstrates their ability to contend with relatively large, incompletely resolved trees. These topology-based methods do not require estimates of relative branch lengths, which should facilitate the analysis of phylogenies, such as supertrees, for which such data are unreliable or unavailable.     

DO FRESHWATER FISHES DIVERSIFY FASTER THAN MARINE FISHES? A TEST USING STATE‐DEPENDENT DIVERSIFICATION ANALYSES AND MOLECULAR PHYLOGENETICS OF NEW WORLD SILVERSIDES (ATHERINOPSIDAE)

Freshwater habitats make up only ～0.01% of available aquatic habitat and yet harbor 40% of all fish species, whereas marine habitats comprise >99% of available aquatic habitat and have only 60% of fish species. One possible explanation for this pattern is that diversification rates are higher in freshwater habitats than in marine habitats. We investigated diversification in marine and freshwater lineages in the New World silverside fish clade Menidiinae (Teleostei, Atherinopsidae). Using a time‐calibrated phylogeny and a state‐dependent speciation–extinction framework, we determined the frequency and timing of habitat transitions in Menidiinae and tested for differences in diversification parameters between marine and freshwater lineages. We found that Menidiinae is an ancestrally marine lineage that independently colonized freshwater habitats four times followed by three reversals to the marine environment. Our state‐dependent diversification analyses showed that freshwater lineages have higher speciation and extinction rates than marine lineages. Net diversification rates were higher (but not significant) in freshwater than marine environments. The marine lineage‐through time (LTT) plot shows constant accumulation, suggesting that ecological limits to clade growth have not slowed diversification in marine lineages. Freshwater lineages exhibited an upturn near the recent in their LTT plot, which is consistent with our estimates of high background extinction rates. All sequence data are currently being archived on Genbank and phylogenetic trees archived on Treebase.     

Phylogenetic relationships in Epidendroideae (Orchidaceae), one of the great flowering plant radiations: progressive specialization and diversification

Background and Aims The largest subfamily of orchids, Epidendroideae, represents one of the most significant diversifications among flowering plants in terms of pollination strategy, vegetative adaptation and number of species. Although many groups in the subfamily have been resolved, significant relationships in the tree remain unclear, limiting conclusions about diversification and creating uncertainty in the classification. This study brings together DNA sequences from nuclear, plastid and mitochrondrial genomes in order to clarify relationships, to test associations of key characters with diversification and to improve the classification.Methods Sequences from seven loci were concatenated in a supermatrix analysis for 312 genera representing most of epidendroid diversity. Maximum-likelihood and parsimony analyses were performed on this matrix and on subsets of the data to generate trees and to investigate the effect of missing values. Statistical character-associated diversification analyses were performed.Key Results Likelihood and parsimony analyses yielded highly resolved trees that are in strong agreement and show significant support for many key clades. Many previously proposed relationships among tribes and subtribes are supported, and some new relationships are revealed. Analyses of subsets of the data suggest that the relatively high number of missing data for the full analysis is not problematic. Diversification analyses show that epiphytism is most strongly associated with diversification among epidendroids, followed by expansion into the New World and anther characters that are involved with pollinator specificity, namely early anther inflexion, cellular pollinium stalks and the superposed pollinium arrangement.Conclusions All tested characters show significant association with speciation in Epidendroideae, suggesting that no single character accounts for the success of this group. Rather, it appears that a succession of key features appeared that have contributed to diversification, sometimes in parallel.     

AUGIST: inferring species trees while accommodating gene tree uncertainty

SUMMARY: AUGIST (accomodating uncertainty in genealogies while inferring species tress) is a new software package for inferring species trees while accommodating uncertainty in gene genealogies. It is written for the Mesquite software system and provides sampling procedures to incorporate uncertainty in gene tree reconstruction while providing confidence estimates for inferred species trees. AVAILABILITY: http://www.lycaenid.org/augist/