Ecological limits and diversification rate: alternative paradigms to explain the variation in species richness among clades and regions

Diversification rate is one of the most important metrics in macroecological and macroevolutionary studies. Here I demonstrate that diversification analyses can be misleading when researchers assume that diversity increases unbounded through time, as is typical in molecular phylogenetic studies. If clade diversity is regulated by ecological factors, then species richness may be independent of clade age and it may not be possible to infer the rate at which diversity arose. This has substantial consequences for the interpretation of many studies that have contrasted rates of diversification among clades and regions. Often, it is possible to estimate the total diversification experienced by a clade but not diversification rate itself. I show that the evidence for ecological limits on diversity in higher taxa is widespread. Finally, I explore the implications of ecological limits for a variety of ecological and evolutionary questions that involve inferences about speciation and extinction rates from phylogenetic data.     

Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates

Time-scales of viral evolution and emergence have been studied widely, but are often poorly understood. Molecular analyses of viral evolutionary time-scales generally rely on estimates of rates of nucleotide substitution, which vary by several orders of magnitude depending on the timeframe of measurement. We analysed data from all major groups of viruses and found a strong negative relationship between estimates of nucleotide substitution rate and evolutionary timescale. Strikingly, this relationship was upheld both within and among diverse groups of viruses. A detailed case study of primate lentiviruses revealed that the combined effects of sequence saturation and purifying selection can explain this time-dependent pattern of rate variation. Therefore, our analyses show that studies of evolutionary time-scales in viruses require a reconsideration of substitution rates as a dynamic, rather than as a static, feature of molecular evolution. Improved modelling of viral evolutionary rates has the potential to change our understanding of virus origins.     

Estimating diversification rates for higher taxa: BAMM can give problematic estimates of rates and rate shifts

Estimates of diversification rates are invaluable for many macroevolutionary studies. Recently, an approach called BAMM (Bayesian Analysis of Macro‐evolutionary Mixtures) has become widely used for estimating diversification rates and rate shifts. At the same time, several articles have concluded that estimates of net diversification rates from the method‐of‐moments (MS) estimators are inaccurate. Yet, no studies have compared the ability of these two methods to accurately estimate clade diversification rates. Here, we use simulations to compare their performance. We found that BAMM yielded relatively weak relationships between true and estimated diversification rates. This occurred because BAMM underestimated the number of rates shifts across each tree, and assigned high rates to small clades with low rates. Errors in both speciation and extinction rates contributed to these errors, showing that using BAMM to estimate only speciation rates is also problematic. In contrast, the MS estimators (particularly using stem group ages), yielded stronger relationships between true and estimated diversification rates, by roughly twofold. Furthermore, the MS approach remained relatively accurate when diversification rates were heterogeneous within clades, despite the widespread assumption that it requires constant rates within clades. Overall, we caution that BAMM may be problematic for estimating diversification rates and rate shifts.     

Nonstochastic variation of species-level diversification rates within angiosperms

Variations in the origination and extinction rates of species over geological time often are linked with a range of factors, including the evolution of key innovations, changes in ecosystem structure, and environmental factors such as shifts in climate and physical geography. Before hypothesizing causality of a single factor, it is critical to demonstrate that the observed variation in diversification is significantly greater than one would expect due to natural stochasticity in the evolutionary branching process. Here, we use a likelihood-ratio test to compare taxonomic rate heterogeneity to a neutral birth-death model, using data on well-supported sister pairs of taxa and their species richness. We test the likelihood that the distribution of extant species among angiosperm genera and families could be the result of constant diversification rates. Results strongly support the conclusion that there is significantly more heterogeneity in diversity at the species level within angiosperms than would be expected due to stochastic processes. This result is consistent in datasets of genus pairs and family pairs and is not affected significantly by degrading pairs to simulate inaccuracy in the assumption of simultaneous origin of sister taxa. When we parse taxon pairs among higher groups of angiosperms, results indicate that a constant rates model is not rejected by rosid and basal eudicot pairs but is rejected by asterid and eumagnoliid pairs. These results provide strong support for the hypothesis that species-level rates of origination and/or extinction have varied nonrandomly within angiosperms and that the magnitude of heterogeneity varies among major groups within angiosperms.     

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.     

Behavioural changes and the adaptive diversification of pigeons and doves

