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.

Molecules, morphology, and ecology indicate a recent, amphibious ancestry for echidnas

The semiaquatic platypus and terrestrial echidnas (spiny anteaters) are the only living egg-laying mammals (monotremes). The fossil record has provided few clues as to their origins and the evolution of their ecological specializations; however, recent reassignment of the Early Cretaceous Teinolophos and Steropodon to the platypus lineage implies that platypuses and echidnas diverged >112.5 million years ago, reinforcing the notion of monotremes as living fossils. This placement is based primarily on characters related to a single feature, the enlarged mandibular canal, which supplies blood vessels and dense electrosensory receptors to the platypus bill. Our reevaluation of the morphological data instead groups platypus and echidnas to the exclusion of Teinolophos and Steropodon and suggests that an enlarged mandibular canal is ancestral for monotremes (partly reversed in echidnas, in association with general mandibular reduction). A multigene evaluation of the echidna–platypus divergence using both a relaxed molecular clock and direct fossil calibrations reveals a recent split of 19–48 million years ago. Platypus-like monotremes (Monotrematum) predate this divergence, indicating that echidnas had aquatically foraging ancestors that reinvaded terrestrial ecosystems. This ecological shift and the associated radiation of echidnas represent a recent expansion of niche space despite potential competition from marsupials. Monotremes might have survived the invasion of marsupials into Australasia by exploiting ecological niches in which marsupials are restricted by their reproductive mode. Morphology, ecology, and molecular biology together indicate that Teinolophos and Steropodon are basal monotremes rather than platypus relatives, and that living monotremes are a relatively recent radiation.

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).

The linking of plate tectonics and evolutionary divergence

It is exciting to be living at a time when the big questions in biology can be investigated using modern genetics and computing [1]. Bauzà-Ribot et al.[2] take on one of the fundamental drivers of biodiversity, the effect of continental drift in the formation of the world’s biota 3 and 4, employing next-generation sequencing of whole mitochondrial genomes and modern Bayesian relaxed molecular clock analysis. Bauzà-Ribot et al.[2] conclude that vicariance via plate tectonics best explains the genetic divergence between subterranean metacrangonyctid amphipods currently found on islands separated by the Atlantic Ocean. This finding is a big deal in biogeography, and science generally [3], as many other presumed biotic tectonic divergences have been explained as probably due to more recent transoceanic dispersal events [4]. However, molecular clocks can be problematic 5 and 6 and we have identified three issues with the analyses of Bauzà-Ribot et al.[2] that cast serious doubt on their results and conclusions. When we reanalyzed their mitochondrial data and attempted to account for problems with calibration 5 and 6, modeling rates across branches 5 and 7 and substitution saturation [5], we inferred a much younger date for their key node. This implies either a later trans-Atlantic dispersal of these crustaceans, or more likely a series of later invasions of freshwaters from a common marine ancestor, but either way probably not ancient tectonic plate movements.

Cloning and characterization of microRNAs from Brassica napus

A library containing approximately 40,000 small RNA sequences was constructed for Brassica napus. Analysis of 3025 sequences obtained from this library resulted in the identification of 11 conserved miRNA families, which were validated by secondary structure prediction using surrounding sequences in the Brassica genome. Two 21 nt small RNA sequences reside within the arm of a pre-miRNA like stem-loop structure, making them likely candidates for novel non-conserved miRNAs in B. napus. Most of the conserved miRNAs were expressed at similar levels in a F1 hybrid B. napus line and its four double haploid progeny that showed marked variations in phenotypes, but many were differentially expressed between B. napus and Arabidopsis. The miR169 family was expressed at high levels in young leaves and stems, but was undetectable in roots and mature leaves, suggesting that miR169 expression is developmentally regulated in B. napus. © 2007 Federation of European Biochemical Societies.

The evolutionary history of termites as inferred from 66 mitochondrial genomes

