Age of cichlids: new dates for ancient lake fish radiations

Timing divergence events allow us to infer the conditions under which biodiversity has evolved and gain important insights into the mechanisms driving evolution. Cichlid fishes are a model system for studying speciation and adaptive radiation, yet, we have lacked reliable timescales for their evolution. Phylogenetic reconstructions are consistent with cichlid origins prior to Gondwanan landmass fragmentation 121-165 MYA, considerably earlier than the first known fossil cichlids (Eocene). We examined the timing of cichlid evolution using a relaxed molecular clock calibrated with geological estimates for the ages of 1) Gondwanan fragmentation and 2) cichlid fossils. Timescales of cichlid evolution derived from fossil-dated phylogenies of other bony fishes most closely matched those suggested by Gondwanan breakup calibrations, suggesting the Eocene origins and marine dispersal implied by the cichlid fossil record may be due to its incompleteness. Using Gondwanan calibrations, we found accumulation of genetic diversity within the radiating lineages of the African Lakes Malawi, Victoria and Barombi Mbo, and Palaeolake Makgadikgadi began around or after the time of lake basin formation. These calibrations also suggest Lake Tanganyika was colonized independently by the major radiating cichlid tribes that then began to accumulate genetic diversity thereafter. These results contrast with the widely accepted theory that diversification into major lineages took place within the Tanganyika basin. Together, this evidence suggests that ancient lake habitats have played a key role in generating and maintaining diversity within radiating lineages and also that lakes may have captured preexisting cichlid diversity from multiple sources from which adaptive radiations have evolved.

The Accumulation of Reproductive Incompatibilities in African Cichlid Fish

The rate at which different components of reproductive isolation accumulate with divergence time between species has only been studied in a limited, but growing, number of species. We measured premating isolation and hybrid inviability at four different ontogenetic stages from zygotes to adults in interspecific hybrids of 26 pairs of African cichlid species, spanning the entire East African haplochromine radiation. We then used multiple relaxed molecular clock calibrations to translate genetic distances into absolute ages to compare evolutionary rates of different components of reproductive isolation. We find that premating isolation accumulates fast initially but then changes little with increasing genetic distance between species. In contrast, postmating isolation between closely related species is negligible but then accumulates rapidly, resulting in complete hybrid inviability after 4.4/8.5/18.4 million years (my). Thus, the rate at which complete intrinsic incompatibilities arise in this system is orders of magnitude lower than rates of speciation within individual lake radiations. Together these results suggest divergent ecological adaptations may prevent populations from interbreeding and help maintain cichlid species diversity, which may be vulnerable to environmental degradation. By quantifying the capacity to produce viable hybrids between allopatric, distantly related lineages our results also provide an upper divergence time limit for the "hybrid swarm origin" model of adaptive radiation.

Host jumps shaped the diversity of extant rust fungi (Pucciniales)

* The aim of this study was to determine the evolutionary time line for rust fungi and date key speciation events using a molecular clock. Evidence is provided that supports a contemporary view for a recent origin of rust fungi, with a common ancestor on a flowering plant. * Divergence times for > 20 genera of rust fungi were studied with Bayesian evolutionary analyses. A relaxed molecular clock was applied to ribosomal and mitochondrial genes, calibrated against estimated divergence times for the hosts of rust fungi, such as Acacia (Fabaceae), angiosperms and the cupressophytes. * Results showed that rust fungi shared a most recent common ancestor with a mean age between 113 and 115 million yr. This dates rust fungi to the Cretaceous period, which is much younger than previous estimations. Host jumps, whether taxonomically large or between host genera in the same family, most probably shaped the diversity of rust genera. Likewise, species diversified by host shifts (through coevolution) or via subsequent host jumps. This is in contrast to strict coevolution with their hosts. * Puccinia psidii was recovered in Sphaerophragmiaceae, a family distinct from Raveneliaceae, which were regarded as confamilial in previous studies.

Accuracy of rate estimation using relaxed-clock models with a critical focus on the early metazoan radiation

