Modelling the ecological-functional diversification of marine Metazoa on geological time scales

The ecological traits and functional capabilities of marine animals have changed significantly since their origin in the late Precambrian. These changes can be analysed quantitatively using multi-dimensional parameter spaces in which the ecological lifestyles of species are represented by particular combinations of parameter values. Here, we present models that describe the filling of this multi-dimensional 'ecospace' by ecological lifestyles during metazoan diversification. These models reflect varying assumptions about the processes that drove ecological diversification; they contrast diffusive expansion with driven expansion and niche conservatism with niche partitioning. Some models highlight the importance of interactions among organisms (ecosystem engineering and predator-prey escalation) in promoting new lifestyles or eliminating existing ones. These models reflect processes that were not mutually exclusive; rigorous analyses will continue to reveal their applicability to episodes in metazoan history.     

Patterns of molecular evolution and diversification in a biodiversity hotspot: the California Floristic Province

The California Floristic Province harbours more endemic plant and animal taxa and more identifiable subspecies than any other area of comparable size in North America. We present evidence that physical historical processes have resulted in congruent patterns of genetic diversity over the past 2-10 million years. Using a molecular clock approach we show that diversification and establishment of spatial genetic structure across six taxonomic groups coincide with the putative age of California's mountain ranges and aridification in the region. Our results demonstrate the importance of geographical barriers and climatological events to species diversification and the overall geographical structure of biodiversity. These results should facilitate conservation efforts in this biodiversity hotspot for taxa whose population genetic structure is still unknown and may suggest the potential utility of this approach in regional conservation planning efforts.     

The origin of Metazoa: a transition from temporal to spatial cell differentiation

For over a century, Haeckel's Gastraea theory remained a dominant theory to explain the origin of multicellular animals. According to this theory, the animal ancestor was a blastula‐like colony of uniform cells that gradually evolved cell differentiation. Today, however, genes that typically control metazoan development, cell differentiation, cell‐to‐cell adhesion, and cell‐to‐matrix adhesion are found in various unicellular relatives of the Metazoa, which suggests the origin of the genetic programs of cell differentiation and adhesion in the root of the Opisthokonta. Multicellular stages occurring in the complex life cycles of opisthokont protists (mesomycetozoeans and choanoflagellates) never resemble a blastula. Here, we discuss a more realistic scenario of transition to multicellularity through integration of pre‐existing transient cell types into the body of an early metazoon, which possessed a complex life cycle with a differentiated sedentary filter‐feeding trophic stage and a non‐feeding blastula‐like larva, the synzoospore. Choanoflagellates are considered as forms with secondarily simplified life cycles.     

The evolution of South American endemic canids: a history of rapid diversification and morphological parallelism

The origin of endemic South American canid fauna has been traditionally linked with the rise of the Isthmus of Panama, suggesting that diversification of the dog fauna on this continent occurred very rapidly. Nevertheless, despite its obvious biogeographic appeal, the tempo of Canid evolution in South America has never been studied thoroughly. This issue can be suitably tackled with the inference of a molecular timescale. In this study, using a relaxed molecular clock method, we estimated that the most recent common ancestor of South American canids lived around 4 Ma, whereas all other splits within the clade occurred after the rise of the Panamanian land bridge. We suggest that the early diversification of the ancestors of the two main lineages of South American canids may have occurred in North America, before the Great American Interchange. Moreover, a concatenated morphological and molecular analysis put some extinct canid species well within the South American radiation, and shows that the dental adaptations to hypercarnivory evolved only once in the South American clade.     

Molecular evidence for the diversification of extant lichens in the late cretaceous and tertiary

A molecular clock based on ITS sequence data from the lichen genera Biatora and Phyllopsora is calibrated with the help of paleoclimatic data and evidence of forest history. The clock indicates that diversification within Biatora started as early as in the Late Cretaceous and took place during periods of climatic cooling, when new types of forest evolved and spread in the Northern Hemisphere. Arctic-alpine species of the genus appear to be of considerable age, dating back to the Late Eocene-Oligocene climatic cooling. By using calibrated phylogenies of epiphytic lichens it may become possible to date many paleoenvironmental events, for which little fossil evidence exists.     