What factors determine the extent of evolutionary diversification remains a major question in evolutionary biology. Behavioural changes have long been suggested to be a major driver of phenotypic diversification by exposing animals to new selective pressures. Nevertheless, the role of behaviour in evolution remains controversial because behavioural changes can also retard evolutionary change by hiding genetic variation from selection. In the present study, we apply recently implemented Ornstein–Uhlenbeck evolutionary models to show that behavioural changes led to associated evolutionary responses in functionally relevant morphological traits of pigeons and doves (Columbiformes). Specifically, changes from terrestrial to arboreal foraging behaviour reconstructed in a set of phylogenies brought associated shorter tarsi and longer tails, consistent with functional predictions. Interestingly, the transition to arboreality accelerated the rates of evolutionary divergence, leading to an increased morphological specialization that seems to have subsequently constrained reversals to terrestrial foraging. Altogether, our results support the view that behaviour may drive evolutionary diversification, but they also highlight that its evolutionary consequences largely depend on the limits imposed by the functional demands of the adaptive zone.     

Exceptional among-lineage variation in diversification rates during the radiation of Australia's most diverse vertebrate clade

The disparity in species richness among groups of organisms is one of the most pervasive features of life on earth. A number of studies have addressed this pattern across higher taxa (e.g. 'beetles'), but we know much less about the generality and causal basis of the variation in diversity within evolutionary radiations at lower taxonomic scales. Here, we address the causes of variation in species richness among major lineages of Australia's most diverse vertebrate radiation, a clade of at least 232 species of scincid lizards. We use new mitochondrial and nuclear intron DNA sequences to test the extent of diversification rate variation in this group. We present an improved likelihood-based method for estimating per-lineage diversification rates from combined phylogenetic and taxonomic (species richness) data, and use the method in a hypothesis-testing framework to localize diversification rate shifts on phylogenetic trees. We soundly reject homogeneity of diversification rates among members of this radiation, and find evidence for a dramatic rate increase in the common ancestor of the genera Ctenotus and Lerista. Our results suggest that the evolution of traits associated with climate tolerance may have had a role in shaping patterns of diversity in this group.     

Biofilm and capsule formation of the diatom Achnanthidium minutissimum are affected by a bacterium

Photoautotrophic biofilms play an important role in various aquatic habitats and are composed of prokaryotic and/or eukaryotic organisms embedded in extracellular polymeric substances (EPS). We have isolated diatoms as well as bacteria from freshwater biofilms to study organismal interactions between representative isolates. We found that bacteria have a strong impact on the biofilm formation of the pennate diatom Achnanthidium minutissimum. This alga produces extracellular capsules of insoluble EPS, mostly carbohydrates (CHO), only in the presence of bacteria (xenic culture). The EPS themselves also have a strong impact on the aggregation and attachment of the algae. In the absence of bacteria (axenic culture), A. minutissimum did not form capsules and the cells grew completely suspended. Fractionation and quantification of CHO revealed that the diatom in axenic culture produces large amounts of soluble CHO, whereas in the xenic culture mainly insoluble CHO were detected. For investigation of biofilm formation by A. minutissimum, a bioassay was established using a diatom satellite Bacteroidetes bacterium that had been shown to induce capsule formation of A. minutissimum. Interestingly, capsule and biofilm induction can be achieved by addition of bacterial spent medium, indicating that soluble hydrophobic molecules produced by the bacterium may mediate the diatom/bacteria interaction. With the designed bioassay, a reliable tool is now available to study the chemical interactions between diatoms and bacteria with consequences for biofilm formation.     

Testing the molecular clock: molecular and paleontological estimates of divergence times in the Echinoidea (Echinodermata)

The phylogenetic relationships of 46 echinoids, with representatives from 13 of the 14 ordinal-level clades and about 70% of extant families commonly recognized, have been established from 3 genes (3,226 alignable bases) and 119 morphological characters. Morphological and molecular estimates are similar enough to be considered suboptimal estimates of one another, and the combined data provide a tree that, when calibrated against the fossil record, provides paleontological estimates of divergence times and completeness of their fossil record. The order of branching on the cladogram largely agrees with the stratigraphic order of first occurrences and implies that their fossil record is more than 85% complete at family level and at a resolution of 5-Myr time intervals. Molecular estimates of divergence times derived from applying both molecular clock and relaxed molecular clock models are concordant with estimates based on the fossil record in up to 70% of cases, with most concordant results obtained using Sanderson's semiparametric penalized likelihood method and a logarithmic-penalty function. There are 3 regions of the tree where molecular and fossil estimates of divergence time consistently disagree. Comparison with results obtained when molecular divergence dates are estimated from the combined (morphology + gene) tree suggests that errors in phylogenetic reconstruction explain only one of these. In another region the error most likely lies with the paleontological estimates because taxa in this region are demonstrated to have a very poor fossil record. In the third case, morphological and paleontological evidence is much stronger, and the topology for this part of the molecular tree differs from that derived from the combined data. Here the cause of the mismatch is unclear but could be methodological, arising from marked inequality of molecular rates. Overall, the level of agreement reached between these different data and methodological approaches leads us to believe that careful application of likelihood and Bayesian methods to molecular data provides realistic divergence time estimates in the majority of cases (almost 80% in this specific example), thus providing a remarkably well-calibrated phylogeny of a character-rich clade of ubiquitous marine benthic invertebrates.     