In recent years, a number of phylogenetic methods have been developed for estimating molecular rates and divergence dates under models that relax the molecular clock constraint by allowing rate change throughout the tree. These methods are being used with increasing frequency, but there have been few studies into their accuracy. We tested the accuracy of several relaxed-clock methods (penalized likelihood and Bayesian inference using various models of rate change) using nucleotide sequences simulated on a nine-taxon tree. When the sequences evolved with a constant rate, the methods were able to infer rates accurately, but estimates were more precise when a molecular clock was assumed. When the sequences evolved under a model of autocorrelated rate change, rates were accurately estimated using penalized likelihood and by Bayesian inference using lognormal and exponential models of rate change, while other models did not perform as well. When the sequences evolved under a model of uncorrelated rate change, only Bayesian inference using an exponential rate model performed well. Collectively, the results provide a strong recommendation for using the exponential model of rate change if a conservative approach to divergence time estimation is required. A case study is presented in which we use a simulation-based approach to examine the hypothesis of elevated rates in the Cambrian period, and it is found that these high rate estimates might be an artifact of the rate estimation method. If this bias is present, then the ages of metazoan divergences would be systematically underestimated. The results of this study have implications for studies of molecular rates and divergence dates.

Relaxed phylogenetics and dating with confidence

In phylogenetics, the unrooted model of phylogeny and the strict molecular clock model are two extremes of a continuum. Despite their dominance in phylogenetic inference, it is evident that both are biologically unrealistic and that the real evolutionary process lies between these two extremes. Fortunately, intermediate models employing relaxed molecular clocks have been described. These models open the gate to a new field of “relaxed phylogenetics.” Here we introduce a new approach to performing relaxed phylogenetic analysis. We describe how it can be used to estimate phylogenies and divergence times in the face of uncertainty in evolutionary rates and calibration times. Our approach also provides a means for measuring the clocklikeness of datasets and comparing this measure between different genes and phylogenies. We find no significant rate autocorrelation among branches in three large datasets, suggesting that autocorrelated models are not necessarily suitable for these data. In addition, we place these datasets on the continuum of clocklikeness between a strict molecular clock and the alternative unrooted extreme. Finally, we present analyses of 102 bacterial, 106 yeast, 61 plant, 99 metazoan, and 500 primate alignments. From these we conclude that our method is phylogenetically more accurate and precise than the traditional unrooted model while adding the ability to infer a timescale to evolution.

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Review ["Telling the Evolutionary Time : Molecular Clocks and the Fossil Record" edited by Philip C. J. Donoghue and M. Paul Smith]

Determining the temporal scale of biological evolution has traditionally been the preserve of paleontology, with the timing of species originations and major diversifications all being read from the fossil record. However, the ages of the earliest (correctly identified) records will underestimate actual origins due to the incomplete nature of the fossil record and the necessity for lineages to have evolved sufficiently divergent morphologies in order to be distinguished. The possibility of inferring divergence times more accurately has been promoted by the idea that the accumulation of genetic change between modern lineages can be used as a molecular clock (Zuckerkandl and Pauling, 1965). In practice, though, molecular dates have often been so old as to be incongruent even with liberal readings of the fossil record. Prominent examples include inferred diversifications of metazoan phyla hundreds of millions of years before their Cambrian fossil record appearances (e.g., Nei et al., 2001) and a basal split between modern birds (Neoaves) that is almost double the age of their earliest recognizable fossils (e.g., Cooper and Penny, 1997).

Accounting for calibration uncertainty in phylogenetic estimation of evolutionary divergence times

The estimation of phylogenetic divergence times from sequence data is an important component of many molecular evolutionary studies. There is now a general appreciation that the procedure of divergence dating is considerably more complex than that initially described in the 1960s by Zuckerkandl and Pauling (1962, 1965). In particular, there has been much critical attention toward the assumption of a global molecular clock, resulting in the development of increasingly sophisticated techniques for inferring divergence times from sequence data. In response to the documentation of widespread departures from clocklike behavior, a variety of local- and relaxed-clock methods have been proposed and implemented. Local-clock methods permit different molecular clocks in different parts of the phylogenetic tree, thereby retaining the advantages of the classical molecular clock while casting off the restrictive assumption of a single, global rate of substitution (Rambaut and Bromham 1998; Yoder and Yang 2000).

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25-50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events. © 2005 Elsevier B.V. All rights reserved.