Origin and rapid diversification of a tropical moss

Molecular sequences rarely evolve at a constant rate. Yet, even in instances where a clock can be assumed or approximated for a particular set of sequences, fossils or clear patterns of vicariance are rarely available to calibrate the clock. Thus, obtaining absolute timing for diversification of natural lineages can prove difficult. Unfortunately, without absolute time we cannot develop a complete understanding of important evolutionary processes, including adaptive radiations and key innovations. In the present study, the coding sequence of the nuclear gene, glyceraldehyde 3-phosphate dehydrogenase (gpd), extracted from the paleotropical moss, Mitthyridium, was found to exhibit clocklike behavior and used to reconstruct the history of 80 distinct molecular lineages that cover the full geographic range of Mitthyridium. Two separate clades endemic to two geographically distinct oceanic archipelagos were revealed by this phylogenetic analysis. This allowed the use of island age (as derived from potassium-argon dating) as a maximum age of origin of each monophyletic group, providing two independent time anchors for the clock found in gpd, the final piece needed to study absolute time. Based on results from both maximum age calibrations, which separately yielded highly consistent estimates, the ancestor of this moss group arose approximately 8 million years ago, and then diversified at the rapid rate of 0.56 +/- 0.004 new lineages per million years. Such a rate is on par with the highest diversification rates reported in the literature including rapidly radiating insular groups like the Hawaiian silversword alliance, a classic example of an adaptive radiation. Using independent sources of data, it was found that neither the age nor diversification estimates were affected by the use of molecular lineages rather than species as the operational taxonomic units. Identifying the cause for this rapid diversification requires further testing, but it appears to be related to a general shift in reproductive strategy from sexual to asexual, which may be a key innovation for this young group.     

Assessing the effects of a sequestered germline on interdomain lateral gene transfer in Metazoa

A sequestered germline in Metazoa has been argued to be an obstacle to lateral gene transfer (LGT), though few studies have specifically assessed this claim. Here, we test the hypothesis that the origin of a sequestered germline reduced LGT events in Bilateria (i.e., triploblast lineages) as compared to early‐diverging Metazoa (i.e., Ctenophora, Cnidaria, Porifera, and Placozoa). We analyze single‐gene phylogenies generated with over 900 species sampled from among Bacteria, Archaea, and Eukaryota to identify well‐supported interdomain LGTs. We focus on ancient interdomain LGT (i.e., those between prokaryotes and multiple lineages of Metazoa) as systematic errors in single‐gene tree reconstruction create uncertainties for interpreting eukaryote‐to‐eukaryote transfer. The breadth of the sampled Metazoa enables us to estimate the timing of LGTs, and to examine the pattern before versus after the evolution of a sequestered germline. We identified 58 LGTs found only in Metazoa and prokaryotes (i.e., bacteria and/or archaea), and seven genes transferred from prokaryotes into Metazoa plus one other eukaryotic clade. Our analyses indicate that more interdomain transfers occurred before the development of a sequestered germline, consistent with the hypothesis that this feature is an obstacle to LGT.     

A general comparison of relaxed molecular clock models

Several models have been proposed to relax the molecular clock in order to estimate divergence times. However, it is unclear which model has the best fit to real data and should therefore be used to perform molecular dating. In particular, we do not know whether rate autocorrelation should be considered or which prior on divergence times should be used. In this work, we propose a general bench mark of alternative relaxed clock models. We have reimplemented most of the already existing models, including the popular lognormal model, as well as various prior choices for divergence times (birth-death, Dirichlet, uniform), in a common Bayesian statistical framework. We also propose a new autocorrelated model, called the "CIR" process, with well-defined stationary properties. We assess the relative fitness of these models and priors, when applied to 3 different protein data sets from eukaryotes, vertebrates, and mammals, by computing Bayes factors using a numerical method called thermodynamic integration. We find that the 2 autocorrelated models, CIR and lognormal, have a similar fit and clearly outperform uncorrelated models on all 3 data sets. In contrast, the optimal choice for the divergence time prior is more dependent on the data investigated. Altogether, our results provide useful guidelines for model choice in the field of molecular dating while opening the way to more extensive model comparisons.     

