The ghosts of Gondwana and Laurasia in modern liverwort distributions

Recent advances in phylogenetics and, in particular, molecular dating, indicate that transoceanic dispersal has played an important role in shaping plant and animal distributions, obscuring any effect of tectonic history. Taxonomic sampling in biogeographic studies is, however, systematically biased towards vertebrates and higher plants and the possibility remains that a much stronger signature of ancient vicariance might be evident among other organisms, particularly among basal land plants. Here, an explicit Bayesian model‐based approach was used to investigate global‐scale biogeographic patterns among liverwort genera and to determine whether the patterns identified are consistent with the expectations of vicariance or dispersal scenarios. The distribution of each genus was mapped onto the phylograms describing the floristic affinities among areas in order to define the synapomorphic transitions supporting the observed groupings. The probabilities of change in a branch were calculated by implementing the Markov model of BayesTraits. The consistent ambiguity in ancestral state reconstructions returned by the unconstrained, two‐rate model indicated that the overall signal in the data was weak, leading us to test the performance of competing, explicit models. The analyses resolved clades of geographic areas that are mostly consistent with the kingdoms traditionally identified for plants and animals, but with strikingly lower rates of endemism. The major split observed in the phylograms is into almost entirely Laurasian and Gondwanan clades. Other patterns recovered by the analyses, including Wallace's line and the South Atlantic Disjunction, have also traditionally been interpreted in terms of vicariance. These observations contrast with the idea that, in spore‐dispersed organisms like bryophytes and pteridophytes, dispersal obscures evidence of vicariance. However, some discrepancies between the liverwort trees and expectations from a continental drift scenario were observed, such as the sister‐group relationship of the Australian and New Zealand floras, which is supported by the co‐occurrence of many genera, often endemic to these two areas. Together with an interpretation of the results within a phylogenetic context, our analyses suggest that patterns, which are at first sight consistent with an ancient vicariance hypothesis, may, in fact, conceal a complex mixture of relictual distributions and more recent, asymmetrical dispersal events. Our results provide a framework for testing specific evolutionary hypotheses concerning the extremely low levels of endemism in bryophytes and in particular, the significance of dispersal and cryptic diversification.     

ANCESTRAL CHARACTER ESTIMATION UNDER THE THRESHOLD MODEL FROM QUANTITATIVE GENETICS

Evolutionary biology is a study of life's history on Earth. In researching this history, biologists are often interested in attempting to reconstruct phenotypes for the long extinct ancestors of living species. Various methods have been developed to do this on a phylogeny from the data for extant taxa. In the present article, I introduce a new approach for ancestral character estimation for discretely valued traits. This approach is based on the threshold model from evolutionary quantitative genetics. Under the threshold model, the value exhibited by an individual or species for a discrete character is determined by an underlying, unobserved continuous trait called “liability.” In this new method for ancestral state reconstruction, I use Bayesian Markov chain Monte Carlo (MCMC) to sample the liabilities of ancestral and tip species, and the relative positions of two or more thresholds, from their joint posterior probability distribution. Using data simulated under the model, I find that the method has very good performance in ancestral character estimation. Use of the threshold model for ancestral state reconstruction relies on a priori specification of the order of the discrete character states along the liability axis. I test the use of a Bayesian MCMC information theoretic criterion based approach to choose among different hypothesized orderings for the discrete character. Finally, I apply the method to the evolution of feeding mode in centrarchid fishes.     

MOLECULAR PHYLOGENETIC ANALYSIS OF LIFE-HISTORY EVOLUTION IN ASTERINID STARFISH