Chromosomal diversification in the flightless Western Mediterranean bushcricket genus Odontura (Orthoptera: Tettigoniidae: Phaneropterinae) inferred from molecular data

We used molecular characters to infer the phylogenetic position of the Western Mediterranean bushcricket genus Odontura and to trace its high karyotype diversity. Analysis of 1391 base pairs of two mitochondrial genes (COI and ND1) and one nuclear sequence (ITS2) was conducted. Phylogenetic topologies were estimated using maximum parsimony, maximum likelihood and likelihood‐based Bayesian inference. The genus Odontura is a phylogenetic outlier in respect of all other European Phaneropterinae genera and has been proposed to have originated from a hitherto unknown ancestor. Our results support the monophyly of the genus Odontura and the recognition of two subgenera: Odontura and Odonturella. We found that both Sicilian taxa of the subgenus Odontura have a completely identical morphology and song patterns. Combining these results, we proposed that both should be treated as subspecies: O. (Odontura) stenoxypha stenoxypha and O. (O.) st. arcuata. Bioacoustic data also proved to support independent markers, with song characteristics reflecting the molecular topology. Mapping the karyotypic characters onto the phylogenetic tree allows a reconstruction of the directions and transitional stages of chromosome differentiation. The number of autosomes within the genus Odontura ranges from 26 to 30. In addition to the ancestral X0 sex determination mechanism, neo‐XY and neo‐X₁X₂Y sex chromosomes have evolved independently.     

Variable evolutionary rates in the molecular evolution of mammalian growth hormones

In mammals pituitary growth hormone (GH) shows a slow basal rate of evolution (0.22 ± 0.03 × 10^–9 substitutions/amino acid site/year) which appears to have increased by at least 25–50-fold on two occasions, during the evolution of primates (to at least 10.8 ± 1.3 X 10^–9 substitutions/amino acid site/year) and artiodactyl ruminants (to at least 5.6 ± 1.3 X 10^–9 substitutions/amino acid site/year). That these rate increases are real, and not due to inadvertent comparison of nonorthologous genes, was established by showing that features of the GH gene sequences that are not expressed as mature hormone do not show corresponding changes in evolutionary rate. Thus, analysis of nonsynonymous substitutions in the coding sequence for the mature protein confirmed the rate increases seen in the primate and ruminant GHs, but analysis of nonsynonymous substitutions in the signal peptide sequence, synonymous substitutions in the coding sequence for signal peptide or mature protein, and 5 and 3 untranslated sequences showed no statistically significant changes in evolutionary rate. Evidence that the increases in evolutionary rate are probably due to positive selection is provided by the observation that in the cases of both ruminant and primate GHs the periods of rapid evolution were followed by a return to a slow rate similar to the basal rate seen in other mammalian GHs. Differences between the biological properties of GHs have been identified which may relate to these periods of rapid adaptive molecular evolution. On the basis of sequence data currently available (but excluding rodent GHs which show an intermediate rate, the basis of which is not clear) for most (90%) of evolutionary time mammalian GHs have been in the slow phase of evolution, with possibly most of the few amino acid substitutions that have occurred being neutral in nature. But most (80%) of the amino acid substitutions that have been introduced into GH during the course of mammalian evolution have been accepted during the rapid phases and were adaptive in nature.     

Interglacial microrefugia and diversification of a cactus species complex: phylogeography and palaeodistributional reconstructions for Pilosocereus aurisetus and allies