Age estimates of Frullania (Frullaniaceae,Porellales) main lineages: another example of rapid and recent diversification in liverwort evolution

Frullania Raddi is an extant genus of liverworts (Bryophytes) widespread around the world. It belongs to the family Frullaniaceae Lorch., with a large number of species distributed into several subgenera, sections and subsections according with different morphological classifications. As shown in previous studies, most Frullania sub-generic classifications are supported by molecular data, indicating that morphological characters appear well suitable to discriminate species. However, deep among-clade relationships remain unclear, probably due to the rapid diversification of main clades, paralleling other plant lineages. In this study, we reconstruct Frullania phylogeny based on available molecular data used in previous studies, and we present time estimates for the origin of its main branches. Results supported the monophyly of most subgenera as demonstrated in previous studies and supported a rapid diversification of these main lineages. Time estimates under a relaxed molecular clock and with integrated fossil evidence and nucleotide mutation rates further suggested a Jurassic origin of the genus and a rapid diversification during Palaeogene and Neogene. This may have been influenced by geographical and climate changes during these periods as predicted for most leafy liverworts.     

Revisiting the phylogeography of Asellus aquaticus in Europe: insights into cryptic diversity and spatiotemporal diversification

相似文献   

Stochastic models of molecular evolution and the estimation of phylogeny and rates of nucleotide substitution in the hominoid primates

We have estimated phylogenetic patterns and rates of nucleotide substitution in the hominoid primates using two different probabilistic models of molecular evolution as applied to three different data sets of nucleic acid sequences. The orang-utan was found to be the out-group of the other hominoids examined. Within the African apes and human clade the sister-group relationship of chimpanzee and human was found to be statistically the best, although the magnitude of the error estimates (a reflection of random statistical fluctuations) makes this conclusion tentative. The ψν-globin data sets were found to be statistically the most consistent and gave estimates of the times of divergence of chimpanzee and human from gorilla and of chimpanzee from human as 7·7 ± 1·5 Ma (Millions of years ago) and 7·4 ± 1·5 Ma respectively, although the speculative nature of these estimates is emphasized. In all cases the calibration point was the assumed divergence of the orang-utan from the remaining hominoids at 14·5 Ma. There was no statistically significant evidence of a slowdown in nucleotide substitution rate for the human lineage, or among the hominoids as a whole with respect to the Old and New World monkeys. We advocate the continued use and development of stochastic models of molecular evolution as a basis for phylogenetic estimation. On this basis one can choose between competing hypotheses of relationship in a statistical manner and can provide estimates of the errors involved in such estimations. The assumptions of all stochastic models are open to test and future refinement.     

Evolutionary origins of germline segregation in Metazoa: evidence for a germ stem cell lineage in the coral Orbicella faveolata (Cnidaria,Anthozoa)

The ability to segregate a committed germ stem cell (GSC) lineage distinct from somatic cell lineages is a characteristic of bilaterian Metazoans. However, the occurrence of GSC lineage specification in basally branching Metazoan phyla, such as Cnidaria, is uncertain. Without an independently segregated GSC lineage, germ cells and their precursors must be specified throughout adulthood from continuously dividing somatic stem cells, generating the risk of propagating somatic mutations within the individual and its gametes. To address the potential for existence of a GSC lineage in Anthozoa, the sister-group to all remaining Cnidaria, we identified moderate- to high-frequency somatic mutations and their potential for gametic transfer in the long-lived coral Orbicella faveolata (Anthozoa, Cnidaria) using a 2b-RAD sequencing approach. Our results demonstrate that somatic mutations can drift to high frequencies (up to 50%) and can also generate substantial intracolonial genetic diversity. However, these somatic mutations are not transferable to gametes, signifying the potential for an independently segregated GSC lineage in O. faveolata. In conjunction with previous research on germ cell development in other basally branching Metazoan species, our results suggest that the GSC system may be a Eumetazoan characteristic that evolved in association with the emergence of greater complexity in animal body plan organization and greater specificity of stem cell functions.     