We analyzed phylogenetic relationships among 12 nominal species of starfish in the genera Patiriella and Asterina (Order Valvatida, Family Asterinidae), based on complete sequences for a mitochondrial protein coding gene (cytochrome oxidase subunit I) and five mitochondrial transfer RNA genes (alanine, leucine, asparagine, glutamine, and proline) (1923 bp total). The resulting phylogeny was used to test a series of hypotheses about the evolution of life-history traits. (1) A complex, feeding, planktonic larva is probably ancestral for these starfish, but this is not the most parsimonious reconstruction of ancestral larval states. (2) The feeding larval form was lost at least four times among these species, and three of these losses occurred among members of a single clade. (3) Small adult size evolved before both cases of hermaphroditism and viviparous brooding, but viviparity was not always preceded by an intermediate form of external brooding. (4) An ordered transformation series from feeding planktonic development to viviparous brooding has been predicted for starfish, but we could not find an example of this transformation series. (5) Viviparity evolved recently (< 2 Mya). (6) Both species selection and transformation of lineages may have contributed to the accumulation of species with nonfeeding development among these starfish. (7) Neither Asterina nor Patiriella are monophyletic genera. Larval forms and life-history traits of these starfish have evolved freely under no obvious constraints, contrary to the widely assumed evolutionary conservatism of early development.     

LOSS AND RE‐EVOLUTION OF COMPLEX LIFE CYCLES IN MARSUPIAL FROGS: DOES ANCESTRAL TRAIT RECONSTRUCTION MISLEAD?

Using phylogeny-based methods to identify evolutionary transitions has become an integral part of evolutionary biology. Here, we demonstrate the potential for these methods to give statistically well-supported but misleading inferences about character evolution. We also show how inferences of character evolution can be informed using GIS-based methods to reconstruct ancestral environmental regimes. We reconstruct a phylogeny for marsupial frogs (Hemiphractidae) using nuclear and mitochondrial DNA sequences and estimate patterns of life-history evolution across the resulting tree. We find that Gastrotheca species with complex life cycles (i.e., egg, tadpole, and adult stages) are phylogenetically nested among species and genera with direct development (i.e., egg and adult stages only). Assuming a single rate for gains and losses in likelihood reconstructions, there is strong statistical support for the hypothesis that the tadpole stage was lost early in the phylogeny but reappeared within Gastrotheca. Assuming different rates of gain and loss, the model with significantly higher statistical support, the tadpole stage seems to have been lost multiple times but never regained. Given that both hypotheses cannot be correct, at least one reconstruction model must be giving well-supported but misleading results. Several lines of evidence (including GIS-based reconstructions of the ancestral climatic regime) suggest that the former hypothesis is correct, and that the tadpole stage has evolved from direct development within Gastrotheca, the only known case of such a reversal in frogs.     

Evolutionary history of vegetative reproduction in Porpidia s.L. (Lichen-forming ascomycota)

The evolutionary history of gains and losses of vegetative reproductive propagules (soredia) in Porpidia s.l., a group of lichen-forming ascomycetes, was clarified using Bayesian Markov chain Monte Carlo (MCMC) approaches to monophyly tests and a combined MCMC and maximum likelihood approach to ancestral character state reconstructions. The MCMC framework provided confidence estimates for the reconstructions of relationships and ancestral character states, which formed the basis for tests of evolutionary hypotheses. Monophyly tests rejected all hypotheses that predicted any clustering of reproductive modes in extant taxa. In addition, a nearest-neighbor statistic could not reject the hypothesis that the vegetative reproductive mode is randomly distributed throughout the group. These results show that transitions between presence and absence of the vegetative reproductive mode within Porpidia s.l. occurred several times and independently of each other. Likelihood reconstructions of ancestral character states at selected nodes suggest that--contrary to previous thought--the ancestor to Porpidia s.l. already possessed the vegetative reproductive mode. Furthermore, transition rates are reconstructed asymmetrically with the vegetative reproductive mode being gained at a much lower rate than it is lost. A cautious note has to be added, because a simulation study showed that the ancestral character state reconstructions were highly dependent on taxon sampling. However, our central conclusions, particularly the higher rate of change from vegetative reproductive mode present to absent than vice versa within Porpidia s.l., were found to be broadly independent of taxon sampling.     

Macroevolutionary synthesis of flowering plant sexual systems