Termites have colonized many habitats and are among the most abundant animals in tropical ecosystems, which they modify considerably through their actions. The timing of their rise in abundance and of the dispersal events that gave rise to modern termite lineages is not well understood. To shed light on termite origins and diversification, we sequenced the mitochondrial genome of 48 termite species and combined them with 18 previously sequenced termite mitochondrial genomes for phylogenetic and molecular clock analyses using multiple fossil calibrations. The 66 genomes represent most major clades of termites. Unlike previous phylogenetic studies based on fewer molecular data, our phylogenetic tree is fully resolved for the lower termites. The phylogenetic positions of Macrotermitinae and Apicotermitinae are also resolved as the basal groups in the higher termites, but in the crown termitid groups, including Termitinae + Syntermitinae + Nasutitermitinae + Cubitermitinae, the position of some nodes remains uncertain. Our molecular clock tree indicates that the lineages leading to termites and Cryptocercus roaches diverged 170 Ma (153-196 Ma 95% confidence interval [CI]), that modern Termitidae arose 54 Ma (46-66 Ma 95% CI), and that the crown termitid group arose 40 Ma (35-49 Ma 95% CI). This indicates that the distribution of basal termite clades was influenced by the final stages of the breakup of Pangaea. Our inference of ancestral geographic ranges shows that the Termitidae, which includes more than 75% of extant termite species, most likely originated in Africa or Asia, and acquired their pantropical distribution after a series of dispersal and subsequent diversification events.

Comment on “Whole-genome analyses resolve early branches in the tree of life of modern birds”

Jarvis et al. (Research Articles, 12 December 2014, p. 1320) presented molecular clock analyses that suggested that most modern bird orders diverged just after the mass extinction event at the Cretaceous-Paleogene boundary (about 66 million years ago). We demonstrate that this conclusion results from the use of a single inappropriate maximum bound, which effectively precludes the Cretaceous diversification overwhelmingly supported by previous molecular studies.

苔藓植物一些东亚特有属的分子系统学研究

苔藓是高等植物（有胚植物或陆地植物）中最原始的一类，但种类却丰富多样，其形态和生长环境的多样化程度高于蕨类和裸子植物，且对极端环境的忍耐力更强，分布范围也更广。“特有”是一个地理概念，它是相对广布而言，当一个类群的分布范围有一定的限制时即为特有现象。“东亚特有”是指分布范围主要局限于中国，朝鲜，日本和蒙古等，向北可及俄罗斯远东地区，少数可分布至中国南部相邻地区的植物类群。东亚地区主要以温带植物区系为主，但也包含一些热带植物区系成分，还因为第四纪以来受冰川活动影响较少，因此植物种类非常丰富。东亚地区也是苔藓植物的多样性中心之一，这里有较多的特有成分。在我国总共分布有苔藓植物东亚特有属35属，其中苔类5属，藓类30属。长期以来，特有成分始终引起人们的极大关注，不仅是因为其在植物地理学上的重要性，还因为特有类群中包含了孓遗类群，往往系统位置比较关键，此外，大部分特有类群对人为干扰比较敏感，对其保护就愈加重要，因为它在这个地区的消失就意味着一个类群的灭绝。 我国对苔藓植物东亚特有类群已有较好的认识，在前人知识积累的基础之上，我们期望通过分子系统学的方法，开展对东亚特有苔藓属的研究，逐步揭开特有属植物的神秘面纱，最终在系统树上找到它们各自应该属于自己的位置。 在本次研究中，我们总共得到十一个苔藓植物东亚特有属的新鲜材料。在实验室中我们对这十一个特有属叶绿体和核的六个基因（叶绿体atpB, rbcL, cp-SSU, cp-LSU 和核18S，26S rDNA）进行了测序，并在此基础之上，构建了来自苔藓植物106个属上述六个基因的联合矩阵，并对它们进行了系统学分析。本文所选十一个特有属中除三个苔类属和一个线齿藓类的属之外，其它七个特有属都属于侧蒴藓类。根据近几年的研究结果，侧蒴藓类中灰藓目被认为是起源自一次快速辐射演化，灰藓目各科之间的关系以及各科的范围都很难确定。即便本实验测序一万多bp，这一支之内的关系仍不能解决。 在以上结果的基础上，本文对线齿藓类的树发藓属（Microdendron）进行了较为详细的研究，我们用最大简约法分析了金发藓目15属，33种的18S, rbcL和trnL-F序列的联合矩阵。对树发藓属的微形态进行了电镜扫描。形态和分子数据的分析结果表明，这个特有属在属级水平是不成立的，它仅是小金发藓属的一个种。此结果支持将这个东亚特有属降为种的等级。此外，本文还对囊绒苔属（Trichocoleopsis）和新绒苔属（Neotrichocolea）的系统位置做了比较详细的研究。我们分别分析了一个苔类植物57属的四基因（cp-SSU, cp-LSU, atpB and rbcL）矩阵和一个苔类植物24属的九基因（cp-SSU, cp-LSU, atpB, psbA, rps4, rbcL, 18S, 26S and nad5）联合矩阵，结果显示囊绒苔属和新绒苔属互为姐妹群关系，而毛叶苔属（Ptilidium）又是它们二者的姐妹群。研究结果支持了囊绒苔属和新绒苔属组成新绒苔科（Neotrichocoleaceae），而不同于前人的观点：将上述两属放置于毛叶苔科（Ptilidiaceae）、绒苔科（Trichocoleaceae）或多囊苔科（Lepidolaenaceae）。另外值得注意的是这两个特有属和毛叶苔属组成的一支位于叶苔类（Leafy liverwort）中“Leafy I”和“Leafy II”两大支之间，但这一支确切的系统位置没有解决，仍有待于进一步研究。 除此之外，本文还利用GenBank中的数据对东亚特有属日鳞苔属（Nipponolejeunea）和耳坠苔属（Ascidiota）（未获得实验材料）进行了初步的系统学分析。结果表明传统上放在细鳞苔科的日鳞苔属与毛耳苔科的毛耳苔属（Jubula）为姐妹群关系，建议将日鳞苔属置于毛耳苔科;耳坠苔属是光萼苔科的成员，属的分类等级是合理的。 最后本文利用罚分似然法，选取多个化石作为标定点，对来自苔藓植物主要类群及其它陆地植物共115个类群5个基因（atpB, rbcL, cp-SSU, cp-LSU, 18S）的矩阵进行了分子钟的分析，初步估算11个东亚特有属的分化时间。