The local-clock permutation test: a simple test to compare rates of molecular evolution on phylogenetic trees

Rates of molecular evolution vary substantially between lineages, and a growing effort is directed at uncovering the causes and consequences of this variation. Comparing local-clocks (rates of molecular evolution estimated from different sets of branches of a phylogenetic tree) is a common tool in this research effort. Here, I show that a commonly used test (the Likelihood Ratio Test, LRT) will not be statistically valid for comparing local-clocks in most cases. Instead, I propose the local-clock permutation test (LCPT), a simple test that can be used to test the significance of differences between local-clocks. The LCPT could also be used to test for differences between any parameter that can be assigned to individual branches on a phylogenetic tree. Using simulated data, I show that the LCPT has good power to detect differences between local-clocks.     

Bayesian predictive information criterion for the evaluation of hierarchical Bayesian and empirical Bayes models

The problem of evaluating the goodness of the predictive distributionsof hierarchical Bayesian and empirical Bayes models is investigated.A Bayesian predictive information criterion is proposed as anestimator of the posterior mean of the expected loglikelihoodof the predictive distribution when the specified family ofprobability distributions does not contain the true distribution.The proposed criterion is developed by correcting the asymptoticbias of the posterior mean of the loglikelihood as an estimatorof its expected loglikelihood. In the evaluation of hierarchicalBayesian models with random effects, regardless of our parametricfocus, the proposed criterion considers the bias correctionof the posterior mean of the marginal loglikelihood becauseit requires a consistent parameter estimator. The use of thebootstrap in model evaluation is also discussed.     

Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes

With about 24,000 extant species, teleosts are the largest group of vertebrates. They constitute more than 99% of the ray-finned fishes (Actinopterygii) that diverged from the lobe-finned fish lineage (Sarcopterygii) about 450 MYA. Although the role of genome duplication in the evolution of vertebrates is now established, its role in structuring the teleost genomes has been controversial. At least two hypotheses have been proposed: a whole-genome duplication in an ancient ray-finned fish and independent gene duplications in different lineages. These hypotheses are, however, based on small data sets and lack adequate statistical and phylogenetic support. In this study, we have made a systematic comparison of the draft genome sequences of Fugu and humans to identify paralogous chromosomal regions ("paralogons") in the Fugu that arose in the ray-finned fish lineage ("fish-specific"). We identified duplicate genes in the Fugu by phylogenetic analyses of the Fugu, human, and invertebrate sequences. Our analyses provide evidence for 425 fish-specific duplicate genes in the Fugu and show that at least 6.6% of the genome is represented by fish-specific paralogons. We estimated the ages of Fugu duplicate genes and paralogons using the molecular clock. Remarkably, the ages of duplicate genes and paralogons are clustered, with a peak around 350 MYA. These data strongly suggest a whole-genome duplication event early during the evolution of ray-finned fishes, probably before the origin of teleosts.     

Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds

We implement a Bayesian Markov chain Monte Carlo algorithm for estimating species divergence times that uses heterogeneous data from multiple gene loci and accommodates multiple fossil calibration nodes. A birth-death process with species sampling is used to specify a prior for divergence times, which allows easy assessment of the effects of that prior on posterior time estimates. We propose a new approach for specifying calibration points on the phylogeny, which allows the use of arbitrary and flexible statistical distributions to describe uncertainties in fossil dates. In particular, we use soft bounds, so that the probability that the true divergence time is outside the bounds is small but nonzero. A strict molecular clock is assumed in the current implementation, although this assumption may be relaxed. We apply our new algorithm to two data sets concerning divergences of several primate species, to examine the effects of the substitution model and of the prior for divergence times on Bayesian time estimation. We also conduct computer simulation to examine the differences between soft and hard bounds. We demonstrate that divergence time estimation is intrinsically hampered by uncertainties in fossil calibrations, and the error in Bayesian time estimates will not go to zero with increased amounts of sequence data. Our analyses of both real and simulated data demonstrate potentially large differences between divergence time estimates obtained using soft versus hard bounds and a general superiority of soft bounds. Our main findings are as follows. (1) When the fossils are consistent with each other and with the molecular data, and the posterior time estimates are well within the prior bounds, soft and hard bounds produce similar results. (2) When the fossils are in conflict with each other or with the molecules, soft and hard bounds behave very differently; soft bounds allow sequence data to correct poor calibrations, while poor hard bounds are impossible to overcome by any amount of data. (3) Soft bounds eliminate the need for "safe" but unrealistically high upper bounds, which may bias posterior time estimates. (4) Soft bounds allow more reliable assessment of estimation errors, while hard bounds generate misleadingly high precisions when fossils and molecules are in conflict.     