Sexual system is a key determinant of genetic variation and reproductive success, affecting evolution within populations and within clades. Much research in plants has focused on evolutionary transitions away from the most common state of hermaphroditism and toward the rare state of dioecy (separate sexes). Rather than transitions predominantly toward greater sexual differentiation, however, evolution may proceed in the direction of lesser sexual differentiation. We analyzed the macroevolutionary dynamics of sexual system in angiosperm genera that contain both dioecious and nondioecious species. Our phylogenetic analyses encompass a total of 2145 species from 40 genera. Overall, we found little evidence that rates of sexual system transitions are greater in any direction. Counting the number of inferred state changes revealed a mild prevalence of transitions away from hermaphroditism and away from dioecy, toward states of intermediate sexual differentiation. We identify genera in which future studies of sexual system evolution might be especially productive, and we discuss how integrating genetic or population‐level studies of sexual system could improve the power of phylogenetic comparative analyses. Our work adds to the evidence that different selective pressures and constraints act in different groups, helping maintain the variety of sexual systems observed among plants.     

Identifying and evaluating causes of alternative community states in wetland plant communities

Ecological communities can vary greatly in species composition. Often this variation is discontinuous, in that abrupt changes in composition occur over small distances in space or short periods of time. A wide range of hypotheses from different subfields of ecology have been proposed to explain these patterns. I suggest a framework to quantitatively evaluate these hypotheses with observational data by characterizing 1) how community composition varies across sites in space, 2) how community composition varies through time, and 3) the possible drivers of this variation. I applied this approach to understand the community composition of producers in temporary and semipermanent wetlands in Michigan, USA. I identified several distinct community states which were variously dominated by particular plant functional groups (submerged, floating or emergent plants) or had no plants throughout a season. Evaluating possible hypotheses to explain this variation, I found that similar communities were not necessarily clustered near each other, suggesting that dispersal was not limited for these plants. Some sites exhibited a great deal of change in plant composition among years, shifting between two community states, but there was relatively little change at sites within a year. Moreover, these shifts did not occur in a particular order to suggest directional change or repeating cycles. Community composition was associated with several environmental variables such as pH, light and depth, and multivariate analyses suggested that species had complex, nonlinear responses to these possible drivers. Alternative stable states and interactions among multiple nonlinear drivers best explained the patterns observed in these wetlands. By formalizing initial data collection in other systems with the framework suggested here, we may gain insight into the causes of alternative community states beyond wetlands and the role of climate change and other anthropogenic forces in precipitating transitions between states.     

Exact and heuristic algorithms for the Indel Maximum Likelihood Problem.

Given a multiple alignment of orthologous DNA sequences and a phylogenetic tree for these sequences, we investigate the problem of reconstructing the most likely scenario of insertions and deletions capable of explaining the gaps observed in the alignment. This problem, that we called the Indel Maximum Likelihood Problem (IMLP), is an important step toward the reconstruction of ancestral genomics sequences, and is important for studying evolutionary processes, genome function, adaptation and convergence. We solve the IMLP using a new type of tree hidden Markov model whose states correspond to single-base evolutionary scenarios and where transitions model dependencies between neighboring columns. The standard Viterbi and Forward-backward algorithms are optimized to produce the most likely ancestral reconstruction and to compute the level of confidence associated to specific regions of the reconstruction. A heuristic is presented to make the method practical for large data sets, while retaining an extremely high degree of accuracy. The methods are illustrated on a 1-Mb alignment of the CFTR regions from 12 mammals.     

Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis

In historical biogeography, model-based inference methods for reconstructing the evolution of geographic ranges on phylogenetic trees are poorly developed relative to the diversity of analogous methods available for inferring character evolution. We attempt to rectify this deficiency by constructing a dispersal-extinction-cladogenesis (DEC) model for geographic range evolution that specifies instantaneous transition rates between discrete states (ranges) along phylogenetic branches and apply it to estimating likelihoods of ancestral states (range inheritance scenarios) at cladogenesis events. Unlike an earlier version of this approach, the present model allows for an analytical solution to probabilities of range transitions as a function of time, enabling free parameters in the model, rates of dispersal, and local extinction to be estimated by maximum likelihood. Simulation results indicate that accurate parameter estimates may be difficult to obtain in practice but also show that ancestral range inheritance scenarios nevertheless can be correctly recovered with high success if rates of range evolution are low relative to the rate of cladogenesis. We apply the DEC model to a previously published, exemplary case study of island biogeography involving Hawaiian endemic angiosperms in Psychotria (Rubiaceae), showing how the DEC model can be iteratively refined from inspecting inferences of range evolution and also how geological constraints involving times of island origin may be imposed on the likelihood function. The DEC model is sufficiently similar to character models that it might serve as a gateway through which many existing comparative methods for characters could be imported into the realm of historical biogeography; moreover, it might also inspire the conceptual expansion of character models toward inclusion of evolutionary change as directly coincident, either as cause or consequence, with cladogenesis events. The DEC model is thus an incremental advance that highlights considerable potential in the nascent field of model-based historical biogeographic inference.     

The polyphyletic origin of endoparasitism in the cyclostome lineages of Braconidae (Hymenoptera)

Abstract. The evolution of endoparasitism in braconid wasps is demonstrably polyphyletic, not only within the family as a whole but even within the rogadine/exothecine assemblage of genera. A cladistic analysis of rogadine and exothecine genera, coupled with an optimization approach to estimating ancestral character states, is used to analyse the evolution of various aspects of endoparasitism within this assemblage. At least two independent origins of endoparasitism are hypothesized-one in the Rodagini and the other in an assemblage including the Gnaptodontini, Opiinae and Alysiinae.
Phylogenetic analyses and hypotheses can be useful not only in interpreting evolutionary trends, but in guiding and providing insights into more general comparative behavioural and physiological studies of parasitism.     

THE DIVERSITY AND EVOLUTION OF BATESIAN MIMICRY IN PAPILIO SWALLOWTAIL BUTTERFLIES

Papilio swallowtail butterflies exhibit a remarkable diversity of Batesian mimicry, manifested in several sex-limited and polymorphic types. There is little understanding of how this diversity is distributed within Papilio , and how different mimicry types have evolved in relation to each other. To answer these questions, I present a graphical model that connects various mimicry types by hypothetical character state changes within a phylogenetic framework. A maximum likelihood analysis of evolution of mimicry types on the Papilio phylogeny showed that sexually monomorphic mimicry and female-limited mimicry have evolved repeatedly but predominantly independently in different clades. However, transitions between these mimicry types are rarely observed. The frequency distribution of character state changes was skewed in favor of the evolution of mimicry, whereas many theoretically plausible character state changes, especially evolutionary loss of mimicry, were not evident. I discuss these findings in relation to studying the tempo of evolutionary change, loss of traits, and directionality and connectivity among character states. The pathway approach and phylogenetic patterns of mimicry demonstrated in Papilio are useful to test novel hypotheses regarding the diversity and evolutionary directionality of Batesian mimicry in other systems.     

Evolution of growth habit, inflorescence architecture, flower size, and fruit type in Rubiaceae: its ecological and evolutionary implications

During angiosperm evolution, innovations in vegetative and reproductive organs have resulted in tremendous morphological diversity, which has played a crucial role in the ecological success of flowering plants. Morindeae (Rubiaceae) display considerable diversity in growth form, inflorescence architecture, flower size, and fruit type. Lianescent habit, head inflorescence, small flower, and multiple fruit are the predominant states, but arborescent habit, non-headed inflorescence, large flower, and simple fruit states occur in various genera. This makes Morindeae an ideal model for exploring the evolutionary appearances and transitions between the states of these characters. We reconstructed ancestral states for these four traits using a bayesian approach and combined nuclear/chloroplast data for 61 Morindeae species. The aim was to test three hypotheses: 1) self-supporting habit is generally ancestral in clades comprising both lianescent and arborescent species; 2) changes from lianescent to arborescent habit are uncommon due to "a high degree of specialization and developmental burden"; 3) head inflorescences and multiple fruits in Morindeae evolved from non-headed inflorescences and simple fruits, respectively. Lianescent habit, head inflorescence, large flower, and multiple fruit are inferred for Morindeae, making arborescent habit, non-headed inflorescence, small flower, and simple fruit derived within the tribe. The rate of change from lianescent to arborescent habit is much higher than the reverse change. Therefore, evolutionary changes between lianescent and arborescent forms can be reversible, and their frequency and trends vary between groups. Moreover, these changes are partly attributed to a scarcity of host trees for climbing plants in more open habitats. Changes from large to small flowers might have been driven by shifts to pollinators with progressively shorter proboscis, which are associated with shifts in breeding systems towards dioecy. A single origin of dioecy from hermaphroditism is supported. Finally, we report evolutionary changes from headed to non-headed inflorescences and multiple to simple fruits.     