Phylogenetic study of complete cytochrome b genes in musk deer (genus Moschus) using museum samples

As an endangered animal group, musk deer (genus Moschus) are not only a great concern of wildlife conservation, but also of special interest to evolutionary studies due to long-standing arguments on the taxonomic and phylogenetic associations in this group. Using museum samples, we sequenced complete mitochondrial cytochrome b genes (1140 bp) of all suggested species of musk deer in order to reconstruct their phylogenetic history through molecular information. Our results showed that the cytochrome b gene tree is rather robust and concurred for all the algorithms employed (parsimony, maximum likelihood, and distance methods). Further, the relative rate test indicated a constant sequence substitution rate among all the species, permitting the dating of divergence events by molecular clock. According to the molecular topology, M. moschiferus branched off the earliest from a common ancestor of musk deer (about 700,000 years ago); then followed the bifurcation forming the M. berezouskii lineage and the lineage clustering M. fuscus, M. chrysogaster, and M. leucogaster (around 370,000 years before present), interestingly the most recent speciation event in musk deer happened rather recently (140,000 years ago), which might have resulted from the diversified habitats and geographic barriers in southwest China caused by gigantic movements of the Qinghai-Tibetan Plateau in history. Combining the data of current distributions, fossil records, and molecular data of this study, we suggest that the historical dispersion of musk deer might be from north to south in China. Additionally, in our further analyses involving other pecora species, musk deer was strongly supported as a monophyletic group and a valid family in Artiodactyla, closely related to Cervidae. (C) 1999 Academic Press.

A multi-locus phylogeny of Nectogalini shrews and influences of the paleoclimate on speciation and evolution

Nectogaline shrews are a major component of the small mammalian fauna of Europe and Asia, and are notable for their diverse ecology, including utilization of aquatic habitats. So far, molecular phylogenetic analyses including nectogaline species have been unable to infer a well-resolved, well-supported phylogeny, thus limiting the power of comparative evolutionary and ecological analyses of the group. Here, we employ Bayesian phylogenetic analyses of eight mitochondrial and three nuclear genes to infer the phylogenetic relationships of nectogaline shrews. We subsequently use this phylogeny to assess the genetic diversity within the genus Episoriculus, and determine whether adaptation to aquatic habitats evolved independently multiple times. Moreover, we both analyze the fossil record and employ Bayesian relaxed clock divergence dating analyses of DNA to assess the impact of historical global climate change on the biogeography of Nectogalini. We infer strong support for the polyphyly of the genus Episoriculus. We also find strong evidence that the ability to heavily utilize aquatic habitats evolved independently in both Neomys and Chimarrogale + Nectogale lineages. Our Bayesian molecular divergence analysis suggests that the early history of Nectogalini is characterized by a rapid radiation at the Miocene/Pliocene boundary, thus potentially explaining the lack of resolution at the base of the tree. Finally, we find evidence that nectogalines once inhabited northern latitudes, but the global cooling and desiccating events at the Miocene/Pliocene and Pliocene/Pleistocene boundaries and Pleistocene glaciation resulted in the migration of most Nectogalini lineages to their present day southern distribution. (C) 2010 Elsevier Inc. All rights reserved.