Molecules, development and fossils in the study of metazoan evolution; Articulata versus Ecdysozoa revisited

Two conflicting hypotheses of protostome relationships, Articulata and Ecdysozoa, are reviewed by evaluating the evidence in favor and against each one of them. Understanding early embryonic development and segmentation in non-arthropod non-annelid protostomes seems crucial to the debate. New ways of coding metazoan matrices, avoiding ground-patterns and higher taxa, and incorporating fossil evidence seems the best way to avoid circular debates. Molecular data served as the catalyzer for the Ecdysozoa hypothesis, although morphological support had been implicitly suggested. Most molecular analyses published so far have shown some support for Ecdysozoa, whereas none has ever supported Articulata. Here, new analyses of up to four nuclear loci, including 18S rRNA, myosin heavy chain II, histone H3 and elongation factor 1- are conducted to test the molecular support for Ecdysozoa, and, at least under some parameter sets, most data sets show a clade formed by the molting animals. In contrast, support for Articulata is not found under any analytical conditions.     

Differentiation trees,a junk DNA molecular clock,and the evolution of neoteny in salamanders

Obligate neotenic salamanders die if forced to metamorphose. We suggest that this can be explained by assuming: 1) their “excess” DNA is “junk” DNA; 2) the “adult” specifying portion of the DNA becomes junk DNA and is available for repeated duplication. This suggests a “new” junk DNA molecular clock. We obtain remarkable agreement in “predicting” the amount of DNA per nucleus in present day non-obligate neotene salamanders from this molecular clock. These observatons are consistent with the idea that the development of these animals is describable in terms of differentiation trees whose branches (gene cascades) corresponding to adult somatic tissues accumulate deleterious mutations over evolutionary time. We show that the amount of DNA per nucleus increases linearly with the phylogenetic age of salamander families. The lack of constraints by natural selection, on unused adult branches, may account for the large amount of so-called “junk DNA” in obligate neotenic salamanders. The effects of this excess DNA, via increased cell size, suggest a positive feedback, ecophysiological explanation for such junk DNA: adaptation to cool water environments is enhanced by the lower metabolism associated with more DNA, larger cells and slower developmental time.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

Biogeography and evolution of the Galapagos: integration of the biological and geological evidence

Biogeographic tracks are mapped for Galapagos endemics representing 25 plant and animal taxa and including organisms with good and poor means of dispersal. These patterns confirm standard biogeographic tracks linking Galapagos with Central America, western North and South America, the Caribbean, Asia and Australasia. Discovery of the Galapagos Gore in the 1970s corroborates the biogeographic prediction for a major tectonic centre associated with the Galapagos. The biogeographic model developed by Croizat in 1958 of Galapagos colonization involving an ancestral biota inhabiting eastern Pacific geosynclinal forelands is congruent with plate tectonic models supporting a Pacific island arc origin for western American terranes. American relatives of Galapagos endemics may have originated within an eastern Pacific paleogeography rather than representing centres of origin for dispersal to the Galapagos. Galapagos colonization by an eastern Pacific biota between late Cretaceous and mid-Tertiary has significant implications for understanding the tempo and mode for both the origins of island biota and general models of evolutionary differentiation. Popular assertions that overwater dispersal represents the only viable origin for the entire Galapagos biota is no longer biogeographically or geologically tenable.