Molecular phylogenetics of squirrelfishes and soldierfishes (Teleostei: Beryciformes: Holocentridae): Reconciling more than 100 years of taxonomic confusion

Squirrelfishes and soldierfishes (Holocentridae) are among the most conspicuous species in the nocturnal reef fish community. However, there is no clear consensus regarding their evolutionary relationships, which is reflected in a complicated taxonomic history. We collected DNA sequence data from multiple single copy nuclear genes and one mitochondrial gene sampled from over fifty percent of the recognized holocentrid species and infer the first species-level phylogeny of the Holocentridae. Our results strongly support the monophyly of the clades Myripristinae (soldierfishes) and Holocentrinae (squirrelfishes). The molecular phylogenies differ with regard to previous hypotheses of relationships within the Myriprisitinae, resolving a clade of cryptic reef associated and deep water non-reef dwelling lineages (Corniger+Plectrypops+Ostichthys) that is the sister lineage to a monophyletic Myripristis. Within Holocentrinae, Neoniphon and Sargocentron are strongly supported as paraphyletic, while Holocentrus is nested within Sargocentron. Using Bayesian ancestral state reconstruction methods, we demonstrate the taxonomically diagnostic characters for Neoniphon and Sargocentron likely represent character states with a complex evolutionary history that is not reflective of shared common ancestry. We propose a new classification for Holocentrinae, recognizing four lineages that are treated as genera: SargocentronFowler, 1904, HolocentrusScopoli, 1777, FlameoJordan and Evermann, 1898, and NeoniphonCastelnau, 1875.     

Historical Biogeography of endemic seed plant genera in the Caribbean: Did GAARlandia play a role?

The Caribbean archipelago is a region with an extremely complex geological history and an outstanding plant diversity with high levels of endemism. The aim of this study was to better understand the historical assembly and evolution of endemic seed plant genera in the Caribbean, by first determining divergence times of endemic genera to test whether the hypothesized Greater Antilles and Aves Ridge (GAARlandia) land bridge played a role in the archipelago colonization and second by testing South America as the main colonization source as expected by the position of landmasses and recent evidence of an asymmetrical biotic interchange. We reconstructed a dated molecular phylogenetic tree for 625 seed plants including 32 Caribbean endemic genera using Bayesian inference and ten calibrations. To estimate the geographic range of the ancestors of endemic genera, we performed a model selection between a null and two complex biogeographic models that included timeframes based on geological information, dispersal probabilities, and directionality among regions. Crown ages for endemic genera ranged from Early Eocene (53.1 Ma) to Late Pliocene (3.4 Ma). Confidence intervals for divergence times (crown and/or stem ages) of 22 endemic genera occurred within the GAARlandia time frame. Contrary to expectations, the Antilles appears as the main ancestral area for endemic seed plant genera and only five genera had a South American origin. In contrast to patterns shown for vertebrates and other organisms and based on our sampling, we conclude that GAARlandia did not act as a colonization route for plants between South America and the Antilles. Further studies on Caribbean plant dispersal at the species and population levels will be required to reveal finer‐scale biogeographic patterns and mechanisms.     

Rate heterogeneity, ancestral character state reconstruction, and the evolution of limb morphology in Lerista (Scincidae, Squamata)