Speciation and Phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: Analysis of the Cytochrome b Gene of mtDNA from Diverse Populations

Gao-Yan Li, Xu-Zhen Wang, Ya-Hui Zhao, Jie Zhang, Chun-Guang Zhang, and Shun-Ping He (2009) Speciation and phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: analysis of the cytochrome b gene of mtDNA from diverse populations. Zoological Studies 48(4): 569-583. The cyprinid fish Opsariichthys bidens Gunther is distributed in all major river systems of continental East Asia, and represents an attractive model for phylogeographic studies among cyprinid species or within a given species. In this study, we investigated the phylogeographic and demographic history of this species, using partial sequences of the cytochrome (cyt) b gene in mitochondrial (mt)DNA. Fish samples were collected from almost all major river systems where O. bidens is distributed in China. Sequence analysis showed remarkably high polymorphism, with 125 haplotypes in the 234 specimens examined, and with 89.8% of haplotypes occurring in only 1 specimen. A neutrality test indicated that some groups were not at mutation-drift equilibrium, suggesting a past population expansion. These results were supported by a mismatch distribution analysis. Based on our analysis, O. bidens consists of 4 groups belonging to 2 clades. The divergence time of the 2 clades was estimated to be 11.06-8.04 my. This value corresponds to the time of the 2nd uplift of the Qinghai-Tibet Plateau, the emergence of the East Asian monsoon, and the Epoch-6 Event. A two species scheme is proposed. http://zoolstud.sinica.edu.tw/Journals/48.4/569.pdf

Bayesian mixed models and divergence time estimation of Chinese cavefishes (Cyprinidae : Sinocyclocheilus)

The genus Sinocyclocheilus is distributed in Yun-Gui Plateau and its surrounding region only, within more than 10 cave species showing different degrees of degeneration of eyes and pigmentation with wonderful adaptations. To present, published morphological and molecular phylogenetic hypotheses of Sinocyclocheilus from prior works are very different and the relationships within the genus are still far from clear. We obtained the sequences of cytochrome b (cyt b) and NADH dehydrogenase subunit 4 (ND4) of 34 species within Sinocyclocheilus, which represent the most dense taxon sampling to date. We performed Bayesian mixed models analyses with this data set. Under this phylogenetic framework, we estimated the divergence times of recovered clades using different methods under relaxed molecular clock. Our phyloegentic results supported the monophyly of Sinocyclocheilus and showed that this genus could be subdivided into 6 major clades. In addition, an earlier finding demonstrating the polyphyletic of cave species and the most basal position of S. jii was corroborated. Relaxed divergence-time estimation suggested that Sinocyclocheilus originated at the late Miocene, about 11 million years ago (Ma), which is older than what have been assumed.

Phylogeny and speciation of the eastern Asian cyprinid genus Sarcocheilichthys