Molecular phylogeny and taxonomy of wood mice (genus Apodemus Kaup, 1829) based on complete mtDNA cytochrome b sequences, with emphasis on Chinese species

Phylogenetic relationships among 15 species of wood mice (genus Apodemus) were reconstructed to explore some long-standing taxonomic problems. The results provided support for the monophyly of the genus Apodemus, but could not reject the hypothesis of paraphyly for this genus. Our data divided the 15 species into four major groups: (1) the Sylvaemus group (A. sylvaticus, A. flavicollis, A. alpicola, and A. uralensis), (2) the Apodemus group (A. peninsulae, A. chevreri, A. agrarius, A. speciosus, A. draco, A. ilex, A. semotus, A. latronum, and A. mystacinus), (3) A. argenteus, and (4) A. gurkha. Our results also suggested that orestes should be a valid subspecies of A. draco rather than an independent species; in contrast, A. ilex from Yunnan may be regarded as a separate species rather than a synonym of orestes or draco. The species level status of A. latronum, tscherga as synonyms of A. uralensis, and A. chevrieri as a valid species and the closest sibling species of A. agrarius were further corroborated by our data. Applying a molecular clock with the divergences of Mus and Rattus set at 12 million years ago (Mya) as a calibration point, it was estimated that five old lineages (A. mystacinus and four major groups above) diverged in the late Miocene (7.82-12.74 Mya). Then the Apodemus group (excluding A. mystacinus) split into two subgroups: agrarius and draco, at about 7.17-9.95 Mya. Four species of the Sylvaemus group were estimated to diverge at about 2.92-5.21 Mya. The Hengduan Mountains Region was hypothesized to have played important roles in Apodemus evolutionary histories since the Pleistocene. (C) 2004 Elsevier Inc. All rights reserved.

Molecular phylogeny of the specialized schizothoracine fishes (Teleostei : Cyprinidae), with their implications for the uplift of the Qinghai-Tibetan Plateau

Molecular phylogeny of three genera containing nine species and subspecies of the specialized schizothoracine fishes are investigated based on the complete nucleotide sequence of mitochondrial cytochrome b gene. Meantime relationships between the main cladogenetic events of the specialized schizothoracine fishes and the stepwise uplift of the Qinghai-Tibetan Plateau are also conducted using the molecular clock, which is calibrated by geological isolated events between the upper reaches of the Yellow River and the Qinghai Lake. Results indicated that the specialized schizothoracine fishes are not a monophyly. Five species and subspecies of Ptychobarbus form a monophyly. But three species of Gymnodiptychus do not form a monophyly. Gd. integrigymnatus is a sister taxon of the highly specialized schizothoracine fishes while Gd. pachycheilus has a close relation with Gd. dybowskii, and both of them are as a sister group of Diptychus maculatus. The specialized schizothoracines fishes might have originated during the Miocene (about 10 MaBP), and then the divergence of three genera happened during late Miocene (about 8 MaBP). Their main specialization occurred during the late Pliocene and Pleistocene (3.54-0.42 MaBP). The main cladogenetic events of the specialized schizothoracine fishes are mostly correlated with the geological tectonic events and intensive climate shift happened at 8, 3.6, 2.5 and 1.7 MaBP of the late Cenozoic. Molecular clock data do not support the hypothesis that the Qinghai-Tibetan Plateau uplifted to near present or even higher elevations during the Oligocene or Miocene, and neither in agreement with the view that the plateau uplifting reached only to an altitude of 2000 in during the late Pliocene (about 2.6 MaBP).

Multiple in-to-Africa dispersals of labeonin fishes (Teleostei: Cyprinidae) revealed by molecular phylogenetic analysis