Rates of phenotypic evolution derive from numerous interrelated processes acting at varying spatial and temporal scales and frequently differ substantially among lineages. Although current models employed in reconstructing ancestral character states permit independent rates for distinct types of transition (forward and reverse transitions and transitions between different states), these rates are typically assumed to be identical for all branches in a phylogeny. In this paper, I present a general model of character evolution enabling rate heterogeneity among branches. This model is employed in assessing the extent to which the assumption of uniform transition rates affects reconstructions of ancestral limb morphology in the scincid lizard clade Lerista and, accordingly, the potential for rate variability to mislead inferences of evolutionary patterns. Permitting rate variation among branches significantly improves model fit for both the manus and the pes. A constrained model in which the rate of digit acquisition is assumed to be effectively zero is strongly supported in each case; when compared with a model assuming unconstrained transition rates, this model provides a substantially better fit for the manus and a nearly identical fit for the pes. Ancestral states reconstructed assuming the constrained model imply patterns of limb evolution differing significantly from those implied by reconstructions for uniform-rate models, particularly for the pes; whereas ancestral states for the uniform-rate models consistently entail the reacquisition of pedal digits, those for the model incorporating among-lineage rate heterogeneity imply repeated, unreversed digit loss. These results indicate that the assumption of identical transition rates for all branches in a phylogeny may be inappropriate in modeling the evolution of phenotypic traits and emphasize the need for careful evaluation of phylogenetic tests of Dollo's law.     

Evolutionary model of protein secondary structure capable of revealing new biological relationships

Ancestral sequence reconstruction has had recent success in decoding the origins and the determinants of complex protein functions. However, phylogenetic analyses of remote homologues must handle extreme amino acid sequence diversity resulting from extended periods of evolutionary change. We exploited the wealth of protein structures to develop an evolutionary model based on protein secondary structure. The approach follows the differences between discrete secondary structure states observed in modern proteins and those hypothesized in their immediate ancestors. We implemented maximum likelihood-based phylogenetic inference to reconstruct ancestral secondary structure. The predictive accuracy from the use of the evolutionary model surpasses that of comparative modeling and sequence-based prediction; the reconstruction extracts information not available from modern structures or the ancestral sequences alone. Based on a phylogenetic analysis of a sequence-diverse protein family, we showed that the model can highlight relationships that are evolutionarily rooted in structure and not evident in amino acid-based analysis.     

Functional divergence in the Arabidopsis beta-1,3-glucanase gene family inferred by phylogenetic reconstruction of expression states

Plant beta-1,3-glucanases (beta-1,3-Gs) (E.C. 3.2.1.39) comprise large, highly complex gene families involved in pathogen defense as well as a wide range of normal developmental processes. In spite of previous phylogenetic analyses that classify beta-1,3-Gs by sequence relatedness, the functional evolution of beta-1,3-Gs remains unclear. Here, expression and phylogenetic analyses have been integrated in order to investigate patterns of functional divergence in the Arabidopsis beta-1,3-G gene family. Fifty beta-1,3-G genes were grouped into expression classes through clustering of microarray data, and functions were inferred based on knowledge of coexpressed genes and existing literature. The resulting expression classes were mapped as discrete states onto a phylogenetic tree and parsimony reconstruction of ancestral expression states was performed, providing a model of expression divergence. Results showed a highly nonrandom distribution of developmental expression states in the phylogeny (P = 0.0002) indicating a significant degree of coupling between sequence and developmental expression divergence. A weaker, yet significant level of coupling was found using stress response data, but not using hormone-response or pathogen-response data. According to the model of developmental expression divergence, the ancestral function was most likely involved in cell division and/or cell wall remodeling. The associated expression state is widely distributed in the phylogeny, is retained by over 25% of gene family members, and is consistent with the known functions of beta-1,3-Gs in distantly related species and gene families. Consistent with previous hypotheses, pathogenesis-related (PR) beta-1,3-Gs appear to have evolved from ancestral developmentally regulated beta-1,3-Gs, acquiring PR function through a number of evolutionary events: divergence from the ancestral expression state, acquisition of pathogen/stress-responsive expression patterns, and loss of the C-terminal region including the glycosylphosphatidylinisotol (GPI)-anchoring site thus allowing for extracellular secretion.