The genus Sarcocheilichthys is a group of small cyprinid fishes comprising 10 species/sub-species widely distributed in East Asia, which represents a valuable model for understanding the speciation of freshwater fishes in East Asia. In the present study, the molecular phylogenetic relationship of the genus Sarcocheilichthys was investigated using a 1140 bp section of the mitochondrial cytochrome b gene. Two different tree-building methods, maximum parsimony (MP) and Bayesian methods, yielded trees with almost the same topology, yielding high bootstrap values or posterior probabilities. The results showed that the genus Sarcocheilichthys consists of two large clades, clades I and II. Clade I contains Sarcocheilichthys lacustris, Sarcocheilichthys sinensis and Sarcocheilichthys parvus, with S. parvus at a basal position. In clade II, Sarcocheilichthys variegatus microoculus is at a basal position; samples of the widespread species, Sarcocheilichthys nigripinnis, form a large subclade containing another valid species Sarcocheilichthys czerskii. Sarcocheilichthys kiangsiensis is retained at an intermediate position. Since S. czerskii is a valid species in the S. nigripinnis clade, remaining samples of S. nigripinnis form a paraphyly. This speciation process is attributed to geographical isolation and special environmental conditions experienced by S. czerskii and stable environments experienced by the other S. nigripinnis populations. This type of speciation process was suggested to be very common. Samples of Sarcocheilichthys sinensis sinensis and Sarcocheilichthys sinensis fukiensis that did not form their own monophyletic groups suggest an early stage of speciation and support their sub-species status. Molecular clock analysis indicates that the two major lineages of the genus Sarcocheilichthys, clades I and II diverged c. 8.89 million years ago (mya). Sarcocheilichthys v. microoculus from Japan probably diverged 4.78 mya from the Chinese group. The northern-southern clades of S. nigripinnis began to diverge c. 2.12 mya, while one lineage of S. nigripinnis evolved into a new species, S. czerski, c. 0.34 mya. (C) 2008 The Authors Journal compilation (C) 2008 The Fisheries Society of the British Isles.

Sox genes in grass carp (Ctenopharyngodon idella) with their implications for genome duplication and evolution

The Sox gene family is found in a broad range of animal taxa and encodes important gene regulatory proteins involved in a variety of developmental processes. We have obtained clones representing the HMG boxes of twelve Sox genes from grass carp (Ctenopharyngodon idella), one of the four major domestic carps in China. The cloned Sox genes belong to group B1, B2 and C. Our analyses show that whereas the human genome contains a single copy of Sox4, Sox11 and Sox14, each of these genes has two co-orthologs in grass carp, and the duplication of Sox4 and Sox11 occurred before the divergence of grass carp and zebrafish, which support the "fish-specific whole-genome duplication" theory. An estimation for the origin of grass carp based on the molecular clock using Sox1, Sox3 and Sox11 genes as markers indicates that grass carp (subfamily Leuciscinae) and zebrafish (subfamily Danioninae) diverged approximately 60 million years ago. The potential uses of Sox genes as markers in revealing the evolutionary history of grass carp are discussed.

Phylogeny and biogeography of Chinese sisorid catfishes re-examined using mitochondrial cytochrome b and 16S rRNA gene sequences

The family Sisoridae is one of the largest and most diverse Asiatic catfish families, most species occurring in the water systems of the Qinhai-Tibetan Plateau and East Himalayas. To date published morphological and molecular phylogenetics hypotheses of sisorid catfishes are part congruent, and there are some areas of significant disagreement with respect to intergeneric relationships. We used mitochondrial cytochrome b and 16S rRNA gene sequences to clarify existing gaps in phylogenetics and to test conflicting vicariant and dispersal biogeographical hypotheses of Chinese sisorids using dispersal-vicariance analysis and weighted ancestral area analysis in combination with palaeogeographical data as well as molecular clock calibration. Our results suggest that: (1) Chinese sisorid catfishes form a monophyletic group with two distinct clades, one represented by (Gagata (Bagarius, Glyptothorax)) and the other by (glyptosternoids, Pseudecheneis); (2) the glyptosternoid is a monophyletic group and Glyptosternum, Glaridoglanis, and Exostoma are three basal species having a primitive position among it; (3) a hypothesis referring to Pseudecheneis as the sister group of the glyptosternoids, based on morphological evidence, is supported; (4) the genus Pareuchiloglanis, as presently defined, is not monophyletic; (5) congruent with previous hypotheses, the uplift of Qinghai-Tibetan Plateau played a primary role in the speciation and radiation of the Chinese sisorids; and (6) an evolutionary scenario combining aspects of both vicariance and dispersal theory is necessary to explain the distribution pattern of the glyptosternoids. In addition, using a cytochrome b substitution rate of 0.91% per million years and 0.23% for 16S rRNA, we tentatively date that the glyptosternoids most possibly originated in Oligocene-Miocene boundary (19-24Myr), and radiated from Miocene to Pleistocene, along with a center of origin in the Irrawaddy-Tsangpo drainages and several rapid speciation in a relatively short time. (c) 2005 Elsevier Inc. All rights reserved.