The role of Pleistocene climate changes in promoting evolutionary diversification in global biota is well documented, but the great majority of data regarding this subject come from North America and Europe, which were greatly affected by glaciation. The effects of Pleistocene changes on cold‐ and/or dry‐adapted species in tropical areas where glaciers were not present remain sparsely investigated. Many such species are restricted to small areas surrounded by unfavourable habitats, which may represent potential interglacial microrefugia. Here, we analysed the phylogeographic structure and diversification history of seven cactus species in the Pilosocereus aurisetus complex that are restricted to rocky areas with high diversity and endemism within the Neotropical savannas of eastern South America. We combined palaeodistributional estimates with standard phylogeographic approaches based on two chloroplast DNA regions (trnT‐trnL and trnS‐trnG), exon 1 of the nuclear gene PhyC and 10 nuclear microsatellite loci. Our analyses revealed a phylogeographic history marked by multiple levels of distributional fragmentation, isolation leading to allopatric differentiation and secondary contact among divergent lineages within the complex. Diversification and demographic events appear to have been affected by the Quaternary climatic cycles as a result of isolation in multiple patches of xerophytic vegetation. These small patches presently harbouring P. aurisetus populations seem to operate as microrefugia, both at present and during Pleistocene interglacial periods; the role of such microrefugia should be explored and analysed in greater detail.     

The molecular evolution of vertebrate growth hormones: A pattern of near-stasis interrupted by sustained bursts of rapid change

It has been demonstrated previously that in mammals the evolution of pituitary growth hormone shows an unusual pattern, with an underlying slow rate and at least two sustained bursts of rapid evolution (in the artiodactyls and primates), during which the rate increased at least 25-fold. It is demonstrated here that a similar pattern applies for growth hormone evolution throughout the vertebrates, with a basal rate similar to that seen in mammals, but bursts of rapid evolution in the amphibia and the elasmobranchs, and several bursts in the teleosts. The placental growth-hormone-like proteins of primates show a similar pattern. It is argued that the bursts of evolution seen for growth hormone are a consequence of selection and that this may reflect changes in the functions of the hormone additional to its basic growth-promoting actions.     

Synonymous and nonsynonymous rate variation in nuclear genes of mammals

A maximum likelihood approach was used to estimate the synonymous and nonsynonymous substitution rates in 48 nuclear genes from primates, artiodactyls, and rodents. A codon-substitution model was assumed, which accounts for the genetic code structure, transition/transversion bias, and base frequency biases at codon positions. Likelihood ratio tests were applied to test the constancy of nonsynonymous to synonymous rate ratios among branches (evolutionary lineages). It is found that at 22 of the 48 nuclear loci examined, the nonsynonymous/synonymous rate ratio varies significantly across branches of the tree. The result provides strong evidence against a strictly neutral model of molecular evolution. Our likelihood estimates of synonymous and nonsynonymous rates differ considerably from previous results obtained from approximate pairwise sequence comparisons. The differences between the methods are explored by detailed analyses of data from several genes. Transition/transversion rate bias and codon frequency biases are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for those factors. The likelihood approach is preferable, even for pairwise sequence comparison, because more-realistic models about the mutation and substitution processes can be incorporated in the analysis. Received: 17 May 1997 / Accepted: 28 September 1997     

Phylogeography of the heavily poached African common pangolin (Pholidota,Manis tricuspis) reveals six cryptic lineages as traceable signatures of Pleistocene diversification

Knowledge on faunal diversification in African rainforests remains scarce. We used phylogeography to assess (i) the role of Pleistocene climatic oscillations in the diversification of the African common pangolin (Manis tricuspis) and (ii) the utility of our multilocus approach for taxonomic delineation and trade tracing of this heavily poached species. We sequenced 101 individuals for two mitochondrial DNA (mtDNA), two nuclear DNA and one Y‐borne gene fragments (totalizing 2602 bp). We used a time‐calibrated, Bayesian inference phylogenetic framework and conducted character‐based, genetic and phylogenetic delineation of species hypotheses within African common pangolins. We identified six geographic lineages partitioned into western Africa, Ghana, the Dahomey Gap, western central Africa, Gabon and central Africa, all diverging during the Middle to Late Pleistocene. MtDNA (cytochrome b + control region) was the sole locus to provide diagnostic characters for each of the six lineages. Tree‐based Bayesian delimitation methods using single‐ and multilocus approaches gave high support for ‘species’ level recognition of the six African common pangolin lineages. Although the diversification of African common pangolins occurred during Pleistocene cyclical glaciations, causative correlation with traditional rainforest refugia and riverine barriers in Africa was not straightforward. We conclude on the existence of six cryptic lineages within African common pangolins, which might be of major relevance for future conservation strategies. The high discriminative power of the mtDNA markers used in this study should allow an efficient molecular tracing of the regional origin of African common pangolin seizures.