An out-of-Africa dispersal route has been proposed for many organisms, including modern man. However, counter examples of in-to-Africa dispersal routes are less common. In the present article, the phylogenetic relationships within the Labeoninae, a subfamily of cyprinid fishes distributed in Asia and Africa, were analyzed to investigate the biogeographic processes governing the modern distribution of these Asian and African cyprinids. The mitochondrial DNA cytochrome b gene was used as a molecular marker. The phylogenetic analysis indicated that the subfamily Labeoninae is a monophyletic group, with some Asian labeonins located at the basal position. Two subclades were found that contained both African and Asian species, which highlighted a need for further biogeographic analysis. Based on this analysis, it is proposed that the centre of origin of the Labeoninae was in East Asia. Molecular clock estimation suggests that the Labeoninae arose by the Early Miocene (similar to 23 MYA) during the period of the second Tibetan uplift. Subsequently, two dispersal events of labeonins from Asia into Africa occured in the Early Miocene (similar to 20 MYA) and Late Miocene (similar to 9 MYA) and serve as examples counter to out-of-Africa dispersal.

Molecular phylogeny of the specialized schizothoracine fishes (Teleostei : Cyprinidae), with their implications for the uplift of the Qinghai-Tibetan Plateau

Molecular phylogeny of three genera containing nine species and subspecies of the specialized schizothoracine fishes are investigated based on the complete nucleotide sequence of mitochondrial cytochrome b gene. Meantime relationships between the main cladogenetic events of the specialized schizothoracine fishes and the stepwise uplift of the Qinghai-Tibetan Plateau are also conducted using the molecular clock, which is calibrated by geological isolated events between the upper reaches of the Yellow River and the Qinghai Lake. Results indicated that the specialized schizothoracine fishes are not a monophyly. Five species and subspecies of Ptychobarbus form a monophyly. But three species of Gymnodiptychus do not form a monophyly. Gd. integrigymnatus is a sister taxon of the highly specialized schizothoracine fishes while Gd. pachycheilus has a close relation with Gd. dybowskii, and both of them are as a sister group of Diptychus maculatus. The specialized schizothoracines fishes might have originated during the Miocene (about 10 MaBP), and then the divergence of three genera happened during late Miocene (about 8 MaBP). Their main specialization occurred during the late Pliocene and Pleistocene (3.54-0.42 MaBP). The main cladogenetic events of the specialized schizothoracine fishes are mostly correlated with the geological tectonic events and intensive climate shift happened at 8, 3.6, 2.5 and 1.7 MaBP of the late Cenozoic. Molecular clock data do not support the hypothesis that the Qinghai-Tibetan Plateau uplifted to near present or even higher elevations during the Oligocene or Miocene, and neither in agreement with the view that the plateau uplifting reached only to an altitude of 2000 in during the late Pliocene (about 2.6 MaBP).

Molecular systematics and biogeography of Crawfurdia, Metagentiana and Tripterospermum (Gentianaceae) based on nuclear ribosomal and plastid DNA sequences

Background and Aims The systematic position of the genus Metagentiana and its phylogenetic relationships with Crawfurdia, Gentiana and Tripterospermum have not been explicitly addressed. These four genera belong to one of two subtribes (Gentianinae) of Gentianeae. The aim of this paper is to examine the systematic position of Crawfurdia, Metagentiana and Tripterospermum and to clarify their phylogenetic affinities more clearly using ITS and trnL intron sequences.Methods Nucleotide sequences from the internal transcribed spacers (ITS) of nuclear ribosomal DNA and the plastid DNA trnL (UAA) intron were analysed phylogenetically. Ten of fourteen Metagentiana species were sampled, together with 40 species of other genera in the subtribe Gentianinae.Key Results The data support several previously published conclusions relating to the separation of Metagentiana from Gentiana and its closer relationships to Crawfurdia and Tripterospermum based on studies of gross morphology, floral anatomy, chromosomes, palynology, embryology and previous molecular data. The molecular clock hypothesis for the tested sequences in subtribe Gentianinae was not supported by the data (P < 0.05), so the clock-independent non-parametric rate smoothing method was used to estimate divergence time. This indicates that the separation of Crawfurdia, Metagentiana and Tripterospermum from Gentiana occurred about 11.4-21.4 Mya (million years ago), and the current species of these three genera diverged at times ranging from 0.4 to 6.2 Mya.Conclusions The molecular analyses revealed that Crawfurdia, Metagentiana and Tripterospermum do not merit status as three separate genera, because sampled species of Crawfurdia and Tripterospermum are embedded within Metagentiana. The speciation and rapid radiation of these three genera is likely to have occurred in western China as a result of upthrust of the Himalayas during the late Miocene and the Pleistocene.